Aromatic Diboronic Acids as Effective KPC/AmpC Inhibitors

Over 30 compounds, including para-, meta-, and ortho-phenylenediboronic acids, ortho-substituted phenylboronic acids, benzenetriboronic acids, di- and triboronated thiophenes, and pyridine derivatives were investigated as potential β-lactamase inhibitors. The highest activity against KPC-type carbapenemases was found for ortho-phenylenediboronic acid 3a, which at the concentration of 8/4 mg/L reduced carbapenems' MICs up to 16/8-fold, respectively. Checkerboard assays revealed strong synergy between carbapenems and 3a with the fractional inhibitory concentrations indices of 0.1-0.32. The nitrocefin hydrolysis test and the whole cell assay with E. coli DH5α transformant carrying blaKPC-3 proved KPC enzyme being its molecular target. para-Phenylenediboronic acids efficiently potentiated carbapenems against KPC-producers and ceftazidime against AmpC-producers, whereas meta-phenylenediboronic acids enhanced only ceftazidime activity against the latter ones. Finally, the statistical analysis confirmed that ortho-phenylenediboronic acids act synergistically with carbapenems significantly stronger than other groups. Since the obtained phenylenediboronic compounds are not toxic to MRC-5 human fibroblasts at the tested concentrations, they can be considered promising scaffolds for the future development of novel KPC/AmpC inhibitors. The complexation of KPC-2 with the most representative isomeric phenylenediboronic acids 1a, 2a, and 3a was modeled by quantum mechanics/molecular mechanics calculations. Compound 3a reached the most effective configuration enabling covalent binding to the catalytic Ser70 residue.

Exploring Covalent Docking Mechanisms of Boron-Based Inhibitors to Class A, C and D β-Lactamases Using Time-dependent Hybrid QM/MM Simulations

Recently, molecular covalent docking has been extensively developed to design new classes of inhibitors that form chemical bonds with their biological targets. This strategy for the design of such inhibitors, in particular boron-based inhibitors, holds great promise for the vast family of β-lactamases produced, inter alia, by Gram-negative antibiotic-resistant bacteria. However, the description of covalent docking processes requires a quantum-mechanical approach, and so far, only a few studies of this type have been presented. This study accurately describes the covalent docking process between two model inhibitors - representing two large families of inhibitors based on boronic-acid and bicyclic boronate scaffolds, and three β-lactamases which belong to the A, C, and D classes. Molecular fragments containing boron can be converted from a neutral, trigonal, planar state with sp2 hybridization to the anionic, tetrahedral sp3 state in a process sometimes referred to as morphing. This study applies multi-scale modeling methods, in particular, the hybrid QM/MM approach which has predictive power reaching well beyond conventional molecular modeling. Time-dependent QM/MM simulations indicated several structural changes and geometric preferences, ultimately leading to covalent docking processes. With current computing technologies, this approach is not computationally expensive, can be used in standard molecular modeling and molecular design works, and can effectively support experimental research which should allow for a detailed understanding of complex processes important to molecular medicine. In particular, it can support the rational design of covalent boron-based inhibitors for β-lactamases as well as for many other enzyme systems of clinical relevance, including SARS-CoV-2 proteins.

Structural characterisation of inhibitory and non-inhibitory MMP-9-TIMP-1 complexes and implications for regulatory mechanisms of MMP-9

MMP-9 plays a number of important physiological functions but is also responsible for many pathological processes, including cancer invasion, metastasis, and angiogenesis. It is, therefore, crucial to understand its enzymatic activity, including activation and inhibition mechanisms. This enzyme may also be partially involved in the "cytokine storm" that is characteristic of COVID-19 disease (SARS-CoV-2), as well as in the molecular mechanisms responsible for lung fibrosis. Due to the variety of processing pathways involving MMP-9 in biological systems and its uniqueness due to the O-glycosylated domain (OGD) and fibronectin-like (FBN) domain, specific interactions with its natural TIMP-1 inhibitor should be carefully studied, because they differ significantly from other homologous systems. In particular, earlier experimental studies have indicated that the newly characterised circular form of a proMMP-9 homotrimer exhibits stronger binding properties to TIMP-1 compared to its monomeric form. However, molecular structures of the complexes and the binding mechanisms remain unknown. The purpose of this study is to fill in the gaps in knowledge. Molecular modelling methods are applied to build the inhibitory and non-inhibitory MMP-9-TIMP-1 complexes, which allows for a detailed description of these structures and should allow for a better understanding of the regulatory processes in which MMP-9 is involved.

The Effect of a Novel c.820C>T (Arg274Trp) Mutation in the Mitofusin 2 Gene on Fibroblast Metabolism and Clinical Manifestation in a Patient

