A structural biology community assessment of AlphaFold2 applications

Comprehensive applications of the artificial intelligence technology in new drug research and development

Purpose

Target-based strategy is a prevalent means of drug research and development (R&D), since targets provide effector molecules of drug action and offer the foundation of pharmacological investigation. Recently, the artificial intelligence (AI) technology has been utilized in various stages of drug R&D, where AI-assisted experimental methods show higher efficiency than sole experimental ones. It is a critical need to give a comprehensive review of AI applications in drug R &D for biopharmaceutical field.

Methods

Relevant literatures about AI-assisted drug R&D were collected from the public databases (Including Google Scholar, Web of Science, PubMed, IEEE Xplore Digital Library, Springer, and ScienceDirect) through a keyword searching strategy with the following terms [("Artificial Intelligence" OR "Knowledge Graph" OR "Machine Learning") AND ("Drug Target Identification" OR "New Drug Development")].

Results

In this review, we first introduced common strategies and novel trends of drug R&D, followed by characteristic description of AI algorithms widely used in drug R&D. Subsequently, we depicted detailed applications of AI algorithms in target identification, lead compound identification and optimization, drug repurposing, and drug analytical platform construction. Finally, we discussed the challenges and prospects of AI-assisted methods for drug discovery.

Conclusion

Collectively, this review provides comprehensive overview of AI applications in drug R&D and presents future perspectives for biopharmaceutical field, which may promote the development of drug industry.

Inhibitory effect of truncated isoforms on GPCR dimerization predicted by combinatorial computational strategy

G protein-coupled receptors (GPCRs) play a pivotal role in fundamental biological processes and disease development. GPCR isoforms, derived from alternative splicing, can exhibit distinct signaling patterns. Some highly-truncated isoforms can impact functional performance of full-length receptors, suggesting their intriguing regulatory roles. However, how these truncated isoforms interact with full-length counterparts remains largely unexplored. Here, we computationally investigated the interaction patterns of three human GPCRs from three different classes, ADORA1 (Class A), mGlu2 (Class C) and SMO (Class F) with their respective truncated isoforms because their homodimer structures have been experimentally determined, and they have truncated isoforms deposited and identified at protein level in Uniprot database. Combining the neural network-based AlphaFold2 and two physics-based protein-protein docking tools, we generated multiple complex structures and assessed the binding affinity in the context of atomistic molecular dynamics simulations. Our computational results suggested all the four studied truncated isoforms showed potent binding to their counterparts and overlapping interfaces with homodimers, indicating their strong potential to block homodimerization of their counterparts. Our study offers insights into functional significance of GPCR truncated isoforms and supports the ubiquity of their regulatory roles.

A review of machine learning-based methods for predicting drug-target interactions

The prediction of drug-target interactions (DTI) is a crucial preliminary stage in drug discovery and development, given the substantial risk of failure and the prolonged validation period associated with in vitro and in vivo experiments. In the contemporary landscape, various machine learning-based methods have emerged as indispensable tools for DTI prediction. This paper begins by placing emphasis on the data representation employed by these methods, delineating five representations for drugs and four for proteins. The methods are then categorized into traditional machine learning-based approaches and deep learning-based ones, with a discussion of representative approaches in each category and the introduction of a novel taxonomy for deep neural network models in DTI prediction. Additionally, we present a synthesis of commonly used datasets and evaluation metrics to facilitate practical implementation. In conclusion, we address current challenges and outline potential future directions in this research field.

AlphaFold 2-based stacking model for protein solubility prediction and its transferability on seed storage proteins

Accurate protein solubility prediction is crucial in screening suitable candidates for food application. Existing models often rely only on sequences, overlooking important structural details. In this study, a regression model for protein solubility was developed using both the sequences and predicted structures of 2983 E. coli proteins. The sequence and structural level properties of the proteins were bioinformatically extracted and subjected to multilayer perceptron (MLP). Moreover, residue level features and contact maps were utilized to construct a graph convolutional network (GCN). The out-of-fold predictions of the two models were combined and fed into multiple meta-regressors to create a stacking model. The stacking model with support vector regressor (SVR) achieved R2 of 0.502 and 0.468 on test and external validation datasets, respectively, displaying higher performance compared to existing regression models. Based on the improved performance compared to its based models, the stacking model effectively captured the strength of its base models as well as the significance of the different features used. Furthermore, the model's transferability was indirectly validated on a dataset of seed storage proteins using Osborne definition as well as on a case study using molecular dynamic simulation, showing potential for application beyond microbial proteins to food and agriculture-related ones.

Prediction of protein structure and AI

AlphaFold, an artificial intelligence (AI)-based tool for predicting the 3D structure of proteins, is now widely recognized for its high accuracy and versatility in the folding of human proteins. AlphaFold is useful for understanding structure-function relationships from protein 3D structure models and can serve as a template or a reference for experimental structural analysis including X-ray crystallography, NMR and cryo-EM analysis. Its use is expanding among researchers, not only in structural biology but also in other research fields. Researchers are currently exploring the full potential of AlphaFold-generated protein models. Predicting disease severity caused by missense mutations is one such application. This article provides an overview of the 3D structural modeling of AlphaFold based on deep learning techniques and highlights the challenges in predicting the pathogenicity of missense mutations.

Plant protein-lipid interfaces studied by molecular dynamics simulations

The delineation of protein-lipid interfaces is essential for understanding the mechanisms of various membrane-associated processes crucial to plant development and growth, including signalling, trafficking, and membrane transport. Due to their highly dynamic nature, the precise characterization of lipid-protein interactions by experimental techniques is challenging. Molecular dynamics simulations provide a powerful computational alternative with a spatial-temporal resolution allowing the atomistic-level description. In this review, we aim to introduce plant scientists to molecular dynamics simulations. We describe different steps of performing molecular dynamics simulations and provide a broad survey of molecular dynamics studies investigating plant protein-lipid interfaces. Our aim is also to illustrate that combining molecular dynamics simulations with artificial intelligence-based protein structure determination opens up unprecedented possibilities for future investigations of dynamic plant protein-lipid interfaces.

The fitness cost of spurious phosphorylation

