Twelfth Novalix Biophysics in Drug Discovery Conference

Last week the Twelfth Novalix Biophysics in Drug Discovery Conference was held for the first time in La Jolla, California. It’s been several years since I wrote about one of these, and I was happy to see that they’ve maintained their reputation for excellent science and convivial conversation. There’s no way to cover the two-dozen talks, but here are a few highlights.

 

One of the things I most enjoy about these meetings is learning about new and emerging methods, and these were well represented. Chris Brosey (AbbVie) discussed time-resolved high-throughput small-angle X-ray scattering (TR-HT-SAXS). As I discussed a couple years ago, the approach can be used to measure the kinetics of protein dimerization in response to fragment-sized ligands.

 

SAXS-based approaches typically require access to a synchrotron, but Takashi Sato (Rigaku) described a related approach, electron density tomography (EDT), using an in-house instrument. Using machine learning, EDT can provide more detailed structural information than standard SAXS, and Takashi provided examples for samples ranging in size from viruses to single-chain variable fragments (scFvs) smaller than 30 kD.

 

Another approach to examine protein complexes is microfluidic diffusional sizing (MDS), described by James Wilkinson of Fluidic Sciences. By assessing the amount of diffusion in disposable chip-based chambers, MDS can determine how the hydrodynamic radius changes in response to ligands. Each chip holds 24 samples, and data can be collected in less than an hour. The minimum observable size change is 5-10% so measuring small molecules directly is unlikely. Still, the technique is useful for observing induced proximity events such as those caused by molecular glues, and it is sufficiently robust that it can be run in pure serum.

 

Among solution-based methods, none have achieved such recent prominence as cryo-EM. Weiru Wang, my colleague at Frontier Medicines, described how this technique was used iteratively to design and characterize bivalent degrader molecules that covalently exploit the E3 ligase DCAF2. (We recently published this work in Structure.)

 

Cryo-EM has revolutionized the types of biological molecules that can be structurally characterized. According to Denis Zeyer (Novalix), the technique accounted for 40% of PDB entries last year. However, despite the “resolution revolution,” most structures are not as detailed as those from X-ray crystallography; Denis noted that only 20% of the new structures were solved to a resolution better than 3 Å, which may have negative implications for machine-learning methods trained on these lower resolution structures.

 

If cryo-EM is the new kid in town, NMR is the grizzled veteran. But proving that it is possible to find new applications for old methods, Matthew Eddy (University of Florida) described a clever ¹⁹F-labeling approach for GPCRs in nanodiscs to quantify the distribution of various states in response to anionic lipids and ligands. This has allowed him to distinguish between antagonists and inverse agonists, which can be difficult using cell-based assays.

 

Turning from solution-based to surface-based methods, SPR has moved into the number two slot for fragment-finding, as we noted in our recent poll. Just as there are new tricks for NMR, the same applies to SPR. Matthew Peterfreund (Bruker Biosensors) described switchSENSE, a fluorescence proximity assay built on an SPR chip that is useful for measuring the binding and kinetics of bifunctional ligands such as PROTACs to two or more proteins. He also introduced the Triceratops SPR#64 instrument, which as its name implies supports 64 sensor spots.

 

Kris Borzilleri (Pfizer) discussed SPR-microscopy (SPRm), which combines an optical microscope with an SPR instrument. This can be used to measure the affinities of ligands binding to receptors in cells grown on SPR chips, and Kris described applications to membrane proteins such as GPCRs and solute carriers. The technique is still quite slow though, at only 10-15 compounds in duplicate per week.

 

Another surface-based approach to screening cells was described by Volker Gatterdam of Lino Biotech. Focal molography relies on changes in diffraction from nanoengineered diffraction gratings, called molograms. Targets, which can include living cells, are immobilized to the molograms, and analyte is flowed over. The instrument contains 64 spots, and assays can be run in complex samples such as tissue lysates.

 

Covalent drug discovery also made an appearance, with talks by Landon Whitby (Lundbeck) and Ben Cravatt (Scripps). Landon provided an overview of chemoproteomics techniques to screen ligands in cells, such as those we wrote about in 2016. Ben continued the theme, including several success stories, and also discussed challenges for finding cryptic ligandable pockets. Despite impressive progress with machine learning, Ben noted that these methods often find only common solutions, while empirical chemoproteomics methods can find rare types of pockets.

 

Of course, as we’ve repeatedly emphasized, biophysics methods are best used in combination, as noted by Daniel Harki (University of Minnesota) and Ann Boriack-Sjodin (Takeda). Daniel presented a screen of 1056 fragments against the cancer target APOBEC3 using NMR and SPR. This yielded just a single validated hit, which interestingly turned out to be the same fragment found against KRAS in a paper we discussed in 2022. And Ann described how biophysics led to multiple clinical compounds against a variety of targets at Epizyme and Accent Therapeutics.

 

I’ll stop here, but please feel free to add your thoughts. And while the date for the next Novalix conference has not yet been scheduled, the location has, with a return to beautiful Strasbourg. Vive la biophysique!

Mapping protein conformations with fragments

Proteins can be remarkably dynamic, and, as we noted recently, different conformational states can reveal different pockets for small molecule ligands. But how can one survey and categorize all the possibilities? In a recent J. Chem. Inf. Model. paper, Doeke Hekstra and colleagues at Harvard University present a new tool for doing so.

 

High-throughput crystallographic fragment screens are becoming faster and more widely accessible, and the researchers wondered whether the information from these screens could be used to map protein conformational landscapes. To do so, they built a Python program called COLAV, short for COnformational LAndscape Visualization. This open-source tool can compile data from hundreds of protein coordinate files and then, for each protein, calculate the dihedral angles between backbone atoms, the pairwise distances between the alpha-carbon atoms, and the strain.

 

To a first approximation, dihedral angles capture local movements, while distances between alpha-carbons capture global movements, such as the distance between the N-terminus and C-terminus. Strain measurements are also local but can reveal particularly important features such as hinge movements. Also, while dihedral and pairwise distances can be calculated for single proteins, strain measurements are calculated after first aligning multiple structures.

 

Having calculated these three parameters for individual protein structures, COLAV can compare them across the selected set of structures using principal component analysis (PCA). These comparisons can reveal clusters with similar dihedral angles, pairwise distances, or strain.

 

The researchers provide two case studies. The first is the metabolic disease target PTP1B, which we recently wrote about here. This enzyme has been pursued intensively for decades, so the researchers were able to draw on 163 individual protein structures deposited in the protein data bank (PDB) as well as 187 structures from a high-throughput crystallographic fragment screen. PTP1B contains two flexible loops, each of which adopts one of two conformations, and COLAV successfully segregated all 350 structures into four clusters. Importantly, these four clusters were found whether the structures were pulled from the PDB (representing experiments conducted across multiple labs and years) or from the fragment screen, suggesting that a single crystallographic fragment screen can identify most or all of the conformational states available to a protein. This is particularly impressive given that most of the fragments bound in allosteric sites while most of the ligands found in the PDB bound in the active site.

 

Next, the researchers turned to the main protease (MPro) of SARS-CoV-2, the subject of intense and successful drug discovery efforts. They used 656 structures from the PDB and 631 structures from high-throughput crystallographic screens to perform COLAV analyses. Unlike PTP1B, discrete conformational clusters were not observed; rather a continuous band was seen, suggesting that the protein can assume myriad conformations. Here too though, the fragment screens were able to sample most of the conformations observed in the PDB.

 

The fact that a single high-throughput crystallographic screen can capture the conformations seen in hundreds of hard-won discrete protein-ligand crystal structures is encouraging, though of course the paper only describes two case studies. Also, as the researchers note, any structure that cannot be crystallized is not sampled. Since COLAV is free to use, it will be fun to see it applied to other proteins.

SAR by TR-HT-SAXS

Well that’s an acronym soup! SAR by NMR was the first practical fragment-finding method, and over the years Practical Fragments has covered lots of other techniques. Small-angle X-ray scattering, or SAXS, has not been among them. As the name suggests, this technique uses X-rays, typically produced at a synchrotron. However, unlike conventional crystallography, it doesn’t require crystalline material. Instead, proteins in solution are analyzed to provide information on their size and shape. The resolution is too low to assess small molecule binding, but suitable for observing dimerization or changes in conformation.

 

