Electronic Binding Features

• Electronic Binding Features are quantitative descriptors that capture electronic contributions, including polarization, protonation, and electron–hole interactions in binding phenomena.
• In excitonic systems, they are measured using techniques like scanning tunneling spectroscopy and photoelectrochemical methods to determine exciton binding energies and related parameters.
• In biomolecular contexts, methods such as alchemical free energy simulations and MBAR/PBSA incorporate corrections for polarization and protonation, improving binding affinity predictions.

Electronic binding features (BFs) denote quantitative descriptors of the electronic contributions to binding free energies in molecular association phenomena. These features, central to both solid-state excitonic systems and biomolecular recognition, reflect underlying mechanisms such as electronic polarization, protonation equilibria, and electron–hole correlation (excitons), producing binding energy terms essential to the accurate modeling and interpretation of electronic structure, spectroscopy, and binding affinity. In alchemical free energy simulations, BFs are defined as explicit free-energy components or corrections capturing the effects of environment-induced electronic polarization and protonation state changes; in condensed matter, excitonic binding energies discern the influence of electron–hole attraction. Experimental and computational determination of electronic BFs enables rigorous ranking, prediction, and engineering of molecular and material properties in both optoelectronic and biochemical contexts (Molina-Mendoza et al., 2015, King et al., 2021).

1. Fundamental Definitions and Context

Electronic binding features comprise the set of contributions to binding that arise from electronic structure effects beyond fixed classical interactions. In molecular association, primary examples include:

• Electronic polarization: Redistribution of atomic charges or multipoles in response to local dielectric environment, altering the electrostatic component of binding free energy.
• Protonation equilibrium: Influence of alternative protonation states on the net binding affinity, especially for titratable groups sensitive to microenvironment pK_a shifts.

In the solid-state domain, excitonic binding energy measures the Coulomb attraction between an electron and a hole, distinguishing the single-particle electronic bandgap ($E_{\mathrm{gap}}^{\mathrm{elec}}$) from the photon energy at the onset of optical absorption ($E_{\mathrm{gap}}^{\mathrm{opt}}$):

$$E_b = E_{\mathrm{gap}}^{\mathrm{elec}} - E_{\mathrm{gap}}^{\mathrm{opt}}$$

(Molina-Mendoza et al., 2015).

In alchemical binding free energy frameworks for biomolecular complexes,

$$\Delta G_\mathrm{bind} = \Delta G_\mathrm{vdw} + \Delta G_\mathrm{elec} + \Delta G_\mathrm{restr} + \Delta G_\mathrm{pol} + \dots$$

where $\Delta G_\mathrm{pol}$ quantifies the polarization correction and $\Delta G_\mathrm{prot}$ quantifies protonation-dependent binding shifts (King et al., 2021).

2. Experimental and Computational Determination

Accurate extraction of electronic BFs necessitates precise experimental and theoretical methodologies.

• Solid-State Excitonic Systems (e.g., TiS$_3$)
  • Electronic bandgap measurement: Scanning tunneling spectroscopy (STS) quantifies $E_{\mathrm{gap}}^{\mathrm{elec}}$ via the bias voltage range yielding zero differential conductance.
  • Optical bandgap measurement: Photoelectrochemical (PEC) methods yield $E_{\mathrm{gap}}^{\mathrm{opt}}$ from Tauc plots of photocurrent as a function of photon energy.
  • Binding energy extraction: The difference yields the exciton binding energy $E_b$. Thermal broadening corrections (e.g., $\Delta E_\mathrm{th} \approx 3.5 k_BT$) ensure fidelity of STS bandgaps (Molina-Mendoza et al., 2015).
• Biomolecular Binding (e.g., protein–ligand complexes)
  • Alchemical free energy protocols: Decouple ligand interactions in both complex and solvent, capturing $\Delta G_\mathrm{bind}$ via molecular dynamics.
  • MBAR/PBSA post-processing: Multistate Bennett acceptance ratio (MBAR) applied to continuum Poisson–Boltzmann (PB) snapshots enables scalable scanning of solute dielectric $\epsilon_\mathrm{in}$, yielding electronic polarization corrections.
  • Protonation enumeration: Explicit sampling and ranking of all relevant ligand and protein (e.g., histidine) protonation microstates quantify $\Delta G_\mathrm{prot}$ shifts (King et al., 2021).

Table: Core Experimental and Computational Measures

Domain	Feature	Measurement / Computation
Excitonic Semiconductors	$E_b$, $E_{\text{gap}}$	STS, PEC, GW/BSE theory
Protein–ligand complexes	$\Delta G_\mathrm{pol}$, $\Delta G_\mathrm{prot}$	Alchemical MD, MBAR/PBSA, protonation-state enumeration

3. Mathematical Formalism

In solid-state and molecular simulation contexts, electronic binding features are formalized by specific equations:

• Exciton binding energy:

$$E_b = E_{\mathrm{gap}}^{\mathrm{elec}} - E_{\mathrm{gap}}^{\mathrm{opt}}$$

with $E_{\mathrm{gap}}^{\mathrm{elec}}$ (electronic bandgap) and $E_{\mathrm{gap}}^{\mathrm{opt}}$ (optical bandgap) as defined above (Molina-Mendoza et al., 2015).

