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

by Corin Wagen · Feb 14, 2025

In a 2018 study (J. Med. Chem. 2018, 61, 1001–1018), Chris Helal and co-workers at Pfizer reported the development of phosphodiesterase 2A (PDE2A) inhibitors as potential treatments for cognitive disorders in schizophrenia. Their core scaffold—pyrazolopyrimidine 1—exhibited good potency but remained relatively lipophilic, leading to high human-liver-microsome clearance and thus an excessively high estimated dose.

The authors investigated peripheral modifications to decrease lipophilicity, but found that the perturbations they tried either compromised activity or led to increased values of MDR BA/AB. (MDR BA/AB is the ratio of basolateral-to-apical flux to apical-to-basolateral flux of a compound in MDCK cells engineered to express P-glycoprotein (MDR1), reflecting P-gp–mediated efflux). Accordingly, the Pfizer team next investigated re-engineering the central heterocyclic core to preserve potency while improving lipophilicity and efflux.

A few different peripheral modifications.

Peripheral modifications to 1 were unsuccessful.

LMP2 Calculations Identify a More Potent Scaffold

Inspired by results from Bayer showing that different heterocyclic isomers could lead to 50x potency improvements for PDE5A, the team explored an imidazotriazine ring to replace the pyrazolopyrimidine scaffold. Counterpoise-corrected LMP2/cc-pVTZ//X3LYP/6-31G** calculations in Jaguar indicated that key hydrogen-bond interactions in the enzyme's active site would be strengthened with the imidazotriazine core. This prediction was borne out experimentally: after extensive optimization, the new scaffold led to the clinical candidate PF-05180999, which demonstrated higher PDE2A affinity and improved brain penetration.

Illustration of the PDE2A binding site.

Illustration of the inhibitor class bound to PDE2A and the key hydrogen bonds made with the protein (PDB 5UZL).

pKBHX Workflow as a Black-Box Alternative

Unfortunately, the LMP2-based calculations run by Helal and co-workers are challenging for non-expert practitioners to perform accurately: proper setup of these calculations requires properly aligning new ligands with experimental or predicted protein–ligand complexes, deleting all but the relevant regions of the protein, and then hoping that the generated structure is both small enough to enable practical MP2 calculations and large enough to capture all relevant interactions.

Rowan's hydrogen-bond-basicity-prediction workflow (pKBHX) offers a more accessible way to estimate site-specific hydrogen-bond acceptor strengths. Whereas LMP2-based workflows require specialized resources and expertise, the pKBHX approach uses only a single density-functional-theory calculation per molecule, automatically handles different conformers, and has been calibrated and validated against a large experimental database. This method identifies the most strongly hydrogen-bonding heteroatoms, enabling medicinal chemists to prioritize promising scaffold hops without advanced quantum-chemistry infrastructure.

Agreement with High-Level QM Calculations

In an independent analysis of the same scaffold hop, pKBHX calculations confirmed that the imidazotriazine ring would generally strengthen the critical hydrogen bond to the five-membered ring (N2 when using the pyrazolopyrimidine numbering) compared to the original pyrazolopyrimidine core: Rowan predicts an increase in pKBHX of 0.88 units (almost an order of magnitude), while Helal and co-workers predict that the hydrogen bond will be 1.4 kcal/mol stronger.

Both the LMP2-based computations and the pKBHX results also agree qualitatively on the overall increase in hydrogen-bond strengths upon scaffold modification. The two methods disagree on the effect of the scaffold hop on N5, suggesting that in this case side-chain and water-mediated interactions can slightly alter outcomes for ligand-only predictions. Nonetheless, the overall conclusion is the same: substituting imidazotriazine yields a net improvement in hydrogen bonding and potency.

ChemDraw comparison of compounds 10 and 11

Comparison of the hydrogen-bond-acceptor strength of the pyrazolipyrimidine 10 and the imidazotriazine 11. More negative LMP2 values imply stronger hydrogen bonding; higher pKBHX values indicate stronger hydrogen-bond acceptors.

