Search

The goal of this paper is predicting the conformational distributions of ligand binding sites using the AlphaFold2 (AF2) protein structure prediction program with stochastic subsampling of the multiple sequence alignment (MSA). We explored the opening of cryptic ligand binding sites in 16 proteins, where the closed and open conformations define the expected extreme points of the conformational variation. Due to the many structures of these proteins in the Protein Data Bank (PDB), we were able to study whether the distribution of X-ray structures affects the distribution of AF2 models. We have found that AF2 generates both a cluster of open and a cluster of closed models for proteins that have comparable numbers of open and closed structures in the PDB and not too many other conformations. This was observed even with default MSA parameters, thus without further subsampling. In contrast, with the exception of a single protein, AF2 did not yield multiple clusters of conformations for proteins that had imbalanced numbers of open and closed structures in the PDB, or had substantial numbers of other structures. Subsampling improved the results only for a single protein, but very shallow MSA led to incorrect structures. The ability of generating both open and closed conformations for six out of the 16 proteins agrees with the success rates of similar studies reported in the literature. However, we showed that this partial success is due to AF2 "remembering" the conformational distributions in the PDB and that the approach fails to predict rarely seen conformations., Competing Interests: Competing interests statement:The authors declare no competing interest.

Cryptic sites can expand the space of druggable proteins, but the potential usefulness of such sites needs to be investigated before any major effort. Given that the binding pockets are not formed, the druggability of such sites is not well understood. The analysis of proteins and their ligands shows that cryptic sites that are formed primarily by the motion of side chains moving out of the pocket to enable ligand binding generally do not bind drug-sized molecules with sufficient potency. By contrast, sites that are formed by loop or hinge motion are potentially valuable drug targets. Arguments are provided to explain the underlying causes in terms of classical enzyme inhibition theory and the kinetics of side chain motion and ligand binding., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Target identification and evaluation is a critical step in the drug discovery process. Although time-intensive and complex, the challenge becomes even more acute in the realm of infectious disease, where the rapid emergence of new viruses, the swift mutation of existing targets, and partial effectiveness of approved antivirals can lead to outbreaks of significant public health concern. The COVID-19 pandemic, caused by the SARS-CoV-2 virus, serves as a prime example of this, where despite the allocation of substantial resources, Paxlovid is currently the only effective treatment. In that case, significant effort pre-pandemic had been expended to evaluate the biological target for the closely related SARS-CoV. In this work, we utilize the computational hot spot mapping method, FTMove, to rapidly identify and rank binding sites for a set of nine SARS-CoV-2 drug/potential drug targets. FTMove takes into account protein flexibility by mapping binding site hot spots across an ensemble of structures for a given target. To assess the applicability of the FTMove approach to a wide range of drug targets for viral pathogens, we also carry out a comprehensive review of the known SARS-CoV-2 ligandable sites. The approach is able to identify the vast majority of all known sites and a few additional sites, which may in fact be yet to be discovered as ligandable. Furthermore, a UMAP analysis of the FTMove features for each identified binding site is largely able to separate predicted sites with experimentally known binders from those without known binders. These results demonstrate the utility of FTMove to rapidly identify actionable sites across a range of targets for a given indication. As such, the approach is expected to be particularly useful for assessing target binding sites for any emerging pathogen, as well as for indications in other disease areas, and providing actionable starting points for structure-based drug design efforts.

Objective: The rapid expansion of the biomedical literature challenges traditional review methods, especially during outbreaks of emerging infectious diseases when quick action is critical. Our study aims to explore the potential of ChatGPT to automate the biomedical literature review for rapid drug discovery., Materials and Methods: We introduce a novel automated pipeline helping to identify drugs for a given virus in response to a potential future global health threat. Our approach can be used to select PubMed articles identifying a drug target for the given virus. We tested our approach on two known pathogens: SARS-CoV-2, where the literature is vast, and Nipah, where the literature is sparse. Specifically, a panel of three experts reviewed a set of PubMed articles and labeled them as either describing a drug target for the given virus or not. The same task was given to the automated pipeline and its performance was based on whether it labeled the articles similarly to the human experts. We applied a number of prompt engineering techniques to improve the performance of ChatGPT., Results: Our best configuration used GPT-4 by OpenAI and achieved an out-of-sample validation performance with accuracy/F1-score/sensitivity/specificity of 92.87%/88.43%/83.38%/97.82% for SARS-CoV-2 and 87.40%/73.90%/74.72%/91.36% for Nipah., Conclusion: These results highlight the utility of ChatGPT in drug discovery and development and reveal their potential to enable rapid drug target identification during a pandemic-level health emergency., Competing Interests: Declaration of Competing Interest The authors declare no competing interests., (Copyright © 2024 Elsevier B.V. All rights reserved.)

