The np * interaction : a rapidly emerging non-covalent interaction

This perspective describes the current status of a recently discovered non-covalent interaction named as the n p* interaction, which is very weak and counterintuitive in nature. In this review, we have provided a brief overview of the widespread presence of this interaction in biomacromolecules, small biomolecules and materials, as well as the physical nature of this interaction explored using various experimental and theoretical techniques. It has been found that this interaction is equally important to other non-covalent interactions for the stability and specific structures of biomolecules and materials. An in-depth understanding of this interaction can help in designing more efficient functional materials as well as drugs. The review also provides a future outlook in terms of exploring the detailed functional role of this interaction in biological processes and its direct spectroscopic evidence, which other commonly known non-covalent interactions (conventional hydrogen bonding, p-hydrogen bonding, p-stacking, etc.) have.


Introduction
Non-covalent interactions are the backbone of biomolecular and material structures, as well as biological processes. [1][2][3][4][5][6][7][8][9][10][11][12][13][14] Thus, an in-depth understanding of these weak interactions may lead to the design of efficient functional materials and drugs. These non-bonding interactions generally have a broad range of classifications depending on their nature, origin and strength. Hydrogen bonding is the most versatile and well understood non-covalent interaction documented in the literature. [1][2][3][4][5] Intermolecular interactions involving p-systems, which are also very popular and well studied, are p-stacking, p-hydrogen bonding and cation-p interactions. [9][10][11][12][13][14][15][16][17][18][19][20][21] However, quite recently, two more non-covalent interactions involving p-electrons, namely anion-p and np* interactions, have been recognized by the scientific community.  Neither of these interactions were appreciated earlier in spite of their widespread presence in biology and materials, because the attractive interaction between an electron rich species and the p-electron cloud is counterintuitive in nature. Although the name ''anion-p interaction'' was conceived by Frontera, Deya, and co-workers 24 in 2002, the interaction between an anion and electron deficient aromatic ring was reported in the literature much earlier. 86 Both np* and anion-p interactions seem to be similar in nature, but the origin of these two attractive interactions is very much different. The anion-p interaction originates mostly from electrostatics (anion-quadrupole) combined with smaller but important contributions of induction (anion-induced polarization) and dispersion. [22][23][24][25][28][29][30][31] On the other hand, it has been found from a few studies that dispersion forces play a significant role in the stability of the np* interaction, but this conclusion demands more detailed investigation in this direction. [32][33][34][35] The cation-p interaction, which is another ion-p interaction, is mainly stabilized by both electrostatics and induction, with a smaller but significant contribution of dispersion. [15][16][17] In some cases with extended p-systems (buckybowls) or alkali and alkaline earth metal cations, the induction effect has a leading contribution in the stabilization of the cation-p interaction. 18,21 The np* interaction is generally classified into two categories: np* Am (amide) and np* Ar (aromatic). 36 The np* Am interaction is widely present in the backbones of proteins and peptoids [37][38][39][40][41][42][43][44][45][46][47][48][49] where the lone pair of electrons on the oxygen atom of a carbonyl group is delocalized into the p* orbital of a neighbouring carbonyl group. [41][42][43][44][45][46][47][48][49] In the case of the np* Ar interaction, lone pair electrons on an oxygen or nitrogen atom are delocalized into the p* orbital of an aromatic ring. [50][51][52][53] Fig. 1 shows examples of np* Am as well as np* Ar interactions and the geometrical parameters used for the characterization of the Burgi-Dunitz trajectory.
The np* Ar interaction is also termed as the lone pair (lp)Á Á Áp interaction. Both types of this np* interaction follow the Burgi-Dunitz trajectory of approach of a nucleophile to an electrophilic system. [54][55][56] To have an np* Am interaction, an OÁ Á ÁCQO distance (d) of r3.2 Å and +OÁ Á ÁCQO (y) of 109 AE 101 are required [see Fig. 1(C)]. It has been noticed that the donor carbonyl distorts the plane of the acceptor carbonyl group, which rises above the plane of its substituents and moves towards the oxygen atom of the donor carbonyl group. Such deviation from planarity can be measured as degree of pyramidalization (D), which is a signature of np* Am interactions in proteins. 43,44 In the case of the np* Ar interaction, a distance (d) of 2.8-3.8 Å between the oxygen or nitrogen atom and aromatic ring centroid as well as an angle (a) of r901 between the plane containing the oxygen or nitrogen atom and aromatic plane are required [see Fig. 1(B)]. The lpÁ Á Áp (np* Ar ) interaction was first discovered in Z-DNA by Egli and co-workers in 1995. 35,53,[58][59][60] Since then, there have been tremendous efforts to recognize the np* interaction in biomolecules, as well as in materials, through extensive search of the X-ray crystal structures in the Cambridge Structural Database (CSD) and Protein Data Bank (PDB). [62][63][64][65][66][67][68] Analysis of the CSD reveals that the occurrence of the lpÁ Á Áp interaction is much more frequent than that of the anion-p interaction. However, the field of the anion-p interaction is much more mature than that of the np* interaction. This fact indeed demands a few reviews on the np* interaction in the literature and this is the motivation for the present perspective.
The most significant and detailed contribution to our current understanding of the np* interaction has come from the recent studies of Raines and co-workers. 44,[69][70][71][72][73][74][75][76] They have investigated the trans/cis ratios in polypeptides, peptoids, and thiopeptoids through NMR spectroscopy and demonstrated that the preferential stability of one of the conformers over the other one is due to the np* interaction. From extensive analysis of the protein crystal structures in the PDB, they have found that the np* interaction, which is present in the backbone of about 45% of amino acid residues of proteins, plays a significant role in the stability of the structures of proteins. 41 The present perspective provides an overview of the recent advancement as well as future outlook on the understanding of this rapidly emerging realm of the np* interaction.
2. np* interaction in biomolecules 2.1. Role of the np* interaction in nucleic acids Egli and co-workers have identified several nucleic acids where significant np* interactions are present. 53,63,65 They have analyzed crystal structures of three forms of Z-DNA: mixed spermine/magnesium, magnesium, and spermine. Z-DNA does not have the proper base stacking that A-DNA and B-DNA have. However, it has been found that the O 4 oxygen of the deoxyribose sugar of the cytidine moiety at each GC site of Z-DNA points directly towards the guanine ring of the guanidium (G) moiety. Fig. 2 exhibits the crystal structure of Z-DNA, showing that the cytidine sugar unit is placed over the guanine base while the stacking interaction between the guanine and cytosine bases is mostly disrupted. The distance between the oxygen atom of the sugar and one of the carbons of the guanine ring in the Professor in the Department of Chemistry and presently, he is an Associate Professor there. His research includes the molecular level understanding of non-covalent interactions through gas phase laser spectroscopy of molecules and complexes relevant to biomolecules and materials using laser desorption as the vaporization source.
Z-DNA has been found to be 2.98 Å, which is less than the sum of the van der Waals radii of oxygen and carbon atoms. The angle between the O 4 atom and the exocyclic C-N bond, i.e. +O 4 Á Á ÁC-N, is 1051. It has been demonstrated from this geometrical arrangement of the oxygen atom of the sugar unit and the guanine ring that an attractive np* interaction between one of the lone pairs of the O 4 atom and the p* orbitals of the guanine ring provides stability in the Z-DNA, in spite of the poor base stacking present there.
The significance of the np* interaction is observed not only in Z-DNA but also in RNA. 63,65 Egli and co-workers have pointed out that there is a possibility of an np* interaction at the U-turns of RNA. U-turns are actually structural motifs found in diverse three dimensional structures of RNA. At a U-turn, the RNA backbone makes a sharp turn of about 1201 from the standard helical conformation between the first and second nucleotides of a UNR sequence, where U is uridine, N is any base and R is usually adenine or guanine. 63 It is interesting to note that U-turns are stabilized by two hydrogen bonds and an np* interaction between an oxygen atom of the phosphate group of the third residue and the face of the nucleobase of the first nucleotide. Such a U-turn is observed in the crystal structure of t-RNA Phe (see Fig. 3A). The GUAA (guanine-uridine-adenine-adenine) tetraloop in hammerhead ribozyme also makes a U-turn where the oxygen atom of the phosphate group of the first adenine is found to be in close contact with the guanine ring. Egli and co-workers have also reported an np* interaction between the oxygen of water molecules and the aromatic ring of unstacked nucleobases in the crystal structure of an RNA pseudoknot obtained from beet western yellow virus (see Fig. 3B). 65 The distance between the oxygen of the water molecule and the ring centroid of cytosine (C8) is 2.93 Å and the dihedral angle between the plane defined by the water molecule and the plane of the ring is 911. It has been found that the np* interaction may have functional importance in the frameshifting activity of the RNA pseudoknot.
In order to understand whether such close contacts between the oxygen of one nucleobase or the phosphates and aromatic ring of other nucleobases can stabilize the structure of Z-DNA, and the U-turns of t-RNA and hammerhead ribozyme, Egli and  co-workers have performed ab initio calculations on model systems and a few of them have been shown in Fig. 4. 63 These model systems are built based on the crystallographic co-ordinates, and single point energies are calculated for each system at the DFT/6-31G* level of theory. The interaction energy for the np* interaction between the oxygen of the ribose sugar of cytosine and the guanine ring in Z-DNA has been found to be insignificant [DE(lp-p)], although it shows stabilization of about À2 kcal mol À1 in the presence of Mg 2+ coordinated to N7 of guanine (Fig. 4A). Protonation of O6 of guanine further increases the interaction energy to À5.6 kcal mol À1 . The interaction between the oxygen of the phosphate group and the ring of uracil in the U-turn of t-RNA has been found to be stable and the interaction energy has been calculated to be À4.8 kcal mol À1 (Fig. 4B). The authors have found that the np* interaction can have a significant contribution towards the stabilization of the structure if the nucleobase is protonated. The np* interactions in the structures of DNA and RNA can be a signature of the polarization or protonation state of a nucleobase. 63 Sankararamakrishnan and co-workers have analyzed high resolution crystal structures of DNA and RNA to find close contacts between the oxygen atom of one base and aromatic centre of another base. 62 They have analyzed 77 crystal structures of DNA having 110 base sequences and 55 crystal structures of RNA containing 76 base sequences. It has been found that the np* interaction is possible in 91 base pairs in DNA and 18 base pairs in RNA. The authors have reported that oxygenaromatic close contacts between the base pairs are very less frequent in RNA compared to in DNA. They have performed ab initio calculations on all the 91 base pairs in DNA showing close contacts between an oxygen atom and the centre of an aromatic ring by taking the co-ordinates of the atoms from the crystal structures. It has been found that such close contacts are indeed attractive in nature.