This case study demonstrates how modeling and simulation tools like the pKBHX workflow can help make complex decisions like scaffold hops more data-driven, accelerating project timelines and decreasing the role of luck in hit-to-lead optimization. While top pharmaceutical companies like Pfizer have long relied on sophisticated computational workflows to help make these high-stakes decisions (as shown by this paper), it's now possible for even small teams to take advantage of these insights.

Rowan is working to democratize computational chemistry and bring advanced design & simulation tools to every scientist, not just computational experts. If you're interested in working with us, reach out to our team at contact@rowansci.com or start using our web-based computational chemistry platform today!

Banner background image

What to Read Next

Partnering with Macrocosmos to Accelerate Next-Generation NNP Development

Partnering with Macrocosmos to Accelerate Next-Generation NNP Development

Starting today, Rowan is teaming up with Macrocosmos to accelerate the development of the next generation of NNPs through Bittensor Subnet 25 - Mainframe.
May 1, 2025 · Ari Wagen
Introducing Egret-1

Introducing Egret-1

trusting computation; speed vs accuracy; Egret-1, Egret-1e, and Egret-1t; benchmarks; speed on CPU and GPU; download Egret-1 or use it through Rowan
Apr 30, 2025 · Eli Mann, Corin Wagen, Jonathon Vandezande, Ari Wagen, and Spencer Schneider
Egret-1: Pretrained Neural Network Potentials For Efficient and Accurate Bioorganic Simulation

Egret-1: Pretrained Neural Network Potentials For Efficient and Accurate Bioorganic Simulation

Here, we present Egret-1, a family of large pre-trained NNPs based on the MACE architecture with general applicability to main-group, organic, and biomolecular chemistry.
Apr 30, 2025 · Elias L. Mann, Corin C. Wagen, Jonathon E. Vandezande, Arien M. Wagen, Spencer C. Schneider
Introducing Egret-1

Introducing Egret-1

Today, we're releasing Egret-1, a family of open-source NNPs for bioorganic simulation.
Apr 30, 2025 · Eli Mann, Corin Wagen, Jonathon Vandezande, Ari Wagen, and Spencer Schneider
Starling: Macroscopic pKa, logD, and Blood–Brain-Barrier Permeability

Starling: Macroscopic pKa, logD, and Blood–Brain-Barrier Permeability

microscopic vs. macroscopic pKa; Uni-pKa and Starling; microstate ensembles; logD and Kp,uu predictions
Apr 25, 2025 · Corin Wagen
Physics-Informed Machine Learning Enables Rapid Macroscopic pKa Prediction

Physics-Informed Machine Learning Enables Rapid Macroscopic pKa Prediction

Here we introduce Starling, a physics-informed neural network based on the Uni-pKa architecture trained to predict per-microstate free energies and compute macroscopic pKa values via thermodynamic ensemble modeling.
Apr 25, 2025 · Corin C. Wagen
Predicting Infrared Spectra and Orb-v3

Predicting Infrared Spectra and Orb-v3

light and its manifold interactions with matter; why IR spectroscopy is useful; predicting IR spectra through Rowan; Orb-v3
Apr 17, 2025 · Ari Wagen, Corin Wagen, and Jonathon Vandezande
What's in a Name?

What's in a Name?

Why our company is named after a tree with no obvious connection to what we do.
Apr 11, 2025 · Corin Wagen and Ari Wagen
2D Scans, More Fukui Options, & Max Concurrency

2D Scans, More Fukui Options, & Max Concurrency

escape from scan-land; making More O'Ferrall–Jencks plots in Rowan; extending Fukui calculations to more levels of theory; managing large numbers of concurrent jobs
Apr 10, 2025 · Corin Wagen, Ari Wagen, Spencer Schneider, and Jonathon Vandezande
Can AI Can Accelerate Scientific Research?

Can AI Can Accelerate Scientific Research?

New research from Aiden Toner-Rodgers at MIT suggests the answer is yes, but with important caveats.
Apr 2, 2025 · Corin Wagen