The fidelity of signal transduction requires the binding of regulatory molecules to their cognate targets. However, the crowded cell interior risks off-target interactions between proteins that are functionally unrelated. How such off-target interactions impact fitness is not generally known. Here, we use Saccharomyces cerevisiae to inducibly express tyrosine kinases. Because yeast lacks bona fide tyrosine kinases, the resulting tyrosine phosphorylation is biologically spurious. We engineered 44 yeast strains each expressing a tyrosine kinase, and quantitatively analysed their phosphoproteomes. This analysis resulted in ~30,000 phosphosites mapping to ~3500 proteins. The number of spurious pY sites generated correlates strongly with decreased growth, and we predict over 1000 pY events to be deleterious. However, we also find that many of the spurious pY sites have a negligible effect on fitness, possibly because of their low stoichiometry. This result is consistent with our evolutionary analyses demonstrating a lack of phosphotyrosine counter-selection in species with tyrosine kinases. Our results suggest that, alongside the risk for toxicity, the cell can tolerate a large degree of non-functional crosstalk as interaction networks evolve.

Artificial Intelligence Learns Protein Prediction

From AlphaGO over StableDiffusion to ChatGPT, the recent decade of exponential advances in artificial intelligence (AI) has been altering life. In parallel, advances in computational biology are beginning to decode the language of life: AlphaFold2 leaped forward in protein structure prediction, and protein language models (pLMs) replaced expertise and evolutionary information from multiple sequence alignments with information learned from reoccurring patterns in databases of billions of proteins without experimental annotations other than the amino acid sequences. None of those tools could have been developed 10 years ago; all will increase the wealth of experimental data and speed up the cycle from idea to proof. AI is affecting molecular and medical biology at giant steps, and the most important might be the leap toward more powerful protein design.

Molecular architecture of coronavirus double-membrane vesicle pore complex

Coronaviruses remodel the intracellular host membranes during replication, forming double-membrane vesicles (DMVs) to accommodate viral RNA synthesis and modifications1,2. SARS-CoV-2 non-structural protein 3 (nsp3) and nsp4 are the minimal viral components required to induce DMV formation and to form a double-membrane-spanning pore, essential for the transport of newly synthesized viral RNAs3-5. The mechanism of DMV pore complex formation remains unknown. Here we describe the molecular architecture of the SARS-CoV-2 nsp3-nsp4 pore complex, as resolved by cryogenic electron tomography and subtomogram averaging in isolated DMVs. The structures uncover an unexpected stoichiometry and topology of the nsp3-nsp4 pore complex comprising 12 copies each of nsp3 and nsp4, organized in 4 concentric stacking hexamer rings, mimicking a miniature nuclear pore complex. The transmembrane domains are interdigitated to create a high local curvature at the double-membrane junction, coupling double-membrane reorganization with pore formation. The ectodomains form extensive contacts in a pseudo-12-fold symmetry, belting the pore complex from the intermembrane space. A central positively charged ring of arginine residues coordinates the putative RNA translocation, essential for virus replication. Our work establishes a framework for understanding DMV pore formation and RNA translocation, providing a structural basis for the development of new antiviral strategies to combat coronavirus infection.

flDPnn2: Accurate and Fast Predictor of Intrinsic Disorder in Proteins

Prediction of the intrinsic disorder in protein sequences is an active research area, with well over 100 predictors that were released to date. These efforts are motivated by the functional importance and high levels of abundance of intrinsic disorder, combined with relatively low amounts of experimental annotations. The disorder predictors are periodically evaluated by independent assessors in the Critical Assessment of protein Intrinsic Disorder prediction (CAID) experiments. The recently completed CAID2 experiment assessed close to 40 state-of-the-art methods demonstrating that some of them produce accurate results. In particular, flDPnn2 method, which is the successor of flDPnn that performed well in the CAID1 experiment, secured the overall most accurate results on the Disorder-NOX dataset in CAID2. flDPnn2 implements a number of improvements when compared to its predecessor including changes to the inputs, increased size of the deep network model that we retrained on a larger training set, and addition of an alignment module. Using results from CAID2, we show that flDPnn2 produces accurate predictions very quickly, modestly improving over the accuracy of flDPnn and reducing the runtime by half, to about 27 s per protein. flDPnn2 is freely available as a convenient web server at http://biomine.cs.vcu.edu/servers/flDPnn2/.

AI-Predicted Protein Deformation Encodes Energy Landscape Perturbation

AI algorithms have proven to be excellent predictors of protein structure, but whether and how much these algorithms can capture the underlying physics remains an open question. Here, we aim to test this question using the Alphafold2 (AF) algorithm: We use AF to predict the subtle structural deformation induced by single mutations, quantified by strain, and compare with experimental datasets of corresponding perturbations in folding free energy ΔΔG. Unexpectedly, we find that physical strain alone-without any additional data or computation-correlates almost as well with ΔΔG as state-of-the-art energy-based and machine-learning predictors. This indicates that the AF-predicted structures alone encode fine details about the energy landscape. In particular, the structures encode significant information on stability, enough to estimate (de-)stabilizing effects of mutations, thus paving the way for the development of novel, structure-based stability predictors for protein design and evolution.

An end-to-end framework for the prediction of protein structure and fitness from single sequence

Significant research progress has been made in the field of protein structure and fitness prediction. Particularly, single-sequence-based structure prediction methods like ESMFold and OmegaFold achieve a balance between inference speed and prediction accuracy, showing promise for many downstream prediction tasks. Here, we propose SPIRED, a single-sequence-based structure prediction model that exhibits comparable performance to the state-of-the-art methods but with approximately 5-fold acceleration in inference and at least one order of magnitude reduction in training consumption. By integrating SPIRED with downstream neural networks, we compose an end-to-end framework named SPIRED-Fitness for the rapid prediction of both protein structure and fitness from single sequence with satisfactory accuracy. Moreover, SPIRED-Stab, the derivative of SPIRED-Fitness, achieves state-of-the-art performance in predicting the mutational effects on protein stability.

AlphaFold2 Reveals Structural Patterns of Seasonal Haplotype Diversification in SARS-CoV-2 Nucleocapsid Protein Variants

The COVID-19 pandemic saw the emergence of various Variants of Concern (VOCs) that took the world by storm, often replacing the ones that preceded them. The characteristic mutant constellations of these VOCs increased viral transmissibility and infectivity. Their origin and evolution remain puzzling. With the help of data mining efforts and the GISAID database, a chronology of 22 haplotypes described viral evolution up until 23 July 2023. Since the three-dimensional atomic structures of proteins corresponding to the identified haplotypes are not available, ab initio methods were here utilized. Regions of intrinsic disorder proved to be important for viral evolution, as evidenced by the targeted change to the nucleocapsid (N) protein at the sequence, structure, and biochemical levels. The linker region of the N-protein, which binds to the RNA genome and self-oligomerizes for efficient genome packaging, was greatly impacted by mutations throughout the pandemic, followed by changes in structure and intrinsic disorder. Remarkably, VOC constellations acted co-operatively to balance the more extreme effects of individual haplotypes. Our strategy of mapping the dynamic evolutionary landscape of genetically linked mutations to the N-protein structure demonstrates the utility of ab initio modeling and deep learning tools for therapeutic intervention.

