How to Stabilize Carbenes in Enzyme Active Sites without Metal Ions

Carbenes are highly reactive compounds with unique value to synthetic chemistry. However, a small number of natural enzymes have been shown to utilize carbene chemistry, and artificial enzymes engineered with directed evolution required transition metal ions to stabilize the carbene intermediates. To facilitate the design of broader classes of enzymes that can take advantage of the rich carbene chemistry, it is thus important to better understand how to stabilize carbene species in enzyme active sites without metal ions. Motivated by our recent studies of the anaerobic ergothioneine biosynthesis enzyme EanB, we examine carbene–protein interaction with both cluster models and QM/MM simulations. The cluster calculations find that an N-heterocyclic carbene interacts strongly with polar and positively charged protein motifs. In particular, the interaction between a guanidinium group and carbene is as strong as ∼30 kcal/mol, making arginine a great choice for the preferential stabilization of carbenes. We also compare the WT EanB and its mutant in which the key tyrosine was replaced by a non-natural analogue (F2Tyr) using DFTB3/MM simulations. The calculations suggest that the carbene intermediate in the F2Tyr mutant is more stable than that in the WT enzyme by ∼3.5 kcal/mol, due to active site rearrangements that enable a nearby arginine to better stabilize the carbene in the mutant. Overall, the current work lays the foundation for the pursuit of enzyme designs that can take advantage of the unique chemistry offered by carbenes without the requirement of metal ions

Hydrogen Abstraction of Camphor Catalyzed by Cytochrome P450_cam: A QM/MM Study

A combined quantum mechanics and molecular mechanics (QM/MM, QM = UB3LYP-D3, MM = AMBER) method is used to study the hydrogen abstraction reaction in P450_cam catalyzed hydroxylation of camphor in the quartet state. Compared to QM/MM calculations in the literature, this study uses larger basis sets for the most important atoms at the active site and QM/MM Hessian harmonic frequency calculations to determine the standard Gibbs free energy of activation and kinetic isotope effect. The QM/MM covalent boundary is treated with a capping hydrogen atom method, which is simple and robust. An energy barrier of 21.3 kcal/mol and a standard free energy of activation of 16.8 kcal/mol are obtained for this hydrogen abstraction reaction. These values are similar to those reported in the literature, suggesting that when a general protocol is followed, QM/MM results are reproducible. It is found that using a sufficiently large basis set is important to minimize basis set errors

Flexibility of Binding Site is Essential to the Ca²⁺ Selectivity in EF-Hand Calcium-Binding Proteins

High binding affinity and selectivity of metal ions are essential to the function of metalloproteins. Thus, understanding the factors that determine these binding characteristics is of major interest for both fundamental mechanistic investigations and guiding of the design of novel metalloproteins. In this work, we perform QM cluster model calculations and quantum mechanics/molecular mechanics (QM/MM) free energy simulations to understand the binding selectivity of Ca2+ and Mg2+ in the wild-type carp parvalbumin and its mutant. While a nonpolarizable MM model (CHARMM36) does not lead to the correct experimental trend, treatment of the metal binding site with the DFTB3 model in a QM/MM framework leads to relative binding free energies (ΔΔGbind) comparable with experimental data. For the wild-type (WT) protein, the calculated ΔΔGbind is ∼6.6 kcal/mol in comparison with the experimental value of 5.6 kcal/mol. The good agreement highlights the value of a QM description of the metal binding site and supports the role of electronic polarization and charge transfer to metal binding selectivity. For the D51A/E101D/F102W mutant, different binding site models lead to considerable variations in computed binding affinities. With a coordination number of seven for Ca2+, which is shown by QM/MM metadynamics simulations to be the dominant coordination number for the mutant, the calculated relative binding affinity is ∼4.8 kcal/mol, in fair agreement with the experimental value of 1.6 kcal/mol. The WT protein is observed to feature a flexible binding site that accommodates a range of coordination numbers for Ca2+, which is essential to the high binding selectivity for Ca2+ over Mg2+. In the mutant, the E101D mutation reduces the flexibility of the binding site and limits the dominant coordination number of Ca2+ to be seven, thereby leading to reduced binding selectivity against Mg2+. Our results highlight that the binding selectivity of metal ions depends on both the structural and dynamical properties of the protein binding site

Flexibility of Binding Site is Essential to the Ca²⁺ Selectivity in EF-Hand Calcium-Binding Proteins

High binding affinity and selectivity of metal ions are essential to the function of metalloproteins. Thus, understanding the factors that determine these binding characteristics is of major interest for both fundamental mechanistic investigations and guiding of the design of novel metalloproteins. In this work, we perform QM cluster model calculations and quantum mechanics/molecular mechanics (QM/MM) free energy simulations to understand the binding selectivity of Ca2+ and Mg2+ in the wild-type carp parvalbumin and its mutant. While a nonpolarizable MM model (CHARMM36) does not lead to the correct experimental trend, treatment of the metal binding site with the DFTB3 model in a QM/MM framework leads to relative binding free energies (ΔΔGbind) comparable with experimental data. For the wild-type (WT) protein, the calculated ΔΔGbind is ∼6.6 kcal/mol in comparison with the experimental value of 5.6 kcal/mol. The good agreement highlights the value of a QM description of the metal binding site and supports the role of electronic polarization and charge transfer to metal binding selectivity. For the D51A/E101D/F102W mutant, different binding site models lead to considerable variations in computed binding affinities. With a coordination number of seven for Ca2+, which is shown by QM/MM metadynamics simulations to be the dominant coordination number for the mutant, the calculated relative binding affinity is ∼4.8 kcal/mol, in fair agreement with the experimental value of 1.6 kcal/mol. The WT protein is observed to feature a flexible binding site that accommodates a range of coordination numbers for Ca2+, which is essential to the high binding selectivity for Ca2+ over Mg2+. In the mutant, the E101D mutation reduces the flexibility of the binding site and limits the dominant coordination number of Ca2+ to be seven, thereby leading to reduced binding selectivity against Mg2+. Our results highlight that the binding selectivity of metal ions depends on both the structural and dynamical properties of the protein binding site

Catalytic Mechanism of Amyloid‑β Peptide Degradation by Insulin Degrading Enzyme: Insights from Quantum Mechanics and Molecular Mechanics Style Møller–Plesset Second Order Perturbation Theory Calculation

Insulin degrading enzyme (IDE), a metalloprotease that degrades amyloid-β (Aβ) peptides and insulin, is associated with Alzheimer’s disease and diabetes. The mechanism of IDE catalyzed degrading of Aβ peptides, which is of fundamental importance in the design of therapeutic methods for Alzheimer’s disease, has not been fully understood. In this work, combined quantum mechanics and molecular mechanics (QM/MM) style Møller–Plesset second order perturbation theory (MP2) geometry optimization calculations are performed to investigate the catalytic mechanism of the Aβ40 Phe19-Phe20 peptide bond cleavage by human IDE. The analyses using QM/MM MP2 optimization suggest that a neutral water molecule is at the active site of the enzyme–substrate (ES) complex. The water molecule is in hydrogen bonding with the nearby anionic Glu111 of IDE but not directly bound to the catalytic Zn ion. This is confirmed by QM/MM DFTB3 molecular dynamics simulation. Our studies also reveal that the hydrolysis of the Aβ40 Phe19-Phe20 peptide bond by IDE consists of four key steps. The neutral water is first activated by moving toward and binding to the Zn ion. A gem–diol intermediate is then formed by the activated neutral water molecule attacking the C atom of the Phe19-Phe20 peptide bond. The next is the protonation of the N atom of Phe19-Phe20 peptide bond to form an intermediate with an elongated C–N bond. The final step is the breaking of the Phe19-Phe20 C–N bond. The final step is the rate-determining step with a calculated Gibbs free energy of activation of 17.34 kcal/mol, in good agreement with the experimental value 16.7 kcal/mol. This mechanism provides the basis for the design of biochemical methods to modulate the activity of IDE in humans

The 8-Silyloxyquinoline Scaffold as a Versatile Platform for the Sensitive Detection of Aqueous Fluoride

Utilizing a novel 8-silyloxyquinoline scaffold, we demonstrate the ability to synthesize fluorogenic probes for the sensitive and selective detection of inorganic fluoride (NaF) in aqueous samples. Our initial probe design (<b>2</b>) is capable of detecting inorganic fluoride at levels as low as 3.8 μM (72 ppb) in aqueous solutions, well below PHS recommended levels for drinking water (0.7–1.2 ppm), placing this probe among the most sensitive fluoride sensors reported to date. Furthermore, our results highlight the utility of the readily modifiable 8-silyloxyquinoline scaffold for the design of tailored fluoride sensing platforms. We demonstrate the ability to rationally tune the fluorescence and physical properties of the 8-silyloxyquinoline scaffold, producing a red-shifted fluoride probe (<b>4</b>) capable of detecting 50 μM (0.95 ppm) NaF in aqueous samples using a straightforward test-strip-based assay format. Taken together this work provides a template for the design of fluoride sensors capable of reporting on relevant concentrations of fluoride in the laboratory and in the field

A Reactive Force Field for Molecular Dynamics Simulations of Glucose in Aqueous Solution

To expand the capabilities of reactive force field (ReaxFF) in simulations of biological processes involving glucose, in this work, using Metropolis Monte Carlo algorithm, new ReaxFF parameters for glucose have been developed to better describe the properties of glucose in water during molecular dynamics (MD) simulations. With the newly trained ReaxFF, the mutarotation of glucose in water can be better described, as suggested by our metadynamics simulations. In addition, the newly trained ReaxFF can better describe the distributions of the three stable conformers along the key dihedral angle of α-anomer and β-anomer. With better descriptions of hydration around glucose, the Raman and Raman optical activity spectra can be more accurately calculated. In addition, the infrared spectra obtained from simulations with the new glucose ReaxFF are more accurate than those obtained with the original ReaxFF. We note that although our trained ReaxFF performs better than the original ReaxFF, it is not generally applicable to all carbohydrates, which require further parametrization. We also find that the absence of explicit water molecules in the training sets may lead to inaccurate descriptions of water–water interactions around the glucose, implicating that it is necessary to optimize the water ReaxFF parameters together with the target molecule. The improved ReaxFF makes it possible to explore interesting biological processes involving glucose more accurately and efficiently

A Reactive Force Field for Molecular Dynamics Simulations of Glucose in Aqueous Solution

To expand the capabilities of reactive force field (ReaxFF) in simulations of biological processes involving glucose, in this work, using Metropolis Monte Carlo algorithm, new ReaxFF parameters for glucose have been developed to better describe the properties of glucose in water during molecular dynamics (MD) simulations. With the newly trained ReaxFF, the mutarotation of glucose in water can be better described, as suggested by our metadynamics simulations. In addition, the newly trained ReaxFF can better describe the distributions of the three stable conformers along the key dihedral angle of α-anomer and β-anomer. With better descriptions of hydration around glucose, the Raman and Raman optical activity spectra can be more accurately calculated. In addition, the infrared spectra obtained from simulations with the new glucose ReaxFF are more accurate than those obtained with the original ReaxFF. We note that although our trained ReaxFF performs better than the original ReaxFF, it is not generally applicable to all carbohydrates, which require further parametrization. We also find that the absence of explicit water molecules in the training sets may lead to inaccurate descriptions of water–water interactions around the glucose, implicating that it is necessary to optimize the water ReaxFF parameters together with the target molecule. The improved ReaxFF makes it possible to explore interesting biological processes involving glucose more accurately and efficiently

DataSheet_1_Causal relationship between rheumatoid arthritis and hypothyroidism or hyperthyroidism: a bidirectional two-sample univariable and multivariable Mendelian randomization study.zip

ObjectiveThe causal relationship between Rheumatoid arthritis (RA) and hypothyroidism/hyperthyroidism remains controversial due to the limitations of conventional observational research, such as confounding variables and reverse causality. We aimed to examine the potential causal relationship between RA and hypothyroidism/hyperthyroidism using Mendelian randomization (MR).MethodWe conducted a bidirectional two-sample univariable analysis to investigate the potential causal relationship between hypothyroidism/hyperthyroidism and RA. Furthermore, we performed a multivariate analysis to account for the impact of body mass index (BMI), smoking quantity, and alcohol intake frequency.ResultsThe univariable analysis indicated that RA has a causative influence on hypothyroidism (odds ratio [OR]=1.07, 95% confidence interval [CI]=1.01–1.14, P=0.02) and hyperthyroidism (OR=1.32, 95% CI=1.15–1.52, P0.05).ConclusionUnivariate and multivariate MR analyses have validated the causal association between RA and hypothyroidism/hyperthyroidism. Hypothyroidism confirmed a causal relationship with RA when employed as an exposure variable, whereas no such relationship was found between hyperthyroidism and RA.</p

Table_1_The causal relationship between 41 inflammatory cytokines and hypothyroidism: bidirectional two-sample Mendelian randomization study.xlsx

ObjectiveInvestigating the association between inflammatory cytokines and hypothyroidism remains challenging due to limitations in traditional observational studies. In this study, we employed Mendelian randomization (MR) to assess the causal relationship between 41 inflammatory cytokines and hypothyroidism.MethodInflammatory cytokines in 30,155 individuals of European ancestry with hypothyroidism and in a GWAS summary containing 8,293 healthy participants were included in the study for bidirectional two-sample MR analysis. We utilized inverse variance weighting (IVW), weighted median (WM), and Mendelian randomization-Egger (MR-Egger) methods. Multiple sensitivity analyses, including MR-Egger intercept test, leave-one-out analysis, funnel plot, scatterplot, and MR-PRESSO, were applied to evaluate assumptions.ResultsWe found evidence of a causal effect of IL-7 and macrophage inflammatory protein-1β (MIP-1β) on the risk of hypothyroidism, and a causal effect of hypothyroidism on several cytokines, including granulocyte colony-stimulating factor (G-CSF), IL-13, IL-16, IL-2rα, IL-6, IL-7, IL-9, interferon-γ-inducible protein 10 (IP10), monokine induced by interferon (IFN)-γ (MIG), macrophage inflammatory protein-1β (MIP-1β), stem cell growth factors-β (SCGF-β), stromal cell derived factor-1α (SDF-1α), and tumor necrosis factor-α (TNF-α).ConclusionOur study suggests that IL-7 and MIP-1β may play a role in the pathogenesis of hypothyroidism, and that hypothyroidism may induce a systemic inflammatory response involving multiple cytokines. These findings may have implications for the prevention and treatment of hypothyroidism and its complications. However, further experimental studies are needed to validate the causal relationships and the potential of these cytokines as drug targets.</p