• Polarization correction in binding free energy:

$$\Delta G_\mathrm{pol}(\epsilon_\mathrm{in}) = \Delta G_\mathrm{decharge}^{\mathrm{PB}}(\epsilon_\mathrm{in}) - \Delta G_\mathrm{decharge}^{\mathrm{PB}}(\epsilon_\mathrm{in} = 1)$$

where $\epsilon_\mathrm{in}$ is the solute dielectric evaluated in the PB model (King et al., 2021).

• Protonation-state shift:

$$\Delta G_\mathrm{prot} = \Delta G_\mathrm{bind}(\mathrm{HID}) - \Delta G_\mathrm{bind}(\mathrm{HIP})$$

capturing the energetic penalty for protonation state transformations of key residues or ligands (King et al., 2021).

Analytic corrections for orientational and positional restraints (e.g., Boresch-style) further refine the absolute free energy estimates.

4. Numerical Results and Comparative Analysis

The magnitude of electronic binding features varies sharply across material and molecular contexts:

• Excitonic Binding Energy (TiS$_3$):
  • Experimental: $E_{\mathrm{gap}}^{\mathrm{elec}} = 1.20 \pm 0.08$ eV, $E_{\mathrm{gap}}^{\mathrm{opt}} = 1.07 \pm 0.01$ eV, $E_b = 0.13$ eV.
  • Theoretical: $E_{\mathrm{gap}}^{\mathrm{elec}} = 1.15$ eV, $E_{\mathrm{gap}}^{\mathrm{opt}} = 1.05$ eV, $E_b = 0.10$ eV.
  • Comparison: $E_b$ of TiS$_3$ ribbons far exceeds classical semiconductors (Si, Ge, GaAs: 1–60 meV); it is comparable to bulk transition-metal dichalcogenides (70 meV) and intermediate between bulk and monolayer TMDCs (up to 900 meV) (Molina-Mendoza et al., 2015).
• Polarization and Protonation BFs (UPA–inhibitor binding):
  • Polarization correction $\Delta G_\mathrm{pol}$: Up to several kcal/mol depending on $\epsilon_\mathrm{in}$, with optimized fits (e.g., all-HID at $\epsilon_\mathrm{in}=1.43$, RMSE = 0.89 kcal/mol) providing significantly improved agreement with experiment (King et al., 2021).
  • Protonation-state dependent $\Delta G_\mathrm{prot}$: Shifts between binding modes (e.g., HID vs. HIP), with both ranking and RMSE depending critically on proper enumeration.

This suggests that proper extraction and inclusion of electronic BFs are critical for both spectroscopy of low-dimensional materials and quantitative binding affinity prediction in molecular recognition.

5. Physical Significance and Applications

Electronic binding features have decisive impact on macroscopic observables and device-level applications:

• Room-temperature excitonic vs. electronic couplings: For semiconductors such as TiS$_3$, $E_b \gg k_BT$ at 300 K assures exciton stability, enabling strong light–matter coupling, narrow absorption, and emission lines—intermediating photocarrier dynamics, trions, and biexcitons in optoelectronics (Molina-Mendoza et al., 2015).
• Protein–ligand binding selectivity and prediction: Inclusion of $\Delta G_\mathrm{pol}$ and $\Delta G_\mathrm{prot}$ as explicit BFs dramatically improves the accuracy, ranking, and interpretability of simulation-based binding energy predictions, especially for highly charged or titratable systems (King et al., 2021).

Potential applications follow accordingly:

• Excitonic devices (photodetectors, flexible solar cells, valley polarization devices)
• Drug design pipelines leveraging BFs as machine learning descriptors, continuum dielectric corrections, and explicit protonation enumeration for improved hit ranking and lead optimization.

6. Methodological Recommendations and Outlook

Findings from recent research advocate adopting the following practices for rigorous inclusion of electronic binding features:

• In optoelectronic materials:
  • Combine experimental STS and PEC measurements with GW/BSE computations for maximal fidelity when extracting $E_{\mathrm{gap}}^{\mathrm{elec}}$, $E_{\mathrm{gap}}^{\mathrm{opt}}$, and $E_b$.
• In biomolecular simulation:
  • Employ explicit MBAR/PBSA processing with varied $\epsilon_\mathrm{in}$ for polarization BFs, ensuring solute radii calibration.
  • Systematically enumerate and include alternative ligand and protein protonation states; carry through and report $\Delta G_\mathrm{prot}$ as a predictive feature.
  • Where computationally feasible, implement fully polarizable force fields (e.g., AMOEBA, Drude oscillators) for highest accuracy, or at minimum apply mean-field dielectric corrections.
  • Include extracted BFs ($\Delta G_\mathrm{pol}$, $\Delta G_\mathrm{prot}$) as descriptors in machine learning-based scoring and lead identification models (King et al., 2021).

A plausible implication is that future electronic binding feature sets will serve as standardized descriptors across fields, linking fundamental quantum chemistry, condensed-matter physics, and computational biophysics in the quantification of electronically driven binding phenomena.