AlphaFold two years on: Validation and impact

Two years on from the initial release of AlphaFold, we have seen its widespread adoption as a structure prediction tool. Here, we discuss some of the latest work based on AlphaFold, with a particular focus on its use within the structural biology community. This encompasses use cases like speeding up structure determination itself, enabling new computational studies, and building new tools and workflows. We also look at the ongoing validation of AlphaFold, as its predictions continue to be compared against large numbers of experimental structures to further delineate the model's capabilities and limitations.

The combination of Asp519Val/Glu665Val and Lys1813Ala mutations in FVIII markedly increases coagulation potential

ABSTRACT

A2 domain dissociation in activated factor VIII (FVIIIa) results in reduced activity. Previous studies demonstrated that some FVIII mutants (D519V/E665V and K1813A) with delayed A2 dissociation enhanced coagulation potential. We speculated, therefore, that FVIII encompassing a combination of these mutations might further enhance coagulant activity. The aim was to assess the D519V/E665V/K1813A-FVIII mutation as a gain of function. The FVIII mutants, D519V/E665V/K1813A, D519V/E665V, and K1813A were expressed in a baby hamster kidney cell system, and global coagulation potential of these mutants was compared with wild-type (WT) FVIII in vitro and in hemophilia A mice in vivo. Kinetic analyses indicated that the apparent Kd for FIXa on the tenase assembly with D519V/E665V and D519V/E665V/K1813A mutants were lower, and that the generated FXa for D519V/E665V/K1813A was significantly greater than WT-FVIII. WT-FVIII activity after thrombin activation increased by ∼12-fold within 5 minutes, and returned to initial levels within 30 minutes. In contrast, The FVIII-related activity of D519V/E665V/K1813A increased further with time after thrombin activation, and showed an ∼25-fold increase at 2 hours. The A2 dissociation rate of D519V/E665V/K1813A was ∼50-fold slower than the WT in a 1-stage clotting assay. Thrombin generation assays demonstrated that D519V/E665V/K1813A (0.125 nM) exhibited coagulation potential comparable with that of the WT (1 nM). In animal studies, rotational thromboelastometry and tail-clip assays showed that the coagulation potential of D519V/E665V/K1813A (0.25 μg/kg) was equal to that of the WT (2 μg/kg). FVIII-D519V/E665V/K1813A mutant could provide an approximately eightfold increase in hemostatic function of WT-FVIII because of increased FVIIIa stability and the association between FVIIIa and FIXa.

Comparative structural studies on Bovine papillomavirus E6 oncoproteins: Novel insights into viral infection and cell transformation from homology modeling and molecular dynamics simulations

Bovine papillomavirus (BPV) infects cattle cells worldwide, leading to hyperproliferative lesions and the potential development of cancer, driven by E5, E6, and E7 oncoproteins along with other cofactors. E6 oncoprotein binds experimentally to various proteins, primarily paxillin and MAML1, as well as hMCM7 and CBP/p300. However, the molecular and structural mechanisms underlying BPV-induced malignant transformation remain unclear. Therefore, we have modeled the E6 oncoprotein structure from non-oncogenic BPV-5 and compared them with oncogenic BPV-1 to assess the relationship between structural features and oncogenic potential. Our analysis elucidated crucial structural aspects of E6, highlighting both conserved elements across genotypes and genotype-specific variations potentially implicated in the oncogenic process, particularly concerning primary target interactions. Additionally, we predicted the location of the hMCM7 binding site on the N-terminal of BPV-5 E6. This study enhances our understanding of the structural characteristics of BPV E6 oncoproteins and their interactions with host proteins, clarifying structural differences and similarities between high and low-risk BPVs. This is important to understand better the mechanisms involved in cell transformation in BPV infection, which could be used as a possible target for therapy.

AlphaFold accelerated discovery of psychotropic agonists targeting the trace amine-associated receptor 1

Artificial intelligence is revolutionizing protein structure prediction, providing unprecedented opportunities for drug design. To assess the potential impact on ligand discovery, we compared virtual screens using protein structures generated by the AlphaFold machine learning method and traditional homology modeling. More than 16 million compounds were docked to models of the trace amine-associated receptor 1 (TAAR1), a G protein-coupled receptor of unknown structure and target for treating neuropsychiatric disorders. Sets of 30 and 32 highly ranked compounds from the AlphaFold and homology model screens, respectively, were experimentally evaluated. Of these, 25 were TAAR1 agonists with potencies ranging from 12 to 0.03 μM. The AlphaFold screen yielded a more than twofold higher hit rate (60%) than the homology model and discovered the most potent agonists. A TAAR1 agonist with a promising selectivity profile and drug-like properties showed physiological and antipsychotic-like effects in wild-type but not in TAAR1 knockout mice. These results demonstrate that AlphaFold structures can accelerate drug discovery.

Multicopper Oxidase from Lactobacillus hilgardii: Mechanism of Degradation of Tyramine and Phenylethylamine in Fermented Food

Elevated levels of biogenic amines (BAs) in fermented food can have negative effects on both the flavor and health. Mining enzymes that degrade BAs is an effective strategy for controlling their content. The study screened a strain of Lactobacillus hilgardii 1614 from fermented food system that can degrade BAs. The multiple copper oxidase genes LHMCO1614 were successfully mined after the whole genome protein sequences of homologous strains were clustered and followed by homology modeling. The enzyme molecules can interact with BAs to stabilize composite structures for catalytic degradation, as shown by molecular docking results. Ingeniously, the kinetic data showed that purified LHMCO1614 was less sensitive to the substrate inhibition of tyramine and phenylethylamine. The degradation rates of tyramine and phenylethylamine in huangjiu (18% vol) after adding LHMCO1614 were 41.35 and 40.21%, respectively. Furthermore, LHMCO1614 demonstrated universality in degrading tyramine and phenylethylamine present in other fermented foods as well. HS-SPME-GC-MS analysis revealed that, except for aldehydes, the addition of enzyme treatment did not significantly alter the levels of major flavor compounds in enzymatically treated fermented foods (p > 0.05). This study presents an enzymatic approach for regulating tyramine and phenylethylamine levels in fermented foods with potential applications both targeted and universal.

The unique allosteric property of crocodilian haemoglobin elucidated by cryo-EM

