Critical Evaluation of Polarizable and Nonpolarizable Force Fields for Proteins Using Experimentally Derived Nitrile Electric Fields.

Molecular dynamics (MD) simulations are frequently carried out for proteins to investigate the role of electrostatics in their biological function. The choice of force field (FF) can significantly alter the MD results, as the simulated local electrostatic interactions lack benchmarking in the absence of appropriate experimental methods. We recently reported that the transition dipole moment (TDM) of the popular nitrile vibrational probe varies linearly with the environmental electric field, overcoming well-known hydrogen bonding (H-bonding) issues for the nitrile frequency and, thus, enabling the unambiguous measurement of electric fields in proteins (J. Am. Chem. Soc. 2022, 144 (17), 7562-7567). Herein, we utilize this new strategy to enable comparisons of experimental and simulated electric fields in protein environments. Specifically, previously determined TDM electric fields exerted onto nitrile-containing o-cyanophenylalanine residues in photoactive yellow protein are compared with MD electric fields from the fixed-charge AMBER FF and the polarizable AMOEBA FF. We observe that the electric field distributions for H-bonding nitriles are substantially affected by the choice of FF. As such, AMBER underestimates electric fields for nitriles experiencing moderate field strengths; in contrast, AMOEBA robustly recapitulates the TDM electric fields. The FF dependence of the electric fields can be partly explained by the presence of additional negative charge density along the nitrile bond axis in AMOEBA, which is due to the inclusion of higher-order multipole parameters; this, in turn, begets more head-on nitrile H-bonds. We conclude by discussing the implications of the FF dependence for the simulation of nitriles and proteins in general.

In Situ Spectroscopic Detection of Large-Scale Reorientations of Transmembrane Helices During Influenza A M2 Channel Opening.

Viroporins are small ion channels in membranes of enveloped viruses that play key roles during viral life cycles. To use viroporins as drug targets against viral infection requires in-depth mechanistic understanding and, with that, methods that enable investigations under in situ conditions. Here, we apply surface-enhanced infrared absorption (SEIRA) spectroscopy to Influenza A M2 reconstituted within a solid-supported membrane, to shed light on the mechanics of its viroporin function. M2 is a paradigm of pH-activated proton channels and controls the proton flux into the viral interior during viral infection. We use SEIRA to track the large-scale reorientation of M2's transmembrane α-helices in situ during pH-activated channel opening. We quantify this event as a helical tilt from 26° to 40° by correlating the experimental results with solid-state nuclear magnetic resonance-informed computational spectroscopy. This mechanical motion is impeded upon addition of the inhibitor rimantadine, giving a direct spectroscopic marker to test antiviral activity. The presented approach provides a spectroscopic tool to quantify large-scale structural changes and to track the function and inhibition of the growing number of viroporins from pathogenic viruses in future studies.

Nitrile Infrared Intensities Characterize Electric Fields and Hydrogen Bonding in Protic, Aprotic, and Protein Environments.

Nitriles are widely used vibrational probes; however, the interpretation of their IR frequencies is complicated by hydrogen bonding (H-bonding) in protic environments. We report a new vibrational Stark effect (VSE) that correlates the electric field projected on the -C≡N bond to the transition dipole moment and, by extension, the nitrile peak area or integrated intensity. This linear VSE applies to both H-bonding and non-H-bonding interactions. It can therefore be generally applied to determine electric fields in all environments. Additionally, it allows for semiempirical extraction of the H-bonding contribution to the blueshift of the nitrile frequency. Nitriles were incorporated at H-bonding and non-H-bonding protein sites using amber suppression, and each nitrile variant was structurally characterized at high resolution. We exploited the combined information available from variations in frequency and integrated intensity and demonstrate that nitriles are a generally useful probe for electric fields.

Protein Electric Fields Enable Faster and Longer-Lasting Covalent Inhibition of ß-Lactamases.

The widespread design of covalent drugs has focused on crafting reactive groups of proper electrophilicity and positioning toward targeted amino-acid nucleophiles. We found that environmental electric fields projected onto a reactive chemical bond, an overlooked design element, play essential roles in the covalent inhibition of TEM-1 ß-lactamase by avibactam. Using the vibrational Stark effect, the magnitudes of the electric fields that are exerted by TEM active sites onto avibactam's reactive C═O were measured and demonstrate an electrostatic gating effect that promotes bond formation yet relatively suppresses the reverse dissociation. These results suggest new principles of covalent drug design and off-target site prediction. Unlike shape and electrostatic complementary which address binding constants, electrostatic catalysis drives reaction rates, essential for covalent inhibition, and deepens our understanding of chemical reactivity, selectivity, and stability in complex systems.

