Demolishing the 'vacuous' argument for the RRW

The Reliable Replacement Warhead program (RRW) is a long-standing and controversial proposal that aims to replace aging plutonium pits and other parts in nukes to modernize the current US nuclear arsenal. Proponents of the RRW say that having nuclear weapons of dubious function and quality will defeat the basic purpose of a deterrent and therefore modernizing the arsenal and replacing worn out parts is essential for its very existence. In addition, Russia has indicated that it is modernizing its nuclear arsenal and this fact has put further pressure on implementation of the RRW. Opponents of the RRW say that by modernizing the arsenal the US will send out the wrong signals to the rest of the world, indicating that nuclear deterrence and weapons development is still an important part of US defense strategy. In any case, very recently the top government advisory group JASON which among others included Freeman Dyson) did a study on the central plutonium 'pits' in the nuclear weapons and concluded that these would last for at least half a decade if not more. Since then several arguments have continued to float around for the RRW.

In the latest issue of The Bulletin, Jeffrey Lewis and Kingston Reif do a neat and clean job in demolishing the latest argument made by a General Chilton who, of all possible reasons, bases his argument on vacuum tubes, the point being that outdated vacuum tubes in nukes necessitate replacement. The last line is priceless and is not exactly BAS-like

Firstly, vacuum tubes are not used in the physics package of a single nuclear weapon design. Vacuum tubes are used only in the radar-fuse, which tells the firing system when the bomb is at the correct altitude for detonation, in some modifications (mods) of one warhead design, the B61 gravity bomb. In total, the B61 bombs that have vacuum tubes in their radar-fuses account for only about one in ten operationally deployed warheads. (Vacuum tubes are used in the radars of three B61 mods: 3, 4, and 7. Mods 10 and 11 have newer radars that use solid-state electronics.) The fuses in these weapons are old, but perfectly functional. To reiterate, vacuum tubes are not in use in any other warhead design, including the W76 warhead, a portion of which would be replaced by the first RRW warhead, the WR1, if it ever were funded and developed.

Secondly, the Energy Department has routinely replaced radars without nuclear testing or redesigning the physics package. In fact, during the 1990s, Sandia National Laboratories scientists developed the MC4033 common radar, which uses solid-state electronics, for planned refurbishments of the B61 and B83 gravity bombs. All B83 bombs now use the common radar, though similar plans to fit a new radar on all B61s have been repeatedly deferred.

Most recently, in 2006, Sandia planned to replace the remaining B61 vacuum tube radars as part of ALT 364/365/366. The National Nuclear Security Administration, which overseas the nuclear weapons complex, canceled these latest ALTs, which would have resulted in the removal of the last vacuum tubes from the U.S. nuclear stockpile, because the U.S. Air Force preferred replacement to life extension. Due to this absurd twist, one could say that vacuum tubes remain in the U.S. nuclear arsenal in part because of the RRW, contrary to Chilton's insistence that the RRW is needed to get rid of them.

The bottom line is that vacuum tubes are used only sparingly in the U.S. nuclear arsenal and can be replaced on short notice if the need arises, independent of whether Congress funds the RRW Program. Of the many reasons that Defense and Energy officials have put forth to justify the RRW Program, the need to replace vacuum tubes is the worst and has no place in the debate about the RRW or modernizing the nuclear stockpile. General Chilton can stick that prop in his, um, pocket.

FDA does, and should, stick to only science

In a welcome reversal of a key politics-driven Bush era mandate, the FDA has approved the Plan B morning-after emergency contraception pill for 17-year olds. Previously the reluctance of FDA to approve the product had led a senior official to rightly resign. Not surprisingly, this decision drew wrath from conservative groups who say that the pill would "encourage promiscuity". This statement is rather typical of conservative statements opposing abortion and promotion of contraceptive measures in school, in spite of the fact that abstinence-only programs have been shown to essentially cause no change or even an increase in "promiscuity".

But here's the thing, and it should be clear all along; the FDA should stick to science and nothing else. Just as the conservative FDA officials during the Bush era were utterly out of line opposing Plan B because of political and religious interests, so should liberals also not applaud the FDA decision as a moral value judgement. The business of the FDA is to determine the efficacy and safety of medical products, period. The moment it starts to pontificate on the moral or political value of its decision its immediately sets itself on a slippery slope.

