skip to main content


Title: Molecular Dynamics Simulations Guide Chimeragenesis and Engineered Control of Chemoselectivity in Diketopiperazine Dimerases
Abstract

In the biosynthesis of the tryptophan‐linked dimeric diketopiperazines (DKPs), cytochromes P450 selectively couple DKP monomers to generate a variety of intricate and isomeric frameworks. To determine the molecular basis for selectivity of these biocatalysts we obtained a high‐resolution crystal structure of selective Csp2−N bond forming dimerase, AspB. Overlay of the AspB structure onto C−C and C−N bond forming homolog NzeB revealed no significant structural variance to explain their divergent chemoselectivities. Molecular dynamics (MD) simulations identified a region of NzeB with increased conformational flexibility relative to AspB, and interchange of this region along with a single active site mutation led to a variant that catalyzes exclusive C−N bond formation. MD simulations also suggest that intermolecular C−C or C−N bond formation results from a change in mechanism, supported experimentally through use of a substrate mimic.

 
more » « less
NSF-PAR ID:
10397637
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Angewandte Chemie
Volume:
135
Issue:
20
ISSN:
0044-8249
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    In the biosynthesis of the tryptophan‐linked dimeric diketopiperazines (DKPs), cytochromes P450 selectively couple DKP monomers to generate a variety of intricate and isomeric frameworks. To determine the molecular basis for selectivity of these biocatalysts we obtained a high‐resolution crystal structure of selective Csp2−N bond forming dimerase, AspB. Overlay of the AspB structure onto C−C and C−N bond forming homolog NzeB revealed no significant structural variance to explain their divergent chemoselectivities. Molecular dynamics (MD) simulations identified a region of NzeB with increased conformational flexibility relative to AspB, and interchange of this region along with a single active site mutation led to a variant that catalyzes exclusive C−N bond formation. MD simulations also suggest that intermolecular C−C or C−N bond formation results from a change in mechanism, supported experimentally through use of a substrate mimic.

     
    more » « less
  2. Abstract

    Infrared (IR) and Raman spectroscopic features of silicate glasses are often interpreted based on the analogy with those of smaller molecules, molecular clusters, or crystalline counterparts; this study tests the accuracy and validity of these widely cited peak assignment schemes by comparing vibrational spectral features with bond parameters of the glass network created by molecular dynamics (MD) simulations. A series of sodium silicate glasses with compositions of [Na2O]x[Al2O3]2[SiO2]98−xwithx = 7, 12, 17, and 22 were synthesized and analyzed with IR and Raman. A silica glass substrate and a crystalline quartz were also analyzed for comparison. Glass structures with the same compositions were generated with MD simulations using three types of potentials: fixed partial charge pairwise (Teter), partial diffuse charge potential (MGFF), and bond order‐based charge transfer potential (ReaxFF). The comparison of simulated and experimental IR spectra showed that, among these three potentials tested, ReaxFF reproduces the concentration dependence of spectral features closest to the experimentally observed trend. Thus, the bond length and angle distributions as well as Si–Qnspecies and ring size distributions of silica and sodium silicate glasses were obtained from ReaxFF‐MD simulations and further compared with the peak assignment or deconvolution schemes—which have been widely used since 1970s and 1980s—(a) correlation between the IR peak position in the Si–O stretch region (1050‐1120 cm−1) and the Si–O–Si bond angle; (b) deconvolution of the Raman bands in the Si–O stretch region with theQnspeciation; and (c) assignment of the Raman bands in the 420‐600 cm−1region to the bending modes of (SiO)nrings with different sizes (typically, n = 3‐6). The comparisons showed that none of these widely used methods is congruent with the bond parameters or structures of silicate glass networks produced via ReaxFF‐MD simulations. This finding invokes that the adequacy of these spectral interpretation methods must be questioned. Alternative interpretations are proposed, which are to be tested independently in future studies.

     
    more » « less
  3. Abstract

    Efficiently hyperpolarizing proton‐dense molecular solids through dynamic nuclear polarization (DNP) solid‐state NMR is still an unmet challenge. Polarizing agents (PAs) developed so far do not perform well on proton‐rich systems, such as organic microcrystals and biomolecular assemblies. Herein we introduce a new PA, cAsymPol‐POK, and report outstanding hyperpolarization efficiency on 12.76 kDa U‐13C,15N‐labeled LecA protein and pharmaceutical drugs at high magnetic fields (up to 18.8 T) and fast magic angle spinning (MAS) frequencies (up to 40 kHz). The performance of cAsymPol‐POK is rationalized by MAS‐DNP simulations combined with electron paramagnetic resonance (EPR), density functional theory (DFT) and molecular dynamics (MD). This work shows that this new biradical is compatible with challenging biomolecular applications and unlocks the rapid acquisition of13C–13C and15N–13C correlations of pharmaceutical drugs at natural isotopic abundance, which are key experiments for structure determination.

     
    more » « less
  4. Abstract

    Efficiently hyperpolarizing proton‐dense molecular solids through dynamic nuclear polarization (DNP) solid‐state NMR is still an unmet challenge. Polarizing agents (PAs) developed so far do not perform well on proton‐rich systems, such as organic microcrystals and biomolecular assemblies. Herein we introduce a new PA, cAsymPol‐POK, and report outstanding hyperpolarization efficiency on 12.76 kDa U‐13C,15N‐labeled LecA protein and pharmaceutical drugs at high magnetic fields (up to 18.8 T) and fast magic angle spinning (MAS) frequencies (up to 40 kHz). The performance of cAsymPol‐POK is rationalized by MAS‐DNP simulations combined with electron paramagnetic resonance (EPR), density functional theory (DFT) and molecular dynamics (MD). This work shows that this new biradical is compatible with challenging biomolecular applications and unlocks the rapid acquisition of13C–13C and15N–13C correlations of pharmaceutical drugs at natural isotopic abundance, which are key experiments for structure determination.

     
    more » « less
  5. Abstract

    Described is a systematic comparison of factors impacting the relative rates and selectivities of C(sp3)−C and C(sp3)−O bond‐forming reactions at high‐valent Ni as a function of oxidation state. Two Ni complexes are compared: a cationic octahedral NiIVcomplex ligated by tris(pyrazolyl)borate and a cationic octahedral NiIIIcomplex ligated by tris(pyrazolyl)methane. Key features of reactivity/selectivity are revealed: 1) C(sp3)−C(sp2) bond‐forming reductive elimination occurs from both centers, but the NiIIIcomplex reacts up to 300‐fold faster than the NiIV, depending on the reaction conditions. The relative reactivity is proposed to derive from ligand dissociation kinetics, which vary as a function of oxidation state and the presence/absence of visible light. 2) Upon the addition of acetate (AcO), the NiIVcomplex exclusively undergoes C(sp3)−OAc bond formation, while the NiIIIanalogue forms the C(sp3)−C(sp2) coupled product selectively. This difference is rationalized based on the electrophilicity of the respective M−C(sp3) bonds, and thus their relative reactivity towards outer‐sphere SN2‐type bond‐forming reactions.

     
    more » « less