Blog Archives

Chemistry Makes the Chemical

The compound shown below in figure 1 is listed by Aldrich as 4,5,6,7-tetrahydroindole, but is it really?


Fig 1. An indole?

To a hardcore organic chemist it is clear that this is not an indole but a pyrrole because  the lack of aromaticity in the fused ring gives this molecule the same reactivity as 2,3-diethyl pyrrole.  If you search the ChemSpider database for ‘tetrahydroindole’ the search returns the following compound with the identical chemical formula C8H11N but with a different hydrogenation pattern: 2,3,3a,4-Tetrahydro-1H-indole


Fig 2. Also listed as an indole

The real indole, upon an electrophilic attack, behaves as a free enamine yielding the product shown in figure 3 in which the substitution occurs in position 3. This compound cannot undergo an Aromatic Electrophilic Susbstitution since that would imply the formation of a sigma complex which would disrupt the aromaticity.


On the contrary, the corresponding pyrrole is substituted in position 2


These differences in reactivity towards electrophiles are easily rationalized when we plot their HOMO orbitals (calculated at the M062X/def2TZVP level of theory):

If we calculate the Fukui indexes at the same level of theory we get the highest value for susceptibility towards an electrophilic attack as follows: 0.20 for C(3) in indole and 0.25 for C(2) in pyrrole, consistent with the previous reaction schemes.

So, why is it listed as an indole? why would anyone search for it under that name? Nobody thinks about cyclohexane as 1,3,5-trihydrobenzene. According to my good friend and colleague Dr. Moisés Romero most names for heterocyles are kept even after such dramatic chemical changes due to historical and mnemonic reasons even when the reactivity is entirely different. This is only a nomenclature issue that we have inherited from the times of Hantzsch more than a century ago. We’ve become used to keeping the trivial (or should I say arbitrary) names and further use them as derivations but this could pose an epistemological problem if students cannot recognize which heterocycle presents which reactivity.

So, in a nutshell:

Chemistry makes the chemical and not the structure.

A thing we all know but sometimes is overlooked for the sake of simplicity.


#CompChem – Can Orbitals Be Directly Observed?


New paper in Tetrahedron #CompChem “Why U don’t React?”

Literature in synthetic chemistry is full of reactions that do occur but very little or no attention is payed to those that do not proceed. The question here is what can we learn from reactions that are not taking place even when our chemical intuition tells us they’re feasible? Is there valuable knowledge that can be acquired by studying the ‘anti-driving force’ that inhibits a reaction? This is the focus of a new manuscript recently published by our research group in Tetrahedron (DOI: 10.1016/j.tet.2016.05.058) which was the basis of Guillermo Caballero’s BSc thesis.



It is well known in organic chemistry that if a molecular structure has the possibility to be aromatic it can somehow undergo an aromatization process to achieve this more stable state. During some experimental efforts Guillermo Caballero found two compounds that could be easily regarded as non-aromatic tautomers of a substituted pyridine but which were not transformed into the aromatic compound by any means explored; whether by treatment with strong bases, or through thermal or photochemical reaction conditions.


These results led us to investigate the causes that inhibits these aromatization reactions to occur and here is where computational chemistry took over. As a first approach we proposed two plausible reaction mechanisms for the aromatization process and evaluated them with DFT transition state calculations at the M05-2x/6-31+G(d,p)//B3LYP/6-31+G(d,p) levels of theory. The results showed that despite the aromatic tautomers are indeed more stable than their corresponding non-aromatic ones, a high activation free energy is needed to reach the transition states. Thus, the barrier heights are the first reason why aromatization is being inhibited; there just isn’t enough thermal energy in the environment for the transformation to occur.