2.2.
Role of the np* interaction in proteins 2.2.1. np* interactions in the backbones of proteins. It has been found that the np* interaction, similarly to the hydrogen bonding interaction, plays a significant role in the stability of the backbones of proteins. The search for the significance of the np* interaction in proteins began after its recognition in the stability of the collagen triple helix in 2002. 38 Collagen consists of three polypeptide chains connected through hydrogen bonds and forms an extended triple helical structure. Each polypeptide chain consists of a Xaa-Yaa-Gly repeating sequence, where Xaa is usually proline (Pro) and Yaa is either proline or 4R-hydroxyproline (4R-Hyp). Although the Xaa iÀ1 -Pro i peptide bond can have possibilities of both trans and cis conformations, the peptide bonds in the collagen triple helix predominantly exist in the trans form. 38 A section of collagen triple helix crystal structure showing the pro-Hyp-gly sequence has been shown in Fig. 5.
Raines and co-workers have performed NMR spectroscopy and computational studies on three model compounds (AcProOMe, Ac-4R-Flp-OMe, and Ac-4S-Flp-OMe) to understand  the role of the np* interaction in the stability of the collagen triple helix. 38 Fig. 6 shows the structures of the trans conformers of the three model compounds. They have found through NMR spectroscopy studies that the equilibrium constant (K trans/cis ) between the trans and cis conformations is greater than 1 for all three of the compounds. The results obtained from theoretical calculations performed on these model compounds correlate well with the experimental results. The trans conformation allows closer approach of the amide oxygen to the carbonyl group of the ester, following the Burgi-Dunitz trajectory. It has been found from the Natural Bond Orbital (NBO) analysis of these model compounds that there is a significant overlap between the lone pair orbital on the amide oxygen and the empty p* orbital of another adjacent carbonyl group in the trans conformation. The second order perturbative energy for this interaction in the trans exo form of Ac-4R-Flp-OMe has been calculated to be 1.38 kcal mol À1 . The crystal structure of Ac-4R-Flp-OMe shown in Fig. 7 also has a good resemblance to the theoretically calculated structure. These results indicate that the trans conformation of the Xaa iÀ1 -Pro i peptide bonds in the collagen triple helix are stabilized by np* interactions between the lone pair of the carbonyl oxygen of position i À 1 and the antibonding orbital of the carbonyl group of position i. But does the np* interaction exist in other secondary structures of proteins too? To answer this question, Raines and co-workers created a computational library of allowed conformations of proteins using a model peptide of four alanine residues (AcAla 4 NHMe) capped with an acetyl group at the N terminus and -NHMe group at the C terminus. 41 They performed NBO analysis on the computationally generated allowed conformations to predict the ones that could have an np* interaction, and those conformations were depicted on a Ramachandran type plot based on the F and C values for the residue pairs. This study was followed by PDB analysis of a high resolution (r1.6 Å) and non-redundant subset of 1731 protein crystal structures to search for close contact between the carbonyl groups as observed in the collagen triple helix. The allowed distance between the oxygen of the carbonyl group of residue i À 1 and the carbon of the carbonyl group of residue i as well as the O iÀ1 Á Á ÁC i QO i angle were set to r3.2 Å and 991-1191, respectively, for the search. It has been found that there are abundant np* interactions in the allowed regions of the Ramachandran plot. Further, the computational search of the possible conformations of proteins which can give rise to an np* interaction correlates well with the actual PDB search. 41 Christian Fufezan followed a similar procedure to determine the significance of the np* interaction in the structural stability of proteins. 66 The author performed PDB analysis of a non-redundant set of 5296 protein crystal structures using python module p3d with similar search criteria for the np* interaction to those applied by Raines and co-workers. 41 About 45.1% of all the residues in the non-redundant dataset satisfied the search criteria. Simulations on model peptides were also performed to search for the possible conformations showing an np* interaction.
Both the studies by Raines and co-workers and those of Christian Fufezan show that the regions for the np* interaction identified through computational methods correspond well to those searched for through PDB analysis. The Ramachandran plots generated in these studies show that the highest propensity for the np* interaction is in left and right handed a-helices. However, the presence of the np* interaction between the i À 1 and i residues in b-strands is insignificant. Nevertheless, the np* interaction can exist in twisted b-strands, which are also known as twisted b-sheets or b-bulges. Raines and co-workers have pointed out two common types of b-bulges, namely G1 and wild types, having a significant amount of np* interaction. 41 This interaction is also observed in 3 10 helices and polyproline II helices. 41 The np* interaction is responsible not only for the stability but also the function of proteins. Raines and co-workers have found a series of four np* interactions between consecutive residues in the selectivity filter of the K + channel from Streptomyces lividans. 41 They have observed that the np* interaction is present only in the region responsible for the filtering of K + ions, while this interaction is nearly absent in that region at high concentrations of K + ions. They pointed out that the np* interactions could stabilize the conformation of the filter region of the protein at low concentrations of K + ions.
2.2.2. np* interactions involving the side chains of proteins. The np* interaction exists not only in the backbone but also in the side chains of proteins. Sankararamakrishnan and co-workers have analyzed high resolution crystal structures of about 500 proteins and reported that the oxygen atom of the backbone carbonyl group interacts frequently with side chain carbonyl groups, as well as aromatic centers. 50,64,68 Thus, both np* Am and np* Ar types of interactions are present in Fig. 6 Structures of the trans conformer of three model compounds designed to understand the role of the np* interaction in the collagen triple helix. 38 Flp stands for fluoroproline. proteins. Four amino acid residues containing aromatic side chains are histidine, tyrosine, phenylalanine and tryptophan. It has been found that there are 249 examples where a close contact of about 3.5 Å between an aromatic center and backbone carbonyl oxygen is present. 64 On the other hand, the occurrence of the interaction between a side chain carbonyl oxygen and aromatic center is very less frequent (37 examples). Fig. 8(A) shows an example of an np* Am interaction between the backbone carbonyl oxygen atom of an alanine residue and the aromatic center of a phenylalanine residue in the protein bacterial luciferase enzyme (PDB ID:1LUC). 64 A similar interaction between the side chain carbonyl oxygen of glutamine and the aromatic center of tyrosine present in the protein Serratia marcescens endonuclease (PDB ID:1QL0) is shown in Fig. 8(B). 64 It has been reported that the aromatic residue that is involved most frequently in this interaction is histidine. 64 Close contact between the carbonyl groups of the backbone and side chain of the same amino acid residue aspartic acid (Asp) has been identified by Sankararamakrishnan and co-workers from an exhaustive search of high resolution protein crystal structures in the PDB. 68 They have found that these intraresidue close contacts between the carbonyl groups are present in about 102 Asp residues. Fig. 9 shows an example of this interaction in an Asp residue of the protein epidermolytic toxin A (PDB ID:1AGJ). 68 The authors have also recognized np* interactions between water oxygen atoms and the p-clouds of aromatic residues in proteins through the detailed analysis of high resolution crystal structures of proteins as well as quantum chemistry calculations of model compounds based on protein structures in the presence of different orientations of water molecules. 50 2.2.3. Co-operativity between hydrogen bonding and np* interactions. NBO analysis performed by Raines and co-workers on a model a-helix peptide shows that the np* interaction is present due to delocalization of the p-rich lone pair of the oxygen atom of the carbonyl group of residue i À 1 over the p* orbital of the carbonyl group of residue i. 41 Interestingly, the s-rich lone pair of the same oxygen is delocalized over the s* antibonding orbital of the N-H bond of residue i + 4 to form a N-HÁ Á ÁOQC hydrogen bond. They have called it kinship of a hydrogen bond and a np* interaction in an a-helix. Raines and co-workers have also performed PDB analysis and found several crystal structures of proteins in which the oxygen atom of the side chain of an asparagine residue is involved in hydrogen bonding with a hydrogen bond donor as well as a np* interaction with a carbonyl group of the backbone of the protein. 73 Fig. 10 shows the involvement of the lone pair orbitals of the oxygen atom of the side-chain of an asparagine residue in the hydrogen bond with the N-H group of residue i + 2 and the np* interaction with the carbonyl group of the main chain. The np* interaction contributes about 0.5 kcal mol À1 in the stabilization of the a-helix. Here the presence of the np* interaction reduces the O iÀ1 Á Á ÁC i QO i distance, which favours the hydrogen bond formation. Thus the existence of the np* interaction increases the strength of the hydrogen bond and this is called positive co-operativity. 73 Very recently, Das and co-workers have studied the interplay between np* Ar and N-HÁ Á ÁN hydrogen bonding interactions in the complexes of 7-azaindole and 2,6-substituted fluoropyridines. 34 They have found that the strength of the hydrogen   bond in these complexes decreases due to the presence of the np* Ar interaction. This negative co-operativity is mostly due to subtle competition between the weak np* Ar (about 0.5 kcal mol À1 ) and strong hydrogen bonding (about 4 kcal mol À1 ) interactions. It is interesting to point out that the weak interaction present here can govern the final geometry of these complexes despite the presence of the strong hydrogen bonding interaction. It is also intriguing to note in this case that the lone pairs on two different nitrogen atoms take part in the hydrogen bonding and np* Ar interactions, while there is participation by two lone pairs on the same oxygen atom in the two non-covalent interactions in an a-helix. Fig. 11 shows the initial and final structures of 7-azaindoleÁ Á Á2,6-difluoropyridine calculated at the M05-2X/cc-pVTZ level of theory. The optimized structure exhibits both N-HÁ Á ÁN hydrogen bonding and np* Ar interactions.

