Research Interests
This research program uses modern methods of electronic structure theory to understand the fundamental nature of interactions between molecules. Chief among the concerns are hydrogen bonds which are so important to structure and function of biomolecules like proteins.
When considering a H-bond, the thoughts of most would ordinarily turn to OH··O or OH··N interactions that pair electronegative atoms like O and N. Indeed, such bonds are quite common in proteins and act as a strong glue, cementing one part of the molecule to another. On the other hand, clues have turned up over the years that even less electronegative atoms like C can participate in H-bonds, albeit weaker ones. Most of these clues have been of the structural variety, wherein a CH group might lie in close proximity to an O atom, and the geometry of the contact might be reminiscent of a H-bond, i.e. the putative bridging H could lie along the direction of one of the lone electron pairs of the oxygen acceptor. Yet structural information of this sort is incapable of distinguishing a "coincidental" contact, where the two groups happen to lie close together for the sake of the structure of the entire molecule, from a true CH··O H-bond. The latter would be accompanied by an attractive interaction, a "fingerprint" electron redistribution pattern, and a number of characteristic spectroscopic features.
Our group makes use of state-of-the-art ab initio quantum chemical methods to probe the underlying nature of the CH··O interaction, paying particular attention to the above properties, extracting information that is inaccessible to experimental measurements. We have learned that under certain conditions, the CH··O interaction is indeed a true H-bond, with an attractive force that can rival the more conventional OH··O bond. And even though the C-H stretching frequency may shift to the blue rather than to lower frequencies, this property is typical of certain weak H-bonds.
The work has demonstrated that these weak H-bonds, commonly overlooked, can play a very important structural role in some of the most important and widely occurring protein structural motifs, such as beta sheets. We have immersed these systems in environments that are characteristic of aqueous solvent or the interior of a protein, and found that many of our findings in the gas phase remain valid. By comparing H-bond strengths in solvents of varying polarizability, it has been possible to estimate the magnitude of the so-called “desolvation penalty” which characterizes the weakening suffered by a H-bond upon transferring from aqueous solvent to the interior of the protein.
And with particular regard to conventional H-bonds, of the NH··O sort that connects a pair of peptide units within a protein, calculations have shown that the strength with which such a bond can form is surprisingly dependent upon the specific conformation of the polypeptide chain on which the bond occurs. For example, the NH··O bond within a stretched polypeptide, as in a β-sheet, is considerably weaker than if the same bond were to occur within a helical portion of the protein.