But this is only the proximal cause, we went then to search for the distal causes (i.e. the reasons behind the high energy of the barriers). The second part of the work was then the calculation of the delocalization energies and frontier molecular orbitals for the non-aromatic tautomers at the HF/cc-pVQZ level of theory to get insights for the large barrier heights. The energies showed a strong electron delocalization of the nitrogen’s lone pair to the oxygen atom in the carbonyl group. Such delocalization promoted the formation of an electron corridor formed with frontier and close-to-frontier molecular orbitals, resembling an extended push-pull effect. The hydrogen atoms that could promote the aromatization process are shown to be chemically inaccessible.


Further calculations for a series of analogous compounds showed that the dimethyl amino moiety plays a crucial role avoiding the aromatization process to occur. When this group was changed for a nitro group, theoretical calculations yielded a decrease in the barrier high, enough for the reaction to proceed. Electronically, the bonding electron corridor is interrupted due to a pull-pull effect that was assessed through the delocalization energies.

The identity of the compounds under study was assessed through 1H, 13C-NMR and 2D NMR experiments HMBC, HMQC so we had to dive head long into experimental techniques to back our calculations.

Fluorescent Chemosensors for Chloride in Water – Sensors and Actuators B: Chemical

A new publication is now available in which we calculated the binding properties of a fluorescent water-soluble chemosensor for halides which is specially sensitive for chloride. Once again, we were working in collaboration with an experimental group who is currently involved in developing all kinds of sustainable chemosensors.

The electronic structure of the chromophore was calculated at the M06-2X/6-311++G(d,p) level of theory under the SMD solvation model (water) at various pH levels which was achieved simply by changing the protonation and charges upon the ligand. Wiberg bond indexes from the Natural Population Analysis showed strong interactions between the chloride ion and the chromophore. Also, Fukui indexes were calculated in order to find the most probable binding sites. A very interesting feature of this compound is its ability to form a cavity without being a macrocycle! I deem it a cavity because of the intramolecular interactions which prevent the entrance of solvent molecules but that can be reversibly disrupted for the inclusion of an anion. In the figure below you can observe the remarkable quenching effect chloride has on the anion.


A quick look to the Frontier Molecular Orbitals (FMO’s) show that the chloride anion acts as an electron donor to the sensor.

Frontier Molecular Orbitals

Frontier Molecular Orbitals

If you are interested in more details please check: Bazany-Rodríguez, I. J., Martínez-Otero, D., Barroso-Flores, J., Yatsimirsky, A. K., & Dorazco-González, A. (2015). Sensitive water-soluble fluorescent chemosensor for chloride based on a bisquinolinium pyridine-dicarboxamide compound. Sensors and Actuators B: Chemical, 221, 1348–1355.

Thanks to Dr. Alejandro Dorazco from CCIQS for asking me to join him in this project which currently includes some other join ventures in the realm of molecular recognition.

Two more students graduated!

It is with great pleasure that I’d like to announce the thesis defense of Guillermo “Memo” Caballero and Howard Diaz who in past days became the second and third students, respectively, to get their B.Sc. degrees with theses completed at our lab. I want to publicly thank them for their hard work which hasn’t only contributed with a thesis to our library but will soon contribute with research papers to our count.

Guillermo “Memo” Caballero worked on the calculation of a reaction mechanism that cannot happen. He started as a synthetic chemist and when he hit a wall at the lab he thought computational chemistry might help him get synthesis on the right direction. He has proven now that the aromatization process of a substituted glutarimide into the corresponding pyridine can only proceed only if substituents with a very strong electron withdrawing effect are used. For two reaction mechanisms proposed, both of them intramolecular rearrangements and only one of them concerted, the calculated energy barriers to reach for the corresponding transition states (QST2 and QST3 methods used) are higher than a pyrolitic decomposition. Memo found also that the delocalization of the pi electron system and its extent goes a long way into the stabilization of the non-aromatic analogue. At first we wanted to treat this problem as a tautomeric equilibrium but since we cannot observe the aromatic tautomer there is no equilibrium and hence no tautomerism. We are still thinking how to name this correspondence between the two compounds when we submit the corresponding paper. It must be said that Guillermo graduated with the highest honors in a most deserved way.

Howard Diaz worked on the design of molecular blockers for the entrance process of the HIV-1 virus into lymphocytes through the GP120 protein. Six known blockers based on phenyl-indoyl-urea were assessed through docking, the binding site of the GP120 protein was described in terms of the interactions formed with each on these compounds and that served as the basis for what in the end came up to be a 36 compound library of blockers, whose structures were first optimized at the B3LYP/6-31G** level of theory. All the 42 blockers were docked in the binding site of the protein and a thorough conformational search was performed. From this set, lead compounds were selected in terms of their binding energies (first calculated heuristically) and further studied at the Density Functional Theory, B97D/cc-pVTZ in order to study the electronic structure of the blocker when interacting with a selection of residues at the binding site. Interaction energies calculated at the quantum level are consistent with the complex formation but since we had to cut the protein to only a few residues little correlation is found with the first calculation; this is fine and still publishable, I just wish we had a more seamless transition between heuristics and quantum chemical calculations. Wiberg indexes were very low, as consistent with a hydrophobic cavity, and delocalization energies calculated with second order perturbation theory analysis on the Natural Bond Orbitals revealed that the two most important interactions are C-H…π and Cl…π, these two were selected as key parameters in our design of new drugs for preventing the HIV-1 virus to bind lymphocytes-T; now we only need to have them synthesized and tested (anyone interested?).

Thank you guys for all your hard work, it has truly payed off. I’m completely certain that no matter what you do and where you go you will be very successful in your careers and I wish you nothing but the very best. This lab’s doors will always remain open for you.


Elements4D – Exploring Chemistry with Augmented Reality

A bit outside the scope of this blog (maybe), but just too cool to overlook. Augmented reality in chemistry education.

Songs | Snaps | Science

This is a guest post from Samantha Morra of, an advertiser on 

Augmented Reality (AR) blurs the line between the physical and digital world. Using cues or triggers, apps and websites can “augment” the physical experience with digital content such as audio, video and simulations. There are many benefits to using AR in education such as giving students opportunities to interact with items in ways that spark inquiry, experimentation, and creativity. There are a quite a few apps and sites working on AR and its application in education.

Elements4D, an AR app from Daqri, allows students explore chemical elements in a fun way while learning about real-life chemistry. To get started, download Elements4D and print the cubes.

There are 6 physical paper cubes printed with different symbols from the periodic table. It takes a while to cut out and put together the cubes, but it…

View original post 475 more words

XIIth Mexican Reunion on Theoretical Physical Chemistry

As every year this month we had the yearly Mexican Reunion on Theoretical Physical Chemistry organized by prominent researchers in the field, such as Dr. Emilio Orgaz (UNAM), Dr. Alberto Vela (CINVESTAV) and many other. Over 150 different works were presented during this edition which took place in Juriquilla, Querétaro at one of the many campuses of the National Autonomous University of Mexico scattered all around the country. Below you can see some pictures from the talks and the first poster session.


This time we contributed with a small poster on a mechanism proposed by Howard Diaz (an undergrad student from UAEM) on the equilibrium transformation of dihydrocinolines into 1-amino-indoles by an intramolecular rearrangement. May this post also serve as the starting point of a -mini-tutorial on how to evaluate a mechanism theoretically using QST3 and IRC in implicitly solvated environments (PCM)


Howard Diaz posing next to his poster

The equilibrium under study and the proposed mechanism  by which it occurs, originally proposed by Frontana-Uribe et al. looks a bit like this:


Dihydrocinolines in equilibrium with 1-aminoindole


Mechanistic proposal by Frontana-Uribe et al.

The energy profile, in which all transition states were calculated with the QST3 method, is presented below, calculated at various levels of theory. Also, the Internal Reaction Coordinate (IRC) connecting both states was calculated and is shown further below in the full poster.

Energy Profile

Energy Profile

From this results we believe that a new mechanistic proposal is needed since the energy barrier for the first step is quite high (~60 kcal/mol) and hence a bit unlikely to occur through that transition state. Nevertheless this is a first approach to elucidating a mechanism and the more knowledge about it the higher the control will be on this chemical transformation.

A full version of the poster is shown below for your convenience (Spanish). See you all at the next RMFQT in Morelia 2014!

Full Poster

Full Poster

#RealTimeChem – Happy birthday DNA!

What a happy coincidence -if indeed it was- that #RealTimeChem week happened to coincide with the sixtieth anniversary of the three seminal papers published in Nature on this day back in 1953, one of which was co-authored by J. Watson and F. Crick; of course I mean the publication for the first time of the structure of deoxyribose nucleic acid, or DNA, as we now call it.

You can get the original Nature papers from 1953 here at: (costs may apply)

Molecular Structure of Nucleic Acids: A Structure for Deoxyribose Nucleic Acid 737


Molecular Structure of Nucleic Acids: Molecular Structure of Deoxypentose Nucleic Acids 738

Molecular Configuration in Sodium Thymonucleate 740

Nature’s podcast released two episodes (called ‘pastcast’) to celebrate DNA’s structure’s birthday, one of them is an interview with Dr. Raymond Gosling who in 1953 worked under Dr. Rosalind Franklin at King’s College London in diffractometry of biological molecules. If you haven’t listened to them you can get them here at Of course, the history around the discovery of DNA’s structure is not without controversy and it has been long argued that the work of Franklin and Gossling didn’t get all deserved credit from Watson and Crick. In their paper W&C acknowledge the contribution of the general nature of DNA from the unpublished results by Franklin’s laboratory but that is as far as they went, they didn’t even mention photo 51 which Crick saw at Wilkins laboratory, who in turn got it from Gossling at Franklin’s suggestion. Still, no one can deny that the helical structure with which we are now familiar is their work, and more importantly the discovery of the specific pairing, which according to Gossling was a stroke of genious that probably couldn’t have happened in his own group, but without Franklin’s diffraction and Gossling’s crystallization  there was little they could do. Details about the process used to crystallize DNA can be heard in the aforementioned podcast, along with an inspiring tale of hard work by Dr. Gossling. Go now and listen to it, its truly inspiring.

For me it was not the story of a helix, that I was familiar with; it was the story of the specific pairing of two hélices
– Dr. Raymond Gosling

Famous Photo 51 by Dr. Rosalind Franklin and Raymond Gosling (Source: Wikipedia)

The iconic Photo 51 by Dr. Rosalind Franklin and Raymond Gosling (Source: Wikipedia)

Above, the iconic Photo 51 taken by Franklin and Gossling (have you ever noticed how most scientists refer to Franklin just as Rosalind but no one refers to Watson as James? Gender bias has a role in this tale too) To a trained crystallographer, the helical symmetry is evident from the diffraction pattern but going from Photo 51 to the representation below was the subject of hard work too.

Modern DNA representation (Source: Wikipedia)

Modern DNA representation (Source: Wikipedia)

There are million of pages written during the last 60 years about DNA’s structure and its role in the chemistry of life; the nature of the pairing and the selectivity of base pairs through hydrogen bond interactions, an interaction found ubiquitously in nature; water itself is a liquid due to the intermolecular hydrogen-bonds, which reminds us about the delicate balance of forces in biochemistry making life a delicate matter. But I digress. Millions of pages have been written and I’m no position of adding a meaningful sentence to them; however, it is a fascinating tale that has shaped the course of mankind, just think of the Human Genome Project and all the possibilities both positive and negative! DNA and its discovery tale will continue to amaze us and inspire us, just like in 2011 it inspired the Genetech company to set a Guiness World Record with the largest human DNA helix.

Genetech SF, Cal. USA (Source

Genetech SF, Cal. USA (Source

Happy birthday, DNA!

%d bloggers like this: