Ligand-Induced Proton Transfer and Low-Barrier Hydrogen Bond Revealed by X-ray Crystallography.

Ligand binding can change the pKa of protein residues and influence enzyme catalysis. Herein, we report three ultrahigh resolution X-ray crystal structures of CTX-M ß-lactamase, directly visualizing protonation state changes along the enzymatic pathway: apo protein at 0.79 Å, precovalent complex with nonelectrophilic ligand at 0.89 Å, and acylation transition state (TS) analogue at 0.84 Å. Binding of the noncovalent ligand induces a proton transfer from the catalytic Ser70 to the negatively charged Glu166, and the formation of a low-barrier hydrogen bond (LBHB) between Ser70 and Lys73, with a length of 2.53 Å and the shared hydrogen equidistant from the heteroatoms. QM/MM reaction path calculations determined the proton transfer barrier to be 1.53 kcal/mol. The LBHB is absent in the other two structures although Glu166 remains neutral in the covalent complex. Our data represents the first X-ray crystallographic example of a hydrogen engaged in an enzymatic LBHB, and demonstrates that desolvation of the active site by ligand binding can provide a protein microenvironment conducive to LBHB formation. It also suggests that LBHBs may contribute to stabilization of the TS in general acid/base catalysis together with other preorganized features of enzyme active sites. These structures reconcile previous experimental results suggesting alternatively Glu166 or Lys73 as the general base for acylation, and underline the importance of considering residue protonation state change when modeling protein-ligand interactions. Additionally, the observation of another LBHB (2.47 Å) between two conserved residues, Asp233 and Asp246, suggests that LBHBs may potentially play a special structural role in proteins.

Can molecular dynamics and QM/MM solve the penicillin binding protein protonation puzzle?

Benzylpenicillin, a member of the ß-lactam antibiotic class, has been widely used to combat bacterial infections since 1947. The general mechanism is well-known: a serine protease enzyme (i.e., DD-peptidase) forms a long lasting intermediate with the lactam ring of the antibiotic known as acylation, effectively preventing biosynthesis of the bacterial cell wall. Despite this overall mechanistic understanding, many details of binding and catalysis are unclear. Specifically, there is ongoing debate about active site protonation states and the role of general acids/bases in the reaction. Herein, a unique combination of MD simulations, QM/MM minimizations, and QM/MM orbital analyses is combined with systematic variation of active site residue protonation states. Critical interactions that maximize the stability of the bound inhibitor are examined and used as metrics. This approach was validated by examining cefoxitin interactions in the CTX-M ß-lactamase from E. coli and compared to an ultra high-resolution (0.88 Å) crystal structure. Upon confirming the approach used, an investigation of the preacylated Streptomyces R61 active site with bound benzylpenicillin was performed, varying the protonation states of His298 and Lys65. We concluded that protonated His298 and deprotonated Lys65 are most likely to exist in the R61 active site.

How does catalase release nitric oxide? A computational structure-activity relationship study.

Hydroxyurea (HU) is the only FDA approved medication for treating sickle cell disease in adults. The primary mechanism of action is pharmacological elevation of nitric oxide (NO) levels which induces propagation of fetal hemoglobin. HU is known to undergo redox reactions with heme based enzymes like hemoglobin and catalase to produce NO. However, specific details about the HU based NO release remain unknown. Experimental studies indicate that interaction of HU with human catalase compound I produces NO. Presently, we combine flexible receptor-flexible substrate induced fit docking (IFD) with energy decomposition analyses to examine the atomic level details of a possible key step in the clinical conversion of HU to NO. Substrate binding modes of nine HU analogs with catalase compound I were investigated to determine the essential properties necessary for effective NO release. Three major binding orientations were found that provide insight into the possible reaction mechanisms for producing NO. Further results show that anion/radical intermediates produced as part of these mechanisms would be stabilized by hydrogen bonding interactions from distal residues His75, Asn148, Gln168, and oxoferryl-heme. These details will ideally contribute to both a clearer mechanistic picture and provide insights for future structure based drug design efforts.

Unlocking the binding and reaction mechanism of hydroxyurea substrates as biological nitric oxide donors.

Hydroxyurea is the only FDA approved treatment of sickle cell disease. It is believed that the primary mechanism of action is associated with the pharmacological elevation of nitric oxide in the blood; however, the exact details of this are still unclear. In the current work, we investigate the atomic level details of this process using a combination of flexible-ligand/flexible-receptor virtual screening coupled with energetic analysis that decomposes interaction energies. Utilizing these methods, we were able to elucidate the previously unknown substrate binding modes of a series of hydroxyurea analogs to hemoglobin and the concomitant structural changes of the enzyme. We identify a backbone carbonyl that forms a hydrogen bond with bound substrates. Our results are consistent with kinetic and electron paramagnetic resonance (EPR) measurements of hydroxyurea-hemoglobin reactions, and a full mechanism is proposed that offers new insights into possibly improving substrate binding and/or reactivity.

Characterizing the mechanism of the double proton transfer in the formamide dimer.

The double proton transfer in the formamide dimer is characterized computationally by combining density functional theory and ab initio methods. The intrinsic reaction coordinate (IRC) is obtained at the B3LYP level of theory. Energies of several points along the IRC are treated by the more rigorous focal point method to test the validity of the B3LYP functional. The reaction mechanism is examined in terms of the energy profile, the reaction force, the chemical potential, and the reaction electronic flux. The energy profile for the activation process of the formamide dimer to the imino ether product obtained with the B3LYP functional is in agreement with the results of the focal point method. Together with the reaction force analysis and the reaction electronic flux a precise assignment of the structural and electronic contributions to the activation barrier becomes possible. The results show that the reaction starts with a structural rearrangement, where the two dimers approach each other, and is followed by electronic changes before the system reaches the transition state. This electronic contribution to the activation barrier steers the activation process. After the transition state is reached, deviations of the B3LYP functional from the more accurate focal point energies become apparent, where the errors may be rationalized in terms of the treatment of exchange. The inconsistency could be assigned to the incapacity of the functional to describe delocalization effects over the whole system.