Time-resolved SAXS, or TR-SAXS, examines SAXS over time in response to a trigger. For example, you can rapidly add a ligand to a protein and watch for changes in conformation. And HT simply means high throughput. A recent Nature Chemical Biology paper from Chris Brosey, John Tainer, and collaborators at the University of Texas MD Anderson Center, Lawrence Berkeley National Laboratory, University of California Santa Cruz, and University of Arkansas for Medical Sciences Little Rock describes structure-activity relationships by time-resolved high throughput small-angle X-ray scattering (TR-HT-SAXS).

 

The researchers were interested in apoptosis-inducing factor (AIF), a mitochondrial protein with potential implications for cancer and other diseases. AIF normally exists as a monomer in complex with an FAD cofactor. Binding of NADH causes reduction of FAD to FADH^- and concomitant dimerization of the protein. Could fragments do the same, allowing dimerization on demand?

 

A library of 2500 fragments purchased from Life Chemicals was screened at 0.75-1.5 mM against the AIF-FAD complex using differential scanning fluorimetry (DSF), and those that raised or lowered the temperature by more than 1.7 ºC were further characterized by microscale thermophoresis (MST). This led to 32 binders and 7 negative controls, or molecules that did not confirm either by DSF or MST. (Side note: although many people discount compounds that give negative thermal shifts, the natural ligand NADH lowers the melting temperature of AIF by a whopping 10.8 ºC.)

 

Next, the fragment binders and negative controls were screened at 0.5-1 mM by TR-SAXS. Intense X-rays cause reduction of the FAD cofactor, but in the absence of NADH or other ligands the AIF protein remains monomeric. However, some fragments did cause dimerization of the protein during TR-SAXS. Interestingly, these fragments were structurally related to one another. Subsequent crystallography revealed that they bind where NADH normally binds and make some of the same interactions to induce protein dimerization. The paper includes much more detailed characterization, including mutagenesis, spectroscopic, and protein crosslinking experiments to further understand the mechanism.

 

TR-SAXS is an interesting addition to our toolbox of biophysical methods suitable for fragment screening. It does have some disadvantages, such as the need for large amounts of protein at high concentrations: 67 µM in this case. Also, the “HT” may be somewhat aspirational, with a current throughput of 100-200 compounds per synchrotron shift. Finally, the technique is probably best suited to well-characterized proteins where SAXS data can be carefully modeled. With these limitations in mind, it will be fun to see how generally TR-SAXS finds fragments that alter the conformation and multimerization of proteins.

Fifth NovAliX Biophysics in Drug Discovery Conference

Last week NovAliX held its biophysics meeting outside of Strasbourg for the first time. Naturally they chose Boston, one of the most European of US cities and a major hub of drug discovery. The event brought together 118 participants from 15 countries, roughly 80% from industry. Although the food and drink could not compare to France, the science and discussion were every bit as satisfying. With 30 talks and 22 posters I won’t attempt to be comprehensive, but as with last year just try to capture a few themes. 

One particularly noteworthy session was devoted to single particle cryo-electron microscopy (cryo-EM), which was recently reviewed in Nat. Rev. Drug Discov. by conference chairman Jean-Paul Renaud and a multinational team of experts. The approach involves flash-freezing a thin film of sample and using transmission electron microscopy to capture two-dimensional “projection” images of your target. If the protein is randomly oriented you can computationally combine thousands of individual images into a three dimensional structure. Although the technique has been around for decades, until recently the resolution was too low to be useful for structure-based drug design. Recent advances in hardware and computation have led to what’s come to be known as the “resolution revolution,” explained Gabe Lander (Scripps).

One advance is the 300 keV Titan Krios – a massive (and massively expensive) instrument that is so widely coveted that Gabe showed pictures of happy scientists hugging newly delivered crates. Indeed, of the ~1000 structures solved to < 4 Å resolution, the vast majority of them were solved on one of more than 130 Krios instruments throughout the world. But Gabe showed that high resolution structures can be obtained with more common 200 keV instruments, including a 2.6 Å resolution structure of aldolase (150 kD), a 2.9 Å structure of hemoglobin (64 kD), and a 2.9 Å resolution structure of alcohol dehydrogenase (81 kD) with bound NAD⁺ cofactor. Although only a handful of sub-2 Å structures have been reported, he thought these would become routine in the next few years.

Bridget Carragher (New York Structural Biology Center) described challenges and how to overcome them. Currently it takes at best eight hours to go from data to structure, but she thought getting this to under one hour would be achievable. Moreover, cryo-EM can be used to characterize different conformational or oligomeric states present in a single sample, as Giovanna Scapin (Merck) demonstrated with insulin binding to its receptor. Indeed, even simple visualization – without fancy computational processing – can provide useful information about protein aggregation, as demonstrated by Wen-ti Liu (NovAliX).

Although primary fragment screening still looks a long way off for cryo-EM, it should start to provide useful structural information for fragments bound to targets less amenable to conventional biophysical techniques, such as membrane proteins – the topic of another session.

Miles Congreve (Heptares) discussed how their stabilized “StaR” GPCRs can provide high-resolution crystal structures suitable for FBDD (see for example here). This has allowed them to discover less lipophilic, more ligand-efficient drug candidates against a variety of targets.

According to Anass Jawhari, it isn’t even necessary to make mutant GPCRs: Calixar has developed proprietary detergents that can stabilize full length adenosine A_2A receptor for a week – more than enough time to perform STD NMR screens of 100 fragments and identify 19 hits, some of which turned out to be functional antagonists. Matthew Eddy (University of Southern California) used two-dimensional NMR on this same protein to reveal dramatic differences in conformational dynamics when bound to agonists vs antagonists.

Indeed, conformational changes and dynamics were a running theme throughout the conference. Keynote speaker and Nobel-laureate Martin Karplus (Harvard) quoted fellow Nobelist Richard Feynman: “everything that living things do can be understood in terms of the jiggling and wiggling of atoms.” (As an aside, Martin’s MCSS method pioneered computational FBDD approaches, predating SAR by NMR.) Göran Dahl (AstraZeneca) described how large scale conformation changes well outside of the active site of PI3Kgamma were responsible for freakishly high selectivity of a class of inhibitors.

But how do you detect conformational changes? We’ve previously mentioned Biodesy’s SHG approach, and Parag Sahasrabudhe (Pfizer) described how this proved useful for classifying ligands for IL-17A. Gerrit Sitters (Lumicks) described a completely different “dynamic single-molecule” (DSM) approach, which involves trapping a single fluorescently labeled protein between DNA strands tethered to two microspheres. Changes in protein conformation caused by ligand binding change the distance between microspheres, and these can be detected to within 1 Å.

Kinetics is intimately linked to dynamics, but the factors responsible for slow binding and dissociation are still poorly understood. Chaohong Sun (AbbVie) examined an archive of 8000 data points and found that on-rates and off-rates each varied by more than five orders of magnitude. There was no correlation with ClogP of the ligands, though larger ligands were more likely to have slower kinetics. There were also significant target effects; on-rates were consistently slow for one target.

As we’ve previously discussed, off-rate screening (ORS) can be used to identify hits in crude reaction mixtures, and Menachem Gunzburg (Monash University) described how this technique is being used in hit-to-lead efforts. Lowering the temperature to 4 °C and adding 5% glycerol further slows dissociation, allowing weaker hits to be discovered.

At the extreme, irreversible inhibitors have an off-rate of 0, and Gregory Craven (Imperial College London) described quantitative irreversible tethering of electrophilic fragments to cysteine residues in proteins using a fluorimetric plate-based assay. As we’ve noted, one challenge with irreversible tethering is deconvoluting intrinsic reactivity from proximity-directed reactivity, which Gregory addresses using a reference thiol such as glutathione.

There is much more to say but in the interest of time I’ll stop here. If you missed the conference you have two chances next year: June 4-7 when it returns to Strasbourg, and November 20-22 when it will be held in Kyoto. And there are still excellent events coming up this year – hope to see you at one!