Blog Archives

A new paper on the Weak Link Approach


Chemically actuating a molecule is a very cool thing to do and the Weak Link Approach (WLA) allows us to do precisely that through the reversible coordination of one or various organometallic centers to a longer ligand that opens or closes a macrocyclic cavity. All this leads to an allosteric effect so important in biological instances available in inorganic molecules. Once again, the Mirkin group at Nortwestern University in Evanston, Illinois, has given me the opportunity to contribute with the calculations to the energetic properties of these actuators as well as their electronic properties for their use as molecular scavengers or selective capsules for various purposes such as drug delivery agents.

As in the previous WLA work (full paper), the NBODel procedure was used at the B97D/LANL2DZ level of theory, only this time the macrocycle consisted of two organometallic centers and for the first time the asymmetric opening of the cavity was achieved, as observed by NMR. With the given fragments, all possibilities shown in scheme 1 were obtained. The calculated bond energies for the Pt – S bonds are around 60 – 70 kcal/mol whereas for the Pt – Cl bonds the values are closer to 90 kcal/mol. This allows for a selective opening of the cavity which can then be closed by removing the chlorine atoms with the help of silver salts.

wla

For the case of complex mixture 4a, 4b, and 4c, the thermochemistry calculations show they are all basically isoenergetic with differences in the thousandths of kcal/mol. The possibilities for the groups in the weakly bonded ligands are enormous; currently, there is work being done about substituting those phenyl rings for calix[4]arenes in order to have a macrucyclic capsule made by macrocylic capusules.

Thanks to Andrea D’Aquino for taking me into her project, for all the stimulating discussions and her great ideas for expanding WLA into new avenues; I’m sure she’ll succeed in surprising us with more possibilities for these allosteric macrocycles.

The full paper is published in Inorganic Chemistry from the ACS (DOI: 10.1021/acs.inorgchem.7b02745). Thanks for reading and -if you made it this far into the post- happy new year!

Advertisements

The Evolution of Photosynthesis


Recently, the journal ACS Central Science asked me to write a viewpoint for their First Reactions section about a research article by Prof. Alán Aspuru-Guzik from Harvard University on the evolution of the Fenna-Matthews-Olson (FMO) complex. It was a very rewarding experience to write this piece since we are very close to having our own work on FMO published as well (stay tuned!). The FMO complex remains a great research opportunity for understanding photosynthesis and thus the origin of life itself.

In said article, Aspuru-Guzik’s team climbed their way up a computationally generated phylogenetic tree for the FMO from different green sulfur bacteria by creating small successive mutations on the protein at a time while also calculating their photochemical properties. The idea is pretty simple and brilliant: perform a series of “educated guesses” on the structure of FMO’s ancestors (there are no fossil records of FMO so this ‘educated guesses’ are the next best thing) and find at what point the photochemistry goes awry. In the end the question is which led the way? did the photochemistry led the way of the evolution of FMO or did the evolution of FMO led to improved photochemistry?

Since both the article and viewpoint are both published as open access by the ACS, I wont take too much space here re-writing the whole thing and will instead exhort you to read them both.

Thanks for doing so!

Collaborations in Inorganic Chemistry


I began my path in computational chemistry while I still was an undergraduate student, working on my thesis under professor Cea at unam, synthesizing main group complexes with sulfur containing ligands. Quite a mouthful, I know. Therefore my first calculations dealt with obtaining Bond indexed for bidentate ligands bonded to tin, antimony and even arsenic; yes! I worked with arsenic once! Happily, I keep a tight bond (pun intended) with inorganic chemists and the recent two papers published with the group of Prof. Mónica Moya are proof of that.

In the first paper, cyclic metallaborates were formed with Ga and Al but when a cycle of a given size formed with one it didn’t with the other (fig 1), so I calculated the relative energies of both analogues while compensating for the change in the number of electrons with the following equation:

Fig 1

Imagen1

Under the same conditions 6-membered rings were formed  with Ga but not with Al and 8-membered rings were obtained for Al but not for Ga. Differences in their covalent radii alone couldn’t account for this fact.

ΔE = E(MnBxOy) – nEM + nEM’ – E(M’nBxOy)                     Eq 1