The principal effect controlling the oxygen affinity of vertebrate haemoglobins (Hbs) is the allosteric switch between R and T forms with relatively high and low oxygen affinity respectively. Uniquely among jawed vertebrates, crocodilians possess Hb that shows a profound drop in oxygen affinity in the presence of bicarbonate ions. This allows them to stay underwater for extended periods by consuming almost all the oxygen present in the blood-stream, as metabolism releases carbon dioxide, whose conversion to bicarbonate and hydrogen ions is catalysed by carbonic anhydrase. Despite the apparent universal utility of bicarbonate as an allosteric regulator of Hb, this property evolved only in crocodilians. We report here the molecular structures of both human and a crocodilian Hb in the deoxy and liganded states, solved by cryo-electron microscopy. We reveal the precise interactions between two bicarbonate ions and the crocodilian protein at symmetry-related sites found only in the T state. No other known effector of vertebrate Hbs binds anywhere near these sites.

I'm Walking into Spiderwebs: Making Sense of Protein-Protein Interaction Data

Protein-protein interactions (PPIs) are at the heart of the molecular landscape permeating life. Proteomics studies can explore this protein interaction landscape using mass spectrometry (MS). Thanks to their high sensitivity, mass spectrometers can easily identify thousands of proteins within a single sample, but that same sensitivity generates tangled spiderwebs of data that hide biologically relevant findings. So, what does a researcher do when she finds herself walking into spiderwebs? In a field focused on discovery, MS data require rigor in their analysis, experimental validation, or a combination of both. In this Review, we provide a brief primer on MS-based experimental methods to identify PPIs. We discuss approaches to analyze the resulting data and remove the proteomic background. We consider the advantages between comprehensive and targeted studies. We also discuss how scoring might be improved through AI-based protein structure information. Women have been essential to the development of proteomics, so we will specifically highlight work by women that has made this field thrive in recent years.

Generation of a high confidence set of domain-domain interface types to guide protein complex structure predictions by AlphaFold

MOTIVATION

While the release of AlphaFold (AF) represented a breakthrough for the prediction of protein complex structures, its sensitivity, especially when using full length protein sequences, still remains limited. Modeling success rates might increase if AF predictions were guided by likely interacting protein fragments. This approach requires available sets of highly confident protein-protein interface types. Computational resources, such as 3did, infer interacting globular domain types from observed contacts in protein structures. Assessing the accuracy of these predicted interface types is difficult because we lack hand-curated reference sets of verified domain-domain interface (DDI) types.

RESULTS

To improve protein complex modeling of DDIs by AF, we manually inspected 80 randomly selected DDI types from the 3did resource to generate a first reference set of DDI types. Identified cases of DDI type nonapproval (40%) primarily resulted from inaccurate Pfam domain matches, crystal contacts, and synthetic protein constructs. Using logistic regression, we predicted a subset of 2411 out of 5724 considered DDI types in 3did to be of high confidence, which we subsequently applied to 53 000 human-protein interactions to predict DDIs followed by AF modeling. We obtained highly confident AF models for 604 out of 1129 predicted DDIs. Of note, for 47% of them no confident AF structural model could be obtained using full length protein sequences.

AVAILABILITY AND IMPLEMENTATION

Code is available at https://github.com/KatjaLuckLab/DDI_manuscript.

Integrating Large-Scale Protein Structure Prediction into Human Genetics Research

The last five years have seen impressive progress in deep learning models applied to protein research. Most notably, sequence-based structure predictions have seen transformative gains in the form of AlphaFold2 and related approaches. Millions of missense protein variants in the human population lack annotations, and these computational methods are a valuable means to prioritize variants for further analysis. Here, we review the recent progress in deep learning models applied to the prediction of protein structure and protein variants, with particular emphasis on their implications for human genetics and health. Improved prediction of protein structures facilitates annotations of the impact of variants on protein stability, protein-protein interaction interfaces, and small-molecule binding pockets. Moreover, it contributes to the study of host-pathogen interactions and the characterization of protein function. As genome sequencing in large cohorts becomes increasingly prevalent, we believe that better integration of state-of-the-art protein informatics technologies into human genetics research is of paramount importance.

The power and pitfalls of AlphaFold2 for structure prediction beyond rigid globular proteins

Artificial intelligence-driven advances in protein structure prediction in recent years have raised the question: has the protein structure-prediction problem been solved? Here, with a focus on nonglobular proteins, we highlight the many strengths and potential weaknesses of DeepMind's AlphaFold2 in the context of its biological and therapeutic applications. We summarize the subtleties associated with evaluation of AlphaFold2 model quality and reliability using the predicted local distance difference test (pLDDT) and predicted aligned error (PAE) values. We highlight various classes of proteins that AlphaFold2 can be applied to and the caveats involved. Concrete examples of how AlphaFold2 models can be integrated with experimental data in the form of small-angle X-ray scattering (SAXS), solution NMR, cryo-electron microscopy (cryo-EM) and X-ray diffraction are discussed. Finally, we highlight the need to move beyond structure prediction of rigid, static structural snapshots toward conformational ensembles and alternate biologically relevant states. The overarching theme is that careful consideration is due when using AlphaFold2-generated models to generate testable hypotheses and structural models, rather than treating predicted models as de facto ground truth structures.

DrugDomain: The evolutionary context of drugs and small molecules bound to domains

Interactions between proteins and small organic compounds play a crucial role in regulating protein functions. These interactions can modulate various aspects of protein behavior, including enzymatic activity, signaling cascades, and structural stability. By binding to specific sites on proteins, small organic compounds can induce conformational changes, alter protein-protein interactions, or directly affect catalytic activity. Therefore, many drugs available on the market today are small molecules (72% of all approved drugs in the last 5 years). Proteins are composed of one or more domains: evolutionary units that convey function or fitness either singly or in concert with others. Understanding which domain(s) of the target protein binds to a drug can lead to additional opportunities for discovering novel targets. The evolutionary classification of protein domains (ECOD) classifies domains into an evolutionary hierarchy that focuses on distant homology. Previously, no structure-based protein domain classification existed that included information about both the interaction between small molecules or drugs and the structural domains of a target protein. This data is especially important for multidomain proteins and large complexes. Here, we present the DrugDomain database that reports the interaction between ECOD of human target proteins and DrugBank molecules and drugs. The pilot version of DrugDomain describes the interaction of 5160 DrugBank molecules associated with 2573 human proteins. It describes domains for all experimentally determined structures of these proteins and incorporates AlphaFold models when such structures are unavailable. The DrugDomain database is available online: http://prodata.swmed.edu/DrugDomain/.

Identification of flavor peptides based on virtual screening and molecular docking from Hypsizygus marmoreuss

Hypsizygus marmoreuss is an under-explored source of flavor peptides that can enhance the flavor of NaCl or MSG, allowing products to be reformulated in line with reduction policies. This study utilized advanced techniques, including UPLC-Q-TOF MS/MS and molecular docking, to identify H. marmoreuss peptides. Sensory evaluations revealed 10 peptides with pronounced umami flavors and seven with dominantly salty tastes. VLPVPQK scored highest for umami intensity (5.2), and EGNPAHQK for salty intensity (6.2). Further investigation influenced by 0.35 % MSG or 0.35 % NaCl exposed peptides with elevated umami and salty thresholds. LDSPATPEK, VVEGEPSLK, and QKLPEKPER had umami-enhancing thresholds of 0.18, 0.18, and 0.35 mM, while LDSPATPEK and VVEGEPSLK had similar thresholds for salt (0.09 mM). Molecular docking revealed that taste receptor proteins interacted with umami peptides through hydrogen, carbon-hydrogen, alkyl, and van der Waals forces. Specific amino acids in the umami receptor T1R1 had roles in bonding with umami peptides through hydrogen and carbon-hydrogen interactions. In conclusion, molecular docking proved to be an effective and efficient method for flavor peptide screening. Further, this study demonstrated that flavor peptides from H. marmoreuss had the capacity to enhance NaCl and MSG flavours and might be useful tools for reformulation, reducing salt and MSG contents.

AlphaFold2 in biomedical research: facilitating the development of diagnostic strategies for disease

Proteins, as the primary executors of physiological activity, serve as a key factor in disease diagnosis and treatment. Research into their structures, functions, and interactions is essential to better understand disease mechanisms and potential therapies. DeepMind's AlphaFold2, a deep-learning protein structure prediction model, has proven to be remarkably accurate, and it is widely employed in various aspects of diagnostic research, such as the study of disease biomarkers, microorganism pathogenicity, antigen-antibody structures, and missense mutations. Thus, AlphaFold2 serves as an exceptional tool to bridge fundamental protein research with breakthroughs in disease diagnosis, developments in diagnostic strategies, and the design of novel therapeutic approaches and enhancements in precision medicine. This review outlines the architecture, highlights, and limitations of AlphaFold2, placing particular emphasis on its applications within diagnostic research grounded in disciplines such as immunology, biochemistry, molecular biology, and microbiology.

Three novel Er blood group system alleles and insights from protein modeling

BACKGROUND

The Er blood group system was recently shown to be defined by PIEZO1. The system consists of high prevalence antigens Era, Er3, ERSA, and ERAMA; and low prevalence antigen Erb. Era/Erb are antithetical with Er(a-b+) defined by the ER*B allele [c.7180G>A p.(Gly2394Ser)]. A nonsense variant c.5289C>G p.(Tyr1763*) is associated with a predicted Ernull phenotype, and a missense variant c.7174G>A p.(Glu2392Lys) in close proximity to p.2394 causes loss of both Era and Erb expression.

STUDY DESIGN AND METHODS

We investigated PIEZO1 in four Er(a-) individuals who presented with anti-Era. Whole genome sequencing (WGS) and Sanger sequencing were performed. The location and structural differences of predicted protein changes were visualized using the predicted 3-D structure of Piezo1 created using AlphaFold2.

RESULTS

One individual was homozygous for the reported ER*B. A second had a novel heterozygous nonsense variant c.3331C>T p.(Gln1111*), but a second allelic variant was not found. In the remaining two individuals, two different heterozygous novel missense variants, c.7184C>T p.(Ala2395Val) or c.7195G>A p.(Gly2399Ser), were in trans to the reported c.7180G>A variant, ER*B. AlphaFold2 protein modeling showed that each of the missense variants is predicted to encode an altered structural conformation near Era and Erb.

CONCLUSIONS

Investigation of archived samples resulted in the identification of three novel PIEZO1 alleles including a predicted Ernull and two missense variants. Structural modeling suggests that the missense changes potentially alter Era/Erb epitope expression with p.2399Ser resulting in a small increase in the negative electrostatic potential.

A pseudoautosomal glycosylation disorder prompts the revision of dolichol biosynthesis

Dolichol is a lipid critical for N-glycosylation as a carrier for activated sugars and nascent oligosaccharides. It is commonly thought to be directly produced from polyprenol by the enzyme SRD5A3. Instead, we found that dolichol synthesis requires a three-step detour involving additional metabolites, where SRD5A3 catalyzes only the second reaction. The first and third steps are performed by DHRSX, whose gene resides on the pseudoautosomal regions of the X and Y chromosomes. Accordingly, we report a pseudoautosomal-recessive disease presenting as a congenital disorder of glycosylation in patients with missense variants in DHRSX (DHRSX-CDG). Of note, DHRSX has a unique dual substrate and cofactor specificity, allowing it to act as a NAD+-dependent dehydrogenase and as a NADPH-dependent reductase in two non-consecutive steps. Thus, our work reveals unexpected complexity in the terminal steps of dolichol biosynthesis. Furthermore, we provide insights into the mechanism by which dolichol metabolism defects contribute to disease.

ProtVar: mapping and contextualizing human missense variation

Genomic variation can impact normal biological function in complex ways and so understanding variant effects requires a broad range of data to be coherently assimilated. Whilst the volume of human variant data and relevant annotations has increased, the corresponding increase in the breadth of participating fields, standards and versioning mean that moving between genomic, coding, protein and structure positions is increasingly complex. In turn this makes investigating variants in diverse formats and assimilating annotations from different resources challenging. ProtVar addresses these issues to facilitate the contextualization and interpretation of human missense variation with unparalleled flexibility and ease of accessibility for use by the broadest range of researchers. By precalculating all possible variants in the human proteome it offers near instantaneous mapping between all relevant data types. It also combines data and analyses from a plethora of resources to bring together genomic, protein sequence and function annotations as well as structural insights and predictions to better understand the likely effect of missense variation in humans. It is offered as an intuitive web server https://www.ebi.ac.uk/protvar where data can be explored and downloaded, and can be accessed programmatically via an API.

Improved protein complex prediction with AlphaFold-multimer by denoising the MSA profile

Structure prediction of protein complexes has improved significantly with AlphaFold2 and AlphaFold-multimer (AFM), but only 60% of dimers are accurately predicted. Here, we learn a bias to the MSA representation that improves the predictions by performing gradient descent through the AFM network. We demonstrate the performance on seven difficult targets from CASP15 and increase the average MMscore to 0.76 compared to 0.63 with AFM. We evaluate the procedure on 487 protein complexes where AFM fails and obtain an increased success rate (MMscore>0.75) of 33% on these difficult targets. Our protocol, AFProfile, provides a way to direct predictions towards a defined target function guided by the MSA. We expect gradient descent over the MSA to be useful for different tasks.

Structure-aware deep learning model for peptide toxicity prediction

Antimicrobial resistance is a critical public health concern, necessitating the exploration of alternative treatments. While antimicrobial peptides (AMPs) show promise, assessing their toxicity using traditional wet lab methods is both time-consuming and costly. We introduce tAMPer, a novel multi-modal deep learning model designed to predict peptide toxicity by integrating the underlying amino acid sequence composition and the three-dimensional structure of peptides. tAMPer adopts a graph-based representation for peptides, encoding ColabFold-predicted structures, where nodes represent amino acids and edges represent spatial interactions. Structural features are extracted using graph neural networks, and recurrent neural networks capture sequential dependencies. tAMPer's performance was assessed on a publicly available protein toxicity benchmark and an AMP hemolysis data we generated. On the latter, tAMPer achieves an F1-score of 68.7%, outperforming the second-best method by 23.4%. On the protein benchmark, tAMPer exhibited an improvement of over 3.0% in the F1-score compared to current state-of-the-art methods. We anticipate tAMPer to accelerate AMP discovery and development by reducing the reliance on laborious toxicity screening experiments.

rAbDesFlow: a novel workflow for computational recombinant antibody design for healthcare engineering

Recombinant antibodies (rAbs) have emerged as a promising solution to tackle antigen specificity, enhancement of immunogenic potential and versatile functionalization to treat human diseases. The development of single chain variable fragments has helped accelerate treatment in cancers and viral infections, due to their favorable pharmacokinetics and human compatibility. However, designing rAbs is traditionally viewed as a genetic engineering problem, with phage display and cell free systems playing a major role in sequence selection for gene synthesis. The process of antibody engineering involves complex and time-consuming laboratory techniques, which demand substantial resources and expertise. The success rate of obtaining desired antibody candidates through experimental approaches can be modest, necessitating iterative cycles of selection and optimization. With ongoing advancements in technology, in silico design of diverse antibody libraries, screening and identification of potential candidates for in vitro validation can be accelerated. To meet this need, we have developed rAbDesFlow, a unified computational workflow for recombinant antibody engineering with open-source programs and tools for ease of implementation. The workflow encompasses five computational modules to perform antigen selection, antibody library generation, antigen and antibody structure modeling, antigen-antibody interaction modeling, structure analysis, and consensus ranking of potential antibody sequences for synthesis and experimental validation. The proposed workflow has been demonstrated through design of rAbs for the ovarian cancer antigen Mucin-16 (CA-125). This approach can serve as a blueprint for designing similar engineered molecules targeting other biomarkers, allowing for a simplified adaptation to different cancer types or disease-specific antigens.

Accelerating Genetic Sensor Development, Scale-up, and Deployment Using Synthetic Biology

Living cells are exquisitely tuned to sense and respond to changes in their environment. Repurposing these systems to create engineered biosensors has seen growing interest in the field of synthetic biology and provides a foundation for many innovative applications spanning environmental monitoring to improved biobased production. In this review, we present a detailed overview of currently available biosensors and the methods that have supported their development, scale-up, and deployment. We focus on genetic sensors in living cells whose outputs affect gene expression. We find that emerging high-throughput experimental assays and evolutionary approaches combined with advanced bioinformatics and machine learning are establishing pipelines to produce genetic sensors for virtually any small molecule, protein, or nucleic acid. However, more complex sensing tasks based on classifying compositions of many stimuli and the reliable deployment of these systems into real-world settings remain challenges. We suggest that recent advances in our ability to precisely modify nonmodel organisms and the integration of proven control engineering principles (e.g., feedback) into the broader design of genetic sensing systems will be necessary to overcome these hurdles and realize the immense potential of the field.

Predicting therapeutic and side effects from drug binding affinities to human proteome structures

Evaluation of the binding affinities of drugs to proteins is a crucial process for identifying drug pharmacological actions, but it requires three dimensional structures of proteins. Herein, we propose novel computational methods to predict the therapeutic indications and side effects of drug candidate compounds from the binding affinities to human protein structures on a proteome-wide scale. Large-scale docking simulations were performed for 7,582 drugs with 19,135 protein structures revealed by AlphaFold (including experimentally unresolved proteins), and machine learning models on the proteome-wide binding affinity score (PBAS) profiles were constructed. We demonstrated the usefulness of the method for predicting the therapeutic indications for 559 diseases and side effects for 285 toxicities. The method enabled to predict drug indications for which the related protein structures had not been experimentally determined and to successfully extract proteins eliciting the side effects. The proposed method will be useful in various applications in drug discovery.

AlphaFold2 structures guide prospective ligand discovery

AlphaFold2 (AF2) models have had wide impact but mixed success in retrospective ligand recognition. We prospectively docked large libraries against unrefined AF2 models of the σ2 and serotonin 2A (5-HT2A) receptors, testing hundreds of new molecules and comparing results with those obtained from docking against the experimental structures. Hit rates were high and similar for the experimental and AF2 structures, as were affinities. Success in docking against the AF2 models was achieved despite differences between orthosteric residue conformations in the AF2 models and the experimental structures. Determination of the cryo-electron microscopy structure for one of the more potent 5-HT2A ligands from the AF2 docking revealed residue accommodations that resembled the AF2 prediction. AF2 models may sample conformations that differ from experimental structures but remain low energy and relevant for ligand discovery, extending the domain of structure-based drug design.

The substrate quality of CK2 target sites has a determinant role on their function and evolution

Most biological processes are regulated by signaling modules that bind to short linear motifs. For protein kinases, substrates may have full or only partial matches to the kinase recognition motif, a property known as "substrate quality." However, it is not clear whether differences in substrate quality represent neutral variation or if they have functional consequences. We examine this question for the kinase CK2, which has many fundamental functions. We show that optimal CK2 sites are phosphorylated at maximal stoichiometries and found in many conditions, whereas minimal substrates are more weakly phosphorylated and have regulatory functions. Optimal CK2 sites tend to be more conserved, and substrate quality is often tuned by selection. For intermediate sites, increases or decreases in substrate quality may be deleterious, as we demonstrate for a CK2 substrate at the kinetochore. The results together suggest a strong role for substrate quality in phosphosite function and evolution. A record of this paper's transparent peer review process is included in the supplemental information.

MoLPC2: improved prediction of large protein complex structures and stoichiometry using Monte Carlo Tree Search and AlphaFold2

MOTIVATION

Today, the prediction of structures of large protein complexes solely from their sequence information requires prior knowledge of the stoichiometry of the complex. To address this challenge, we have enhanced the Monte Carlo Tree Search algorithms in MoLPC to enable the assembly of protein complexes while simultaneously predicting their stoichiometry.

RESULTS

In MoLPC2, we have improved the predictions by allowing sampling alternative AlphaFold predictions. Using MoLPC2, we accurately predicted the structures of 50 out of 175 nonredundant protein complexes (TM-score ≥ 0.8) without knowing the stoichiometry. MoLPC2 provides new opportunities for predicting protein complex structures without stoichiometry information.

AVAILABILITY AND IMPLEMENTATION

MoLPC2 is freely available at https://github.com/hychim/molpc2. A notebook is also available from the repository for easy use.

Overview of AlphaFold2 and breakthroughs in overcoming its limitations

Predicting three-dimensional (3D) protein structures has been challenging for decades. The emergence of AlphaFold2 (AF2), a deep learning-based machine learning method developed by DeepMind, became a game changer in the protein folding community. AF2 can predict a protein's three-dimensional structure with high confidence based on its amino acid sequence. Accurate prediction of protein structures can dramatically accelerate our understanding of biological mechanisms and provide a solid foundation for reliable drug design. Although AF2 breaks through the barriers in predicting protein structures, many rooms remain to be further studied. This review provides a brief historical overview of the development of protein structure prediction, covering template-based, template-free, and machine learning-based methods. In addition to reviewing the potential benefits (Pros) and considerations (Cons) of using AF2, this review summarizes the diverse applications, including protein structure predictions, dynamic changes, point mutation, integration of language model and experimental data, protein complex, and protein-peptide interaction. It underscores recent advancements in efficiency, reliability, and broad application of AF2. This comprehensive review offers valuable insights into the applications of AF2 and AF2-inspired AI methods in structural biology and its potential for clinically significant drug target discovery.

Random, de novo, and conserved proteins: How structure and disorder predictors perform differently

Understanding the emergence and structural characteristics of de novo and random proteins is crucial for unraveling protein evolution and designing novel enzymes. However, experimental determination of their structures remains challenging. Recent advancements in protein structure prediction, particularly with AlphaFold2 (AF2), have expanded our knowledge of protein structures, but their applicability to de novo and random proteins is unclear. In this study, we investigate the structural predictions and confidence scores of AF2 and protein language model-based predictor ESMFold for de novo and conserved proteins from Drosophila and a dataset of comparable random proteins. We find that the structural predictions for de novo and random proteins differ significantly from conserved proteins. Interestingly, a positive correlation between disorder and confidence scores (pLDDT) is observed for de novo and random proteins, in contrast to the negative correlation observed for conserved proteins. Furthermore, the performance of structure predictors for de novo and random proteins is hampered by the lack of sequence identity. We also observe fluctuating median predicted disorder among different sequence length quartiles for random proteins, suggesting an influence of sequence length on disorder predictions. In conclusion, while structure predictors provide initial insights into the structural composition of de novo and random proteins, their accuracy and applicability to such proteins remain limited. Experimental determination of their structures is necessary for a comprehensive understanding. The positive correlation between disorder and pLDDT could imply a potential for conditional folding and transient binding interactions of de novo and random proteins.

The synergy between the PscC subunits for electron transfer to the P840 special pair in Chlorobaculum tepidum

The primary photochemical reaction of photosynthesis in green sulfur bacteria occurs in the homodimer PscA core proteins by a special chlorophyll pair. The light induced excited state of the special pair producing P840+ is rapidly reduced by electron transfer from one of the two PscC subunits. Molecular dynamics (MD) simulations are combined with bioinformatic tools herein to provide structural and dynamic insight into the complex between the two PscA core proteins and the two PscC subunits. The microscopic dynamic model involves extensive sampling at atomic resolution and at a cumulative time-scale of 22µs and reveals well defined protein-protein interactions. The membrane complex is composed of the two PscA and the two PscC subunits and macroscopic connections are revealed within a putative electron transfer pathway from the PscC subunit to the special pair P840 located within the PscA subunits. Our results provide a structural basis for understanding the electron transport to the homodimer RC of the green sulfur bacteria. The MD based approach can provide the basis to further probe the PscA-PscC complex dynamics and observe electron transfer therein at the quantum level. Furthermore, the transmembrane helices of the different PscC subunits exert distinct dynamics in the complex.

Decoding the molecular interplay in the central dogma: An overview of mass spectrometry-based methods to investigate protein-metabolite interactions

With the emergence of next-generation nucleotide sequencing and mass spectrometry-based proteomics and metabolomics tools, we have comprehensive and scalable methods to analyze the genes, transcripts, proteins, and metabolites of a multitude of biological systems. Despite the fascinating new molecular insights at the genome, transcriptome, proteome and metabolome scale, we are still far from fully understanding cellular organization, cell cycles and biology at the molecular level. Significant advances in sensitivity and depth for both sequencing as well as mass spectrometry-based methods allow the analysis at the single cell and single molecule level. At the same time, new tools are emerging that enable the investigation of molecular interactions throughout the central dogma of molecular biology. In this review, we provide an overview of established and recently developed mass spectrometry-based tools to probe metabolite-protein interactions-from individual interaction pairs to interactions at the proteome-metabolome scale.

Imbalance in Unc80 RNA Editing Disrupts Dynamic Neuronal Activity and Olfactory Perception

A-to-I RNA editing, catalyzed by the ADAR protein family, significantly contributes to the diversity and adaptability of mammalian RNA signatures, aligning with developmental and physiological needs. Yet, the functions of many editing sites are still to be defined. The Unc80 gene stands out in this context due to its brain-specific expression and the evolutionary conservation of its codon-altering editing event. The precise biological functions of Unc80 and its editing, however, are still largely undefined. In this study, we first demonstrated that Unc80 editing occurs in an ADAR2-dependent manner and is exclusive to the brain. By employing the CRISPR/Cas9 system to generate Unc80 knock-in mouse models that replicate the natural editing variations, our findings revealed that mice with the "gain-of-editing" variant (Unc80G/G) exhibit heightened basal neuronal activity in critical olfactory regions, compared to the "loss-of-editing" (Unc80S/S) counterparts. Moreover, an increase in glutamate levels was observed in the olfactory bulbs of Unc80G/G mice, indicating altered neurotransmitter dynamics. Behavioral analysis of odor detection revealed distinctive responses to novel odors-both Unc80 deficient (Unc80+/-) and Unc80S/S mice demonstrated prolonged exploration times and heightened dishabituation responses. Further elucidating the olfactory connection of Unc80 editing, transcriptomic analysis of the olfactory bulb identified significant alterations in gene expression that corroborate the behavioral and physiological findings. Collectively, our research advances the understanding of Unc80's neurophysiological functions and the impact of its editing on the olfactory sensory system, shedding light on the intricate molecular underpinnings of olfactory perception and neuronal activity.

GTExome: Modeling commonly expressed missense mutations in the human genome

A web application, GTExome, is described that quickly identifies, classifies, and models missense mutations in commonly expressed human proteins. GTExome can be used to categorize genomic mutation data with tissue specific expression data from the Genotype-Tissue Expression (GTEx) project. Commonly expressed missense mutations in proteins from a wide range of tissue types can be selected and assessed for modeling suitability. Information about the consequences of each mutation is provided to the user including if disulfide bonds, hydrogen bonds, or salt bridges are broken, buried prolines introduced, buried charges are created or lost, charge is swapped, a buried glycine is replaced, or if the residue that would be removed is a proline in the cis configuration. Also, if the mutation site is in a binding pocket the number of pockets and their volumes are reported. The user can assess this information and then select from available experimental or computationally predicted structures of native proteins to create, visualize, and download a model of the mutated protein using Fast and Accurate Side-chain Protein Repacking (FASPR). For AlphaFold modeled proteins, confidence scores for native proteins are provided. Using this tool, we explored a set of 9,666 common missense mutations from a variety of tissues from GTEx and show that most mutations can be modeled using this tool to facilitate studies of protein-protein and protein-drug interactions. The open-source tool is freely available at https://pharmacogenomics.clas.ucdenver.edu/gtexome/.

Mapping adipocyte interactome networks by HaloTag-enrichment-mass spectrometry

Mapping protein interaction complexes in their natural state in vivo is arguably the Holy Grail of protein network analysis. Detection of protein interaction stoichiometry has been an important technical challenge, as few studies have focused on this. This may, however, be solved by artificial intelligence (AI) and proteomics. Here, we describe the development of HaloTag-based affinity purification mass spectrometry (HaloMS), a high-throughput HaloMS assay for protein interaction discovery. The approach enables the rapid capture of newly expressed proteins, eliminating tedious conventional one-by-one assays. As a proof-of-principle, we used HaloMS to evaluate the protein complex interactions of 17 regulatory proteins in human adipocytes. The adipocyte interactome network was validated using an in vitro pull-down assay and AI-based prediction tools. Applying HaloMS to probe adipocyte differentiation facilitated the identification of previously unknown transcription factor (TF)-protein complexes, revealing proteome-wide human adipocyte TF networks and shedding light on how different pathways are integrated.

Predicted mechanistic impacts of human protein missense variants

Genome sequencing efforts have led to the discovery of tens of millions of protein missense variants found in the human population with the majority of these having no annotated role and some likely contributing to trait variation and disease. Sequence-based artificial intelligence approaches have become highly accurate at predicting variants that are detrimental to the function of proteins but they do not inform on mechanisms of disruption. Here we combined sequence and structure-based methods to perform proteome-wide prediction of deleterious variants with information on their impact on protein stability, protein-protein interactions and small-molecule binding pockets. AlphaFold2 structures were used to predict approximately 100,000 small-molecule binding pockets and stability changes for over 200 million variants. To inform on protein-protein interfaces we used AlphaFold2 to predict structures for nearly 500,000 protein complexes. We illustrate the value of mechanism-aware variant effect predictions to study the relation between protein stability and abundance and the structural properties of interfaces underlying trans protein quantitative trait loci (pQTLs). We characterised the distribution of mechanistic impacts of protein variants found in patients and experimentally studied example disease linked variants in FGFR1.

Expression landscape of cancer-FOXP3 and its prognostic value in pancreatic adenocarcinoma

FOXP3, a key identifier of Treg, has also been identified in tumor cells, which is referred to as cancer-FOXP3 (c-FOXP3). Human c-FOXP3 undergoes multiple alternative splicing events, generating several isoforms, like c-FOXP3FL and c-FOXP3Δ3. Previous research on c-FOXP3 often ignore its cellular source (immune or tumor cells) and isoform expression patterns, which may obscure our understanding of its clinical significance. Our immunohistochemistry investigations which conducted across 18 tumors using validated c-FOXP3 antibodies revealed distinct expression landscapes for c-FOXP3 and its variants, with the majority of tumors exhibited a predominantly expression of c-FOXP3Δ3. In pancreatic ductal adenocarcinoma (PDAC), we further discovered a potential link between nuclear c-FOXP3Δ3 in tumor cells and poor prognosis. Overexpression of c-FOXP3Δ3 in tumor cells was associated with metastasis. This work elucidates the expression pattern of c-FOXP3 in pan-cancer and indicates its potential as a prognostic biomarker in clinical settings, offering new perspectives for its clinical application.

Proton-coupled transport mechanism of the efflux pump NorA

Efflux pump antiporters confer drug resistance to bacteria by coupling proton import with the expulsion of antibiotics from the cytoplasm. Despite efforts there remains a lack of understanding as to how acid/base chemistry drives drug efflux. Here, we uncover the proton-coupling mechanism of the Staphylococcus aureus efflux pump NorA by elucidating structures in various protonation states of two essential acidic residues using cryo-EM. Protonation of Glu222 and Asp307 within the C-terminal domain stabilized the inward-occluded conformation by forming hydrogen bonds between the acidic residues and a single helix within the N-terminal domain responsible for occluding the substrate binding pocket. Remarkably, deprotonation of both Glu222 and Asp307 is needed to release interdomain tethering interactions, leading to opening of the pocket for antibiotic entry. Hence, the two acidic residues serve as a "belt and suspenders" protection mechanism to prevent simultaneous binding of protons and drug that enforce NorA coupling stoichiometry and confer antibiotic resistance.

Discovering predisposing genes for hereditary breast cancer using deep learning

Breast cancer (BC) is the most common malignancy affecting Western women today. It is estimated that as many as 10% of BC cases can be attributed to germline variants. However, the genetic basis of the majority of familial BC cases has yet to be identified. Discovering predisposing genes contributing to familial BC is challenging due to their presumed rarity, low penetrance, and complex biological mechanisms. Here, we focused on an analysis of rare missense variants in a cohort of 12 families of Middle Eastern origins characterized by a high incidence of BC cases. We devised a novel, high-throughput, variant analysis pipeline adapted for family studies, which aims to analyze variants at the protein level by employing state-of-the-art machine learning models and three-dimensional protein structural analysis. Using our pipeline, we analyzed 1218 rare missense variants that are shared between affected family members and classified 80 genes as candidate pathogenic. Among these genes, we found significant functional enrichment in peroxisomal and mitochondrial biological pathways which segregated across seven families in the study and covered diverse ethnic groups. We present multiple evidence that peroxisomal and mitochondrial pathways play an important, yet underappreciated, role in both germline BC predisposition and BC survival.