- Mini review
- Open Access
Towards EMIC rational design: setting the molecular simulation toolbox for enantiopure molecularly imprinted catalysts
© The Author(s) 2016
- Received: 3 June 2016
- Accepted: 20 October 2016
- Published: 25 October 2016
- Racemic mixtures
- Prochiral substrates
- Transition states
- Ab-initio simulations
- Molecular dynamics
- Monte Carlo
Nature as a whole is a chiral system, many of the molecules that constitute living organisms are chiral and, in the vast majority of cases, preference is shown for one of the enantiomers. For example, proteins are formed exclusively of the L form of amino acids. Meanwhile, saccharide units of the D form singularly constitute carbohydrates; in the same manner, enantiomeric forms in the building blocks of DNA and RNA (d-ribose or d-deoxyribose) have been observed .
Enantiomers have the same physical properties with the exception that they interact differently with polarised light. Regarding the chemical properties, both enantiomers solely differ in their reactivity with other chiral molecules. Hence, a chiral molecule only manifests itself as such by the influence of polarised light or other chiral molecules.
Biological systems, such as proteins and enzymes that catalyse life’s essential reactions have a three-dimensional structure and establish preferences to interact with one of the enantiomers of other molecules. The effect of these interactions is the basis for the study of chiral drugs. As a result of their chirality, racemic drugs can have different effects on our bodies. There are chiral drugs in which each one of the enantiomers could produce opposite effects in the organism, in other cases, the effect is similar, but one of the enantiomers is more active than the other (eutomer and distomer, respectively). While in some cases, one enantiomer is active and the other is inactive and also can occur that one enantiomer has a beneficial effect meanwhile the other is toxic.
Through an evolutionary pathway, nature has become stereoselective, being capable of synthesising the best for a purpose of the enantiomers. A practically endless list of chiral compounds provided by nature can be compiled. The tobacco leaves only produce the levorotatory S-nicotine. The coca only makes S-cocaine. The sugarcane generates d-sucrose exclusively. Limonene is an interesting case, which implies that genetic information drives the biosynthesis of the enantiomers, the dextrorotatory d-limonene is found in the orange or lemon peel. Meanwhile, in the mint, it is found as the levorotatory l-limonene and in the turpentine (derived from pines) as the racemic mixture (±)-limonene .
Chiral drugs dominate the modern pharmaceutical landscape, making up to 40–50% of the market in 2013 with 9 of the top 10 bestseller drugs being chiral . These drugs are sold as racemic mixtures or as a single enantiomer. Currently, there is a significant trend in the pharmaceutical industry to produce what is called “chiral switches”: chiral drugs already commercialised as racemates that could be developed as a single enantiomer . The idea behind these chiral switches is the fact that the enantiomers exhibit different behaviour when they are exposed to the chiral environment that is the human body. This discrimination between enantiomers—or chiral recognition—depends on the degree of interaction that each enantiomer exhibits with the chiral binding site in the body.
Pointing out the enantioselective action of chiral drugs at the beginning of modern pharmacology was regarded as vain within the global profile of drug activity. Nowadays, this is no longer the case. At this very moment, most of the existing patents for drugs consisting of racemic mixtures are coming to an end and the race to obtain new ones for enantiopure production has already begun .
Therefore, there is a need for systematic studies to enhance the understanding of eutomers and to guide their stereoselective synthesis. This work introduces the most relevant molecular simulation methods to help in the design of enantiopure molecularly imprinted catalysts, EMICs. A well-designed EMIC would create a considerable impact in the way the synthesis of enantiopure drugs is performed. An EMIC could circumvent the effort involved in separating racemic mixtures and enable direct access to the eutomer, which in turn reduces the necessary dosage and the chronic side effects of the racemate as well as simplifying dosage-effect studies.
The description of different molecular simulation techniques for the study and development of these efficient catalysts throughout their synthesis stages is the principal purpose of this contribution and the central pillar of our on-going research efforts, translating the principles of enzyme catalysis to the design of EMICs from a molecular perspective.
Currents paths from racemate to enantiopure drugs
The separation of racemates into their enantiomers is a difficult task, e.g. distillation cannot be employed, as both enantiomers will have the same bubble point. To achieve an enantiopure separation the technique used must discriminate based on the stereo orientation of the enantiomer. The most relevant categories of chiral drugs and the current ways to obtain the enantiopure ones follow.
Most chiral drugs have one key bioactive enantiomer. In this case, one of the enantiomers, the eutomer, is much more active and efficient than the other. The distomer can either be less active, toxic or produce undesirable effects. Drugs that fall under this category will often benefit from the synthesis of an enantiopure drug, e.g. ethambutol, whereas the (S,S)-(+)-enantiomer is used to treat tuberculosis, the (R,R)-(−)-ethambutol causes blindness .
Some chiral drugs have equally bioactive enantiomers. Here, the two enantiomers would have the same activity and identical pharmacodynamic properties. There are only a few chiral medications that may fall under this category, but none has been confirmed .
Finally, some chiral drugs can undergo chiral inversion in the body. These drugs have the unique property that the eutomer or distomer can be converted into the other by our body. For these drugs, it can be unnecessary to develop a single enantiomer drug. For example, in the case of ibuprofen, while the racemic mixture is 50/50 when administered, some distomers are converted into eutomers in the body, ultimately making the drug more potent .
Synthesis of diastereomeric salts by treatment with an enantiomer. The salts of the two enantiomers have different solubilities, allowing them to be separated from each other.
Utilising the various reaction rates of the two enantiomers with the addition of a different enantiopure compound. Up to 50% of the enantiomer that reacts more slowly can be recovered from the racemate.
Other resolution of racemates also takes advantage of the differing reaction rates to separate the mixture. However, the unrecovered enantiomer is converted back into a racemic mixture. This process is then repeated until a higher yield of the eutomer is recovered.
Some approaches take advantage of naturally occurring enantiopure compounds. The natural enantiopure compound is modified to create the desired enantiopure drug. This method is extremely useful when the product you want has a similar chemical structure to the naturally occurring enantiomer and is used in such cases.
Synthesis of the enantiopure compound from prochiral substrates by the introduction of a chiral auxiliary to the racemic substrate mixture to separate the two enantiomers. The auxiliary is then removed post-separation. This method is effective, but the auxiliary is required in a stoichiometric quantity. Because of this, the auxiliary must be cheap and easy to produce.
Selective adsorption, a stereoselective adsorbent is used to remove only one enantiomer thoroughly from the racemate. This method has a large advantage over using an auxiliary because sub-stoichiometric amounts of the adsorbent can be used (and re-used) for an adequate separation. Current separation techniques include the use of enzymes and homogeneous chiral metal containing complexes.
Here, we propose EMICs as a seventh alternative to obtaining the eutomer avoiding the racemic separation. Such a catalyst could be done by exploiting on the field of molecularly imprinted polymers (MIPs) , the basic concepts of molecular imprinting can be adopted to create a catalytic polymer network that will promote the transition state (TS) of a particular reaction in a lock and key fashion.
Following nature’s example
Natural enzymes possess an arrangement of functional groups responsible for their specificity [11, 12]. The substrate-enzyme binding interactions are rather complex and consist of a combination of electrostatic interactions, hydrogen bonds, hydrophobic interactions, and other contributions. Then, some prerequisites have to be fulfilled for the preparation of a material showing enzyme-like catalytic activity towards the eutomer, i.e. to construct an EMIC.
First, a cavity has to be made with a defined shape. This shape can correspond to the substrate or, even better, to the TS of the reaction. Due to the TS instability a transition state analogue, TSA, must be found. The cavity can also adopt the shape of the eutomer. Functional groups have to be introduced to act as binding sites within the cavity in a defined stereochemistry. These requirements were introduced with the imprinting protocol conceptualised by Dickey  and implemented by Wulff  and Mosbach .
Following the technique described above, recent years have seen remarkable progress in the design of molecularly imprinted catalysts . Numerous reviews on the molecular imprinting procedure have been published [21–27]. However, a comparison between these catalysts and enzymes , shows that enzymes are still in every case several orders of magnitude catalytically more efficient, but in a few cases, the efforts have reached the activity of catalytic antibodies, e.g., in the hydrolysis of carbamates .
What do modelling techniques have to offer for the design of better MICs and EMICs?
The binding site homogeneity in enzymes is high, whereas MICs, have a broad distribution of activity, and there is no method available at the moment to reduce this broadness significantly. Though some progress has been made in the preparation of MICs, for modest use in industry and wider application in research, refining the experimental imprinting procedures with insights gained from statistical mechanics tools could make further developments.
Improvement of the mass transfer in the imprinted networks, reduction of the polyclonality of cavities, an increase of available active sites (in particular with the frequent noncovalent interaction) and development of further suitable groupings for catalysis are just some problems at the forefront of investigations [20, 22]. Some researchers have concluded that a larger extent of self-assembly can result in a higher specificity. Others have claimed that the shape of the imprinted cavity is the main aspect of molecular recognition and that a change in the form will result in a lower level of identification. More recently, more and more researchers tend to support a modest extent of self-assembly as the condition for the strongest molecular recognition. Undoubtedly, the design of MICs is attracting an extensive research effort .
The idea behind EMICs is clear and straightforward, but the huge pool of variables in its synthesis and characterization requires some rational screening strategies. We believe these strategies could evolve from simultaneous and synergic use of modelling tools with experimental work for the sound design of EMICs in silico. All the variables involved in the synthesis can be independently controlled, and their impact systematically assessed to prepare better catalysts. With molecular models, we seek to understand how imprinted materials are created and what happens to TSAs and substrate molecules in these imprinted cavities to become eutomers. This will help to elucidate the different contributions of each parameter to the overall catalytic effect with the use of proper control systems, for the ultimate developments of better MICs for various reactions and EMICs for specific diseases. Figure 1c, obtained through our molecular dynamics (MD) methodology , shows the complexation, polymerization and cavity rebinding points for a pyridine-selective polymer.
Some detailed atomistic simulations have been employed for the computational design of imprinted polymers [30–34]. Recently, a more general approach derives a set of design principles and backs up the possibility of efficiently imprinting drugs , although very few specific examples of molecular modelling efforts could be found for MICs design . While some interesting insights have been gained [30–34], most of these efforts suffer from two significant drawbacks. First, they focus on a single cavity (neglecting issues related to the heterogeneity of binding sites and porosity). Second, the material optimisation is reduced to a simplified scoring function based on the internal energy of complexation, rather than on proper adsorption or rebinding isotherms or reaction yield as measured in experiments.
We should aim to develop models and methodologies that feature a sufficient level of realism and detail, specifically based on accurate force fields, and that reflect some underlying principles behind the materials formation and function. These protocols should imitate the actual process of MICs formation, characterization and applications within four stages of development:
Stage 1 involves a mixture of TSA (or substrate or product), functional monomers and cross-linkers. Ab initio calculations are envisaged to identify the plausible TSA of the desired reaction, obtain the partial charge distributions of these structures, and describe the complex TSA cavity-regarding binding site energy. The equilibrium properties of the mixture (TSA, functional monomer, cross-linker and solvent) could be obtained by molecular dynamics  or Monte Carlo  approaches to mimic the synthesis conditions in NPT and NVT ensembles.
In Stage 2, the polymerization of functional monomers and cross-linkers should be modelled from the equilibrated mixture structure (a direct outcome of Stage 1). The idea is to focus on the generation of the functional/selective catalytic cavity-ignoring, for a while, the details of the network formation. However, the explicit account of the new bonds formed during the polymerization step can be attained by kinetic Monte Carlo .
In Stage 3, the imprint TSA is removed. The model can be further extended to imitate some post-formational modifications, such as cavity shrinking and introduction of defects. This resulting structure would serve as a porous matrix for both the structural characterization and applications. The resulting model EMIC could be used as the material structure for adsorption studies. The grand canonical Monte Carlo (GCMC) method  is appropriate to describe the re-binding behaviour of the EMIC under study.
In Stage 4, before the reaction, the TSA is bound to the EMIC in a pre-equilibrium step. The bound TSA is converted under catalysis of the TSA-EMIC to the product and is then released. At this final stage, the different reaction kinetics (regarding rates of reactions of different orders of magnitude expected for the diffusion and binding of the template to the polymer) can be investigated by using the probability-weighted dynamic Monte Carlo method . Complex molecular geometries may require the employment of advanced techniques such as configurational bias Monte Carlo  and cavity/energy bias Monte Carlo  to efficiently explore the binding sites.
As briefly described above, required modelling tools for EMIC’s rational design protocol (models and methods) are available, but so far these tools remain unrelated to the field. The compilation of such a computational toolbox would encompass linking various pieces of research together in a consistent workflow, standardising inputs and outputs between the stages and methods. Many aspects of the sketched protocol are challenging and require substantial expertise in the areas of molecular modelling, programming and statistical mechanics. However, the expected outcomes in the understanding of these systems are worth the effort.
The largest gain of the proposed theoretical approach to these systems is that it will allow going beyond current knowledge and exploring these novel formulations. However, the validation of any computer simulation strategy requires the comparison with nature; i.e. the model must be able to reproduce the essential properties of a system that has been already explored experimentally. An excellent source for test cases is the first book in the area of MICs with its substantial amount of experimental work and applications .
This toolbox could be very useful in improving the scope and applicability of MICs for more advanced catalysis, as the EMICs proposed here, (i.e. the selective catalysis of enantiopure drugs). The fundamental efforts of the described tasks would help to ask, and hopefully to answer, the “what if” questions for a range of possible catalytic systems, focusing on the in silico performance rank of the candidate materials, bypassing the economic constraints of such search through real experiments. As a result, EMICs could be synthesised to corroborate, or dispute, the predictions and guide the ultimately necessary experimental work.
For instance, one can compare the re-binding affinity of a synthesised EMIC using the TSA, against the theoretical adsorption affinity of that EMIC but imprinted with the TS. Such a comparison would serve as a characterization approach and pre-screening of plausible EMIC formulations. This type of study will be experimentally inconclusive, due to the instability associated with the TS. However, this useful exercise will be setting a theoretical limit to help identify the best TSA for a specified system. The presented computational techniques allow us to fulfil both characterization processes (i.e. the structural and the energetic ones) using an entirely controlled framework.
MICs are easy to prepare and handle, and the EMICs will inherit these qualities. MICs can be prepared in large quantities by suspension polymerization, and stable particles of uniform diameter can be easily obtained . In addition to beads or broken particles, MICs can also be prepared in other very different forms, such as monoliths, microcapsules, membranes or surfaces . MICs have both excellent mechanical and thermal stability. Frequently, they can be used for a long time in a continuous process, or they can be reused many times. As a result of their insolubility, they can be easily filtered off after a reaction, or they can be placed in a flow reactor. Whereas enzymes and antibodies degrade under harsh conditions such as high temperature, chemically aggressive media, and high and low pH, MICs show better behaviour in most cases, and they can be applied directly in chemical processes since they are rather stable materials.
The authors researched and wrote this review to the same extent. All authors read and approved the final manuscript.
The authors thank Prof. Paul R. Raithby and Dr Julian Rose from the EPSRC Directed Assembly Network for the invitation to write this contribution.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Cox MM, Doudna JA, O’Donnell M (2012) Molecular biology principle and practice, 1st edn. W.H. Freeman and Company, New YorkGoogle Scholar
- Bhat SV, Nagasampagi BA, Sivakumar M (2006) Chemistry of natural products, 2nd edn. Springer Narosa, ChennaiGoogle Scholar
- Sekhon BS (2013) Exploiting the power of stereochemistry in drugs: an overview of racemic and enantiopure drugs. J Mod Med Chem 1:10–36Google Scholar
- Somogyi A, Bochner F, Foster D (2004) Inside the isomers: the tale of chiral switches. Aust Prescr 27:47–49View ArticleGoogle Scholar
- Tucker GT (2000) Chiral switches. Lancet 355:1085–1087View ArticleGoogle Scholar
- Nguyen LA, He H, Pham-Huy C (2006) Chiral drugs: an overview. Int J Biomed Sci 2:85–100Google Scholar
- Jimenez-Lucho VE, del Busto R, Odel J (1987) Isoniazid and ethambutol as a cause of optic neuropathy. Eur J Respir Dis 71:42–45Google Scholar
- Baillie TA, Adams WJ, Kaiser DG, Olanoff LS, Halstead GW, Harpootlian H, Van Giessen G (1988) Use of deuterium labelling in mechanistic studies of the metabolic chiral inversion of ibuprofen. Int J Radiat Appl Instrum Part A Appl Radiat Isot 39(6):548View ArticleGoogle Scholar
- Arroniz C, Escolano C (2012) 7. Strategies for the synthesis of enantiopure compounds focused on organocatalysis. In: Munoz-Torrero D, Haro D, Valles J (eds) Recent advances in pharmaceutical sciences II. Transworld Research Network, Kerala, pp 116–119Google Scholar
- Sellergren B (2003) Molecularly imprinted polymers man-made mimics of antibodies and their application in analytical chemistry, 2nd edn. Elsevier, AmsterdamGoogle Scholar
- Schwyzer R (1970) Organization and read-out of biological information in polypeptides. Proc Fourth Int Congr Pharmacol 5:196–209Google Scholar
- Wulff G, Sarhan A, Zabrocki K (1973) Enzyme-analogue built polymers and their use for the resolution of racemates. Tetrahedron Lett 44:4329–4332View ArticleGoogle Scholar
- Dickey FH (1949) The preparation of specific adsorbents. Proc Natl Acad Sci 35(5):227–229View ArticleGoogle Scholar
- Wulff G (1995) Molecular imprinting in cross-linked materials with the aid of molecular templates—a way towards artificial antibodies. Angew Chem Int Ed Engl 34(17):1812–1832View ArticleGoogle Scholar
- Mosbach K, Ramstrom O (1996) The emerging technique if molecular imprinting and its future impact on biotechnology. Biotechnology 14(2):163–170View ArticleGoogle Scholar
- Pauling L (1946) Molecular architecture and biological reactions. Chem Eng News 24:1375–1377View ArticleGoogle Scholar
- Jencks WP (1987) Catalysis in chemistry and enzymology, 1st edn. Dover, New YorkGoogle Scholar
- Lerner RA, Benkovic SJ, Schulz PG (1991) At the crossroads of chemistry and immunology: catalytic antibodies. Science 252:659–667View ArticleGoogle Scholar
- Schultz PG (1989) Antikörper als Katalysatoren. Angew Chem 101:1336–1348View ArticleGoogle Scholar
- Li S, Cao S, Piletsky SA, Turner APF (2015) Molecularly imprinted catalysts: principles, syntheses, and applications, 1st edn. Elsevier, OxfordGoogle Scholar
- Davis ME (1997) Catalytic materials via molecular imprinting. Cattech 1:19–26Google Scholar
- Wulff G (2002) Enzyme-like catalysis by molecularly imprinted polymers. Chem Rev 102(1):1–28View ArticleGoogle Scholar
- Brady PA, Sanders JKM (1997) Selection approaches to catalytic systems. Chem Soc Rev 26:327–336View ArticleGoogle Scholar
- Vidyasankar S, Arnold FH (1995) Molecular imprinting: selective materials for separations, sensors and catalysis. Curr Opin Biotechnol 6:218–224View ArticleGoogle Scholar
- Shea KJ (1994) Molecular imprinting of synthetic network polymers: the de novo synthesis of macromolecular binding and catalytic sites. Trends Polym Sci 2:166–173Google Scholar
- Takeuchi T, Matsui J (1996) Molecular imprinting: an approach to “tailor-made” synthetic polymers with biomimetic functions. Acta Polym 47:471–480View ArticleGoogle Scholar
- Alexander C, Smith CR, Whitcombe MJ, Vulfson EN (1999) Imprinted polymers as protecting groups for regioselective modification of polyfunctional substrates. J Am Chem Soc 121:6640–6651View ArticleGoogle Scholar
- Strikovsky AG, Kasper D, Grün M, Green BS, Hradil J, Wulff G (2000) Catalytic molecularly imprinted polymers using conventional bulk polymerization or suspension polymerization: selective hydrolysis of diphenyl carbonate and diphenyl carbamate. J Am Chem Soc 122:6295–6296View ArticleGoogle Scholar
- Herdes C, Sarkisov L (2009) Computer simulation of volatile organic compound adsorption in atomistic models of molecularly imprinted polymers. Langmuir 25(9):5352–5359View ArticleGoogle Scholar
- Chianella I, Lotierzo M, Piletsky SA, Tothill IE, Chen BN, Karim K, Turner APF (2002) Rational design of a polymer specific for microcystin-LR using a computational approach. Anal Chem 74:1288–1293View ArticleGoogle Scholar
- Karim K, Breton F, Rouillon R, Piletska EV, Guerreiro A, Chianella I, Piletsky SA (2005) How to find effective functional monomers for effective molecularly imprinted polymers? Adv Drug Deliv Rev 57:1795–1808View ArticleGoogle Scholar
- Chianella I, Karim K, Piletska EV, Preston C, Piletsky SA (2006) Computational design and synthesis of molecularly imprinted polymers with high binding capacity for pharmaceutical applications-model case: adsorbent for abacavir. Anal Chim Acta 559:73–78View ArticleGoogle Scholar
- Pavel D, Lagowski J (2005) Computationally designed monomers and polymers for molecular imprinting of theophylline and its derivatives. Part I. Polymer 46:7528–7542View ArticleGoogle Scholar
- Pavel D, Lagowski J, Lepage CJ (2006) Computationally designed monomers for molecular imprinting of chemical warfare agents—Part V. Polymer 47:8389–8399View ArticleGoogle Scholar
- Curk T, Dobnikar J, Frenkel D (2016) Rational design of molecularly imprinted polymers. Soft Matter 12:35–44View ArticleGoogle Scholar
- Schumacher C, Gonzalez J, Wright PA, Seaton NA (2006) Generation of atomistic models of periodic mesoporous silica by kinetic Monte Carlo simulation of the synthesis of the material. J Phys Chem B 110:319–333View ArticleGoogle Scholar
- Adams DJ (1974) Chemical potential of hard-sphere fluids by Monte Carlo methods. Mol Phys 28(5):1241–1252View ArticleGoogle Scholar
- Resat H, Wiley HS, Dixon DA (2001) Probability-weighted dynamic Monte Carlo method for reaction kinetics simulations. J Phys Chem B 105:11026–11034View ArticleGoogle Scholar
- Siepmann JI, Frenkel D (1992) Configurational bias Monte Carlo: a new sampling scheme for flexible chains. Mol Phys 75(1):59–70View ArticleGoogle Scholar
- Brennan JK (2005) Cavity-bias sampling in reaction ensemble Monte Carlo simulations. Mol Phys 103(19):2647–2654View ArticleGoogle Scholar