Significance of the np* interaction in materials
The np* interaction plays a significant role in the structures of not only biomacromolecules but also materials. [78][79][80][81][82][83][84][85] There are quite extensive reports in the literature which show the existence of the np* interaction as well as its role in the supramolecular self assembly of molecules to give rise to specific structures of materials. Reedijk and co-workers have named this interaction as a new supramolecular bond. 67 They have synthesized a [Zn 4 (oxodentriz)Cl 8 ](CH 3 CN) 2 (H 2 O) complex in acetonitrile at room temperature and determined the crystal structure, which shows that the nitrogen atom of acetonitrile is in close contact with the aromatic ring of triazine (see Fig. 12). 78 The distance between the nitrogen atom of the acetonitrile molecule and the center of the triazine ring is 3.25 Å and the angle of approach of the nitrogen atom of the acetonitrile molecule towards the center of the triazine ring is 75.21. These favourable geometrical parameters indicate the presence of an lpÁ Á Áp interaction between the nitrogen atom of acetonitrile and heteroaromatic ring triazine.
Further, Reedijk and co-workers have designed two triazinebased compounds which self assemble to form structures stabilized by both C-HÁ Á Áp and lpÁ Á Áp interactions. 79 The crystal packing of the self assembled structure of one of the compounds has been shown in Fig. 13, which shows that each molecular unit interacts with four neighbouring molecules by means of two ClÁ Á Áp and two weak C-HÁ Á Áp interactions. Later, the same group reported the crystal structure of a co-ordination complex, [Cu(Hspet)(NO 3 ) 2 ](NO 3 )Á0.5H 2 O, synthesized from the reaction of the ligand 2,4,6-tris-N-[N-2-(methylsulfanyl)-N-(pyridine-2-ylmethyl)ethanamine]-1,3,5-triazine (spet) with copper II nitrate. 81 This crystal structure is depicted in Fig. 14, and it shows that the third substituent nitrogen atom (N19) of the 1,3,5-triazine ring is in close contact with the triazine ring of the adjacent co-ordination molecule, and the N19 atom of the latter co-ordination molecule is also in close contact with the centroid of the triazine ring of the former co-ordination molecule. Thus, a double lpÁ Á Áp interaction holds the two complexes together in their crystal packing. 81 Mukhopadhyay and co-workers have designed two co-ordination complexes of nickel II malonate, namely (C 5 H 7 N 2 ) 4     Caracelli and co-workers have reviewed several crystal structures where antimony (Sb) and bismuth (Bi) participate in the lpÁ Á Áp interaction. 83 They have performed CSD analysis to find such structures. It has been found that the metal lpÁ Á Áp interaction is independent of other supramolecular binding motifs like hydrogen bonding, pÁ Á Áp stacking, etc. Thus, the structure can be stabilized by only SbÁ Á Áp interactions or BiÁ Á Áp interactions. A representative structure having SbÁ Á Áp interactions obtained from the CSD search is shown in Fig. 16 This structure shows that the lpÁ Á Áp interaction between the metal atom (Sb) and aromatic ring can lead to the formation of small onedimensional or large three-dimensional structure. The authors have also found the existence of the Sb(lp)Á Á Áp interaction in biomolecules. Sb is found in the FAD (flavine adenine dinucleotide) active site of chains A and B of trypanothione reductase. 83 Another type of lpÁ Á Áp interaction involving a sulfur atom and an aromatic ring, i.e. sulphurÁ Á Áp interaction, has been found to be very important for the stabilization of protein structures. Sherrill and co-workers have studied the benzeneÁ Á ÁH 2 S complex as a prototype for the sulphurÁ Á Áp interaction and reported that the configuration with the sulphur atom facing towards the p-cloud is favorable. 87,88 The presence of sulphurÁ Á Áp interactions in solid state crystal structures has also been reported recently. 89