A seamless substitution would imply ΔE = 0 when changing from M to M’

Imagen2.jpg

Hipothetical compounds optimized at the B3LYP/6-31G(d,p) level of theory

The calculated ΔE were: ΔE(3/3′) = -81.38 kcal/mol; ΔE(4/4′) = 40.61 kcal/mol; ΔE(5/5′) = 70.98 kcal/mol

In all, the increased stability and higher covalent character of the Ga-O-Ga unit compared to that of the Al analogue favors the formation of different sized rings.

Additionally, a free energy change analysis was performed to assess the relative stability between compounds. Changes in free energy can be obtained easily from the thermochemistry section in the FREQ calculation from Gaussian.

This paper is published in Inorganic Chemistry under the following citation: Erandi Bernabé-Pablo, Vojtech Jancik, Diego Martínez-Otero, Joaquín Barroso-Flores, and Mónica Moya-Cabrera* “Molecular Group 13 Metallaborates Derived from M−O−M Cleavage Promoted by BH3” Inorg. Chem. 2017, 56, 7890−7899

The second paper deals with heavier atoms and the bonds the formed around Yttrium complexes with triazoles, for which we calculated a more detailed distribution of the electronic density and concluded that the coordination of Cp to Y involves a high component of ionic character.

This paper is published in Ana Cristina García-Álvarez, Erandi Bernabé-Pablo, Joaquín Barroso-Flores, Vojtech Jancik, Diego Martínez-Otero, T. Jesús Morales-Juárez, Mónica Moya-Cabrera* “Multinuclear rare-earth metal complexes supported by chalcogen-based 1,2,3-triazole” Polyhedron 135 (2017) 10-16

We keep working on other projects and I hope we keep on doing so for the foreseeable future because those main group metals have been in my blood all this century. Thanks and a big shoutout to Dr. Monica Moya for keeping me in her highly productive and competitive team of researchers; here is to many more years of joint work.

New paper in Computational and Theoretical Chemistry


I always get very happy to have a new paper out there! I find it exciting but most of all liberating since it makes you feel like your work is going somewhere but most of all that it is making its way ‘out there’; there is a strong feeling of validation, I guess.

Two very different families of calix[n]arenes (Fig 1) were tested as drug carriers for a very small molecule with a huge potential as a chemotherapeutic agent against Leukemia, namely, 3-phenyl-1H-[1]benzofuro[3,2-c]pyrazole a.k.a. GTP which has proven to be an effective in vitro Tyrosine Kinase III inhibitor. Having such a low molecular weight it is expected to have a very high excretion rate therefore the use of a carrier could increase its retention time and hence its activity. This time we considered n = 4, 5, 6 and 8 for the size of the cavities and R = -SO3H and -OEt as functional groups on the upper rim as to evaluate only a polar coordinating group and a non-polar non-coordinating one since GTP offers two H-bond acceptor sites and one H-bond donor one along the π electron density that could form π – π stacking interactions between the aromatic groups on GTP and the walls of the calixarene.

Fig 1. Calixarenes under study and their complexes with GTP

Fig 1. Calixarenes under study and their complexes with GTP

Once again calculations were carried out at the B97D/6-31G(d,p) level of theory along with molecular dynamics simulations for over 100 ns of production runs. NBO Deletion interaction energies were computed in order to discern which hosts formed the most stable complexes.

NBO Del interaction energies B97D/6-31G(d,p)

NBO Del interaction energies B97D/6-31G(d,p)

You may find a link to the ScienceDirect website for downloading the paper from this link. Last, but certainly not least, I’d like to thank all coauthors for their contributions and patience in getting this study published: Dr. Rodrigo Galindo-Murillo; Alberto Olmedo-Romero; Eduardo Cruz-Flores; Dr. Petronela M. Petrar and Prof. Dr. Kunsági-Máté Sándor. Thanks a lot for everything!

fig8

Donor and acceptor H-bond sites increases the probability of keeping the drug in place for a higher retention rate

Donor and acceptor H-bond sites increases the probability of keeping the drug in place for a higher retention rate

%d bloggers like this: