Long-Range Electrostatics-Induced Two-Proton Transfer Captured by Neutron Crystallography in an Enzyme Catalytic Site

"Blade of Polarized Water Molecule" Is the Key to Hydrolase Catalysis Regulation

Hydrolysis catalyzed by aspartic proteases is a crucial reaction in many biological processes. However, anchoring water molecules and unifying multiple catalytic pathways remain significant challenges. Consequently, molecular design often compromises by focusing on enhancing substrate specificity. Using our self-developed polarizable point charge (PPC) force field, we determined the significant role of polarization in the hydrolase of pepsin for the first time. To be stably anchored in the active site, the water should be intensely polarized with a charge higher than -0.94e. Induced by this polarization, the pepsin was shown to support three general base/general acid pathways, with a preference for the gemdiol-intermediate-based pathway. Consequently, we proposed the "Blade of Polarized Water Molecule" model for rational enzyme design, highlighting that the polarization of both the attacking water and the attacked carbonyl is crucial for enhancing hydrolysis. Mutants D290Q and S172P showed activity enhancements of 191.23% and 324.70%, respectively. The improved polarization of water, carbonyl, and relevant nucleophilic attack distances in the mutants reaffirmed the crucial role of polarization in improving hydrolysis. This study provides a new perspective on hydrolase analysis and modification.

Electronic structure theory on modeling short-range noncovalent interactions between amino acids

While short-range noncovalent interactions (NCIs) are proving to be of importance in many chemical and biological systems, these atypical bindings happen within the so-called van der Waals envelope and pose an enormous challenge for current computational methods. We introduce SNCIAA, a database of 723 benchmark interaction energies of short-range noncovalent interactions between neutral/charged amino acids originated from protein x-ray crystal structures at the "gold standard" coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory with a mean absolute binding uncertainty less than 0.1 kcal/mol. Subsequently, a systematic assessment of commonly used computational methods, such as the second-order Møller-Plesset theory (MP2), density functional theory (DFT), symmetry-adapted perturbation theory (SAPT), composite electronic-structure methods, semiempirical approaches, and the physical-based potentials with machine learning (IPML) on SNCIAA is carried out. It is shown that the inclusion of dispersion corrections is essential even though these dimers are dominated by electrostatics, such as hydrogen bondings and salt bridges. Overall, MP2, ωB97M-V, and B3LYP+D4 turned out to be the most reliable methods for the description of short-range NCIs even in strongly attractive/repulsive complexes. SAPT is also recommended in describing short-range NCIs only if the δMP2 correction has been included. The good performance of IPML for dimers at close-equilibrium and long-range conditions is not transferable to the short-range. We expect that SNCIAA will assist the development/improvement/validation of computational methods, such as DFT, force-fields, and ML models, in describing NCIs across entire potential energy surfaces (short-, intermediate-, and long-range NCIs) on the same footing.

Cooperative protein-solvent tuning of proton transfer energetics: carbonic anhydrase as a case study

We investigate the coupling between the proton transfer (PT) energetics and the protein-solvent dynamics using the intra-molecular PT in wild type (wt) human carbonic anhydrase II and its ten-fold faster mutant Y7F/N67Q as a test case. We calculate the energy variation upon PT, and from that we also calculate the PT reaction free energy, making use of a hybrid quantum mechanics/molecular dynamics approach. In agreement with the experimental data, we obtain that the reaction free energy is basically the same in the two systems. Yet, we show that the instantaneous PT energy is on average lower in the mutant possibly contributing to the faster PT rate. Analysis of the contribution to the PT energetics of the solvent and of each protein residue, also not in the vicinity of the active site, provides evidence for electrostatic tuning of the PT energy arising from the combined effect of the solvent and the protein environment. These findings open up a way to the more general task of the rational design of mutants with either enhanced or reduced PT rate.

On the Case of the Misplaced Hydrogens

