Search papers, labs, and topics across Lattice.
This paper introduces a novel unbiased Monte Carlo method for estimating Gaussian molecular shape overlap, addressing the limitations of existing analytical methods that tend to overestimate overlap due to first-order approximations. The new stochastic estimator not only reproduces the analytic overlap accurately but also computes the exact union volume of Gaussian mixtures with linear computational cost per sample. On drug-like molecules, it achieves a mean relative error of just 0.07% compared to high-resolution methods, while significantly reducing sampling needs and providing confidence bounds for virtual screening.
A new Monte Carlo estimator for Gaussian shape overlap achieves unprecedented accuracy and efficiency, cutting sampling needs by 94% while delivering reliable uncertainty estimates.
Gaussian descriptions of molecular shape underpin 3D shape-based virtual screening, but existing methods evaluate Gaussian overlap analytically. The widely used first-order approximation is fast but systematically overestimates overlap, whereas the exact molecular volume requires a combinatorial inclusion-exclusion expansion. We introduce the first stochastic estimator of Gaussian shape overlap: an unbiased Monte Carlo method that importance-samples directly from a molecule's Gaussian mixture. The estimator reproduces analytic overlap without bias and extends to the exact union volume of all inclusion-exclusion orders with O(N) cost per sample. On drug-like molecules, the union estimator matches high-resolution grid quadrature with a mean relative error of 0.07 percent, while the first-order approximation overestimates the true union volume by 3.4x on average. The estimator provides analytic standard errors, enabling confidence-bounded screening that reduces sampling by 94 percent while preserving ranking. Implemented in JAX, it is fully differentiable and supports gradient-based rigid alignment on CPU, GPU, and TPU. On the DUD-E and LIT-PCBA benchmarks, the method achieves shape-only enrichment comparable to existing single-conformer approaches while additionally providing unbiased absolute volumes and uncertainty estimates.
