Jeremy G. Vinter

Substituent effects on aromatic stacking interactions

Synthetic supramolecular zipper complexes have been used to quantify substituent effects on the free energies of aromatic stacking interactions. The conformational properties of the complexes have been characterised using NMR spectroscopy in CDCl(3), and by comparison with the solid state structures of model compounds. The structural similarity of the complexes makes it possible to apply the double mutant cycle method to evaluate the magnitudes of 24 different aromatic stacking interactions. The major trends in the interaction energy can be rationalised using a simple model based on electrostatic interactions between the pi-faces of the two aromatic rings. However, electrostatic interactions between the substituents of one ring and the pi-face of the other make an additional contribution, due to the slight offset in the stacking geometry. This property makes aromatic stacking interactions particularly sensitive to changes in orientation as well as the nature and location of substituents.

Substituent effects on cation–π interactions: A quantitative study

A synthetic supramolecular complex has been adapted to quantify cation–π interactions in chloroform by using chemical double-mutant cycles. The interaction of a pyridinium cation with the π-face of an aromatic ring is found to be very sensitive to the π-electron density. Electron-donating substituents lead to a strong attractive interaction (−8 kJ/mol−1), but electron-withdrawing groups lead to a repulsive interaction (+2 kJ/mol−1).

FieldScreen: virtual screening using molecular fields. Application to the DUD data set.

FieldScreen, a ligand-based Virtual Screening (VS) method, is described. Its use of 3D molecular fields makes it particularly suitable for scaffold hopping, and we have rigorously validated it for this purpose using a clustered version of the Directory of Useful Decoys (DUD). Using thirteen pharmaceutically relevant targets, we demonstrate that FieldScreen produces superior early chemotype enrichments, compared to DOCK. Additionally, hits retrieved by FieldScreen are consistently lower in molecular weight than those retrieved by docking. Where no X-ray protein structures are available, FieldScreen searches are more robust than docking into homology models or apo structures.

Noncovalent Assembly of [2]Rotaxane Architectures

Reversible zinc-pyridine coordination and hydrogen-bonding interactions have been used to assemble a [2]rotaxane from three components. Cooperativity in the macrocyclization process that results in the porphyrin dimer makes the system exceptionally stable. However, the kinetic lability of the zinc-porphyrin interaction means the dimer is in dynamic equilibrium with its monomer, and this has been exploited in the construction of a [2]rotaxane.

An evaluation of force-field treatments of aromatic interactions

Experimental measurements of edge-to-face aromatic interactions have been used to test a series of molecular mechanics force fields. The experimental data were determined for a range of differently substituted aromatic rings using chemical double mutant cycles on hydrogen-bonded zipper complexes. These complexes were truncated for the purposes of the molecular mechanics calculations so that problems of conformational searching and the optimisation of large structures could be avoided. Double-mutant cycles were then carried out in silico using these truncated systems. Comparison of the experimental aromatic interaction energies and the X-ray crystal structures of these truncated complexes with the calculated data show that conventional molecular mechanics force fields (MM2, MM3, AMBER and OPLS) do not perform well. However, the XED force field which explicitly represents electron anisotropy as an expansion of point charges around each atom reproduces the trends in interaction energy and the three-dimensional structures exceedingly well. Collapsing the XED charges onto atom centres or the use of semi-empirical atom-centred charges within the XED force field gives poor results. Thus the success of XED is not related to the methods used to assign the atomic charge distribution but can be directly attributed to the use of off-atom centre charges.

Quantitative determination of intermolecular interactions with fluorinated aromatic rings.

The chemical double mutant cycle approach has been used to investigate substituent effects on intermolecular interactions between aromatic rings and pentafluorophenyl pi-systems. The complexes have been characterised using 1H and 19F NMR titrations, X-ray crystal structures of model compounds and molecular mechanics calculations. In the molecular zipper system used for these experiments, H-bonds and the geometries of the interacting surfaces favour the approach of the edge of the aromatic ring with the face of the pentafluorophenyl pi-system. The interactions are generally repulsive and this repulsion increases with more electron-withdrawing substituents up to a limit of +2.2 kJ mol(-1), when the complex distorts to minimise the unfavourable interaction. Strongly electron-donating groups cause a change in the geometry of the aromatic interaction and attractive stacking interactions are found (-1.6 kJ mol(-1) for NMe2). These results are generally consistent with an electrostatic model: the polarisation of the pentafluorophenyl ring leads to a partial positive charge located at the centre and this leads to repulsive interactions with the positive charges on the protons on the edge of the aromatic ring; when the aromatic ring has a high pi-electron density there is a large electrostatic driving force in favour of the stacked geometry which places this pi-electron density over the centre of the positive charge on the pentafluorophenyl group.

Preferential Solvation and Hydrogen Bonding in Mixed Solvents

The molecular-recognition events that underpin many processes in chemistry and biology are intimately linked to accompanying changes in solvation. The subtle interplay of the different factors that contribute to these changes in solvation gives rise to phenomena that are difficult to interpret at the molecular level. However, a detailed understanding of these effects is essential for the development of quantitative approaches to the analysis of noncovalent molecular interactions in solution. We recently introduced the 1:1 complex formed between tri-n-butylphosphine oxide (1) and perfluoro-tert-butanol (2 ; Scheme 1) as a

Rationalizing the activities of diverse cholecystokinin 2 receptor antagonists using molecular field points.

Cholecystokinin 2 receptor antagonists encompass a wide range of structures. This makes them unsuitable candidates for existing 3D-QSAR methods and has led us to develop an alternative approach to account for their observed biological activities. A diverse set of 21 antagonists was subjected to a novel molecular field-based similarity analysis. The hypothesis is that compounds with similar field patterns will bind at the same target site regardless of their underlying structure. This initial report demonstrates a linear correlation between ligand similarity and biological activity for this challenging data set. A model generated with three molecules was used to predict the activity of 18 test compounds, with different chemotypes, with a root-mean-square error of 0.68 pKB units. The ability to automatically derive a molecular alignment without knowledge of the protein structure represents an improvement over existing pharmacophore methods and makes the method particularly suitable for scaffold-hopping.

Influence of solvent polarity on preferential solvation of molecular recognition probes in solvent mixtures.

The association constants for formation of 1:1 complexes between a H-bond acceptor, tri-n-butylphosphine oxide, and a H-bond donor, 4-phenylazophenol, have been measured in a range of different solvent mixtures. Binary mixtures of n-octane and a more polar solvent (ether, ester, ketone, nitrile, sulfoxide, tertiary amide, and halogenated and aromatic solvents) have been investigated. Similar behavior was observed in all cases. When the concentration of the more polar solvent is low, the association constant is identical to that observed in pure n-octane. Once a threshold concentration of the more polar solvent in reached, the logarithm of the association constant decreases in direct proportion to the logarithm of the concentration of the more polar solvent. This indicates that one of the two solutes is preferentially solvated by the more polar solvent, and it is competition with this solvation equilibrium that determines the observed association constant. The concentration of the more polar solvent at which the onset of preferential solvation takes place depends on solvent polarity: 700 mM for toluene, 60 mM for 1,1,2,2-tetrachloroethane, 20 mM for the ether, ester, ketone, and nitrile, 0.2 mM for the tertiary amide, and 0.1 mM for the sulfoxide solvents. The results can be explained by a simple model that considers only pairwise interactions between specific sites on the surfaces of the solutes and solvents, which implies that the bulk properties of the solvent have little impact on solvation thermodynamics.

The role of functional group concentration in solvation thermodynamics

High throughput experiments using a molecular recognition probe reveal a simple relationship between solvent functional group concentration and selective solvation. The 1 : 1 association constant for the H-bonding interaction between tri-n-butylphosphine oxide and 4-phenylazophenol was measured in 1088 different alkane–ether mixtures using a UV-Vis plate reader. Although the stability of the complex decreased with increasing concentration of the more polar ether cosolvent as expected, the results show that it is the functional group composition rather than the constitution of the solvent molecules or the properties of the bulk liquid that determines the solvation thermodynamics. Thus the solvent properties of a simple ether can be reproduced by an appropriate mixture of a polyether and an alkane that has the same net concentration of ether oxygen functional groups. The results suggest that solvation may be understood at the molecular level simply by considering the polarities and the concentrations of the functional groups present in the solvent, because these are the parameters that affect local solvation interactions with the solutes.