So just like the NAS and the NCSE should stick to demonstrating the evidence for evolution and lack of evidence for ID/creationism and not pass judgement on whether science and religion are compatible, so should the FDA stick to the science behind the approval of medical products. Not making political or religious statements, either conservative or liberal, would be in the safe and best interests of both the FDA and society.

Timeless Classics

As someone who loves to read more than anything else, I have long also been addicted to classic textbooks. Some of the most memorable moments of my student life involved walking into ghostly libraries looking like medieval castles and dusting off inches of dust collected on tomes which I regarded as treasures, volumes of great works that had not been checked out in 25 years and were languishing in anonymity, begging to be touched and read.

Sadly very few seem to bother about these anymore and regard them as outdated. I can bet that no modern undergraduate that I can meet has browsed Pauling's classic "The Nature of the Chemical Bond", a book that is regarded by many as one of the most influential works in chemistry of all time. In my opinion these students confuse outdated with poorly written. But many of the basics of chemistry don't change, and many of these old works provide crystal clear treatments of basics that are lacking in more modern books. As far as fundamentals go these books have stood the test of time and several are still in print, although some regrettably are not. A comment by Srini about Morrison & Boyd brought back fond memories of favourite classics...

Morrison & Boyd: I have already mentioned it before. Crystal clear treatments of mechanism with an especially outstanding chapter on electrophilic aromatic substitution. For organic chemistry, I have to admit that the new book by Clayden et al. is probably the single-best book on the subject I have seen, but the elegance of explanation in M & B is still hard to beat. I also remember the book by Roberts and Caserio also being pretty good. For what it's worth, the book which compresses the most number of insights in the fewest number of words is a slim volume by Peter Sykes whose clarity in explaining mechanistic concepts in short, crisp paragraphs is unprecedented.

The Great Linus Pauling: I first saw "The Nature of the Chemical Bond" as a freshman. While I then perceived it as boring and too detailed, it was only later that I recognized its monumental significance. Many famous scientists including Max Perutz and Francis Crick have learnt chemistry from it and Pauling's other book. The number of ground-breaking concepts that Pauling invented and put into this book is staggering. Especially check out the chapters on hybridization and partial ionic character of bonds. I have browsed all three editions, and the second one is probably the best-written, although the third edition is the most up-to-date and still in print. A measure of the book's significance in the history of science can be gained from the simple fact that after publication of the first edition the volume was cited no less than 16,000 times in the next 10 years. One constantly keeps on finding new papers in journals like Science and Nature that still cite it. Pauling's "The Nature..." did for modern chemistry what the Principia did for natural philosophy; it infused its subject with logic and tied together disparate threads to formulate a comprehensive and lasting science.

As if one work were not enough, Pauling also authored "General Chemistry". Again, it's a model of simplicity and clarity (note for instance how he explains the source of the difference in the three pKa values of phosphoric acid) although its emphasis on more descriptive chemistry makes it look a little quaint. The text is still widely read and in print as a Dover reprint edition; I have a copy on my shelf for a while now and recently saw one in the Barnes & Noble@GeorgiaTech.

Finally, "Introduction to Quantum Mechanics with Applications to Chemistry" co-authored with E Bright Wilson at Harvard was the first book to explain quantum mechanics to chemists. I will admit the book is not easy to read, but with effort one can find many gems in the first few chapters, especially the treatment of the hydrogen atom. Again, a Dover reprint is available and cheap.

With these three books, the prodigious Pauling secured his place in history not only as the greatest chemist of all time, but one of the most successful and greatest scientific writers of the century.

Glasstone: Samuel Glasstone was a remarkably prolific and versatile technical writer. It was in high-school that I came across a compendium of nuclear science that he had written, "Sourcebook on Atomic Energy". It is hands down the single-best example of technical science writing that I have come across, and wherever I have been since then, I have always had a copy on my bookshelf. The book is a model for comprehensive, all-inclusive writing that is clear as water from a virgin glacier. It also satisfies the difficult condition of being extremely valuable for both laymen and scientists. No praise for the book matches the praise that then AEC chairman and Nobel Laureate Glenn Seaborg penned in the foreword. Seaborg quipped that this was "technical writing at its very best" and that Glasstone was a man who had fortuitously come along at the right time to fulfill the needs of science and technology.

The range of Glasstone's writing is amazing; multi-volume works on physical chemistry and treatments of thermodynamics, electrochemistry, nuclear reactor engineering and even a sourcebook on space sciences. An exhaustively detailed and yet comprehensible book on the effects of nuclear weapons served as a standard declassified guide for years and is still in print. Although his PChem books are now really out-of-date, they still educated me in the basics when I first found them, and I learnt a lot from his book on thermodynamics. But again, nothing beats "Sourcebook on Atomic Energy" a book against which I think every other technical work should be measured.

"Valence" by Charles Coulson: Very few people in history have had the capacity to be both fine scientists and excellent writers. Charles's Coulson's book was the first book, even before Pauling and Wilson, to make quantum chemistry comprehensible to students. When it comes to pedagogical explanation it's hard to beat the British, and this is the finest example of that. I was fortunate to secure a copy at the famous Powell's bookstore in Portland, OR.

Classic books are like old wine. They should be cherished, preserved, and sampled one concept at a time.

The Pope of Cosmology 'very ill'

For a 67 year old man with ALS who has already defied medical science, this is not good news at all. Remember what happened to Christopher Reeve. When you are in a condition like this, even otherwise normal ailments may become life-threatening.

I have been recently reading a lot about Hawking in Leonard Susskind's splendid book "The Black Hole War: My Battle with Stephen Hawking to Make the World Safe for Quantum Mechanics". The rather grandiose title of the book obscures a perfectly entertaining and informative romp through the world of black holes; this is about as close as possible to black hole thermodynamics, string theory and quantum mechanics that we laymen can get without being drowned in a whirlpool of math. Susskind who is a professor at Stanford tells the story of the paradox of information falling into a black hole and supposedly disappearing with lots of verve, hilarious personal anecdotes and tributes to famous physicists. Being a prime participant in the debate with Hawking on the other side, he is in a unique position to tell the story. His recounting of the way the physicist Jacob Bekenstein used high-school math to derive the formula for the entropy of black holes is astounding; very rarely has someone used such simple physics and mathematics to discover such profound relationships and the act reminded me of Bell's Theorem, another spectacular twentieth-century physics result that can essentially be derived using high-school mathematics.

But more than anyone else, it is Hawking's figure that looms large in the book. Susskind describes how his physical disability, his strange disembodied computer voice and his astonishingly brilliant and creative mind guarantees the kind of reverence and silence wherever he appears that otherwise only seems to be reserved for the Pope. Susskind vividly describes a typical Q & A session after a Hawking lecture; Hawking's physical condition means that he can compose even a "yes/no" answer only after several minutes, and what's striking is that during such times Susskind has witnessed audiences of thousands maintain stand-still silence with not a whisper spoken for sometimes fifteen minutes while the great man painfully communicates himself. Hawking may be the only living scientist whose presence provokes utter and rapt silence and attention that one would observe only during religious prayer. No wonder Hawking is compared to God by many, a comparison which only makes him uncomfortable. Susskind describes a particular time in a restaurant where a passerby went to his knees and virtually kissed Hawking's feet. Needless to say Hawking was embarrassed and galled.

In any case, we can only hope that Hawking feels better. However in one way we can rest assured; Stephen Hawking's name has been etched in the annals of science forever. That's the power of ideas. Their timelessness assures us that they remain youthful and vibrant, irrespective of the age and condition of their source. But let's all hope Hawking springs back from this illness to his mischievous, witty self.

Assessing the known and unknown unknowns: WYSI(N)WYG

Ken Dill and David Mobley from UCSF have a really nice review in Structure on computational modeling of protein-drug interactions and the problems inherent in the process. I would strongly recommend anyone interested in the challenges of calculating protein-drug binding to read the review, if not for anything else for the copious references provided. The holy grail of most such modeling is to accurately calculate the free energy of binding. For doing this we frequently start with a known structure of a protein-ligand complex. The main point that the authors emphasize is that when we are looking at a single protein-ligand complex, deduced either through crystallography or NMR, we are missing a lot of important things.

Perhaps the most important factor is entropy which is not at all obvious in a single structure. Typically both the protein and the ligand will populate several different conformations in solution. Both will have to pay complex entropic penalties to bind one another. The ligand strain energy (usually estimated at 2-3 kcal/mol for most ligands) also plays an important role. The desolvation cost for the ligand also can prominently figure. In addition both protein and ligand will have some residual entropy even in the bound state. As if this were not enough of a problem, much of the binding energy can come from the entropic gain that the release of water molecules from active sites engenders. Calculating all these entropies for protein, ligand and solvent is important for accurately calculating the free energy of protein-ligand binding. But there are few methods that can accomplish this complex task.

Among the methods reviewed in the article are most of the important methods used currently. Usually the tradeoff for each method is between cost and accuracy. Methods like docking are fast but inaccurate although they can work well on relatively rigid and well-parameterized systems. Docking also typically does not take protein motion and induced-fit effects into account. Slightly better methods are MM-PBSA or MM-GBSA which as the names indicate, combine docking poses with an implicit solvent model (PBSA or GBSA). Entropy and especially protein entropy is largely ignored, but since we are usually comparing similar ligands, such errors are expected to cancel. Going to more advanced techniques, relative free-energy calculations use molecular dynamics (MD) to try to map the detailed potential energy surfaces for both protein and ligand. Absolute free-energy perturbation calculations are perhaps the gold standard in calculating free energies but are hideously expensive. They work best for ligands that are simple.

There is clearly a long way to go before calculation of ∆Gs becomes a practical endeavor in the pharmaceutical industry. There are essentially two factors that contribute to the recalcitrance of the problem. The first factor as indicated is the sheer complexity of the problem; assessing the thermodynamic features of protein, ligand and solvent in multiple configurational and conformational states. The second problem is a problem inherent in nature; the sensitivity of the binding constant to the free energy. As iterated before, the all-holy relation ∆G = -RT ln K ensures that an error of even 1 kcal/mol in calculation will translate to a large error in the binding constant. The myriad complex factors noted above ensure that errors of 2-3 kcal/mol already constitute the limit of what the best methods can give us. Recall that an error of 3 kcal/mol means that you are dead and buried.

But we push on. One equal temper of heroic hearts. Made weak by time and fate, but strong in will. To strive, to seek, to find, and not to yield. At some point we will reach 1 kcal/mol. And then we will sail.

Reference:
Mobley, D., & Dill, K. (2009). Binding of Small-Molecule Ligands to Proteins: “What You See” Is Not Always “What You Get” Structure, 17 (4), 489-498 DOI: 10.1016/j.str.2009.02.010

Reading C P Snow and The Two Cultures

Over at The Intersection blog which I often read and comment on, Chris Mooney (author of "The Republican War on Science") has initiated an informal reading of C. P. Snow's "The Two Cultures". Anyone who is interested is more than welcome to read the influential and very short lecture and blog or comment on it. The schedule is listed in the post. The recommended edition is the Canto edition, with a very readable introduction by Stefan Collini. Incidentally it was this version that I read many years ago. Time now for a re-reading.

I have encountered Snow in two other interesting books. The first one- "The Physicists: A Generation that Changed the World"- was authored by him and contains clear and abundant photographs as well as recollections and insights on some of the most famous physicists of the century whom he closely knew. In this for instance I read his generous assessment of Enrico Fermi that captures the supreme greatness of the man's talents and achievements

"If Fermi had been born twenty years earlier, it is possible to envisage him first discovering Rutherford's nucleus and then discovering Bohr's atom. If this sounds like hyperbole, anything about Fermi is likely to sound like hyperbole"

Snow also thought that Robert Oppenheimer's real tragedy was not his sidelining or victimization during the 1950s witch hunts but the fact that he would have thrown away all his fame, brilliance and glory if he had the privilege to make one timeless discovery like Pauli's exclusion principle.

Another book with Snow in it is a fascinating piece of "scientific fiction" written by John Casti. "The Cambridge Quintet: A Work Of Scientific Speculation" features four famous scientists and intellectuals- Ludwig Wittgenstein, Erwin Schrödinger, J B S Haldane and Alan Turing- being invited over to Snow's house for a multi-course dinner. As the dinner unfolds, so do the conversations between these stalwarts. The topic is artificial intelligence, and the participants hold forth in myriad and fascinating ways on the subject with excursions that not surprisingly take them into avenues like the philosophy of mind and language, epistemology and metaphysical questions. Very much worth reading.

In any case, I am looking forward to reading The Two Cultures again and writing about it. Anyone who is interested is more than welcome. The entire lecture is 50 pages and could be read in a few hours of thoughtful contemplation. The topic is as relevant today as it was then, which explains the lecture's enduring appeal. The consequences though could be vastly more pronounced.

The rest is all noise: errors in R values, and the greatness of Carl Friedrich Gauss reiterated

One of the questions seldom asked when building a model or assessing experimental data is "What's the error in that?". Unless we know the errors in the measurement of a variable, fitting predicted to experimental values may be a flawed endeavor. For instance when one expects a linear relationship between calculated and experimental values and does not see it, it could either mean that there is a flaw in the underlying expectation or calculation (commonly deduced) or that there is a problem with the errors in the measurements (not always discussed).

Unfortunately it's not easy to find out the distribution of errors in experimental values. The law of large numbers and central limit theorem often thwart us here; most of the times the values are not adequate enough to get a handle on the type of error. But in the absence of such concrete error estimation, nature has provided us with a wonderful measure to handle error; assume that the errors are normally distributed. The Gaussian or normal distribution of quantities in nature is an observation and assumption that is remarkably consistent and spectacularly universal in its application. You can apply it to people's heights, car accidents, length of nails, frequency of sex, number of photos emitted by a source and virtually any other variable. While the Gaussian distribution is not always followed (and strictly speaking it applies only when the central limit theorem holds), I personally regard it to be as much of a view into the "mind of God" as anything else.

In any case, it's thus important to calculate the distribution of errors in a dataset, Gaussian or otherwise. In biological assays where compounds are tested, this becomes especially key. An illustration of the importance in error estimation in these common assays is provided by this recent analysis of model performance by Philip Hajduk and Steve Muchmore's group at Abbott. Essentially what they do is to estimate the standard deviations or errors in a set of in-house measurements on compound activities and look at the effect of those errors on the resulting R values during comparison of calculated activities with these experimental ones. The R value or correlation coefficient is a time-tested and standard measure of fit between two datasets. The authors apply the error they have obtained in the form of "Gaussian noise" to a hypothetical set of calculated vs predicted activity plots with 4, 10 and 20 points. They find that after applying the error, the R-value itself adopts a Gaussian distribution that varies from 0.7 to 0.9 in case of the 20 point measurement. This immediately tells us that any such measurement in the real world that gives us, say, a R value of 0.95 is suspicious since the probability of such a value arising is very low (0.1%), given the errors in its distribution.

You know what should come next. The authors apply this methodology and look at several cases of calculated R values for various calculated vs measured biological activities in leading journals during 2006-2007. As they themselves say,

It is our opinion that the majority of R-values obtained from this (small) literature sample are unsubstantiated given the properties of the underlying data.

Following this analysis they then apply similar noise to measurements for High-Throughput Screening (HTS) and Lead Optimization (LO). Unlike HTS, LO usually deals with molecules sequentially synthesized by medicinal chemists that are separated by small changes in activity. To investigate the effect of such errors, enrichment factors (EFs) are calculated for both scenarios. The EF denotes the percentage of active molecules found or "enriched" in the top fraction of screened molecules relative to random screening, with a higher EF corresponding to better performance. The observation for HTS is that small errors give large EFs, but what is interesting is that even large errors in measurement can give modest enrichment, thus obscuring the presence of such error. For LO the dependence of enrichment on error is less, reflecting the relatively small changes in activity engendered by structure optimization.

The take home message from all this is of course that one needs to be aware of errors and needs to apply those errors in quantifying measures of model assessment. God is in the details and in this case his name is Carl Friedrich Gauss, who must be constantly beaming from his Hanover grave.

Reference:
Brown, S., Muchmore, S., & Hajduk, P. (2009). Healthy skepticism: assessing realistic model performance Drug Discovery Today, 14 (7-8), 420-427 DOI: 10.1016/j.drudis.2009.01.012

Much ado about protein dynamics

Let me alert you, in case you haven't noticed, to the latest issue of Science which is a special issue on protein dynamics.

There is much of merit here, but this article is especially relevant to drug discovery. It talks about the interaction of small molecules and how it reshapes the energy landscape of protein conformational motion. One of the most useful ways of thinking about small molecule-protein interactions is to visualize a protein that fluctuates between several conformational states which are in equilibrium. A small molecule can inhibit the protein by preferentially stabilizing one of these states.

The article illustrates this concept with several examples, most notably inhibition of kinases. Many kinases exist in an inactive and active state and kinase inhibitors stabilize and block one of these states. Such conformational trapping can also reduce mobility of the protein. The article also describes how certain kinase inhibitors such as imatinib and nilotinib trap the kinase in the inactive state while others such as dasatinib trap it in an active state. Although all three of these are classified as ATP-competitive inhibitors, dasatinib blocks the Abl-Bcr kinase by effecting an allosteric movement of a particular loop. Allosteric inhibitors of kinases are of value since they won't target the highly conserved ATP-binding site, thus reducing problems with selectivity. But allosteric targeting is difficult since many times it involves targeting shallow, poorly defined sites including those involved in protein-protein interactions. HTS campaigns aimed at disrupting P-P interactions usually give very poor results. However, recent tools and especially NMR with labeled residues may improve the detection of weakly binding molecules that may be missed in assays (where the limit is usually 30 µM). HSQC spectra are generally taken of the protein, with and without the inhibitor, and changes in residue resonances can give an indication of conformational changes.

In any case, this article and the others are worth reading. Basically it seems that the remodeling of energy landscapes of proteins by either small molecules or other signals is a concept acquiring central traction. Such a concept could essentially tie together the dual problems of protein folding and inhibiting proteins with small molecules.

Reference:
Lee, G., & Craik, C. (2009). Trapping Moving Targets with Small Molecules Science, 324 (5924), 213-215 DOI: 10.1126/science.1169378

New ligands for everyone's favorite protein

A landmark event in structural biology and pharmacology occurred in 2007 when the structure of the ß2-adrenergic receptor was solved using xray crystallography by Brian Kobilka's and Raymond Stevens's groups at Stanford and Scripps respectively. The structure was co-crystallized with the inverse agonist carazolol. Until then the only GPCR structure available was that of rhodopsin and all homology models of GPCR were based on this structure. The availability of this new high resolution structure opened new avenues for structure-based GPCR ligand discovery.

The ß2 binding pocket is especially suited for drug design since it is tight, narrow and lined with mostly hydrophobic residues with polar residues well-separated. Two crucial residues, an Asp and a Ser bind to the ubiquitous charged amino nitrogen present in most catecholamines and the aromatic section of the molecule docks deep into the hydrophobic pocket. These particular features also make computational docking more facile; a mix of polar and non-polar features with bridging waters can make docking and scoring more challenging.

Since the ß2 structure has been published, attempts are being made to use it as a template to build homology models of other GPCRs. A couple of months back I described an interesting proof-of-principle paper by Stefano Costanzi that sought to investigate how well a homology model based on the ß2 would perform. In that study carazolol itself was used as a ligand for docking into the homology model. Comparison with the original crystal structure revealed that while the ligand docked more or less satisfactorily, an important deviation in its orientation could be explained by a counterintuitive orientation of a Phe residue in the binding site. The study indicated that the devil is in the details when one is considering homology models.

However, finding ligands for the ß2 itself is also an important and interesting endeavor. Virtual screening could help in such studies. To this end Brian Shoichet, Brian Kobilka and their group have used the DOCK program to virtually screen one million lead-like ligands from their ZINC database against the ß2. Out of the 1 million ranked poses, they chose and clustered the top 500 compounds (0.05% of the database) into 25 unique chemotypes, a choice also guided by visual inspection of the protein-ligand interactions and commercial availability. They then tested these 25 compounds against the ß2 and found 6 compounds with IC50s better than 4 µM. One of these compounds with an IC50 of 9 nM is perhaps the most potent inverse agonist of the ß2 known. The binding poses revealed substantial overlap of similar functional groups with the carazolol structure. Two compounds turned out to have novel chemotypes and bore very little similarity with known ß2 ligands. A negative test was also run where a known predicted binder was chemical modified so that it would not bind.

Interestingly all the compounds found were inverse agonists. The ZINC library is somewhat biased against aminergic ligands as is most of chemical space. The catecholamine scaffold is one of the favourite scaffolds in medicinal chemistry. However, subtle difference in protein structure can sometimes turn an inverse agonist into an agonist. In this case, small changes in the orientation of the crucial Ser residue near the mouth of the binding pocket. In a past study for instance, slightly changing the rotameric features of the Ser residue thus resulting in a different orientation of the hydroxyl was sufficient to retrieve agonists.

The study thus shows the value of virtual screening in the discovery of new ß2 ligands and indicates the effect of library bias and protein structure on such ligand discovery. Many factors can contribute to the success or failure of such a search; nature is a multi-armed demon.

Reference:
Kolb, P., Rosenbaum, D., Irwin, J., Fung, J., Kobilka, B., & Shoichet, B. (2009). Structure-based discovery of ß2-adrenergic receptor ligands Proceedings of the National Academy of Sciences DOI: 10.1073/pnas.0812657106

Post-docking as a post-doc, and some fragment docking

I am now ready to post-doc. I am also now ready to post-dock, that is, engage in activities beyond docking. Sorry, I could not resist cracking that terrible joke. It's been a long journey and I have enjoyed every most moments of it. Thanks to everyone in the chemistry blogworld who regaled, informed, provoked and entertained on this blog. I am now ready to move on to the freakingly chilly Northeast. Location not disclosed for now, but maybe later.



Speaking of docking, here is a nice paper from the Shoichet group in which they use fragment docking to divine hits from a large library for a beta-lactamase. Fragment docking can often be tricky compared to "normal" docking since fragments being small usually demonstrate promiscuity, low-affinity and non-selectivity in binding. Fragment docking thus is not yet a completely validated technique.

In their study, the present authors screen their ZINC library for fragments binding to the ß lactamase CTX-M by docking using the program DOCK. They also screen a lead-like library for larger molecules. The top hits from the fragment docking results were assayed and showed micromolar inhibition against the lactamase. These included tetrazole scaffolds not seen before. Importantly, five of these hits could be crystallized and the high-res crystal structures validated the docking modes.

What was interesting was that the same tetrazole scaffolds in the larger lead-like library were ranked very low (>900) and would not have ever been selected had their tetrazole fragments not showed up at the top in the fragment docking results. These compounds, when assayed showed sub-milimolar to micromolar activity against the lactamase. Thus, the protocol essentially demonstrated that fragment docking can reveal hits that can be missed by docking larger lead-like molecules. One of the reasons DOCK succeeds in this capacity is because of its use of a physics-based scoring function that has no bias against fragments. It also helps that the active site of CTM-X is relatively rigid with little protein motion.

The fragments were also assayed against another lactamase for Amp C. Usually, hits for CTM-X and Amp C are mutually exclusive. What was seen was that the higher the potency of the fragments for CTX-M, the higher the specificity for CTX-M, not surprising considering that increased potency translates to a much better complementary fit of the fragments for CTX-M.

Fragment docking can be messy since fragments can bind non-selectively and haphazardly to many different parts of many different proteins. But this study indicates that fragment docking is not an uninteresting strategy to possibly find hits from other lead-like libraries that may be otherwise concealed.

The potencies of the compounds found may look pretty weak, but because there are extremely few molecules inhibiting these medicinally important lactamases, such advances are welcome. Lactamases are of course an important target for overcoming resistance in antibiotic treatment.

Reference:
Chen, Y., & Shoichet, B. (2009). Molecular docking and ligand specificity in fragment-based inhibitor discovery Nature Chemical Biology DOI: 10.1038/nchembio.155