Local Electric Field Changes during the Photoconversion of the Bathy Phytochrome Agp2.

Phytochromes switch between a physiologically inactive and active state via a light-induced reaction cascade, which is initiated by isomerization of the tetrapyrrole chromophore and leads to the functionally relevant secondary structure transition of a protein segment (tongue). Although details of the underlying cause-effect relationships are not known, electrostatic fields are likely to play a crucial role in coupling chromophores and protein structural changes. Here, we studied local electric field changes during the photoconversion of the dark state Pfr to the photoactivated state Pr of the bathy phytochrome Agp2. Substituting Tyr165 and Phe192 in the chromophore pocket by para-cyanophenylalanine (pCNF), we monitored the respective nitrile stretching modes in the various states of photoconversion (vibrational Stark effect). Resonance Raman and IR spectroscopic analyses revealed that both pCNF-substituted variants undergo the same photoinduced structural changes as wild-type Agp2. Based on a structural model for the Pfr state of F192pCNF, a molecular mechanical-quantum mechanical approach was employed to calculate the electric field at the nitrile group and the respective stretching frequency, in excellent agreement with the experiment. These calculations serve as a reference for determining the electric field changes in the photoinduced states of F192pCNF. Unlike F192pCNF, the nitrile group in Y165pCNF is strongly hydrogen bonded such that the theoretical approach is not applicable. However, in both variants, the largest changes of the nitrile stretching modes occur in the last step of the photoconversion, supporting the view that the proton-coupled restructuring of the tongue is accompanied by a change of the electric field.

The C-Terminal VPRTES Tail of LL-37 Influences the Mode of Attachment to a Lipid Bilayer and Antimicrobial Activity.

Cathelicidins are a family of host defense antimicrobial peptides in mammalian species. Among them, LL-37 is the only peptide of this family found in humans. Although LL-37 has been intensively investigated in the past, the mode of exerting its bactericidal activity through the specific interactions with bacterial membranes remains elusive. In this work, we combined microbiological and computational approaches with a tool box of experimental biophysical techniques, including conventional and surface-enhanced infrared absorption spectroscopy as well as fluorescence spectroscopy to characterize the structural and dynamic properties of LL-37 and shorter variants adsorbed on POPC/POPG (9:1) lipid bilayers as mimics of bacterial membranes. First, microbiological assays demonstrate that, while LL-32 and, in a lesser degree, LL-37 show hemolysis and antimicrobial activity, LL-20 remains practically inactive. Second, by comparing experimental and computational data of LL-37 with LL-20, we explained the bactericidal activity of the active peptide core as a consequence of an increased flexibility of the peptide structure, leading to reactive dangling charged side chains. Third, permeabilization assays showed a concentration-dependent membrane disruption activity of LL-37 and LL-32: at high peptide concentrations, LL-32 shows higher activity than LL-37, while, at low peptide concentrations, both peptides show similar activities. Responsible for this behavior is the C-terminal VPRTES tail (Ct-VPRTES tail), which, according to atomistic simulations, is able to promote the insertion of the peptide in the membrane and plays an essential role in controlling ordered peptide oligomerization on the surface of the membrane.

On the pH-Modulated Ru-Based Prodrug Activation Mechanism.

The RuIII-based prodrug AziRu efficiently binds to proteins, but the mechanism of its release is still disputed. Herein, in order to test the hypothesis of a reduction-mediated Ru release from proteins, a Raman-assisted crystallographic study on AziRu binding to a model protein (hen egg white lysozyme), in two different oxidation states, RuII and RuIII, was carried out. Our results indicate Ru reduction, but the Ru release upon reduction is dependent on the reducing agent. To better understand this process, a pH-dependent, spectroelectrochemical surface-enhanced Raman scattering (SERS) study was performed also on AziRu-functionalized Au electrodes as a surrogate and simplest model system of RuII- and RuIII-based drugs. This SERS study provided a p Ka of 6.0 ± 0.4 for aquated AziRu in the RuIII state, which falls in the watershed range of pH values separating most cancer environments from their physiological counterparts. These experiments also indicate a dramatic shift of the redox potential E0 by >600 mV of aquated AziRu toward more positive potentials upon acidification, suggesting a selective AziRu reduction in cancer lumen but not in healthy ones. It is expected that the nature of the ligands (e.g., pyridine vs imidazole, present in well-known RuIII complex NAMI-A) will modulate the p Ka and E0, without affecting the underlying reaction mechanism.

Monitoring the Orientational Changes of Alamethicin during Incorporation into Bilayer Lipid Membranes.

Antimicrobial peptides (AMPs) are the first line of defense after contact of an infectious invader, for example, bacterium or virus, with a host and an integral part of the innate immune system of humans. Their broad spectrum of biological functions ranges from cell membrane disruption over facilitation of chemotaxis to interaction with membrane-bound or intracellular receptors, thus providing novel strategies to overcome bacterial resistances. Especially, the clarification of the mechanisms and dynamics of AMP incorporation into bacterial membranes is of high interest, and different mechanistic models are still under discussion. In this work, we studied the incorporation of the peptaibol alamethicin (ALM) into tethered bilayer lipid membranes on electrodes in combination with surface-enhanced infrared absorption (SEIRA) spectroscopy. This approach allows monitoring the spontaneous and potential-induced ion channel formation of ALM in situ. The complex incorporation kinetics revealed a multistep mechanism that points to peptide-peptide interactions prior to penetrating the membrane and adopting the transmembrane configuration. On the basis of the anisotropy of the backbone amide I and II infrared absorptions determined by density functional theory calculations, we employed a mathematical model to evaluate ALM reorientations monitored by SEIRA spectroscopy. Accordingly, ALM was found to adopt inclination angles of ca. 69°-78° and 21° in its interfacially adsorbed and transmembrane incorporated states, respectively. These orientations can be stabilized efficiently by the dipolar interaction with lipid head groups or by the application of a potential gradient. The presented potential-controlled mechanistic study suggests an N-terminal integration of ALM into membranes as monomers or parallel oligomers to form ion channels composed of parallel-oriented helices, whereas antiparallel oligomers are barred from intrusion.

Catalytic Activity and Proton Translocation of Reconstituted Respiratory Complex I Monitored by Surface-Enhanced Infrared Absorption Spectroscopy.

Respiratory complex I (CpI) is a key player in the way organisms obtain energy, being an energy transducer, which couples nicotinamide adenine dinucleotide (NADH)/quinone oxidoreduction with proton translocation by a mechanism that remains elusive so far. In this work, we monitored the function of CpI in a biomimetic, supported lipid membrane system assembled on a 4-aminothiophenol (4-ATP) self-assembled monolayer by surface-enhanced infrared absorption spectroscopy. 4-ATP serves not only as a linker molecule to a nanostructured gold surface but also as pH sensor, as indicated by concomitant density functional theory calculations. In this way, we were able to monitor NADH/quinone oxidoreduction-induced transmembrane proton translocation via the protonation state of 4-ATP, depending on the net orientation of CpI molecules induced by two complementary approaches. An associated change of the amide I/amide II band intensity ratio indicates conformational modifications upon catalysis which may involve movements of transmembrane helices or other secondary structural elements, as suggested in the literature [ Di Luca , Proc. Natl. Acad. Sci. U.S.A. , 2017 , 114 , E6314 - E6321 ].

Voltage-dependent structural changes of the membrane-bound anion channel hVDAC1 probed by SEIRA and electrochemical impedance spectroscopy.

The voltage-dependent anion channel (VDAC) is a transmembrane protein that regulates the transfer of metabolites between the cytosol and the mitochondrium. Opening and partial closing of the channel is known to be driven by the transmembrane potential via a mechanism that is not fully understood. In this work, we employed a spectroelectrochemical approach to probe the voltage-induced molecular structure changes of human VDAC1 (hVDAC1) embedded in a tethered bilayer lipid membrane on a nanostructured Au electrode. The model membrane consisted of a mixed self-assembled monolayer of 6-mercaptohexanol and (cholesterylpolyethylenoxy)thiol, followed by the deposition of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine vesicles including hVDAC1. The stepwise assembly of the model membrane and the incorporation of hVDAC1 were monitored by surface enhanced infrared absorption and electrochemical impedance spectroscopy. Difference spectra allowed for identifying the spectral changes which may be associated with the transition from the open to the "closed" states by shifting the potential above or below the transmembrane potential determined to be ca. 0.0 V vs. the open circuit potential. These spectral changes were interpreted on the basis of the orientation- and distance-dependent IR enhancement and indicate alterations of the inclination angle of the ß-strands as crucial molecular events, reflecting an expansion or contraction of the ß-barrel pore. These protein structural changes that do not confirm nor exclude the reorientation of the α-helix are either directly induced by the electric field or a consequence of a potential-dependent repulsion or attraction of the bilayer.

