In the search for breakthrough medicines in biomedical research, the rapid preparation of new, complex molecules for biological evaluation is of paramount importance, but substrates for the synthesis of such compounds are invariably in short supply. In later stages of chemistry development, substrates are abundant, and state-of-the-art microvial (8
) or microfluidic
tools can be effective in “turning on” reactions that were otherwise unsuccessful by exploring combinations of catalysts
, reagents, and other key reaction variables (12
). Such studies require at least milligram (micromole) quantities of substrate per reaction—a prohibitively large amount in early drug discovery, where new molecules are prepared for the first time. Consequently, the set of desirable compounds designed to test a biological hypothesis is often winnowed to the much smaller subset of compounds that can be successfully synthesized using a single set of reaction conditions, with little opportunity to study and improve unsuccessful syntheses (14, 15).
Miniaturizing chemistry to the nanomole scale is a potential solution to this problem that has heretofore met with substantial engineering problems, such as accurately bringing together extremely small charges of materials that are often heterogeneous, effectively agitating reaction mixtures, preventing loss of volatile solvents, and incorporating general analytical approaches to assay reaction outcomes. We present the first results of a study aimed at developing general nanomole reaction screening capabilities to support the rapid synthesis of complex, highly functionalized drug leads.
Fig. 1. Nanomole-scale reaction screening.
) High-complexity electrophilic
from the Merck compound collection. (B
) High-polarity nucleophilic
building blocks 9
) Route scouting to identify successful reaction conditions for previously failed chemistry can be accomplished using nanomole-scale reaction screening with minimal consumption of precious substrate (0.02 mg per reaction). Successful reactions run on nanoscale (0.02 mg, 50 nmol) can be repeated on a 25-mg scale (factor of 1000 scale-up). A 96-member parallel coupling library using best reaction conditions located by screening model nucleophiles (table S3) gave 54 products (white squares) and 42 failed reactions (blue squares) (table S4). After nanomole-scale screening, productive conditions were located for an additional 21 of the 96 products (yellow squares) (tables S6 and S8). (D
) Reaction mixtures that are fully soluble in high-boiling solvents can be dosed with high-precision nanoliter robotics. Room-temperature reactions require no stirring and minimal sealing. (E
) High-performance LC-MS MISER
analysis allows for rapid label-free analysis of reactions with diverse products. (F
) Reaction discovery: Organic superbases 24
promote room-temperature cross-coupling in DMSO
with second- and third-generation biaryl palladium precatalysts 30
) Ninety-six combinations of bases and ligands were investigated in triplicate in glass microvials with stirring and reproduced in 1536-well plastic microtiter plates without stirring.
Panels A and B give a list of eight electrophilic cores from Merck collection and 12 nucleophilic building blocks.
Panel F gives a list of six bases and 16 Palladium catalysts.
The ultimate purpose
Panel C shows that out of 96 (eight electrophiles times 12 nucleophiles) reaction combinations, blue blocks (42) represent failed reactions, but 21 (yellow blocks) of these succeed after implementing the “nanomole route scouting” step.
Panel D shows that nanoliter robotics comes in handy to prepare reaction mixtures.
Panel E shows that HPLC-MISER can be used for rapid analysis of products resulting from these reactions.
Validating the nanomole concept
Panel G shows that 96 (six times 16) base/catalyst combinations with a single electrophile/nucleophile (bromide/amine) combination at room temperature with DMSO as solvent, when evaluated at 1 µL (100 nmol) scale in microtiter plates gave the same reaction profiles as reactions evaluated at 25 µL scale in glass microvials.
Soluble superbase coupling conditions for miniaturization from micromole to nanomole scale
Using the catalysis cross-coupling space defined in Fig. 1, A and B
, as a test case, this work evaluates the idea that miniaturized chemistry experimentation to inform preparative scale synthesis (Fig. 1C) could be run on the same scale as contemporary biological assays using the high-precision nanoliter robotics commonly used in biochemistry labs (Fig. 1D) in conjunction with high-throughput mass spectrometric analysis
). The engineering complexities of this task could potentially be overcome by using ambient-temperature chemistry
that operates in a solubilizing, low-volatility, plastic-compatible reaction medium such as dimethyl sulfoxide (DMSO
) or N
). We chose DMSO
for the initial proof of concept because it is an environmentally friendly solvent and represents the potential to directly link chemistry with biochemical experimentation. Very few efficient ambient-temperature palladium-catalyzed cross-couplings in DMSO or other polar solvents have been reported (21)
, in part because of the strong solvent-metal coordination
of such solvents and their incompatibility with the strong bases typically used in low-temperature C-N couplings (NaOt
Bu, LiN(SiMe3)2, Zn(N(SiMe3)2)2 etc.) (22
). It seemed possible that recently reported highly hindered, electron-rich ligands (23) may effectively protect Pd from DMSO coordination and that non-nucleophilic organic superbases (24), which are soluble in most solvents and should be strongly basic enough to promote Pd C-N couplings at room temperature, could be compatible with these solvents.
Initially working in glass microvials, which represent the current validated lower limits of scalable batch miniaturization (8
), the coupling of amine 21 with bromide 22 (2.5 μmol in 25 μl of DMSO) was evaluated in a 96-reaction array of superbases 24 to 29 with catalysts 30 to45
[10 mole percent (mol %)] (Fig. 1F
), and several highly productive room-temperature coupling conditions were identified
. All solutions were homogeneous but were stirred nonetheless. This same 96-reaction array was then evaluated on a 1.0-μl scale (100 nmol of substrate) in a plastic 1536-well plate without stirring and with reagent dosing via a nanoliter liquid-handling robot (TTP LabTech Mosquito HTS) inside a nitrogen-filled glovebox
, demonstrating the same hits and reaction performance
). The Mosquito robot allows sequential aspiration from different wells of a source plate into a single pipette tip to create multicomponent mixtures that are dosed as single reaction drops into a 1536-well plate. This ensures proper mixing of the reaction components and obviates the need for stirring. This pipetting mechanism, in conjunction with a material-sparing source plate with minimal dead volume, ensures that almost no material is wasted, and the overage can be recovered if desired. Also, with this batch reaction approach, some degree of sample heterogeneity can be tolerated, unlike for microfluidic technologies where channel clogging can lead to systemic failure.
The coupling of bromide 22 was then evaluated with 16 nitrogen, oxygen, carbon, phosphorus, and sulfur nucleophiles under the same 96 combinations of 24 to 45 in a single miniaturized 1536-reaction experiment (Fig. 2A and table S3). In this experiment, dosing required 30 min (four components per well, 6144 total reagent doses) followed by 22 hours of reaction time, 1 hour for automated reaction quenching and sampling, and 52 hours for analysis by ultra performance liquid chromatography (UPLC)
. This experiment uncovered extremely mild, soluble coupling conditions for amines, alcohols, amides, sulfonamides, carbamates, amidines, aryl boronates, alkynes, and malonates, most of which have not been reported at ambient temperature with non-nucleophilic bases (table S3). Because the reactions were by design homogeneous and required no heating or cooling, translation to a larger scale was straightforward, delivering 68 to 98% yield when scaled up by a factor of 3000, to 320 μmol of22 (50 mg), now running at 5 mol % catalyst loading and 0.2 M substrate concentration.
Fig. 2A. high-throughput nanomole-scale chemistry evaluation using rapid MISER and mixed-MISER LC-MS.
) Heat map
showing data from 1536 nanomole-scale reactions analyzed by UPLC
(total analysis time ~52 hours). Sixteen diverse N/C/O/P/S nucleophiles were screened against 96 catalyst-base combinations to find room-temperature coupling conditions to 3-bromopyridine 22
. Best conditions for each nucleophile were repeated on a 50-mg scale, giving 68 to 95% isolated yields (table S3). Reaction information about best catalyst-base combination was matched to the polar nucleophile used in the 96-member parallel coupling library of 1
, initial synthesis). (B
) Heat map showing data from 1536 nanomole-scale reactions analyzed by MISER
. Thirty-two electrophile-nucleophile combinations that failed under previous best conditions were investigated under 48 catalyst-base combinations. Best conditions for each combination were repeated on a 50-μmol scale and assigned a “PASS
” (Fig. 1C
, after screening) if a >95% purity sample was obtained. (C
chromatographic output (MISER
gram) for 48 catalyst-base combinations that were screened with 0.02 mg of 2
per reaction, using 9
as nucleophile, and depiction of pooled MISER
) Nanomole-scale DOE
: three-factorial, four-level DOE
surface model produced by screening 6
against reagent charges of 10
, and 42
. *Conditions identified in 0.02 mg scale reactions informed practical 25-mg and 1-g scale-ups.
1536 Stage 1 – Simple substrates
Panel A shows that 16 nucleophiles were coupled with one electrophile (3-bromopyridine) and 96 base-catalyst combinations and these reactions were evaluated in a 1536 plate.
Reaction time was about 52 hours and ultra performance liquid chromatography (UPLC) was used for product analysis.
Best conditions for each nucleophile gave 68% to 95% product yields at 50-mg scale.
1536 Stage 2 – Complex substrates
32 electrophile-nucleophile combinations (from initial 96-member coupling library (one to eight times nine to 20) that had failed under previous best conditions) were coupled with 48 catalyst-base combinations and reactions analyzed by MISER (typical chromatographic output shown in Panel C).
Reaction time was about 2.5 hours. 21 combinations were confirmed as hits and scaled up to 25-mg scale. Pure product was obtained for 16 of these reactions. Panel C just shows a sample chromatographic output from MISER for product identification from each reaction.
Design of experiments (DOE) approach and scale up
Panel D shows how three factorial, four level statistical approach to DOE can be run at nanomole scale to identify optimal reaction conditions on practical scale such as 25 mg and then 1 g.
For example, electrophile No. 6 combined with nucleophile No. 10 and base/catalyst numbers 29 and 42 was identified as a successful reaction combination at 0.02-mg scale and could be scaled up to 25-mg scale and 1-g scale with 79% and 76% yields, respectively.
Implementation of nanomole screening for complex molecule synthesis
To apply nanomole-scale reaction screening to this problem, we subjected 32 different electrophile-nucleophile substrate combinations, each of which produced little or no material from the 96-reaction array, to 48 reaction conditions (8 organic superbases × 6 ligands; Fig. 2B). The substrate concentration was reduced to 0.05 M and the catalyst loading was doubled to 20 mol % so that running all 48 reactions would consume less than 1 mg of each substrate (50 nmol, ~0.02 mg per reaction) (27). In addition, the reaction time was reduced to 2 hours and a faster liquid chromatography–mass spectrometry (LC-MS) approach was pursued, with a goal of running and analyzing 1536 reactions in less than 24 hours. Of the available label-free techniques (28), MISER (multiple injections in a single experimental run) LC-MS was selected for its ease of data acquisition and analysis
). Thus, 48 reactions that had the same potential product along a row of the 1536-well plate were analyzed by multiple injections (22 s per sample) in a single isocratic chromatographic run with mass detection settings observing the molecular ion for the desired product (Fig. 2C
). In this way, the 1536-well plate could be analyzed in ~9 hours. To reduce the time required for analysis even further, we combined four rows of reaction samples that could give different product masses and simultaneously monitored the four desired products’ m/z
values (mass/charge ratios) in these pooled samples (Fig. 2C
). By pooling four mass-encoded wells, the analysis of all 1536 reactions could be performed in ~2.5 hours. In addition, a more rigorous approach to this same analysis was pursued using a parallel two-channel LC
system connected to a triple-quadrupole mass spectrometer. This instrumentation used a fast gradient to provide improved resolution and additional MS
structure confirmation data and required ~2.5 hours for analysis of all 1536 reactions. The MISER analyses located the expected product masses for 21 of 32 substrate combinations, and the best hits were then confirmed with UPLC-MS. These 21 hits were scaled up by a factor of 1000 (25-mg scale), and a pure product sample was obtained for 16 of these reactions
Nanomole reaction optimization to inform gram-scale synthesis
An important consideration for this nanomole-scale approach toward discovery synthesis is the optimization of reactions for application to larger-scale synthesis. Compounds that are deemed interesting in initial biological evaluation or are synthetic intermediates in a multistep route must rapidly be prepared in larger quantities once conditions are identified. Statistical tools such as DOE
(design of experiments) are powerful tools for synthesis (30
), yet they typically focus on performing minimal numbers of experiments, as large numbers of experiments are often resource- and time-intensive to conduct. However, because large numbers of experiments are readily feasible with this nanomole-scale chemistry platform, we were able to construct a three-factorial, four-level response surface modeling experiment to study the loading of catalyst against varying stoichiometries of base and nucleophile in the reaction of chloride 6 with amine 10. In this DOE experiment, each condition was repeated twice, resulting in 128 total reactions with <3 mg of 6
. Indeed, a high-quality response surface model was generated with the nanomole-scale chemistry approach (Fig. 2D),
which helped to define the critical charges of nucleophile and base for optimal reaction performance. The optimized conditions used 15 mol % 42 at 0.05 M concentration; by translating to more practical conditions of 5 mol % 42 and 0.24 M concentration, we obtained full conversion and a 79% isolated yield of 46 on a 25-mg scale, which was reproduced to obtain a 76% isolated yield on a 1-g scale (Fig. 2D). This result shows that advanced statistical reaction analysis, which is typically reserved for chemistry opportunities where material is plentiful, can be applied to reactions in the material-limited front lines of drug discovery or natural product synthesis.
In biomedical research, chemical synthesis should not limit access to any molecule that is designed to answer a biological question. This work demonstrates an example of how conditions for complex Pd-catalyzed C-O, C-N, and C-C cross-coupling reactions can be evolved into a powerful, substrate-focused approach to chemistry miniaturization to overcome limited access to complex products. With innovative research, other high-value modern chemistry reactions could be similarly designed into this paradigm to improve synthesis in material-limited environments by evolution of catalysts
and reagents to perform in DMSO
, or other high-boiling solvents at ambient temperature.
Materials and Methods
Figs. S1 to S32
Tables S1 to S11
Data Files S1 to S5
References and Notes
M. R. Friedfeld et al., Science 342, 1076–1080 (2013).
D. A. DiRocco et al., Angew. Chem. Int. Ed. 53, 4802–4806 (2014).
D. W. Robbins, J. F. Hartwig, Science 333, 1423–1427 (2011).
A. McNally, C. K. Prier, D. W. C. MacMillan, Science 334, 1114–1117 (2011).
K. D. Collins, T. Gensch, F. Glorius, Nat. Chem. 6, 859–871 (2014).
R. Moreira, M. Havranek, D. Sames, J. Am. Chem. Soc. 123, 3927–3921 (2001).
S. M. Preshlock et al., J. Am. Chem. Soc. 135, 7572–7582 (2013).
A. Bellomo et al., Angew. Chem. Int. Ed. 51, 6912–6915 (2012).
J. R. Schmink, A. Bellomo, S. Berritt, Aldrichim. Acta 46, 71–80 (2013).
M. Peplow, Nature 512, 20–22 (2014).
T. Rodrigues, P. Schneider, G. Schneider, Angew. Chem. Int. Ed. 53, 5750–5758 (2014).
S. Monfette, J. M. Blacquiere, D. E. Fogg, Organometallics 30, 36–42 (2011).
P. M. Murray, S. N. G. Tyler, J. D. Moseley, Org. Process Res. Dev. 17, 40–46 (2013).
S. D. Roughley, A. M. Jordan, J. Med. Chem. 54, 3451–3479 (2011).
">T. W. J. Cooper, I. B. Campbell, S. J. F. Macdonald, Angew. Chem. Int. Ed. 49, 8082–8091 (2010).
A. Nadin, C. Hattotuwagama, I. Churcher, Angew. Chem. Int. Ed. 51, 1114–1122 (2012).
Merck internal study of electronic notebooks.
M. M. Hann, G. M. Keserü, Nat. Rev. Drug Discov. 11, 355–365 (2012).
F. Lovering, J. Bikker, C. Humblet, J. Med. Chem. 52, 6752–6756 (2009).
H. A. Malik et al., Chem. Sci. 5, 2352–2361 (2014).
R. E. Tundel, K. W. Anderson, S. L. Buchwald, J. Org. Chem. 71, 430–433 (2006).
D. S. Surry, S. L. Buchwald, Chem. Sci. 2, 27−50 (2011).
N. C. Bruno, M. T. Tudge, S. L. Buchwald, Chem. Sci. 4, 916–920 (2013).
T. Ishikawa, Y. Kondo, H. Kotsuki, T. Kumamoto, D. Margetic, K. Nagasawa, W. Nakanishi, in Superbases for Organic Synthesis: Guanidines, Amidines, Phosphazenes and Related Organocatalysts, T. Ishikawa, Ed. (Wiley, West Sussex, UK, ed. 1, 2009), pp. 1–326.
Compounds were generally purified by MS-directed purification. Isolated yields ranged from 1 to 100%, but we made no attempt to maximize the isolated yields in these reactions and instead focused on obtaining high-purity compounds as quickly as possible, which is typical in most medicinal chemistry campaigns. Some reactions showed product formation by UPLC-MS analysis but were either insufficiently pure or too low in yield for purification.
M. Liu et al., ACS Comb. Sci. 14, 51–59 (2012).
Some electrophiles were not fully soluble in DMSO, so NMP was used instead. Even though three stock solutions in NMP still displayed mild insolubility, the TTP Mosquito operates on positive-displacement pipetting, so viscous solutions or suspensions of small particulates are easily transferred.
W. Schafer, X. Bu, X. Gong, L. A. Joyce, C. J. Welch, in Comprehensive Organic Synthesis, C. J. Welch, Ed. (Elsevier, Oxford, ed. 2, 2014), vol. 9, pp. 28−53.
C. J. Welch et al., Tetrahedron Asymmetry 21, 1674–1681 (2010).
J. C. Ianni, V. Annamalai, P.-W. Phuan, M. Panda, M. C. Kozlowski, Angew. Chem. Int. Ed. 45, 5502–5505 (2006).
S. E. Denmark, C. R. Butler, J. Am. Chem. Soc. 130, 3690–3704 (2008).
K. C. Harper, M. S. Sigman, Science 333, 1875–1878 (2011).
Acknowledgements We thank S. Krska, M. Tudge, G. Hughes, and E. Parmee for helpful discussions; M. Liu, E. Streckfuss, T. Meng, N. Pissarniski, and W. Li for assistance in purification of compounds; M. Christensen and J. Voigt for experimental assistance; and S. M. O’Brien and M. McColgan for graphic design. S.B. was supported by an NSF GOALI Grant associated with the University of Pennsylvania. Supported by the MRL Postdoctoral Research Fellows Program (A.B.S. and E.L.R.).