Significance of the np* interaction in small molecules
The significance of the np* interaction is observed even in small molecules like N-acyl homoserine lactones, aspirin, poly-(lactic acid), g-aminobutyric acid (GABA) and peptoids. 36,44,70,71,77 N-Acyl homoserine lactone (see Fig. 17A) acts as a signalling molecule and is released in response to the population density of a bacterial colony, a phenomenon which is known as quorum sensing. 44 It binds to LuxR-type receptor proteins and causes changes in growth, virulence, antibiotic resistance and other phenotypes. N-Acyl homoserine lactones have two carbonyl groups adjacent to each other, as in proline residue (see Fig. 17). 44 Raines and co-workers have argued that N-acyl homoserine lactones, similarly to the proline residue, can exhibit an np* interaction between the proximal carbonyl groups. 44 They were able to crystallize N-trimethyl homoserine lactone, which is not a natural N-acyl homoserine lactone, but can mimic the natural one. Natural N-acyl homoserine lactones cannot be crystallized easily due to the presence of a large alkyl group. The crystal structure of N-trimethyl homoserine lactone shows that the distance between the N-acyl carbonyl oxygen and carbon of the lactone carbonyl group is 2.73 Å and the angle of approach of the donor oxygen towards the acceptor carbonyl group is 90.61. These geometrical parameters are very close to the Burgi-Dunitz trajectory for the attack of a nucleophile on an electrophile. This indicates that there is a possibility of an np* interaction between the donor oxygen of the N-acyl carbonyl group and the acceptor p* orbital of the lactone carbonyl group. The authors have further performed DFT calculations to obtain the most stable geometry of N-trimethyl homoserine lactone and carried out NBO analysis of this structure. NBO analysis shows significant overlap between the lone pair orbital (n) of the donor oxygen and the acceptor p* orbital of the lactone carbonyl group with an interaction energy of 0.64 kcal mol À1 . They have examined the crystal structure of bound N-acyl homoserine lactones with their LuxR receptors and found that the bound-state conformation of N-acyl homoserine lactones changes due to the formation of a hydrogen bond between the N-acyl carbonyl oxygen and the tyrosine residue of the receptor prohibiting an np* interaction. 44 Thus, there is a competition between the hydrogen bonding and np* interaction at the time of binding of the N-acyl homoserine lactone with its receptor. Any factor which can attenuate the strength of the np* interaction will increase the binding affinity of N-acyl homoserine lactones with their receptors. Thus, the np* interaction has valuable significance in governing the biological activity of N-acyl homoserine lactones.
g-Amino butyric acid (GABA) is an inhibitory neurotransmitter in the brain and spinal cord. Lopez and co-workers have identified the most stable conformation of GABA using laser ablation (LA) molecular beam Fourier transform microwave spectroscopy (LA-FTMW) combined with ab initio calculations. 77 Microwave spectroscopic data show evidence of the nine conformers of GABA in the gas phase, which differ in the folding and orientation of the functional groups present in the molecule. Out of the nine conformers, five conformers have folded form with intramolecular non-covalent interactions between the functional groups, while the remaining four structures with no intramolecular non-covalent interactions are called extended conformations. Among the five folded conformations, three show the possibility of intramolecular hydrogen bonding and two show the possibility of an np* interaction between the donor lone pair (n) orbital of nitrogen and acceptor (p*) orbital of the carbonyl group. Thus, the np* interaction also contributes to the conformation of an important inhibitory neurotransmitter.
Raines and co-workers have found the presence of the np* interaction in poly-(lactic acid) (PLA), a versatile polyester. 70 They have examined the crystal structure of a-PLA and found that its backbone torsional angles, i.e. F and C, are similar to those of the strands of the collagen triple helix. Fig. 18 shows the structure of a-PLA. NBO analysis of the crystal structure of a-PLA at the B3LYP/6-311++G(2d,p) level of theory shows that the energy for each np* interaction is 0.44 kcal mol À1 . This indicates that the np* interaction has significance in stabilizing the structure of polyesters, which are similar to polyamides. The authors have pointed out that poly-(lactic acid) is similar to polyalanine and the only difference is that the amide is replaced by an ester bond, resulting in the absence of a N-H hydrogen bond donor in poly-(lactic acid). The existence of the np* interaction even in the absence of hydrogen bonding indicates that this interaction can exist independently and it is not just a short contact imposed by the geometric constraint due to the hydrogen bonding patterns. Thus, the np* interaction can exist in proteins and its presence is independent of the hydrogen bonding patterns.
The np* interaction is also found to regulate the conformations of peptoids. Blackwell and co-workers have designed small peptoids to study these interactions. 36 Fig. 19A shows the general structure of a peptoid, while structures of peptoids representing np* Am and np* Ar interactions have been depicted in Fig. 19B and C, respectively. In the case of the peptoid in Fig. 19B, the trans-amide is more stable than the cis-amide, while the peptoid in Fig. 19C shows a propensity for the cis-amide compared to the trans-amide. The authors have determined the preferential stability of one of the conformers of these peptoids by measuring K cis/trans using NMR spectroscopy and this conformational preference could be only due to the presence of the np* interaction. Raines and co-workers have found the significance of the np* interaction in the structure as well as chemical and biological functions of the well known analgesic medicine aspirin. 71 Aspirin is an O-acetylated form of salicylic acid. The authors have found that the crystal structure of aspirin shows close contact between the oxygen atom of the alcohol group of the carboxylic acid and the carbon atom of the carbonyl group of the ester (see Fig. 20A) in a Burgi-Dunitz trajectory fashion. The crystal structure also shows significant pyramidalization of the acceptor carbonyl group, which is the key signature of the np* interaction. The pyramidalization parameter (D) has been found to be 0.011 Å. They have found six low energy conformations of aspirin by performing ab initio calculations. Interestingly, the lowest energy conformer shows an np* interaction between the oxygen atom of the carbonyl group of the carboxylic acid and the carbonyl group of the ester. The interaction energy of the np* interaction calculated through NBO analysis of this conformer has been found to be 1.73 kcal mol À1 . The second lowest energy conformer resembles the crystal structure of aspirin where the donor oxygen atom is of the alcohol group of the carboxylic acid. The NBO analysis of this conformer also shows an np* interaction of 0.95 kcal mol À1 in energy. Fig. 20B shows these two structures that have an np* interaction. It has been found that aspirin needs to pass through a paraffin-like channel without any hydration to perform its pharmacological activity. 71 The authors have explained that the intramolecular np* interaction in aspirin may shield the donor oxygen from interaction with solvent and help its passage through the hydrophobic channel. Thus, it is interesting to point out that the np* interaction could facilitate the chemical and pharmacological functions of aspirin.

The np* Ar or lpÁ Á Áp interaction
Among the various ion-p interactions, in general, the cation-pi interaction is quite obvious and easy to realize but neither anion-p nor lpÁ Á Áp interactions are straightforward. These two interactions (anion-p and lpÁ Á Áp) are generally found to be favorable when the or aromatic ring is electron deficient or p-acidic. However, it has been reported recently that both lpÁ Á Áp and anion-p interactions can be favorable in the case of electron rich p-clouds. 27,33,35,91,92 5.1.1. The lpÁ Á Áp interaction in complexes of electron deficient aromatic rings. The very first theoretical study on the lpÁ Á Áp interaction was performed on the C 6 F 6 Á Á ÁH 2 O   complex by Dougherty and co-workers in 1999. 57 They found that water binds to hexafluorobenzene through the oxygen atom pointing towards the face of the p system. On the other hand, benzene interacts with water through p-hydrogen bonding between a hydrogen atom of the water and the p-cloud of the benzene ring. The structures of C 6 H 6 Á Á ÁH 2 O and C 6 F 6 Á Á ÁH 2 O calculated at the cp-MP2/6-31G** level of theory have been provided in Fig. 21. Benzene and hexafluorobenzene have a quadrupole moment that is equal in magnitude but opposite in sign. The electrostatic potential above and below the ring is negative in benzene while the same is positive in hexafluorobenzene. As the lone pair electrons on the oxygen atom of water are attracted towards the electron deficient aromatic ring in the C 6 F 6 Á Á ÁH 2 O complex, the authors have argued that the lpÁ Á Áp interaction in the complex originates from electrostatics. Besnard and co-workers have also reported the same geometry of the C 6 F 6 Á Á ÁH 2 O complex from ab initio calculations. 58 The binding energy of the C 6 F 6 Á Á ÁH 2 O complex has been calculated to be about À2 kcal mol À1 at the MP2 level of theory. Alkorta and co-workers have performed an ab initio study on complexes of hexafluorobenzene with small electron rich moieties (FH, HLi, CH 2 , NCH, and CNH) at the HF, MP2 and B3LYP levels of theory. 59 In the case of all of these complexes, geometrical minima have been found when the electron rich atoms of these small molecules point directly towards the face of the hexafluorobenzene. It seems from these studies on the complexes of hexafluorobenzene that lone pair electrons on the electronegative atoms can interact only with electron deficient aromatic rings through electrostatic interaction.
5.1.2. The lpÁ Á Áp interaction in complexes of electron rich aromatic rings. The question now should be asked whether the lpÁ Á Áp interaction can occur between a lone pair containing atom and an electron rich aromatic ring. If this interaction is possible with an electron rich aromatic ring, electrostatics cannot be the major driving force. There are a few theoretical calculations reported in the literature to address this issue. It is indeed interesting that the lpÁ Á Áp interaction is favorable between the oxygen atom of H 2 O and polycyclic aromatic hydrocarbons (PAHs, i.e. anthracene, triphenylene, coronene, circumcoronene, etc.) although the same interaction is repulsive in the benzeneÁ Á ÁH 2 O complex. 91,92 It has been also reported that the lpÁ Á Áp interaction in PAHÁ Á ÁH 2 S is much stronger than that in PAHÁ Á ÁH 2 O. 91 Interestingly, this interaction is attractive even in the benzeneÁ Á ÁH 2 S complex. 91 Very recently, it has been demonstrated that substituted benzenes with electron donating groups have attractive lpÁ Á Áp interactions with dimethyl ether (Me 2 O) and trimethyl ammonia (NMe 3 ) while similar interaction is absent with water and ammonia. 33,35 All of these results reveal that the increase in the dispersion in these systems switches on the lpÁ Á Áp interaction by overcoming the repulsive electrostatic interaction. Thus, it seems that a dispersion interaction is a key to having a lpÁ Á Áp interaction in the system. However, it is very difficult to conclude from these data that dispersion is the main contribution in the stability of the lpÁ Á Áp interaction.
In general, the introduction of bulky groups in the systems increases the dispersion contribution in the total interaction energy for all types of non-covalent interactions. Recently, Hoja et al. studied several dimers starting from a water dimer to a tertiary-butyl alcohol dimer by increasing the size of the alkyl substituent. 93 They have reported that the dispersion becomes very important in the stabilization of the hydrogen bonded (O-HÁ Á ÁO) complexes as the substituents become bulkier. Cabaleiro-Lago and co-workers have shown that the anion-p interaction is also stable in electron rich aromatic systems, i.e. extended PAHs (buckybowls). 27 The dispersion interaction also becomes significant in the case of the cation-p interaction when the bulkiness of the cation increases. 19,20 In Section 2.2.3, the competition between the lpÁ Á Áp interaction and hydrogen bonding in the complexes of 7-azaindole and several 2,6-substituted fluoropyridines studied by Das and co-workers has been discussed. 34 They found that the binding energies of the 7-azaindoleÁ Á Á2,6-substituted fluoropyridines calculated at the B3LYP level of theory are very much less than those obtained at the dispersion-corrected DFT and MP2 levels. Fig. 22 shows the plots of the binding energies of these complexes calculated at different levels of theory. 34 The discrepancy  in the results obtained from the B3LYP and dispersion-corrected DFT levels indicates that the dispersion interaction has a significant role in the stability of these complexes.
The most significant finding from the study of 7-azaindoleÁ Á Á2,6substituted fluoropyridines has come out through decomposition of the total interaction energy. It has been found that there is a significant increase in the dispersion component of the total interaction energy in 7-azaindoleÁ Á Á2,6-substituted fluoropyridines with the increase in fluorine substitution on the pyridine ring. Fig. 23 shows the different components of the total interaction energy of these complexes, calculated at the M05-2X/cc-pVDZ level of theory. It is intriguing to note that the increase in the dispersion force compared to the electrostatic one is greater with the increase in fluorine substitution, although fluorine has a smaller size and less polarizability. Thus, the point is that the increase in fluorine substitution in the system cannot increase the dispersion contribution, unlike the case in benzeneÁ Á ÁNMe 3 or benzeneÁ Á ÁMe 2 O complexes where the introduction of bulky methyl groups increases the dispersion. Indeed, it has been reported that there is hardly any change in the dispersion contribution with the increase in fluorine substitution in the case of the complexes of 7-azaindole and fluorosubstituted pyridines, which have planar double hydrogen bonded structures. The different components of the total interaction energy in the complexes of 7-azaindole with 2,3-difluoropyridine, 2,3,4-trifluoropyridine and 2,3,4,5tetrafluoropyridine are depicted in Fig. 24. This observation clearly points out that the increase in the dispersion force with the increase in fluorine substitution in 7-azaindoleÁ Á Á2,6fluoropyridines is due to the lpÁ Á Áp interaction only. Thus, it could be stated that dispersion plays a major role in the stability of the lpÁ Á Áp interaction.
The lpÁ Á Áp interaction in complexes of electron rich aromatic rings has also been explored through experiments by Gung and co-workers. They have designed triptycene scaffolds (see Fig. 25) which have anti and syn conformations. 51 The syn conformation of triptycene allows close contact between the benzyl ring, attached to the C(9) carbon atom and oxygen atom of the OMe/RCOO (R = H, Me, Et, i-Pr, CF 3 ) group attached to the C(1) carbon atom. Both electron withdrawing groups (NO 2 , CN, CF 3 , F, Br, Cl) and electron donating groups (Me, MeO, Me 2 N) have been used as a substituent (X) on the benzyl ring, attached to the C(9) carbon atom, keeping the substituents (MeO or ester group) at the C(1) carbon atom fixed. The authors have measured the NMR spectra of these triptycene scaffolds at low temperature to determine the ratio of syn and anti conformers in each of these cases. The syn to anti ratio reflects the degree of interaction of the arene ring with the oxygen atom of the OMe or ester group attached to the C(1) carbon atom. The ratio has been found to be 2 : 1 when no interaction is present between the C(1) and C(9) substituents. Therefore a ratio greater than 2 : 1 will indicate the preference of the syn over anti conformation due to the presence of the attractive arene-oxygen interaction.
Gung and co-workers have observed that the syn/anti ratio with OMe as the substituent at C1 is greater than 2 (B2.4) even if the aromatic ring is electron rich, i.e. X = Me, MeO. The highest syn/anti ratio (3.5-5.5) is observed when X is an electron withdrawing group like CF 3 , Cl, etc. 51 They have initially concluded that the origin of the attraction could be electrostatic as the interaction between the arene and oxygen atom increases with the decrease in electron density of the aromatic ring. To explain the observed value of the syn/anti ratio greater than 2 in the presence of the electron donating substituent, they have pointed out that the oxygen atom of the methoxy group is close to the edge of the aromatic ring, which has a belt of positive electrostatic potential. Therefore the syn conformation, even in the case of the arene ring with an electron donating substituent, is preferentially stable over the anti conformation.
In a later study, Gung and co-workers replaced the OMe group in triptycene with a methylmethoxy group (MOM) to place the oxygen atom of the methoxy group close to the centre of the aromatic ring in the syn conformation, and studied the system with low temperature NMR spectroscopy (see Fig. 26). 32 Keeping the MOM group fixed at C1, the substituent X of the arene ring has been varied from electron withdrawing to electron donating groups. They have observed a syn/anti ratio of 7-8 in the case of the electron withdrawing groups (X) and a ratio of 3-4 even in the case of the electron donating group. The syn/anti ratio has been found to be 25 and 6.8 when the arene ring is hexafluorobenzene and benzene, respectively. Here the authors have pointed out that the syn conformation of trypticene with electron donating substituents on the arene ring is still preferred over the anti conformation in spite of the unanticipated close proximity between the oxygen atom and electron rich aromatic ring in the former. These results indicate that there is an attractive interaction between the oxygen atom of the MOM group and the electron rich aromatic ring. The above observation cannot be rationalized on the basis of the electrostatic potential of the ring. Thus, eventually the authors have concluded that the major driving force for such an interaction cannot be electrostatic.

The np* Am interaction
To unearth the origin of the np* Am interaction, Raines and co-workers have studied the nature of the carbonyl-carbonyl interaction in proteins using the model system shown in Fig. 27. 37 The model system contains a pyrrolidine ring with acetyl substitution at the C2 position and acetyl or thioacetyl substitution at the N1 position. This model system can exist in both trans and cis conformations. However, close contact between the adjacent carbonyl groups in the Burgi-Dunitz trajectory fashion can occur only in the trans conformation. Further, this model system mimics the proline amino acid residue in proteins. When the acetyl group is substituted at the N1 position, the trans/cis ratio has been observed to be 4.6, i.e. the trans conformer is more stable. Here the authors have asked the question whether the nature of the carbonyl-carbonyl interaction is due to charge-charge, dipole-dipole or np* electronic delocalization (see Fig. 28). In order to understand the source of this interaction, they have replaced the oxygen donor (O iÀ1 ) with a sulphur donor (S À1 ) in their model system and measured the equilibrium constant (K trans/cis ) between trans and cis conformations through NMR studies. There will be a decrease in the trans/cis ratio if a charge-charge interaction is responsible for the close contact of the two carbonyl groups since sulphur is less electronegative than oxygen. On the other hand, an increase in the trans/cis ratio will indicate that the origin of the carbonyl-carbonyl interaction is due to dipoledipole or np* electronic delocalization because CQS has a larger dipole moment than CQO and sulphur, being a softer base than oxygen, is more efficient in donating an electron pair than oxygen.
Raines and co-workers have observed an increase in the trans/cis ratio when the oxygen donor has been replaced by sulphur. 37 This result indicates that the preferential stability of the trans conformer over cis is not due to a charge-charge interaction and is rather due to a dipole-dipole interaction or np* electronic delocalization. Delocalization of electrons is possible when there is an appreciable overlap between the donor and acceptor orbitals. To have proper orbital overlap, the distance between the donor and acceptor atoms should be less than the sum of the van der Waals radii of both the atoms.   As the np* interaction is more sensitive to the distance between the donor and acceptor atoms in comparison to the dipole-dipole interaction, an increase in the distance between these two interacting atoms will affect the strength of the np* interaction. Generally, an electron withdrawing group substituted at the C4 carbon atom (where the C4 carbon atom is in the S-configuration) of the pyrrolidine ring (see Fig. 27), influences the pyrrolidine ring to be in C g -endo puckered conformation. In this conformation, the distance between the O iÀ1 or S iÀ1 donor and C i QO i acceptor increases and it is greater than the sum of the van der Waals radii of oxygen/ sulphur and carbon. This indicates that proper orbital overlap is not possible in the 4S configuration but dipole-dipole interaction is possible. Thus, the authors argued that an increase in the trans/cis ratio upon replacement of the oxygen donor by sulphur in the C g -endo puckered configuration of the pyrrolidine ring will indicate that the dipole-dipole interaction is the driving force for the carbonyl-carbonyl interaction. They have observed that the value of K trans/cis for the endo puckered conformation of the pyrrolidine ring with a 4S EWG is similar with both oxygen and sulphur as donors. This indicates that the trans conformation is not stabilized by dipole-dipole but rather the np* interaction is the origin of the stabilization.
To have a favourable dipole-dipole interaction between two carbonyl groups, two interacting dipoles should have proper orientation, i.e. they should be orthogonal to each other. 72 In other words, it restricts all the four atoms to a fixed orientation. On the other hand, close contact between two carbonyl groups for the np* interaction demands proper orbital overlap between O iÀ1 and C i QO i and this does not restrict the positions of all the four atoms. In the case of the np* interaction, only the donor atom (O iÀ1 ) and acceptor carbonyl group (C i QO i ) are required to be in favourable orientation and the positive pole of the donor carbonyl group, i.e. C iÀ1 , is not restricted. It has been found that the two adjacent carbonyl dipoles in an a-helix are in a repulsive orientation but there is still close contact between them. Raines and co-workers have demonstrated that the secondary structures of proteins exhibit close contact between the adjacent carbonyl groups and there is an orientation preference on the O iÀ1 Á Á ÁC i QO i angle but not the CQOÁ Á ÁCQO dihedral angle. 72 This indicates that the positive pole of the donor carbonyl group is not restricted in terms of its orientation. In order to further clarify whether a dipole is required for the intimate interaction between two carbonyl groups, the close contacts between a monopole (halide ion, X À1 ) and carbonyl group in the crystal structures of small molecules containing a halide ion and carbonyl group have been analyzed. 72 It was observed that the distance between the halide ion and the carbon of the carbonyl group, as well as the X À1 Á Á ÁCQO angle, fall mainly within the Burgi-Dunitz trajectory. This suggests that a dipole is not required for close contact between two carbonyl groups.
All the theoretical as well as experimental studies described above indicate that the nature of the np* Am interaction is neither charge-charge nor dipole-dipole, but is rather the delocalization of lone pair electrons on the donor atom into the antibonding orbital (p*) of the acceptor carbonyl group. Further, Raines and co-workers have studied the trans/cis ratio of the amide bond in N-formylproline phenylester to understand the energetics of the np* interaction. 42 Fig. 29 shows the cis and trans conformers of N-formylproline phenylester. The trans conformer, having an np* interaction, is more stable than the cis when X (substitution) is an electron withdrawing group, while the cis conformer is more stable than the trans when X is an electron donating group. It has been found that the substituents modulate the electrophilicity of C i QO i , which affects the strength of the np* interaction. This suggests that the electron withdrawing group decreases the energy of the p* orbital of C i QO i , making it more electrophilic, and thus allows closer overlap of the lone pair orbital of the oxygen with the p* orbital of the C i QO i carbonyl group. Thus, the strength of the np* interaction is affected substantially due to the electronic effect of the substituents attached close to the acceptor C i QO i group.
Raines and co-workers have also studied the effect on the strength of the np* interaction on changing the substituent attached to the donor carbonyl group (C iÀ1 QO iÀ1 ) and the nature of the donor atom. 37,69 The donor oxygen (O iÀ1 ) of N-acetylproline methylester was replaced by sulphur (S iÀ1 ) to give N-thioacetylproline methylester. NMR experiments show that the trans/cis ratio of N-thioacetylproline methylester is 2 times greater than that measured for N-acetylproline methylester. NBO analysis also shows that the interaction energy of the np* interaction is more when the donor atom is sulphur compared to when it is an oxygen atom. Sulphur, being a better electron pair donor, enhances the np* interaction. Thus, stronger electron pair donor atoms can enhance the np* interaction and favor one of the conformers over the other one.