The knowledge of ligand binding hot spots and of the important interactions within such hot spots is crucial for the design of lead compounds in the early stages of structure-based drug discovery. The computational solvent mapping server FTMap can reliably identify binding hot spots as consensus clusters, free energy minima that bind a variety of organic probe molecules. However, in its current implementation, FTMap provides limited information on regions within the hot spots that tend to interact with specific pharmacophoric features of potential ligands. E-FTMap is a new server that expands on the original FTMap protocol. E-FTMap uses 119 organic probes, rather than the 16 in the original FTMap, to exhaustively map binding sites, and identifies pharmacophore features as atomic consensus sites where similar chemical groups bind. We validate E-FTMap against a set of 109 experimentally derived structures of fragment-lead pairs, finding that highly ranked pharmacophore features overlap with the corresponding atoms in both fragments and lead compounds. Additionally, comparisons of mapping results to ensembles of bound ligands reveal that pharmacophores generated with E-FTMap tend to sample highly conserved protein-ligand interactions. E-FTMap is available as a web server at https://eftmap.bu.edu.

Introduction: A fundamental challenge in computational vaccinology is that most B-cell epitopes are conformational and therefore hard to predict from sequence alone. Another significant challenge is that a great deal of the amino acid sequence of a viral surface protein might not in fact be antigenic. Thus, identifying the regions of a protein that are most promising for vaccine design based on the degree of surface exposure may not lead to a clinically relevant immune response. Methods: Linear peptides selected by phage display experiments that have high affinity to the monoclonal antibody of interest ("mimotopes") usually have similar physicochemical properties to the antigen epitope corresponding to that antibody. The sequences of these linear peptides can be used to find possible epitopes on the surface of the antigen structure or a homology model of the antigen in the absence of an antigen-antibody complex structure. Results and Discussion: Herein we describe two novel methods for mapping mimotopes to epitopes. The first is a novel algorithm named MimoTree that allows for gaps in the mimotopes and epitopes on the antigen. More specifically, a mimotope may have a gap that does not match to the epitope to allow it to adopt a conformation relevant for binding to an antibody, and residues may similarly be discontinuous in conformational epitopes. MimoTree is a fully automated epitope detection algorithm suitable for the identification of conformational as well as linear epitopes. The second is an ensemble approach, which combines the prediction results from MimoTree and two existing methods., Competing Interests: Authors MS and AC were employed by Fractal Therapeutics Inc. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Li, Wilderotter, Stoddard, Van Egeren, Chakravarty and Joseph-McCarthy.)

The neural network-based program AlphaFold2 (AF2) provides high accuracy structure prediction for a large fraction of globular proteins. An important question is whether these models are accurate enough for reliably docking small ligands. Several recent papers and the results of CASP15 reveal that local conformational errors reduce the success rates of direct ligand docking. Here, we focus on the ability of the models to conserve the location of binding hot spots, regions on the protein surface that significantly contribute to the binding free energy of the protein-ligand interaction. Clusters of hot spots predict the location and even the druggability of binding sites, and hence are important for computational drug discovery. The hot spots are determined by protein mapping that is based on the distribution of small fragment-sized probes on the protein surface and is less sensitive to local conformation than docking. Mapping models taken from the AlphaFold Protein Structure Database show that identifying binding sites is more reliable than docking, but the success rates are still 5% to 10% lower than based on mapping X-ray structures. The drop in accuracy is particularly large for models of multidomain proteins. However, both the model binding sites and the mapping results can be substantially improved by generating AF2 models for the ligand binding domains of interest rather than the entire proteins and even more if using forced sampling with multiple initial seeds. The mapping of such models tends to reach the accuracy of results obtained by mapping the X-ray structures.

SARS-CoV-2, the virus that causes COVID-19, led to a global health emergency that claimed the lives of millions. Despite the widespread availability of vaccines, the virus continues to exist in the population in an endemic state which allows for the continued emergence of new variants. Most of the current vaccines target the spike glycoprotein interface of SARS-CoV-2, creating a selection pressure favoring viral immune evasion. Antivirals targeting other molecular interactions of SARS-CoV-2 can help slow viral evolution by providing orthogonal selection pressures on the virus. GRP78 is a host auxiliary factor that mediates binding of the SARS-CoV-2 spike protein to human cellular ACE2, the primary pathway of cell infection. As GRP78 forms a ternary complex with SARS-CoV-2 spike protein and ACE2, disrupting the formation of this complex is expected to hinder viral entry into host cells. Here, we developed a model of the GRP78-Spike RBD-ACE2 complex. We then used that model together with hot spot mapping of the GRP78 structure to identify the putative binding site for spike protein on GRP78. Next, we performed structure-based virtual screening of known drug/candidate drug libraries to identify binders to GRP78 that are expected to disrupt spike protein binding to the GRP78, and thereby preventing viral entry to the host cell. A subset of these compounds has previously been shown to have some activity against SARS-CoV-2. The identified hits are starting points for the further development of novel SARS-CoV-2 therapeutics, potentially serving as proof-of-concept for GRP78 as a potential drug target for other viruses., (© 2024. The Author(s), under exclusive licence to Springer Nature Switzerland AG.)

