1
|
Reinhardt CR, Manetsch MT, Li WL, Román-Leshkov Y, Head-Gordon T, Kulik HJ. Computational Screening of Putative Catalyst Transition Metal Complexes as Guests in a Ga 4L 612- Nanocage. Inorg Chem 2024; 63:14609-14622. [PMID: 39049593 DOI: 10.1021/acs.inorgchem.4c02113] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/27/2024]
Abstract
Metal-organic cages form well-defined microenvironments that can enhance the catalytic proficiency of encapsulated transition metal complexes (TMCs). We introduce a screening protocol to efficiently identify TMCs that are promising candidates for encapsulation in the Ga4L612- nanocage. We obtain TMCs from the Cambridge Structural Database with geometric and electronic characteristics amenable to encapsulation and mine the text of associated manuscripts to curate TMCs with documented catalytic functionality. By docking candidate TMCs inside the nanocage cavity and carrying out electronic structure calculations, we identify a subset of successfully optimized candidates (TMC-34) and observe that encapsulated guests occupy an average of 60% of the cavity volume, in line with previous observations. Notably, some guests occupy as much as 72% of the cavity as a result of linker rotation. Encapsulation has a universal effect on the electrostatic potential (ESP), systematically decreasing the ESP at the metal center of each TMC in the TMC-34 data set, while minimally altering TMC metal partial charges. Collectively these observations support geometry-based screening of potential guests and suggest that encapsulation in Ga4L612- cages could electrostatically stabilize diverse cationic or electropositive intermediates. We highlight candidate guests with associated known reactivity and solubility most amenable for encapsulation in experimental follow-up studies.
Collapse
Affiliation(s)
- Clorice R Reinhardt
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Melissa T Manetsch
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Wan-Lu Li
- Kenneth S. Pitzer Center for Theoretical Chemistry, Berkeley, California 94720, United States
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Yuriy Román-Leshkov
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Teresa Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Berkeley, California 94720, United States
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, United States
- Department of Bioengineering, University of California, Berkeley, California 94720, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
2
|
Boaler P, Piskorz TK, Bickerton LE, Wang J, Duarte F, Lloyd-Jones GC, Lusby PJ. Origins of High-Activity Cage-Catalyzed Michael Addition. J Am Chem Soc 2024; 146:19317-19326. [PMID: 38976816 PMCID: PMC11258793 DOI: 10.1021/jacs.4c05160] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2024] [Revised: 05/23/2024] [Accepted: 05/24/2024] [Indexed: 07/10/2024]
Abstract
Cage catalysis continues to create significant interest, yet catalyst function remains poorly understood. Herein, we report mechanistic insights into coordination-cage-catalyzed Michael addition using kinetic and computational methods. The study has been enabled by the detection of identifiable catalyst intermediates, which allow the evolution of different cage species to be monitored and modeled alongside reactants and products. The investigations show that the overall acceleration results from two distinct effects. First, the cage reaction shows a thousand-fold increase in the rate constant for the turnover-limiting C-C bond-forming step compared to a reference state. Computational modeling and experimental analysis of activation parameters indicate that this stems from a significant reduction in entropy, suggesting substrate coencapsulation. Second, the cage markedly acidifies the bound pronucleophile, shifting this equilibrium by up to 6 orders of magnitude. The combination of these two factors results in accelerations up to 109 relative to bulk-phase reference reactions. We also show that the catalyst can fundamentally alter the reaction mechanism, leading to intermediates and products that are not observable outside of the cage. Collectively, the results show that cage catalysis can proceed with very high activity and unique selectivity by harnessing a series of individually weak noncovalent interactions.
Collapse
Affiliation(s)
- Patrick
J. Boaler
- EaStCHEM
School of Chemistry, University of Edinburgh, Joseph Black Building, David Brewster
Road, Edinburgh, Scotland EH9 3FJ, U.K.
| | - Tomasz K. Piskorz
- Chemistry
Research Laboratory, University of Oxford, Oxford OX1 3TA, U.K.
| | - Laura E. Bickerton
- EaStCHEM
School of Chemistry, University of Edinburgh, Joseph Black Building, David Brewster
Road, Edinburgh, Scotland EH9 3FJ, U.K.
| | - Jianzhu Wang
- EaStCHEM
School of Chemistry, University of Edinburgh, Joseph Black Building, David Brewster
Road, Edinburgh, Scotland EH9 3FJ, U.K.
| | - Fernanda Duarte
- Chemistry
Research Laboratory, University of Oxford, Oxford OX1 3TA, U.K.
| | - Guy C. Lloyd-Jones
- EaStCHEM
School of Chemistry, University of Edinburgh, Joseph Black Building, David Brewster
Road, Edinburgh, Scotland EH9 3FJ, U.K.
| | - Paul J. Lusby
- EaStCHEM
School of Chemistry, University of Edinburgh, Joseph Black Building, David Brewster
Road, Edinburgh, Scotland EH9 3FJ, U.K.
| |
Collapse
|
3
|
Piskorz TK, Martí-Centelles V, Spicer RL, Duarte F, Lusby PJ. Picking the lock of coordination cage catalysis. Chem Sci 2023; 14:11300-11331. [PMID: 37886081 PMCID: PMC10599471 DOI: 10.1039/d3sc02586a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2023] [Accepted: 08/29/2023] [Indexed: 10/28/2023] Open
Abstract
The design principles of metallo-organic assembly reactions have facilitated access to hundreds of coordination cages of varying size and shape. Many of these assemblies possess a well-defined cavity capable of hosting a guest, pictorially mimicking the action of a substrate binding to the active site of an enzyme. While there are now a growing collection of coordination cages that show highly proficient catalysis, exhibiting both excellent activity and efficient turnover, this number is still small compared to the vast library of metal-organic structures that are known. In this review, we will attempt to unpick and discuss the key features that make an effective coordination cage catalyst, linking structure to activity (and selectivity) using lessons learnt from both experimental and computational analysis of the most notable exemplars. We will also provide an outlook for this area, reasoning why coordination cages have the potential to become the gold-standard in (synthetic) non-covalent catalysis.
Collapse
Affiliation(s)
- Tomasz K Piskorz
- Chemistry Research Laboratory, University of Oxford Oxford OX1 3TA UK
| | - Vicente Martí-Centelles
- Instituto Interuniversitario de Investigación de Reconocimiento Molecular y Desarrollo Tecnológico (IDM), Universitat Politècnica de València, Universitat de València Camino de Vera, s/n 46022 Valencia Spain
| | - Rebecca L Spicer
- Department of Chemistry, Lancaster University Lancaster LA14YB UK
| | - Fernanda Duarte
- Chemistry Research Laboratory, University of Oxford Oxford OX1 3TA UK
| | - Paul J Lusby
- EaStCHEM School of Chemistry, University of Edinburgh Edinburgh Scotland EH9 3FJ UK
| |
Collapse
|
4
|
Blyth MT, Coote ML. Manipulation of N-heterocyclic carbene reactivity with practical oriented electric fields. Phys Chem Chem Phys 2023; 25:375-383. [DOI: 10.1039/d2cp04507a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Electric fields can be used to tune the nucleophilicity and electrophilicity of N-heterocyclic carbenes and enhance their catalytic activity.
Collapse
Affiliation(s)
- Mitchell T. Blyth
- Research School of Chemistry, Australian National University, Canberra, Australian Capital Territory 2601, Australia
| | - Michelle L. Coote
- Institute for Nanoscale Science and Technology, College of Science and Engineering, Flinders University, Bedford Park, South Australia 5042, Australia
| |
Collapse
|
5
|
Norjmaa G, Himo F, Maréchal J, Ujaque G. Catalysis by [Ga 4 L 6 ] 12- Metallocage on the Nazarov Cyclization: The Basicity of Complexed Alcohol is Key. Chemistry 2022; 28:e202201792. [PMID: 35859038 PMCID: PMC9804567 DOI: 10.1002/chem.202201792] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2022] [Indexed: 01/05/2023]
Abstract
The Nazarov cyclization is investigated in solution and within K12 [Ga4 L6 ] supramolecular organometallic cage by means of computational methods. The reaction needs acidic condition in solution but works at neutral pH in the presence of the metallocage. The reaction steps for the process are analogous in both media: (a) protonation of the alcohol group, (b) water loss and (c) cyclization. The relative Gibbs energies of all the steps are affected by changing the environment from solvent to the metallocage. The first step in the mechanism, the alcohol protonation, turns out to be the most critical one for the acceleration of the reaction inside the metallocage. In order to calculate the relative stability of protonated alcohol inside the cavity, we propose a computational scheme for the calculation of basicity for species inside cavities and can be of general use. These results are in excellent agreement with the experiments, identifying key steps of catalysis and providing an in-depth understanding of the impact of the metallocage on all the reaction steps.
Collapse
Affiliation(s)
- Gantulga Norjmaa
- Departament de Química and Centro de Innovación en Química Avanzada (ORFEO-CINQA)Universitat Autònoma de Barcelona08193Cerdanyola del VallesBarcelona, CataloniaSpain
| | - Fahmi Himo
- Department of Organic ChemistryArrhenius LaboratoryStockholm University10691StockholmSweden
| | - Jean‐Didier Maréchal
- Departament de Química and Centro de Innovación en Química Avanzada (ORFEO-CINQA)Universitat Autònoma de Barcelona08193Cerdanyola del VallesBarcelona, CataloniaSpain
| | - Gregori Ujaque
- Departament de Química and Centro de Innovación en Química Avanzada (ORFEO-CINQA)Universitat Autònoma de Barcelona08193Cerdanyola del VallesBarcelona, CataloniaSpain
| |
Collapse
|
6
|
Nguyen QNN, Xia KT, Zhang Y, Chen N, Morimoto M, Pei X, Ha Y, Guo J, Yang W, Wang LP, Bergman RG, Raymond KN, Toste FD, Tantillo DJ. Source of Rate Acceleration for Carbocation Cyclization in Biomimetic Supramolecular Cages. J Am Chem Soc 2022; 144:11413-11424. [PMID: 35699585 DOI: 10.1021/jacs.2c04179] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
The results of quantum chemical and molecular dynamics calculations reveal that polyanionic gallium-based cages accelerate cyclization reactions of pentadienyl alcohols as a result of substrate cage interactions, preferential binding of reactive conformations of substrate/H3O+ pairs, and increased substrate basicity. However, the increase in basicity dominates. Experimental structure-activity relationship studies in which the metal vertices and overall charge of the cage are varied confirm the model derived via calculations.
Collapse
Affiliation(s)
- Quynh Nhu N Nguyen
- Department of Chemistry, University of California, Davis, 1 Shields Avenue, Davis, California 95616, United States
| | - Kay T Xia
- Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Yue Zhang
- Department of Chemistry, University of California, Davis, 1 Shields Avenue, Davis, California 95616, United States
| | - Nanhao Chen
- Department of Chemistry, University of California, Davis, 1 Shields Avenue, Davis, California 95616, United States
| | - Mariko Morimoto
- Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Xiaokun Pei
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Yang Ha
- Advanced Light Source, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Jinghua Guo
- Advanced Light Source, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Wanli Yang
- Advanced Light Source, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Lee-Ping Wang
- Department of Chemistry, University of California, Davis, 1 Shields Avenue, Davis, California 95616, United States
| | - Robert G Bergman
- Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Kenneth N Raymond
- Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - F Dean Toste
- Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Dean J Tantillo
- Department of Chemistry, University of California, Davis, 1 Shields Avenue, Davis, California 95616, United States
| |
Collapse
|
7
|
Piskorz TK, Martí-Centelles V, Young TA, Lusby PJ, Duarte F. Computational Modeling of Supramolecular Metallo-organic Cages-Challenges and Opportunities. ACS Catal 2022; 12:5806-5826. [PMID: 35633896 PMCID: PMC9127791 DOI: 10.1021/acscatal.2c00837] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2022] [Revised: 04/09/2022] [Indexed: 01/18/2023]
Abstract
![]()
Self-assembled
metallo-organic
cages have emerged as promising
biomimetic platforms that can encapsulate whole substrates akin to
an enzyme active site. Extensive experimental work has enabled access
to a variety of structures, with a few notable examples showing catalytic
behavior. However, computational investigations of metallo-organic
cages are scarce, not least due to the challenges associated with
their modeling and the lack of accurate and efficient protocols to
evaluate these systems. In this review, we discuss key molecular principles
governing the design of functional metallo-organic cages, from the
assembly of building blocks through binding and catalysis. For each
of these processes, computational protocols will be reviewed, considering
their inherent strengths and weaknesses. We will demonstrate that
while each approach may have its own specific pitfalls, they can be
a powerful tool for rationalizing experimental observables and to
guide synthetic efforts. To illustrate this point, we present several
examples where modeling has helped to elucidate fundamental principles
behind molecular recognition and reactivity. We highlight the importance
of combining computational and experimental efforts to speed up supramolecular
catalyst design while reducing time and resources.
Collapse
Affiliation(s)
- Tomasz K. Piskorz
- Chemistry Research Laboratory, University of Oxford, Mansfield Road, Oxford OX1 3TA, United Kingdom
| | - Vicente Martí-Centelles
- Instituto Interuniversitario de Investigación de Reconocimiento Molecular y Desarrollo Tecnológico (IDM), Universitat Politècnica de València, Universitat de València, Valencia 46022, Spain
| | - Tom A. Young
- Chemistry Research Laboratory, University of Oxford, Mansfield Road, Oxford OX1 3TA, United Kingdom
| | - Paul J. Lusby
- EaStCHEM School of Chemistry, University of Edinburgh, Joseph Black Building, David Brewster Road, Edinburgh, Scotland EH9 3FJ, United Kingdom
| | - Fernanda Duarte
- Chemistry Research Laboratory, University of Oxford, Mansfield Road, Oxford OX1 3TA, United Kingdom
| |
Collapse
|
8
|
Tarzia A, Jelfs KE. Unlocking the computational design of metal-organic cages. Chem Commun (Camb) 2022; 58:3717-3730. [PMID: 35229861 PMCID: PMC8932387 DOI: 10.1039/d2cc00532h] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2022] [Accepted: 02/22/2022] [Indexed: 12/11/2022]
Abstract
Metal-organic cages are macrocyclic structures that can possess an intrinsic void that can hold molecules for encapsulation, adsorption, sensing, and catalysis applications. As metal-organic cages may be comprised from nearly any combination of organic and metal-containing components, cages can form with diverse shapes and sizes, allowing for tuning toward targeted properties. Therefore, their near-infinite design space is almost impossible to explore through experimentation alone and computational design can play a crucial role in exploring new systems. Although high-throughput computational design and screening workflows have long been known as powerful tools in drug and materials discovery, their application in exploring metal-organic cages is more recent. We show examples of structure prediction and host-guest/catalytic property evaluation of metal-organic cages. These examples are facilitated by advances in methods that handle metal-containing systems with improved accuracy and are the beginning of the development of automated cage design workflows. We finally outline a scope for how high-throughput computational methods can assist and drive experimental decisions as the field pushes toward functional and complex metal-organic cages. In particular, we highlight the importance of considering realistic, flexible systems.
Collapse
Affiliation(s)
- Andrew Tarzia
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, White City Campus, Wood Lane, London, W12 0BZ, UK.
| | - Kim E Jelfs
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, White City Campus, Wood Lane, London, W12 0BZ, UK.
| |
Collapse
|
9
|
Norjmaa G, Maréchal JD, Ujaque G. Origin of the Rate Acceleration in the C-C Reductive Elimination from Pt(IV)-complex in a [Ga 4 L 6 ] 12- Supramolecular Metallocage. Chemistry 2021; 27:15973-15980. [PMID: 34545974 PMCID: PMC9293218 DOI: 10.1002/chem.202102250] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2021] [Indexed: 12/29/2022]
Abstract
The reductive elimination on [(Me3 P)2 Pt(MeOH)(CH3 )3 ]+ , 2P, complex performed in MeOH solution and inside a [Ga4 L6 ]12- metallocage are computationally analysed by mean of QM and MD simulations and compared with the mechanism of gold parent systems previously reported [Et3 PAu(MeOH)(CH3 )2 ]+ , 2Au. The comparative analysis between the encapsulated Au(III) and Pt(IV)-counterparts shows that there are no additional solvent MeOH molecules inside the cavity of the metallocage for both systems. The Gibbs energy barriers for the 2P reductive elimination calculated at DFT level are in good agreement with the experimental values for both environments. The effect of microsolvation and encapsulation on the rate acceleration are evaluated and shows that the latter is far more relevant, conversely to 2Au. Energy decomposition analysis indicates that the encapsulation is the main responsible for most of the energy barrier reduction. Microsolvation and encapsulation effects are not equally contributing for both metal systems and consequently, the reasons of the rate acceleration are not the same for both metallic systems despite the similarity between them.
Collapse
Affiliation(s)
- Gantulga Norjmaa
- Departament de Química and Centro de Innovación en Química Avanzada (ORFEO-CINQA), Universitat Autònoma de Barcelona, 08193, Cerdanyola del Valles, Barcelona, Catalonia, Spain
| | - Jean-Didier Maréchal
- Departament de Química and Centro de Innovación en Química Avanzada (ORFEO-CINQA), Universitat Autònoma de Barcelona, 08193, Cerdanyola del Valles, Barcelona, Catalonia, Spain
| | - Gregori Ujaque
- Departament de Química and Centro de Innovación en Química Avanzada (ORFEO-CINQA), Universitat Autònoma de Barcelona, 08193, Cerdanyola del Valles, Barcelona, Catalonia, Spain
| |
Collapse
|
10
|
Norjmaa G, Vidossich P, Maréchal JD, Ujaque G. Modeling Kinetics and Thermodynamics of Guest Encapsulation into the [M 4L 6] 12- Supramolecular Organometallic Cage. J Chem Inf Model 2021; 61:4370-4381. [PMID: 34505774 PMCID: PMC8479806 DOI: 10.1021/acs.jcim.1c00348] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
The encapsulation
of molecular guests into supramolecular hosts
is a complex molecular recognition process in which the guest displaces
the solvent from the host cavity, while the host deforms to let the
guest in. An atomistic description of the association would provide
valuable insights on the physicochemical properties that guide it.
This understanding may be used to design novel host assemblies with
improved properties (i.e., affinities) toward a given class of guests.
Molecular simulations may be conveniently used to model the association
processes. It is thus of interest to establish efficient protocols
to trace the encapsulation process and to predict the associated magnitudes
ΔGbind and ΔGbind⧧. Here, we report the calculation of the Gibbs energy barrier and
Gibbs binding energy by means of explicit solvent molecular simulations
for the [Ga4L6]12– metallocage
encapsulating a series of cationic molecules. The ΔGbind⧧ for encapsulation was estimated by means of umbrella sampling simulations.
The steps involved were identified, including ion-pair formation and
naphthalene rotation (from L ligands of the metallocage) during the
guest’s entrance. The ΔGbind values were computed using the attach–pull–release
method. The results reveal the sensitivity of the estimates on the
force field parameters, in particular on atomic charges, showing that
higher accuracy is obtained when charges are derived from implicit
solvent quantum chemical calculations. Correlation analysis identified
some indicators for the binding affinity trends. All computed magnitudes
are in very good agreement with experimental observations. This work
provides, on one side, a benchmarked way to computationally model
a highly charged metallocage encapsulation process. This includes
a nonstandard parameterization and charge derivation procedure. On
the other hand, it gives specific mechanistic information on the binding
processes of [Ga4L6]12– at
the molecular level where key motions are depicted. Taken together,
the study provides an interesting option for the future design of
metal–organic cages.
Collapse
Affiliation(s)
- Gantulga Norjmaa
- Departament de Química and Centro de Innovación en Química Avanzada (ORFEO-CINQA), Universitat Autònoma de Barcelona, Cerdanyola del Valles, 08193 Cerdanyola del Vallès, Catalonia, Spain
| | - Pietro Vidossich
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy
| | - Jean-Didier Maréchal
- Departament de Química and Centro de Innovación en Química Avanzada (ORFEO-CINQA), Universitat Autònoma de Barcelona, Cerdanyola del Valles, 08193 Cerdanyola del Vallès, Catalonia, Spain
| | - Gregori Ujaque
- Departament de Química and Centro de Innovación en Química Avanzada (ORFEO-CINQA), Universitat Autònoma de Barcelona, Cerdanyola del Valles, 08193 Cerdanyola del Vallès, Catalonia, Spain
| |
Collapse
|
11
|
Léonard NG, Dhaoui R, Chantarojsiri T, Yang JY. Electric Fields in Catalysis: From Enzymes to Molecular Catalysts. ACS Catal 2021; 11:10923-10932. [DOI: 10.1021/acscatal.1c02084] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Nadia G. Léonard
- Department of Chemistry, University of California, Irvine, California 92697, United States
| | - Rakia Dhaoui
- Department of Chemistry, University of California, Irvine, California 92697, United States
| | - Teera Chantarojsiri
- Department of Chemistry and Center of Excellence for Innovation in Chemistry, Faculty of Science, Mahidol University, Bangkok, 10400, Thailand
| | - Jenny Y. Yang
- Department of Chemistry, University of California, Irvine, California 92697, United States
| |
Collapse
|
12
|
Daver H, Rebek J, Himo F. Modeling the Reaction of Carboxylic Acids and Isonitriles in a Self-Assembled Capsule. Chemistry 2020; 26:10861-10870. [PMID: 32428333 PMCID: PMC7522688 DOI: 10.1002/chem.202001735] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2020] [Revised: 05/11/2020] [Indexed: 11/10/2022]
Abstract
Quantum chemical calculations were used to study the reaction of carboxylic acids with isonitriles inside a resorcinarene-based self-assembled capsule. Experimentally, it has been shown that the reactions between p-tolylacetic acid and n-butyl isonitrile or isopropyl isonitrile behave differently in the presence of the capsule compared both with each other and also with their solution counterparts. Herein, the reasons for these divergent behaviors are addressed by comparing the detailed energy profiles for the reactions of the two isonitriles inside and outside the capsule. An energy decomposition analysis was conducted to quantify the different factors affecting the reactivity. The calculations reproduce the experimental findings very well. Thus, encapsulation leads to lowering of the energy barrier for the first step of the reaction, the concerted α-addition and proton transfer, which in solution is rate-determining, and this explains the rate acceleration observed in the presence of the capsule. The barrier for the final step of the reaction, the 1,3 O→N acyl transfer, is calculated to be higher with the isopropyl substituent inside the capsule compared with n-butyl. With the isopropyl substituent, the transition state and the product of this step are significantly shorter than the preceding intermediate, and this results in energetically unfavorable empty spaces inside the capsule, which cause a higher barrier. With the n-butyl substituent, on the other hand, the carbon chain can untwine and hence uphold an appropriate guest length.
Collapse
Affiliation(s)
- Henrik Daver
- Department of Organic ChemistryArrhenius LaboratoryStockholm University106 91StockholmSweden
- Present address: Department of Drug Design and PharmacologyUniversity of CopenhagenUniversitetsparken 22100CopenhagenDenmark
| | - Julius Rebek
- The Skaggs Institute for Chemical Biology and Department of ChemistryThe Scripps Research Institute10550 North Torrey Pines RoadLa JollaCalifornia92037USA
- Center for Supramolecular Chemistry and CatalysisShanghai UniversityShanghai200444P.R. China
| | - Fahmi Himo
- Department of Organic ChemistryArrhenius LaboratoryStockholm University106 91StockholmSweden
| |
Collapse
|
13
|
Norjmaa G, Maréchal J, Ujaque G. Reaction Rate Inside the Cavity of [Ga
4
L
6
]
12−
Supramolecular Metallocage is Regulated by the Encapsulated Solvent. Chemistry 2020; 26:6988-6992. [PMID: 32125031 DOI: 10.1002/chem.201905608] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2019] [Indexed: 11/05/2022]
Affiliation(s)
- Gantulga Norjmaa
- Departament de Química Universitat Autònoma de Barcelona 08193 Cerdanyola del Valles, Barcelona Catalonia, and Centro de Innovaciónen Química Avanzada (ORFEO-CINQA Spain
| | - Jean‐Didier Maréchal
- Departament de Química Universitat Autònoma de Barcelona 08193 Cerdanyola del Valles, Barcelona Catalonia, and Centro de Innovaciónen Química Avanzada (ORFEO-CINQA Spain
| | - Gregori Ujaque
- Departament de Química Universitat Autònoma de Barcelona 08193 Cerdanyola del Valles, Barcelona Catalonia, and Centro de Innovaciónen Química Avanzada (ORFEO-CINQA Spain
| |
Collapse
|
14
|
Nijamudheen A, Datta A. Gold-Catalyzed Cross-Coupling Reactions: An Overview of Design Strategies, Mechanistic Studies, and Applications. Chemistry 2019; 26:1442-1487. [PMID: 31657487 DOI: 10.1002/chem.201903377] [Citation(s) in RCA: 101] [Impact Index Per Article: 20.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2019] [Revised: 10/28/2019] [Indexed: 12/14/2022]
Abstract
Transition-metal-catalyzed cross-coupling reactions are central to many organic synthesis methodologies. Traditionally, Pd, Ni, Cu, and Fe catalysts are used to promote these reactions. Recently, many studies have showed that both homogeneous and heterogeneous Au catalysts can be used for activating selective cross-coupling reactions. Here, an overview of the past studies, current trends, and future directions in the field of gold-catalyzed coupling reactions is presented. Design strategies to accomplish selective homocoupling and cross-coupling reactions under both homogeneous and heterogeneous conditions, computational and experimental mechanistic studies, and their applications in diverse fields are critically reviewed. Specific topics covered are: oxidant-assisted and oxidant-free reactions; strain-assisted reactions; dual Au and photoredox catalysis; bimetallic synergistic reactions; mechanisms of reductive elimination processes; enzyme-mimicking Au chemistry; cluster and surface reactions; and plasmonic catalysis. In the relevant sections, theoretical and computational studies of AuI /AuIII chemistry are discussed and the predictions from the calculations are compared with the experimental observations to derive useful design strategies.
Collapse
Affiliation(s)
- A Nijamudheen
- School of Chemical Sciences, Indian Association for the, Cultivation of Sciences, 2A & 2B Raja S C Mullick Road, Kolkata, 700032, India.,Department of Chemical & Biomedical Engineering, Florida A&M University-Florida State University, Joint College of Engineering, 2525 Pottsdamer Street, Tallahassee, FL, 32310, USA
| | - Ayan Datta
- School of Chemical Sciences, Indian Association for the, Cultivation of Sciences, 2A & 2B Raja S C Mullick Road, Kolkata, 700032, India
| |
Collapse
|
15
|
Norjmaa G, Maréchal JD, Ujaque G. Microsolvation and Encapsulation Effects on Supramolecular Catalysis: C-C Reductive Elimination inside [Ga 4L 6] 12- Metallocage. J Am Chem Soc 2019; 141:13114-13123. [PMID: 31390202 DOI: 10.1021/jacs.9b04909] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The host effect of the supramolecular [Ga4L6]12- tetrahedral metallocage on reductive elimination of substrate by encapsulated Au(III) complexes is investigated by means of computational methods. The behavior of the reactants in solution and within the metallocage is initially evaluated by means of classical molecular dynamics simulations. These results guided the selection of proper computational models to describe the reaction in solution and inside the metallocage at the DFT level. The calculated Gibbs energy barriers are in very good agreement with experiment both in solution and inside the metallocage. The analysis in solution revealed that microsolvation around the Au(III) complex increases the Gibbs energy barrier. The analysis within the metallocage shows that its encapsulation favors the reaction. The process can be formally described as removing explicit microsolvation around the gold complex and encapsulating the metal complex inside the metallocage. Both processes are important for the reaction, but the removal of the solvent molecules surrounding the Au(III) metal complex is fundamental for the reduction of the reaction barrier. The energy decomposition analysis of the barrier among strain, interaction, and thermal terms shows that strain term is very low whereas the contribution of thermal (entropic) effects is moderate. Interestingly, the key term responsible for reducing the Gibbs energy barrier is the interaction. This term can be mainly associated with electrostatic interactions in agreement with previous examples in the literature.
Collapse
Affiliation(s)
- Gantulga Norjmaa
- Departament de Química and Centro de Innovación en Química Avanzada (ORFEO-CINQA) , Universitat Autònoma de Barcelona , Cerdanyola del Valles , 08193 Barcelona , Catalonia , Spain
| | - Jean-Didier Maréchal
- Departament de Química and Centro de Innovación en Química Avanzada (ORFEO-CINQA) , Universitat Autònoma de Barcelona , Cerdanyola del Valles , 08193 Barcelona , Catalonia , Spain
| | - Gregori Ujaque
- Departament de Química and Centro de Innovación en Química Avanzada (ORFEO-CINQA) , Universitat Autònoma de Barcelona , Cerdanyola del Valles , 08193 Barcelona , Catalonia , Spain
| |
Collapse
|
16
|
Nothling MD, Xiao Z, Bhaskaran A, Blyth MT, Bennett CW, Coote ML, Connal LA. Synthetic Catalysts Inspired by Hydrolytic Enzymes. ACS Catal 2018. [DOI: 10.1021/acscatal.8b03326] [Citation(s) in RCA: 61] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Mitchell D. Nothling
- Department of Chemical and Biomolecular Engineering, The University of Melbourne, Parkville, Victoria 3010, Australia
| | - Zeyun Xiao
- Chongqing Institute of Green and Intelligent Technology, Chinese Academy of Sciences, Chongqing 400714, P. R. China
| | - Ayana Bhaskaran
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia
| | - Mitchell T. Blyth
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia
| | - Christopher W. Bennett
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia
| | - Michelle L. Coote
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia
| | - Luke A. Connal
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia
| |
Collapse
|
17
|
|
18
|
Vaissier Welborn V, Head-Gordon T. Electrostatics Generated by a Supramolecular Capsule Stabilizes the Transition State for Carbon-Carbon Reductive Elimination from Gold(III) Complex. J Phys Chem Lett 2018; 9:3814-3818. [PMID: 29939756 DOI: 10.1021/acs.jpclett.8b01710] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Tetrahedral assemblies of stoichiometry M4L6 have been proven to catalyze a range of chemical reactions including the carbon-carbon reductive elimination reaction from transition metals such as gold. Here, we perform quantum chemical calculations of Gold(III) transition metal complexes in vacuum, and encapsulated in Ga4L612- or Si4L68- assemblies within both a reaction field continuum solvent and in an aqueous molecular environment with counterions, to rationalize the rate enhancements observed experimentally for the reductive elimination reaction. We find that the Ga4L612- assembly lowers the energy barrier of the reaction compared to Si4L68-, which is consistent with kinetic trends observed experimentally. We have determined that the primary factor for catalytic rate acceleration stems from the electrostatic environment emanating from the Ga4L612- capsule as opposed to the water or counterions.
Collapse
Affiliation(s)
- Valerie Vaissier Welborn
- Chemical Sciences Division , Lawrence Berkeley National Laboratories , Berkeley , California 94720 , United States
| | - Teresa Head-Gordon
- Chemical Sciences Division , Lawrence Berkeley National Laboratories , Berkeley , California 94720 , United States
| |
Collapse
|
19
|
Hong CM, Morimoto M, Kapustin EA, Alzakhem N, Bergman RG, Raymond KN, Toste FD. Deconvoluting the Role of Charge in a Supramolecular Catalyst. J Am Chem Soc 2018; 140:6591-6595. [PMID: 29767972 DOI: 10.1021/jacs.8b01701] [Citation(s) in RCA: 69] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We have demonstrated that the microenvironment of a highly anionic supramolecular catalyst can mimic the active sites of enzymes and impart rate accelerations of a million-fold or more. However, these microenvironments can be challenging to study, especially in the context of understanding which specific features of the catalyst are responsible for its high performance. We report here the development of an experimental mechanistic probe consisting of two isostructural catalysts. When examined in parallel transformations, the behavior of these catalysts provides insight relevant to the importance of anionic host charge on reactivity. These two catalysts exhibit similar host-substrate interactions, but feature a significant difference in overall anionic charge (12- and 8-). Within these systems, we compare the effect of constrictive binding in a net neutral aza-Cope rearrangement. We then demonstrate how the magnitude of anionic host charge has an exceptional influence on the reaction rates for a Nazarov cyclization, evidenced by an impressive 680-fold change in reaction rate as a consequence of a 33% reduction in catalyst charge.
Collapse
Affiliation(s)
- Cynthia M Hong
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States.,Department of Chemistry , University of California , Berkeley , California 94720 , United States
| | - Mariko Morimoto
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States.,Department of Chemistry , University of California , Berkeley , California 94720 , United States
| | - Eugene A Kapustin
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States.,Department of Chemistry , University of California , Berkeley , California 94720 , United States
| | - Nicola Alzakhem
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States.,Department of Chemistry , University of California , Berkeley , California 94720 , United States
| | - Robert G Bergman
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States.,Department of Chemistry , University of California , Berkeley , California 94720 , United States
| | - Kenneth N Raymond
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States.,Department of Chemistry , University of California , Berkeley , California 94720 , United States
| | - F Dean Toste
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States.,Department of Chemistry , University of California , Berkeley , California 94720 , United States
| |
Collapse
|
20
|
Goehry C, Besora M, Maseras F. Computational Description of a Huisgen Cycloaddition Inside a Self-Assembled Nanocapsule. European J Org Chem 2018. [DOI: 10.1002/ejoc.201800179] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Charles Goehry
- Institute of Chemical Research of Catalonia (ICIQ); The Barcelona Institute of Science and Technology; Avgda. Paı̈sos Catalans 16 43007 Tarragona Catalonia Spain
| | - Maria Besora
- Institute of Chemical Research of Catalonia (ICIQ); The Barcelona Institute of Science and Technology; Avgda. Paı̈sos Catalans 16 43007 Tarragona Catalonia Spain
| | - Feliu Maseras
- Institute of Chemical Research of Catalonia (ICIQ); The Barcelona Institute of Science and Technology; Avgda. Paı̈sos Catalans 16 43007 Tarragona Catalonia Spain
- Departament de Química; Universitat Autònoma de Barcelona; 08193 Bellaterra Catalonia Spain
| |
Collapse
|
21
|
Daver H, Harvey JN, Rebek J, Himo F. Quantum Chemical Modeling of Cycloaddition Reaction in a Self-Assembled Capsule. J Am Chem Soc 2017; 139:15494-15503. [PMID: 29019655 DOI: 10.1021/jacs.7b09102] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Dispersion-corrected density functional theory is used to study the cycloaddition reaction between phenyl acetylene and phenyl azide inside a synthetic, self-assembled capsule. The capsule is first characterized computationally and a previously unrecognized structure is identified as being the most stable. Next, an examination of the free energies of host-guest complexes is conducted, considering all possible reagent, solvent, and solvent impurity combinations as guests. The experimentally observed relative stabilities of host-guest complexes are quite well reproduced, when the experimental concentrations are taken into account. Experimentally, the presence of the host capsule has been shown to accelerate the cycloaddition reaction and to yield exclusively the 1,4-regioisomer product. Both these observations are reproduced by the calculations. A detailed energy decomposition analysis shows that reduction of the entropic cost of bringing together the reactants along with a geometric destabilization of the reactant supercomplex are the major contributors to the rate acceleration compared to the background reaction. Finally, a sensitivity analysis is conducted to assess the stability of the results with respect to the choice of methodology.
Collapse
Affiliation(s)
- Henrik Daver
- Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University , SE-106 91 Stockholm, Sweden
| | - Jeremy N Harvey
- Department of Chemistry, KU Leuven , Celestijnenlaan 200F, B-3001 Heverlee, Belgium
| | - Julius Rebek
- The Skaggs Institute for Chemical Biology and Department of Chemistry, The Scripps Research Institute , 10550 North Torrey Pines Road, La Jolla, California 92037, United States.,Center for Supramolecular Chemistry and Catalysis, Shanghai University , Shanghai 200444, China
| | - Fahmi Himo
- Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University , SE-106 91 Stockholm, Sweden
| |
Collapse
|
22
|
Wang X, Zhang JZH, He X. Quantum mechanical calculation of electric fields and vibrational Stark shifts at active site of human aldose reductase. J Chem Phys 2016; 143:184111. [PMID: 26567650 DOI: 10.1063/1.4935176] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Abstract
Recent advance in biophysics has made it possible to directly measure site-specific electric field at internal sites of proteins using molecular probes with C = O or C≡N groups in the context of vibrational Stark effect. These measurements directly probe changes of electric field at specific protein sites due to, e.g., mutation and are very useful in protein design. Computational simulation of the Stark effect based on force fields such as AMBER and OPLS, while providing good insight, shows large errors in comparison to experimental measurement due to inherent difficulties associated with point charge based representation of force fields. In this study, quantum mechanical calculation of protein's internal electrostatic properties and vibrational Stark shifts was carried out by using electrostatically embedded generalized molecular fractionation with conjugate caps method. Quantum calculated change of mutation-induced electric field and vibrational Stark shift is reported at the internal probing site of enzyme human aldose reductase. The quantum result is in much better agreement with experimental data than those predicted by force fields, underscoring the deficiency of traditional point charge models describing intra-protein electrostatic properties.
Collapse
Affiliation(s)
- Xianwei Wang
- Center for Optics and Optoelectronics Research, College of Science, Zhejiang University of Technology, Hangzhou, Zhejiang 310023, China
| | - John Z H Zhang
- State Key Laboratory of Precision Spectroscopy, Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China
| | - Xiao He
- State Key Laboratory of Precision Spectroscopy, Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China
| |
Collapse
|
23
|
Frushicheva MP, Mills MJL, Schopf P, Singh MK, Warshel A. Computer aided enzyme design and catalytic concepts. Curr Opin Chem Biol 2014; 21:56-62. [PMID: 24814389 PMCID: PMC4149935 DOI: 10.1016/j.cbpa.2014.03.022] [Citation(s) in RCA: 71] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2014] [Revised: 03/19/2014] [Accepted: 03/26/2014] [Indexed: 11/21/2022]
Abstract
Gaining a deeper understanding of enzyme catalysis is of great practical and fundamental importance. Over the years it has become clear that despite advances made in experimental mutational studies, a quantitative understanding of enzyme catalysis will not be possible without the use of computer modeling approaches. While we believe that electrostatic preorganization is by far the most important catalytic factor, convincing the wider scientific community of this may require the demonstration of effective rational enzyme design. Here we make the point that the main current advances in enzyme design are basically advances in directed evolution and that computer aided enzyme design must involve approaches that can reproduce catalysis in well-defined test cases. Such an approach is provided by the empirical valence bond method.
Collapse
Affiliation(s)
- Maria P. Frushicheva
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02142, USA
| | - Matthew J. L. Mills
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089, USA
| | - Patrick Schopf
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089, USA
| | - Manoj K. Singh
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089, USA
| | - Arieh Warshel
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089, USA
| |
Collapse
|
24
|
Raynal M, Ballester P, Vidal-Ferran A, van Leeuwen PWNM. Supramolecular catalysis. Part 2: artificial enzyme mimics. Chem Soc Rev 2013; 43:1734-87. [PMID: 24365792 DOI: 10.1039/c3cs60037h] [Citation(s) in RCA: 665] [Impact Index Per Article: 60.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The design of artificial catalysts able to compete with the catalytic proficiency of enzymes is an intense subject of research. Non-covalent interactions are thought to be involved in several properties of enzymatic catalysis, notably (i) the confinement of the substrates and the active site within a catalytic pocket, (ii) the creation of a hydrophobic pocket in water, (iii) self-replication properties and (iv) allosteric properties. The origins of the enhanced rates and high catalytic selectivities associated with these properties are still a matter of debate. Stabilisation of the transition state and favourable conformations of the active site and the product(s) are probably part of the answer. We present here artificial catalysts and biomacromolecule hybrid catalysts which constitute good models towards the development of truly competitive artificial enzymes.
Collapse
Affiliation(s)
- Matthieu Raynal
- Institute of Chemical Research of Catalonia (ICIQ), Av. Països Catalans 16, 43007 Tarragona, Spain.
| | | | | | | |
Collapse
|