Conclusions and outlook
In comparison to the commonly known non-covalent interactions in the literature, the np* interaction is quite new to the scientific community and it demands extensive exploration from experimental as well as theoretical points of view. The np* interaction, being very weak and counterintuitive in nature, was overlooked by scientists for a long time. We have reviewed here the recent discovery of this special non-covalent interaction and the rapid growth in its understanding through experiment and theory. It has been found that the np* interaction is widely present in biomolecules, i.e. proteins and nucleic acids, as well as in materials. This review describes the detailed literature studies on the propensity of the np* interaction in biomolecules and materials from PDB and CSD searches, respectively, the physical nature of this interaction, the interplay of this interaction with other non-covalent interactions, and experimental as well as theoretical evidences of the np* interaction.
It has been found that the np* interaction is as important as other non-covalent interactions for the stability of the structures of biomolecules and materials, although the former is generally very weak (B0.5-1 kcal mol À1 ) in strength. It is revealed from recent studies that the np* interaction can even modulate the overall structural motifs in the presence of strong hydrogen bonding interactions. Till now, experimental approaches to determine the presence of the np* interaction were limited to only NMR spectroscopy and X-ray crystallographic studies. NMR spectroscopy experiments on small model compounds show preferential population of the conformer having an np* interaction over the conformer without an np* interaction. X-ray crystallographic studies confirm the presence of the np* interaction in biomolecules and materials in terms of the Burgi-Dunitz trajectory, but there is no direct spectroscopic evidence for the np* interaction through probing the vibrational frequencies of the functional groups involved in the interaction. Generally, a shift in the X-H stretching frequency is a signature for the X-HÁ Á ÁY (both X and Y are generally electronegative atoms) hydrogen bond formation. In the future, it will be interesting to probe the stretching frequency of the carbonyl group in peptides or peptoids using IR spectroscopy to obtain direct spectroscopic proof for the np* interaction. As this non-covalent interaction is very weak compared to other relatively strong noncovalent interactions (i.e. hydrogen bonding), isolated gas phase spectroscopy in supersonic jet or matrix isolation infrared spectroscopy will be ideal to determine the strength of this interaction accurately. Although there are thorough literature studies on conventional hydrogen bonding, p-stacking and p-hydrogen bonding interactions in isolated gas phase, similar studies on the np* interaction are scarce in the literature, which is definitely due to unfamiliarity and the weak nature of this interaction. In fact, gas phase spectroscopy of small model compounds is highly desirable to check the results available on the np* interaction from quantum chemistry calculations.
There are very limited theoretical studies in the literature to determine the physical origin of the np* interaction. It has been demonstrated that an increase of dispersion in the system favors the np* interaction, overcoming repulsive electrostatic interactions. However, it cannot be stated with certainty that dispersion is the main source of stabilization for this interaction. Further theoretical studies, including energy decomposition analysis, are required for this purpose. In the future, theoretical studies should be aimed at molecular systems without bulky groups (i.e. alkyl) or extended conjugation, which contribute to dispersion. One recent work indeed shows that the dispersion contribution in an np* bound complex increases with increasing fluorine substitution, which itself does not introduce dispersion in the system. 34 This result indicates that the enhancement of the dispersion contribution could be due to the np* interaction.
Although the widespread presence of the np* interaction in biomolecules and materials is well documented in the literature, the functional role of this interaction in biological processes, i.e. protein-DNA interactions, protein-protein interactions, protein-ligand binding, and drug-DNA interactions, has not been explored yet. Thus, future studies should include detailed analysis of the crystal structures of protein-DNA, protein-protein complexes, etc. to trace the role of the np* interaction in these processes. It has been already pointed out that the np* interaction governs the chemical and biological activity of the drug aspirin, the regio-selective phosphorylation of anhydroribonucleoside, and the pharmacological activity of N-acyl homoserine lactones by inhibiting hydrolysis. 44,71,76 In a nutshell, detailed understanding of this interaction will be helpful in crystal engineering as well as drug design to obtain effective functional materials and drugs. The requirement of the development of force fields that will account for this counterintuitive interaction, leading to the refinement of the theoretical results of protein modelling, has been already recommended by Raines and co-workers. 41