Intramolecular Hydrogen Bonds for Brain-Penetrant Therapeutics

by Corin Wagen · Aug 6, 2024

The blood–brain barrier (BBB) is a semi-permeable membrane that separates the central nervous system (CNS) from the bloodstream. This serves to protect the CNS from toxins and other environmental harms, but also makes it difficult to deliver therapeutics to the CNS: it has been estimated that over 98% of small-molecule drugs are not BBB-penetrant.

Modifications to Increase BBB Penetrance

Extensive empirical studies have found a number of physicochemical properties which typically predict whether a given molecule will be BBB penetrant or not. A study of top CNS-active drugs found that virtually all of them had a low molecular weight, low total polar surface area, high lipophilicity, and—perhaps most crucially—at most one hydrogen-bond donor. These constraints make it extremely difficult to develop potent CNS therapeutics, as polar groups and hydrogen-bond donors are often key design elements that lead to high binding affinity and low IC50 values.

To help make potent drugs more CNS-penetrant without unduly sacrificing potency, hydrogen-bond donors can be "masked" by adding a weak internal hydrogen-bond acceptor to form an intramolecular hydrogen bond. This often increases the permeability of the molecule while maintaining binding affinity: the intramolecular hydrogen bond is robust enough to allow molecules to diffuse through the BBB, but not so strong that the molecule cannot adopt the preferred binding pose and engage a hydrogen-bond acceptor on the target.

Identifying a BBB-Penetrant IP6K1 Inhibitor

A recent paper from Barrow and co-workers at Johns Hopkins used this strategy to design an IP6K1 inhibitor capable of passing through the blood–brain barrier. IP6K1 is a kinase that's key to the inositol phosphate signaling pathway and has attracted considerable interest as a target for many different indications. Sadly, existing IP6K1 inhibitors all contain key carboxylic acids, which are ionized at physiological pH and thus prevent BBB penetration.

Replacing the carboxylic acid with an amide doesn't improve BBB penetrance

Replacing the carboxylic acid with an amide doesn't improve BBB penetrance

The authors tried a ton of different modifications to turn their existing IP6K1 inhibitor into a CNS-compatible IP6K1 inhibitor. Replacing the carboxylic acid with a carboxamide residue (which is neutral at pH 7.4) preserved most of the potency but led to even lower BBB penetrance, which the authors attribute to the additional hydrogen-bond donor moiety introduced.

Adding ortho-fluorine atoms dramatically improves permeability

Adding ortho-fluorine atoms dramatically improves permeability

The authors went on to explore introduction of an internal hydrogen bond to mask the carboxamide: while introducing a pyridine group improved permeability at the expense of potency, adding an ortho-fluorine preserved potency and dramatically increased BBB permeability for both the acid and the amide. Ultimately, they chose to move forward with the amide and found that it had a 173x greater brain/plasma ratio than the starting compound and good pharmacokinetic properties.

Modeling Intramolecular Hydrogen Bonding

This paper serves as a great case study of how introducing intramolecular hydrogen bonding can improve BBB penetrance. To understand the effects described in this paper in more detail, we can use Rowan to quickly compute the potential energy surface for the hydrogen-bond acceptor. Here's a comparison of scans for the parent phenyl compound, the pyridyl compound, and the ortho-fluoro compound.

A comparison of different hydrogen-bond acceptors; click to view the full interactive page

We can see that in the absence of a hydrogen-bond acceptor the s-cis carboxylic acid conformation is strongly preferred. Addition of a strong hydrogen-bond acceptor (pyridine) shifts the ground state towards the s-trans carboxylic acid conformation, consistent with the increased BBB penetrance and decreased potency observed. The ortho-fluorine serves as a weak hydrogen-bond acceptor, making the intramolecular hydrogen-bonded conformation energetically accessible but maintaining the same ground state.

These calculations suggest that there's a "Goldilocks" zone of hydrogen-bond-acceptor strength for intramolecular bonds: acceptors need to be strong enough to mask hydrogen-bond donors in a membrane environment, but not so strong that they compete in the binding site. This chameleonicity has been remarked upon as a key feature of macrocycles, enabling them to change their physicochemical properties in an environment-specific manner.

Prospective Design of BBB-Penetrant Therapeutics

This study, while retrospective, illustrates how calculations can be used for the rational design of CNS-compatible small molecules. The calculations run here took only minutes through Rowan's cloud platform, and the figures were generated automatically through Rowan's molecule viewer.

Rowan provides an optimized workflow for running scans that's ideally suited to running calculations like the ones shown here, where a quick estimate of the relative energies of two states is needed. Rowan also contains cutting-edge machine-learned interatomic potentials that can generate torsional profiles with chemical accuracy in a fraction of the time demanded by traditional quantum chemical methods, meaning that your scans can run in minutes instead of days.

Banner background image

What to Read Next

Hydrogen-Bond-Basicity Predictions for Scaffold Hopping in PDE2A Inhibitors

Hydrogen-Bond-Basicity Predictions for Scaffold Hopping in PDE2A Inhibitors

How new computational workflows can make it possible to design complex modifications to heterocyclic cores.
Feb 14, 2025 · Corin Wagen
Intrinsic Reaction Coordinates

Intrinsic Reaction Coordinates

expanding Rowan's reaction modeling toolkit; verifying transition states; reaction mechanism insights
Feb 6, 2025 · Jonathon Vandezande and Ari Wagen
Reactions from the Bottom Up

Reactions from the Bottom Up

Building up an understanding of how energy barriers and the potential energy surface affect the rate of a reaction.
Feb 4, 2025 · Jonathon Vandezande
A New RDKit-Native API

A New RDKit-Native API

cultural barriers in science; integrating RDKit with quantum chemistry; Rowan's new API; changes to billing
Jan 31, 2025 · Corin Wagen and Spencer Schneider
Hydrogen-Bond Basicity Prediction Made Easy

Hydrogen-Bond Basicity Prediction Made Easy

not all hydrogen-bond donors are created equal; the pKBHX scale; predicting pKBHX in Rowan; case studies & a preprint
Jan 24, 2025 · Corin Wagen
Efficient Black-Box Prediction of Hydrogen-Bond-Acceptor Strength

Efficient Black-Box Prediction of Hydrogen-Bond-Acceptor Strength

Here, we report a robust black-box workflow for predicting site-specific pKBHX values in organic molecules with minimal computational cost.
Jan 24, 2025 · Corin C. Wagen
Benchmarking NNPs, Orb-v2, and MACE-MP-0

Benchmarking NNPs, Orb-v2, and MACE-MP-0

benchmarking as driver of systematic methodological improvement; our new benchmarking website; new NNPs on Rowan; GPU-based inference coming to more users
Jan 17, 2025 · Ari Wagen
Density-Functional-Theory Functionals Quiz

Density-Functional-Theory Functionals Quiz

Ready to test your knowledge of density-functional-theory functionals in a multiple-choice game of "real or fake"?
Jan 10, 2025 · Jonathon Vandezande
Wiggle150: Benchmarking Density Functionals And Neural Network Potentials On Highly Strained Conformers

Wiggle150: Benchmarking Density Functionals And Neural Network Potentials On Highly Strained Conformers

We introduce Wiggle150, a benchmark comprising 150 highly strained conformations of adenosine, benzylpenicillin, and efavirenz, to validate computational protocols involving non-equilibrium systems and guide the development of new density functionals and neural network potentials.
Jan 8, 2025 · Joseph Gair, Corin Wagen, et al., ChemRxiv
The "Charlotte's Web" of Density-Functional Theory

The "Charlotte's Web" of Density-Functional Theory

A layman's guide to cutting your way through the web of DFT functionals, explaining GGAs, mGGAs, hybrids, range-separated hybrids, double hybrids, and dispersion corrections.
Dec 20, 2024 · Jonathon Vandezande