Cation Dependence of Enniatin B/Membrane-Interactions Assessed Using Surface-Enhanced Infrared Absorption (SEIRA) Spectroscopy.

Enniatins are mycotoxins with well-known antibacterial, antifungal, antihelmintic and antiviral activity, which have recently come to attention as potential mitochondriotoxic anticancer agents. The cytotoxicity of enniatins is traced back to ionophoric properties, in which the cyclodepsipeptidic structure results in enniatin:cation-complexes of various stoichiometries proposed as membrane-active species. In this work, we employed a combination of surface-enhanced infrared absorption (SEIRA) spectroscopy, tethered bilayer lipid membranes (tBLMs) and density functional theory (DFT)-based computational spectroscopy to monitor the cation-dependence (Mz+=Na+, K+, Cs+, Li+, Mg2+, Ca2+) on the mechanism of enniatin B (EB) incorporation into membranes and identify the functionally relevant EBn : Mz+ complexes formed. We find that Na+ promotes a cooperative incorporation, modelled via an autocatalytic mechanism and mediated by a distorted 2 : 1-EB2 : Na+ complex. K+ (and Cs+) leads to a direct but less efficient insertion into membranes due to the adoption of "ideal" EB2 : K+ sandwich complexes. In contrast, the presence of Li+, Mg2+, and Ca2+ causes a (partial) extraction of EB from the membrane via the formation of "belted" 1 : 1-EB : Mz+ complexes, which screen the cationic charge less efficiently. Our results point to a relevance of the cation dependence for the transport into the malignant cells where the mitochondriotoxic anticancer activity is exerted.

Resonance Raman characterization of the ammonia-generated oxo intermediate of cytochrome c oxidase from Paracoccus denitrificans.

A novel oxo state of cytochrome c oxidase from Paracoccus denitrificans generated by successive addition of excess H2O2 and ammonia was investigated using resonance Raman (RR) spectroscopy. Addition of ammonia to the H2O2-generated artificial F state resulted in an upshift of the oxoferryl stretching vibration from 790 to 796 cm(-1), indicating that ammonia influences ligation of the heme-bound oxygen in the binuclear center. Concomitantly performed RR measurements in the high-frequency region between 1300 and 1700 cm(-1) showed a high-spin to low-spin transition of heme a3 upon generation of the F state that was not altered by addition of ammonia. Removal of H2O2 by addition of catalase resulted in the disappearance of the oxoferryl stretching vibration and major back transformation of heme a3 into the high-spin state. The ratio of high-spin to low-spin states was identical for intermediates created with and without ammonia, leading to the conclusion that ammonia does not interact directly with heme a3. Only for the ammonia-created state was a band at 612 nm observed in the UV-visible difference spectrum that was shifted to 608 nm after addition of catalase. Our results support the hypothesis by von der Hocht et al. [von der Hocht, I., et al. (2011) Proc. Natl. Acad. Sci. U.S.A. 108, 3964-3969] that addition of ammonia creates a novel oxo intermediate state called PN where ammonia binds to CuB once the oxo intermediate F state has been formed.

Potential-dependent surface-enhanced resonance Raman spectroscopy at nanostructured TiO2 : a case study on cytochrome b5.

Nanostructured titanium dioxide (TiO2 ) electrodes, prepared by anodization of titanium, are employed to probe the electron-transfer process of cytochrome b5 (cyt b5 ) by surface-enhanced resonance Raman (SERR) spectroscopy. Concomitant with the increased nanoscopic surface roughness of TiO2 , achieved by raising the anodization voltage from 10 to 20 V, the enhancement factor increases from 2.4 to 8.6, which is rationalized by calculations of the electric field enhancement. Cyt b5 is immobilized on TiO2 under preservation of its native structure but it displays a non-ideal redox behavior due to the limited conductivity of the electrode material. The electron-transfer efficiency which depends on the crystalline phase of TiO2 has to be improved by appropriate doping for applications in bioelectrochemistry.

Potential Distribution across Model Membranes.

Membrane models assembled on electrodes are widely used tools to study potential-dependent molecular processes at or in membranes. However, the relationship between the electrode potential and the potential across the membrane is not known. Here we studied lipid bilayers immobilized on mixed self-assembled monolayers (SAM) on Au electrodes. The mixed SAM was composed of thiol derivatives of different chain lengths such that between the islands of the short one, mercaptobenzonitrile (MBN), and the tethered lipid bilayer an aqueous compartment was formed. The nitrile function of MBN, which served as a reporter group for the vibrational Stark effect (VSE), was probed by surface-enhanced infrared absorption spectroscopy to determine the local electric field as a function of the electrode potential for pure MBN, mixed SAM, and the bilayer system. In parallel, we calculated electric fields at the VSE probe by molecular dynamics (MD) simulations for different charge densities on the metal, thereby mimicking electrode potential changes. The agreement with the experiments was very good for the calculations of the pure MBN SAM and only slightly worse for the mixed SAM. The comparison with the experiments also guided the design of the bilayer system in the MD setups, which were selected to calculate the electrode potential dependence of the transmembrane potential, a quantity that is not directly accessible by the experiments. The results agree very well with estimates in previous studies and thus demonstrate that the present combined experimental-theoretical approach is a promising tool for describing potential-dependent processes at biomimetic interfaces.

A two-directional vibrational probe reveals different electric field orientations in solution and an enzyme active site.

The catalytic power of an electric field depends on its magnitude and orientation with respect to the reactive chemical species. Understanding and designing new catalysts for electrostatic catalysis thus requires methods to measure the electric field orientation and magnitude at the molecular scale. We demonstrate that electric field orientations can be extracted using a two-directional vibrational probe by exploiting the vibrational Stark effect of both the C=O and C-D stretches of a deuterated aldehyde. Combining spectroscopy with molecular dynamics and electronic structure partitioning methods, we demonstrate that, despite distinct polarities, solvents act similarly in their preference for electrostatically stabilizing large bond dipoles at the expense of destabilizing small ones. In contrast, we find that for an active-site aldehyde inhibitor of liver alcohol dehydrogenase, the electric field orientation deviates markedly from that found in solvents, which provides direct evidence for the fundamental difference between the electrostatic environment of solvents and that of a preorganized enzyme active site.

The Photoreaction of the Proton-Pumping Rhodopsin 1 From the Maize Pathogenic Basidiomycete Ustilago maydis.

Microbial rhodopsins have recently been discovered in pathogenic fungi and have been postulated to be involved in signaling during the course of an infection. Here, we report on the spectroscopic characterization of a light-driven proton pump rhodopsin (UmRh1) from the smut pathogen Ustilago maydis, the causative agent of tumors in maize plants. Electrophysiology, time-resolved UV/Vis and vibrational spectroscopy indicate a pH-dependent photocycle. We also characterized the impact of the auxin hormone indole-3-acetic acid that was shown to influence the pump activity of UmRh1 on individual photocycle intermediates. A facile pumping activity test was established of UmRh1 expressed in Pichia pastoris cells, for probing proton pumping out of the living yeast cells during illumination. We show similarities and distinct differences to the well-known bacteriorhodopsin from archaea and discuss the putative role of UmRh1 in pathogenesis.

The Interplay of Electrostatics and Chemical Positioning in the Evolution of Antibiotic Resistance in TEM ß-Lactamases.

The interplay of enzyme active site electrostatics and chemical positioning is important for understanding the origin(s) of enzyme catalysis and the design of novel catalysts. We reconstruct the evolutionary trajectory of TEM-1 ß-lactamase to TEM-52 toward extended-spectrum activity to better understand the emergence of antibiotic resistance and to provide insights into the structure-function paradigm and noncovalent interactions involved in catalysis. Utilizing a detailed kinetic analysis and the vibrational Stark effect, we quantify the changes in rates and electric fields in the Michaelis and acyl-enzyme complexes for penicillin G and cefotaxime to ascertain the evolutionary role of electric fields to modulate function. These data are combined with MD simulations to interpret and quantify the substrate-dependent structural changes during evolution. We observe that this evolutionary trajectory utilizes a large preorganized electric field and substrate-dependent chemical positioning to facilitate catalysis. This governs the evolvability, substrate promiscuity, and protein fitness landscape in TEM ß-lactamase antibiotic resistance.