The rapid emergence of immune-evading viral variants of SARS-CoV-2 calls into question the practicality of a vaccine-only public-health strategy for managing the ongoing COVID-19 pandemic. It has been suggested that widespread vaccination is necessary to prevent the emergence of future immune-evading mutants. Here, we examined that proposition using stochastic computational models of viral transmission and mutation. Specifically, we looked at the likelihood of emergence of immune escape variants requiring multiple mutations and the impact of vaccination on this process. Our results suggest that the transmission rate of intermediate SARS-CoV-2 mutants will impact the rate at which novel immune-evading variants appear. While vaccination can lower the rate at which new variants appear, other interventions that reduce transmission can also have the same effect. Crucially, relying solely on widespread and repeated vaccination (vaccinating the entire population multiple times a year) is not sufficient to prevent the emergence of novel immune-evading strains, if transmission rates remain high within the population. Thus, vaccines alone are incapable of slowing the pace of evolution of immune evasion, and vaccinal protection against severe and fatal outcomes for COVID-19 patients is therefore not assured.

At the outset of an emergent viral respiratory pandemic, sequence data is among the first molecular information available. As viral attachment machinery is a key target for therapeutic and prophylactic interventions, rapid identification of viral "spike" proteins from sequence can significantly accelerate the development of medical countermeasures. For six families of respiratory viruses, covering the vast majority of airborne and droplet-transmitted diseases, host cell entry is mediated by the binding of viral surface glycoproteins that interact with a host cell receptor. In this report it is shown that sequence data for an unknown virus belonging to one of the six families above provides sufficient information to identify the protein(s) responsible for viral attachment. Random forest models that take as input a set of respiratory viral sequences can classify the protein as "spike" vs. non-spike based on predicted secondary structure elements alone (with 97.3% correctly classified) or in combination with N-glycosylation related features (with 97.0% correctly classified). Models were validated through 10-fold cross-validation, bootstrapping on a class-balanced set, and an out-of-sample extra-familial validation set. Surprisingly, we showed that secondary structural elements and N-glycosylation features were sufficient for model generation. The ability to rapidly identify viral attachment machinery directly from sequence data holds the potential to accelerate the design of medical countermeasures for future pandemics. Furthermore, this approach may be extendable for the identification of other potential viral targets and for viral sequence annotation in general in the future., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Walker et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

In the face of a long-running pandemic, understanding the drivers of ongoing SARS-CoV-2 transmission is crucial for the rational management of COVID-19 disease burden. Keeping schools open has emerged as a vital societal imperative during the pandemic, but in-school transmission of SARS-CoV-2 can contribute to further prolonging the pandemic. In this context, the role of schools in driving SARS-CoV-2 transmission acquires critical importance. Here we model in-school transmission from first principles to investigate the effectiveness of layered mitigation strategies on limiting in-school spread. We examined the effect of masks and air quality (ventilation, filtration and ionizers) on steady-state viral load in classrooms, as well as on the number of particles inhaled by an uninfected person. The effectiveness of these measures in limiting viral transmission was assessed for variants with different levels of mean viral load (ancestral, Delta, Omicron). Our results suggest that a layered mitigation strategy can be used effectively to limit in-school transmission, with certain limitations. First, poorly designed strategies (insufficient ventilation, no masks, staying open under high levels of community transmission) will permit in-school spread even if some level of mitigation is present. Second, for viral variants that are sufficiently contagious, it may be difficult to construct any set of interventions capable of blocking transmission once an infected individual is present, underscoring the importance of other measures. Our findings provide practical recommendations; in particular, the use of a layered mitigation strategy that is designed to limit transmission, with other measures such as frequent surveillance testing and smaller class sizes (such as by offering remote schooling options to those who prefer it) as needed., Competing Interests: Authors AC and MS are employed by Fractal Therapeutics. Author RN is employed by Halozyme Therapeutics. Authors AC, MS, RN, DV, and DJ-M are shareholders in Fractal Therapeutics. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Van Egeren, Stoddard, Malakar, Ghosh, Acharya, Mainuddin, Majumdar, Luo, Nolan, Joseph-McCarthy, White, Hochberg, Basu and Chakravarty.)