Charcot-Marie-Tooth disease type 2A (CMT2A) is an autosomal dominant axonal peripheral neuropathy caused by mutations in the mitofusin 2 gene (MFN2). Mitofusin 2 is a GTPase protein present in the outer mitochondrial membrane and responsible for regulation of mitochondrial network architecture via the fusion of mitochondria. As that fusion process is known to be strongly dependent on the GTPase activity of mitofusin 2, it is postulated that the MFN2 mutation within the GTPase domain may lead to impaired GTPase activity, and in turn to mitochondrial dysfunction. The work described here has therefore sought to verify the effects of MFN2 mutation within its GTPase domain on mitochondrial and endoplasmic reticulum morphology, as well as the mtDNA content in a cultured primary fibroblast obtained from a CMT2A patient harboring a de novo Arg274Trp mutation. In fact, all the parameters studied were affected significantly by the presence of the mutant MFN2 protein. However, using the stable model for mitofusin 2 obtained by us, we were next able to determine that the Arg274Trp mutation does not impact directly upon GTP binding. Such results were also confirmed for GTP-hydrolysis activity of MFN2 protein in patient fibroblast. We therefore suggest that the biological malfunctions observable with the disease are not consequences of impaired GTPase activity, but rather reflect an impaired contribution of the GTPase domain to other MFN2 activities involving that region, for example protein-protein interactions.

A QM-MD simulation approach to the analysis of FRET processes in (bio)molecular systems. A case study: complexes of E. coli purine nucleoside phosphorylase and its mutants with formycin A

Predicting FRET pathways in proteins using computer simulation techniques is very important for reliable interpretation of experimental data. A novel and relatively simple methodology has been developed and applied to purine nucleoside phosphorylase (PNP) complexed with a fluorescent ligand — formycin A (FA). FRET occurs between an excited Tyr residue (D*) and FA (A). This study aims to interpret experimental data that, among others, suggests the absence of FRET for the PNPF159A mutant in complex with FA, based on novel theoretical methodology. MD simulations for the protein molecule containing D*, and complexed with A, are carried out. Interactions of D* with its molecular environment are accounted by including changes of the ESP charges in S₁, compared to S₀, and computed at the SCF-CI level. FRET probability W_F depends on the inverse six-power of the D*-A distance, R_da. The orientational factor 0 < k² < 4 between D* and A is computed and included in the analysis. Finally W_F is time-averaged over the MD trajectories resulting in its mean value. The red-shift of the tyrosinate anion emission and thus lack of spectral overlap integral and thermal energy dissipation are the reasons for the FRET absence in the studied mutants at pH 7 and above. The presence of the tyrosinate anion results in a competitive energy dissipation channel and red-shifted emission, thus in consequence in the absence of FRET. These studies also indicate an important role of the phenyl ring of Phe159 for FRET in the wild-type PNP, which does not exist in the Ala159 mutant, and for the effective association of PNP with FA. In a more general context, our observations point out very interesting and biologically important properties of the tyrosine residue in its excited state, which may undergo spontaneous deprotonation in the biomolecular systems, resulting further in unexpected physical and/or biological phenomena. Until now, this observation has not been widely discussed in the literature.

Two novel presenilin 1 gene mutations connected with frontotemporal dementia-like clinical phenotype: Genetic and bioinformatic assessment

Mutations in the amyloid precursor protein (APP), presenilin 1 (PSEN1) and presenilin 2 (PSEN2) genes are associated with early-onset familial Alzheimer's disease (EOAD). There are several reports describing mutations in PSEN1 in cases with frontotemporal dementia (FTD). We identified two novel mutations in the PSEN1 gene: L226F and L424H. The first mutation was detected in a patient with a clinical diagnosis of FTD and a post-mortem diagnosis of AD. The second mutation is connected with a clinical phenotype of variant AD with strong FTD signs. In silico modeling revealed that the mutations, as well as mutations used for comparison (F177L and L424R), change the local structure, stability and/or properties of the transmembrane regions of the presenilin 1 protein (PS1). In contrast, a silent non-synonymous substitution F175S is eclipsed by external residues and has no influence on PS1 interfacial surface. We suggest that in silico analysis of PS1 substitutions can be used to characterize novel PSEN1 mutations, to discriminate between silent polymorphisms and a potential disease-causing mutation. We also propose that PSEN1 mutations should be considered in FTD patients with no MAPT mutations.

Pulling single bacteriorhodopsin out of a membrane: Comparison of simulation and experiment

Mechanical unfolding of single bacteriorhodopsins from a membrane bilayer is studied using molecular dynamics simulations. The initial conformation of the lipid membrane is determined through all-atom simulations and then its coarse-grained representation is used in the studies of stretching. A Go-like model with a realistic contact map and with Lennard-Jones contact interactions is applied to model the protein-membrane system. The model qualitatively reproduces the experimentally observed differences between force-extension patterns obtained on bacteriorhodopsin at different temperatures and predicts a lack of symmetry in the choice of the terminus to pull by. It also illustrates the decisive role of the interactions of the protein with the membrane in determining the force pattern and thus the stability of transmembrane proteins.