Noncovalent interactions of a benzo[a]pyrene diol epoxide with DNA base pairs: insight into the formation of adducts of (+)-BaP DE-2 with DNA.

Noncovalent complexes of a tumorigenic benzo[a]pyrene diol epoxide with the guanine-cytosine (GC) and adenine-thymine (AT) base pairs have been examined computationally. (+)-BaP DE-2 forms covalent adducts with DNA via nucleophilic attack on the (+)-BaP DE-2 epoxide. Computational results predict five thermodynamically accessible complexes of AT with (+)-BaP DE-2 that are compatible with intact DNA. Among these, two are expected to lead to adenine adducts. In the lowest energy AT...(+)-BaP DE-2 complex, which has a gas-phase interaction energy of -20.9 kcal mol(-1), the exocyclic NH(2) of adenine is positioned for backside epoxide attack and formation of a trans adduct. The most energetically favorable complex leading to formation of a cis ring-opened adduct lies only 0.6 kcal mol(-1) higher in energy. For GC...(+)-BaP DE-2, there are only two thermodynamically accessible complexes. The higher-lying complex, bound in the gas phase by 24.4 kcal mol(-1) relative to separated GC and (+)-BaP DE-2, would lead to a trans ring-opened N(2)-guanine adduct. In the global minimum energy GC...(+)-BaP DE-2 complex, bound by 27.3 kcal mol(-1), the exocyclic NH(2) group of cytosine is positioned for cis epoxide addition. However, adducts of (+)-BaP DE-2 with cytosine are rarely observed experimentally. The paucity of cytosine adducts, despite the predicted thermodynamic stability of this GC...(+)-BaP DE-2 complex, is attributed to the electrostatic destabilization of the benzylic cation intermediate thought to precede cis addition.

Short intramolecular hydrogen bonds: derivatives of malonaldehyde with symmetrical substituents.

A systematic study of various derivatives of malonaldehyde has been carried out to explore very short hydrogen bonds (r(OO) < 2.450 A). Various electron-withdrawing groups, including CN, NO(2), and BH(2), have been attached to the central carbon atom, C(2). To C(1) and C(3), strong electron donors and/or sterically hindered substituents were used to strengthen the intramolecular hydrogen bond, including but not limited to NH(2), N(CH(3))(2), and C(CH(3))(3). Seven molecules (Figure 2 ) were found to have extremely short intramolecular hydrogen bonds. The chemical systems investigated are intriguing due to their low energetic barriers for the intramolecular hydrogen atom transfers. Classical barriers were predicted using correlated methods including second-order perturbation theory and coupled cluster theory in conjunction with the Dunning hierarchy of correlation consistent basis sets, cc-pVXZ (X = D, T, Q, 5). Focal point analyses allowed for the barriers to be evaluated at the CBS limit including core correlation and zero-point vibrational energy corrections. B3LYP energies are benchmarked against highly accurate correlated energies for intramolecular hydrogen bonded systems. The focal point extrapolated method, including coupled cluster full triple excitation contributions, gives a hydrogen transfer barrier for malonaldehyde of approximately 4 kcal mol(-1). We describe two compounds with extremely low classical barriers, nitromalonamide (0.43 kcal mol(-1)) and 2-borylmalonamide (0.60 kcal mol(-1)). An empirical relationship was drawn between the B3LYP energetic barriers and the predicted coupled cluster barriers at the CBS limit. By relating these two quantities, barrier heights may be estimated for systems too large to presently use highly correlated electronic structure methods.

Identification and Characterization of Noncovalent Interactions That Drive Binding and Specificity in DD-Peptidases and ß-Lactamases.

Bacterial resistance to standard (i.e., ß-lactam-based) antibiotics has become a global pandemic. Simultaneously, research into the underlying causes of resistance has slowed substantially, although its importance is universally recognized. Key to unraveling critical details is characterization of the noncovalent interactions that govern binding and specificity (DD-peptidases, antibiotic targets, versus ß-lactamases, the evolutionarily derived enzymes that play a major role in resistance) and ultimately resistance as a whole. Herein, we describe a detailed investigation that elicits new chemical insights into these underlying intermolecular interactions. Benzylpenicillin and a novel ß-lactam peptidomimetic complexed to the Stremptomyces R61 peptidase are examined using an arsenal of computational techniques: MD simulations, QM/MM calculations, charge perturbation analysis, QM/MM orbital analysis, bioinformatics, flexible receptor/flexible ligand docking, and computational ADME predictions. Several key molecular level interactions are identified that not only shed light onto fundamental resistance mechanisms, but also offer explanations for observed specificity. Specifically, an extended π-π network is elucidated that suggests antibacterial resistance has evolved, in part, due to stabilizing aromatic interactions. Additionally, interactions between the protein and peptidomimetic substrate are identified and characterized. Of particular interest is a water-mediated salt bridge between Asp217 and the positively charged N-terminus of the peptidomimetic, revealing an interaction that may significantly contribute to ß-lactam specificity. Finally, interaction information is used to suggest modifications to current ß-lactam compounds that should both improve binding and specificity in DD-peptidases and their physiochemical properties.