Testing the Limitations of MD-Based Local Electric Fields Using the Vibrational Stark Effect in Solution: Penicillin G as a Test Case.

Noncovalent interactions underlie nearly all molecular processes in the condensed phase from solvation to catalysis. Their quantification within a physically consistent framework remains challenging. Experimental vibrational Stark effect (VSE)-based solvatochromism can be combined with molecular dynamics (MD) simulations to quantify the electrostatic forces in solute-solvent interactions for small rigid molecules and, by extension, when these solutes bind in enzyme active sites. While generalizing this approach toward more complex (bio)molecules, such as the conformationally flexible and charged penicillin G (PenG), we were surprised to observe inconsistencies in MD-based electric fields. Combining synthesis, VSE spectroscopy, and computational methods, we provide an intimate view on the origins of these discrepancies. We observe that the electric fields are correlated to conformation-dependent effects of the flexible PenG side chain, including both the local solvation structure and solute conformational sampling in MD. Additionally, we identified that MD-based electric fields are consistently overestimated in three-point water models in the vicinity of charged groups; this cannot be entirely ameliorated using polarizable force fields (AMOEBA) or advanced water models. This work demonstrates the value of the VSE as a direct method for experiment-guided refinements of MD force fields and establishes a general reductionist approach to calibrating vibrational probes for complex (bio)molecules.

Quantification of Local Electric Field Changes at the Active Site of Cytochrome c Oxidase by Fourier Transform Infrared Spectroelectrochemical Titrations.

Cytochrome c oxidase (CcO) is a transmembrane protein complex that reduces molecular oxygen to water while translocating protons across the mitochondrial membrane. Changes in the redox states of its cofactors trigger both O2 reduction and vectorial proton transfer, which includes a proton-loading site, yet unidentified. In this work, we exploited carbon monoxide (CO) as a vibrational Stark effect (VSE) probe at the binuclear center of CcO from Rhodobacter sphaeroides. The CO stretching frequency was monitored as a function of the electrical potential, using Fourier transform infrared (FTIR) absorption spectroelectrochemistry. We observed three different redox states (R4CO, R2CO, and O), determined their midpoint potential, and compared the resulting electric field to electrostatic calculations. A change in the local electric field strength of +2.9 MV/cm was derived, which was induced by the redox transition from R4CO to R2CO. We performed potential jump experiments to accumulate the R2CO and R4CO species and studied the FTIR difference spectra in the protein fingerprint region. The comparison of the experimental and computational results reveals that the key glutamic acid residue E286 is protonated in the observed states, and that its hydrogen-bonding environment is disturbed upon the redox transition of heme a3. Our experiments also suggest propionate A of heme a3 changing its protonation state in concert with the redox state of a second cofactor, heme a. This supports the role of propionic acid side chains as part of the proton-loading site.

Biosynthetic Incorporation of Site-Specific Isotopes in ß-Lactam Antibiotics Enables Biophysical Studies.

A biophysical understanding of the mechanistic, chemical, and physical origins underlying antibiotic action and resistance is vital to the discovery of novel therapeutics and the development of strategies to combat the growing emergence of antibiotic resistance. The site-specific introduction of stable-isotope labels into chemically complex natural products is particularly important for techniques such as NMR, IR, mass spectrometry, imaging, and kinetic isotope effects. Toward this goal, we developed a biosynthetic strategy for the site-specific incorporation of 13C labels into the canonical ß-lactam carbonyl of penicillin G and cefotaxime, the latter via cephalosporin C. This was achieved through sulfur-replacement with 1-13C-l-cysteine, resulting in high isotope incorporations and milligram-scale yields. Using 13C NMR and isotope-edited IR difference spectroscopy, we illustrate how these molecules can be used to interrogate interactions with their protein targets, e.g., TEM-1 ß-lactamase. This method provides a feasible route to isotopically labeled penicillin and cephalosporin precursors for future biophysical studies.