Few other elements play a more central role in biology than hydrogen. The interactions, bonding and movement of hydrogen atoms are central to biological catalysis, structure and function. Yet owing to the elusive nature of a single hydrogen atom few experimental and computational techniques can precisely determine its location. This is exemplified in short hydrogen bonds (SHBs) where the location of the hydrogen atom is indicative of the underlying strength of the bonds, which can vary from 1-5 kcal/mol in canonical hydrogen bonds, to an almost covalent nature in single-well hydrogen bonds. Owing to the often-times inferred position of hydrogen, the role of SHBs in biology has remained highly contested and debated. This has also led to discrepancies in computational, biochemical and structural studies of proteins thought to use SHBs in performing chemistry and stabilizing interactions. Herein, we discuss in detail two distinct examples, namely the conserved catalytic triad and the photoreceptor, photoactive yellow protein, where studies of these SHB-containing systems have permitted contextualization of the role these unique hydrogen bonds play in biology.

Opportunities and challenges in neutron crystallography

Neutron and X-ray crystallography are complementary to each other. While X-ray scattering is directly proportional to the number of electrons of an atom, neutrons interact with the atomic nuclei themselves. Neutron crystallography therefore provides an excellent alternative in determining the positions of hydrogens in a biological molecule. In particular, since highly polarized hydrogen atoms (H+) do not have electrons, they cannot be observed by X-rays. Neutron crystallography has its own limitations, mainly due to inherent low flux of neutrons sources, and as a consequence, the need for much larger crystals and for different data collection and analysis strategies. These technical challenges can however be overcome to yield crucial structural insights about protonation states in enzyme catalysis, ligand recognition, as well as the presence of unusual hydrogen bonds in proteins.

Direct Molecular Evidence of Proton Transfer and Mass Dynamics at the Electrode-Electrolyte Interface

Proton transfer has been widely regarded as a key step in many electrochemical and biological processes. However, direct molecular evidence has long been lacking. In this work, we chose the electrochemical oxidation of acetaminophen (APAP) as a model system and utilized in situ liquid time-of-flight secondary ion mass spectroscopy (ToF-SIMS) to molecularly examine proton solvation and transfer in this process. In addition, we successfully captured and identified the transient radical intermediate, providing solid molecular evidence to resolve an important debate in electron transfer-proton transfer oxidation mechanism of APAP. Moreover, the potential-dependent behaviors of both inert ions and electroactive species during the dynamic potential scanning were chemically monitored in real time and the mass diffusion mechanism regarding the electroactive and nonelectroactive species was revealed under polarized conditions. The results are consistent with our computer simulations. The observations in this work greatly improved our understanding of proton transfer and mass dynamics occurring at the electrode-electrolyte interface in complex electrochemical processes.

Proton in a Confined Space: Structural Studies of H+ ⊂Crypt-111 Iodide and Some Halogen-Bonded Derivatives

Experimental observations and modeling data are reported on the solid-state structural features of crypt- 111⋅HI (1) and the three-component co-crystals that 1 forms with α,ω-diiodoperfluoroalkanes 2 a-d. X-ray analyses indicate that, in all five systems and at low temperature, the caged proton is covalently bonded to a single nitrogen atom and is involved in a network of intramolecular hydrogen bonds. In contrast, room-temperature, solid-state ¹⁵ N NMR spectroscopy suggests magnetic equivalency of the two N atoms of crypt-111 in both 1 and co-crystals of 1 with diiodoperfluoroalkanes. Computational modelling confirms that the acidic hydrogen inside the cavity preferentially sits along the internitrogen axis and is covalently bonded to one nitrogen. The computed energy barriers suggest that the hopping of the encapsulated proton between the two N atoms of the cage can occur in the halogen-bonded co-crystals of 1⋅2, but it is hardly possible in the pure H⁺ ⊂crypt-111 iodide 1. These different pictures of the proton position and dynamics obtained by using different techniques and conditions confirm the unique characteristics of the confined space within the cavity of crypr-111 and the distinctive features of processes occurring therein.