Targeting the cell wall of Mycobacterium tuberculosis: a molecular modeling investigation of the interaction of imipenem and meropenem with L,D-transpeptidase 2

Role of UDP-N-acetylmuramic acid in the regulation of MurA activity revealed by molecular dynamics simulations

The peptidoglycan biosynthesis pathway plays a vital role in bacterial cells, and facilitates peptidoglycan layer formation, a fundamental structural component of the bacterial cell wall. The enzymes in this pathway are candidates for antibiotic development, as most do not have mammalian homologues. The UDP-N-acetylglucosamine (UNAG) enolpyruvyl transferase enzyme (MurA) in the peptidoglycan pathway cytoplasmic step is responsible for the phosphoenolpyruvate (PEP)-UNAG catalytic reaction, forming UNAG enolpyruvate and inorganic phosphate. Reportedly, UDP-N-acetylmuramic acid (UNAM) binds tightly to MurA forming a dormant UNAM-PEP-MurA complex and acting as a MurA feedback inhibitor. MurA inhibitors are complex, owing to competitive binding interactions with PEP, UNAM, and UNAG at the MurA active site. We used computational methods to explore UNAM and UNAG binding. UNAM showed stronger hydrogen-bond interactions with the Arg120 and Arg91 residues, which help to stabilize the closed conformation of MurA, than UNAG. Binding free energy calculations using end-point computational methods showed that UNAM has a higher binding affinity than UNAG, when PEP is attached to Cys115. The unbinding process, simulated using τ-random acceleration molecular dynamics, showed that UNAM has a longer relative residence time than UNAG, which is related to several complex dissociation pathways, each with multiple intermediate metastable states. This prevents the loop from opening and exposing the Arg120 residue to accommodate UNAG and potential new ligands. Moreover, we demonstrate the importance of Cys115-linked PEP in closed-state loop stabilization. We provide a basis for evaluating novel UNAM analogues as potential MurA inhibitors. PUBLIC SIGNIFICANCE: MurA is a critical enzyme involved in bacterial cell wall biosynthesis and is involved in antibiotic resistance development. UNAM can remain in the target protein's active site for an extended time compared to its natural substrate, UNAG. The prolonged interaction of this highly stable complex known as the 'dormant complex' comprises UNAM-PEP-MurA and offers insights into antibiotic development, providing potential options against drug-resistant bacteria and advancing our understanding of microbial biology.

Anxiolytic and anticonvulsant effect of Ibuprofen derivative through GABAergic neuromodulation in adult Zebrafish

Anxiety and epilepsy affect millions of people worldwide, and the treatment of these pathologies involves the use of Benzodiazepines, drugs that have serious adverse effects such as dependence and sedation, so the discovery of new anxiolytic and antiepileptic drugs are necessary. Many routes for synthesizing ibuprofen derivatives have been developed, and these derivatives have shown promising pharmacological effects. Therefore, this study aims to evaluate its anxiolytic and anticonvulsant effect against the adult Zebrafish animal model of Ibuprofen (IBUACT) and its interaction with the GABAergic receptor through in silico studies. The light/dark preference test (Scototaxis test) was used to evaluate the anxiolytic behavior of adult Zebrafish acutely treated with IBUACT and Diazepam, and their anticonvulsant effects were investigated through the pentylenetetrazol (PTZ)-induced seizure model. Animals treated with IBUACT showed anxiolytic behavior similar to Diazepam, and pretreatment with flumazenil reversed this behavior. PTZ-induced seizures were delayed by IBUACT in all three stages and were shown to bind strongly in the Diazepam region of GABAA. In addition, this work presents evidence of new pharmacological applications of ibuprofen derivative in pathologies of the central nervous system (CNS), opening the horizon for new studies.Communicated by Ramaswamy H. Sarma.

Structure and Function of L,D- and D,D-Transpeptidase Family Enzymes from Mycobacterium tuberculosis

Peptidoglycan, the exoskeleton of bacterial cell and an essential barrier that protects the cell, is synthesized by a pathway where the final steps are catalysed by transpeptidases. Knowledge of the structure and function of these vital enzymes that generate this macromolecule in M. tuberculosis could facilitate the development of potent lead compounds against tuberculosis. This review summarizes the experimental and computational studies to date on these aspects of transpeptidases in M. tuberculosis that have been identified and validated. The reported structures of L,D- and D,D-transpeptidases, as well as their functionalities, are reviewed and the proposed enzymatic mechanisms for L,D-transpeptidases are summarized. In addition, we provide bioactivities of known tuberculosis drugs against these enzymes based on both experimental and computational approaches. Advancing knowledge about these prominent targets supports the development of new drugs with novel inhibition mechanisms overcoming the current need for new drugs against tuberculosis.

Competition Between Phenothiazines and BH3 Peptide for the Binding Site of the Antiapoptotic BCL-2 Protein

The study of proteins and mechanisms involved in the apoptosis and new knowledge about cancer's biology are essential for planning new drugs. Tumor cells develop several strategies to gain proliferative advantages, including molecular alterations to evade from apoptosis. Failures in apoptosis could contribute to cancer pathogenesis, since these defects can cause the accumulation of dividing cells and do not remove genetic variants that have malignant potential. The apoptosis mechanism is composed by proteins that are members of BCL-2 and cysteine-protease families. BH3-only peptides are the “natural” intracellular ligands of BCL-2 family proteins. On the other hand, studies have proved that phenothiazine compounds influence the induction of cellular death. To understand the characteristics of phenothiazines and their effects on tumoral cells and organelles involved in the apoptosis, as well as evaluating their pharmacologic potential, we have carried out computational simulation with the purpose of relating the structures of the phenothiazines with their biological activity. Since the tridimensional (3D) structure of the target protein is known, we have employed the molecular docking approach to study the interactions between compounds and the protein's active site. Hereafter, the molecular dynamics technique was used to verify the temporal evolution of the BCL-2 complexes with phenothiazinic compounds and the BH3 peptide, the stability and the mobility of these molecules in the BCL-2 binding site. From these results, the calculation of binding free energy between the compounds and the biological target was carried out. Thus, it was possible to verify that thioridazine and trifluoperazine tend to increase the stability of the BCL-2 protein and can compete for the binding site with the BH3 peptide.

Identification of potent L,D-transpeptidase 5 inhibitors for Mycobacterium tuberculosis as potential anti-TB leads: virtual screening and molecular dynamics simulations

Virtual screening is a useful in silico approach to identify potential leads against various targets. It is known that carbapenems (doripenem and faropenem) do not show any reasonable inhibitory activities against L,D-transpeptidase 5 (Ldt_Mt5) and also an adduct of meropenem exhibited slow acylation. Since these drugs are active against L,D-transpeptidase 2 (Ldt_Mt2), understanding the differences between these two enzymes is essential. In this study, a ligand-based virtual screening of 12,766 compounds followed by molecular dynamics (MD) simulations was applied to identify potential leads against Ldt_Mt5. To further validate the obtained virtual screening ranking for Ldt_Mt5, we screened the same libraries of compounds against Ldt_Mt2 which had more experimetal and calculated binding energies reported. The observed consistency between the binding affinities of Ldt_Mt2 validates the obtained virtual screening binding scores for Ldt_Mt5. We subjected 37 compounds with docking scores ranging from - 7.2 to - 9.9 kcal mol^-1 obtained from virtual screening for further MD analysis. A set of compounds (n = 12) from four antibiotic classes with ≤ - 30 kcal mol^-1 molecular mechanics/generalized born surface area (MM-GBSA) binding free energies (ΔG_bind) was characterized. A final set of that, all β-lactams (n = 4), was considered. The outcome of this study provides insight into the design of potential novel leads for Ldt_Mt5. Graphical abstract.

The 1-β-methyl group confers a lower affinity of l,d-transpeptidase LdtMt2 for ertapenem than for imipenem

L,D-transpeptidases, widely distributed in bacteria and even in the difficult-to-treat ESKAPE pathogens, can confer antibacterial resistance against the traditional β-lactam antibiotics through bypass of the 4 → 3 transpeptide linkage. Ldt_Mt2, a l,d-transpeptidase in Mycobacteria tuberculosis, is essential for bacterial virulence and is considered as a potential anti-tuberculosis target inhibited by carbapenems. Diverse interaction modes between carbapenems and Ldt_Mt2 have been reported, there are only limited evidences to validate those interaction modes. Herein, we identified the stable binding states of two carbapenems, imipenem and ertapenem, via crystallographic and biochemical studies, discovered that they adopt similar binding conformations. We further demonstrate the absence of the 1-β-methyl group in imipenem and the presence of both Y308 and Y318 residues in Ldt_Mt2 synergistically resulted in one order of magnitude higher affinity for imipenem than ertapenem. Our study provides a structural basis for the rational drug design and evolvement of novel carbapenems against bacterial L,D-transpeptidases.

Inhibition mechanism of L,D-transpeptidase 5 in presence of the β-lactams using ONIOM method

Tuberculosis (TB) is one of the world's deadliest diseases resulting from infection by the bacterium, Mycobacterium tuberculosis (M.tb). The L,D-transpeptidase enzymes catalyze the synthesis of 3 → 3 transpeptide linkages which are predominant in the peptidoglycan of the M.tb cell wall. Carbapenems is class of β-lactams that inactivate L,D-transpeptidases by acylation, although differences in antibiotic side chains modulate drug binding and acylation rates. Herein, we used a two-layered our Own N-layer integrated Molecular Mechanics ONIOM method to investigate the catalytic mechanism of L,D-transpeptidase 5 (Ldt_Mt5) by β-lactam derivatives. Ldt_Mt5 complexes with six β-lactams, ZINC03788344 (1), ZINC02462884 (2), ZINC03791246 (3), ZINC03808351 (4), ZINC03784242 (5) and ZINC02475683 (6) were simulated. The QM region (high-level) comprises the β-lactam, one water molecule and the Cys360 catalytic residue, while the rest of the Ldt_Mt5 residues were treated with AMBER force field. The activation energies (ΔG^#) were calculated with B3LYP, M06-2X and ωB97X density functionals with 6-311++G(2d, 2p) basis set. The ΔG^# for the acylation of Ldt_Mt5 by the selected β-lactams were obtained as 13.67, 20.90, 22.88, 24.29, 27.86 and 28.26 kcal mol^-1respectively. Several of the compounds showed an improved ΔG^# when compared to the previously calculated energies for imipenem and meropenem for the acylation step for Ldt_Mt5. This model provides further validation of the catalytic inhibition mechanism of LDTs with atomistic detail.

Molecular insight on the non-covalent interactions between carbapenems and L,D-transpeptidase 2 from Mycobacterium tuberculosis: ONIOM study

Tuberculosis remains a dreadful disease that has claimed many human lives worldwide and elimination of the causative agent Mycobacterium tuberculosis also remains elusive. Multidrug-resistant TB is rapidly increasing worldwide; therefore, there is an urgent need for improving the current antibiotics and novel drug targets to successfully curb the TB burden. L,D-Transpeptidase 2 is an essential protein in Mtb that is responsible for virulence and growth during the chronic stage of the disease. Both D,D- and L,D-transpeptidases are inhibited concurrently to eradicate the bacterium. It was recently discovered that classic penicillins only inhibit D,D-transpeptidases, while L,D-transpeptidases are blocked by carbapenems. This has contributed to drug resistance and persistence of tuberculosis. Herein, a hybrid two-layered ONIOM (B3LYP/6-31G+(d): AMBER) model was used to extensively investigate the binding interactions of Ldt_Mt2 complexed with four carbapenems (biapenem, imipenem, meropenem, and tebipenem) to ascertain molecular insight of the drug-enzyme complexation event. In the studied complexes, the carbapenems together with catalytic triad active site residues of Ldt_Mt2 (His187, Ser188 and Cys205) were treated at with QM [B3LYP/6-31+G(d)], while the remaining part of the complexes were treated at MM level (AMBER force field). The resulting Gibbs free energy (ΔG), enthalpy (ΔH) and entropy (ΔS) for all complexes showed that the carbapenems exhibit reasonable binding interactions towards Ldt_Mt2. Increasing the number of amino acid residues that form hydrogen bond interactions in the QM layer showed significant impact in binding interaction energy differences and the stabilities of the carbapenems inside the active pocket of Ldt_Mt2. The theoretical binding free energies obtained in this study reflect the same trend of the experimental  observations. The electrostatic, hydrogen bonding and Van der Waals interactions between the carbapenems and Ldt_Mt2 were also assessed. To further examine the nature of intermolecular interactions for carbapenem-Ldt_Mt2 complexes, AIM and NBO analysis were performed for the QM region (carbapenems and the active residues of Ldt_Mt2) of the complexes. These analyses revealed that the hydrogen bond interactions and charge transfer from the bonding to anti-bonding orbitals between catalytic residues of the enzyme and selected ligands enhances the binding and stability of carbapenem-Ldt_Mt2 complexes. The two-layered ONIOM (B3LYP/6-31+G(d): Amber) model was used to evaluate the efficacy of FDA approved carbapenems antibiotics towards Ldt_Mt2.

Computational analyses of interactions between ALK-5 and bioactive ligands: insights for the design of potential anticancer agents

Activin Receptor-Like Kinase 5 (ALK-5) is related to some types of cancer, such as breast, lung, and pancreas. In this study, we have used molecular docking, molecular dynamics simulations, and free energy calculations in order to explore key interactions between ALK-5 and six bioactive ligands with different ranges of biological activity. The motivation of this work is the lack of crystal structure for inhibitor-protein complexes for this set of ligands. The understanding of the molecular structure and the protein-ligand interaction could give support for the development of new drugs against cancer. The results show that the calculated binding free energy using MM-GBSA, MM-PBSA, and SIE is correlated with experimental data with r² = 0.88, 0.80, and 0.94, respectively, which indicates that the calculated binding free energy is in excellent agreement with experimental data. In addition, the results demonstrate that H bonds with Lys232, Glu245, Tyr249, His283, Asp351, and one structural water molecule play an important role for the inhibition of ALK-5. Overall, we discussed the main interactions between ALK-5 and six inhibitors that may be used as starting points for designing new molecules to the treatment of cancer.

Differential flap dynamics in l,d-transpeptidase2 from mycobacterium tuberculosis revealed by molecular dynamics

Despite the advances in tuberculosis treatment, TB is still one the most deadly infectious diseases and remains a major global health quandary. Mycobacterium tuberculosis (Mtb) is the only known mycobacterium with a high content of 3→3 crosslinks in the biosynthesis of peptidoglycan, which is negligible in most bacterial species. An Mtb lacking Ldt_Mt2 leads to alteration of the colony morphology and loss of virulence which makes this enzyme an attractive target. Regardless of the vital role of Ldt_Mt2 for cell wall survival, the impact of ligand binding on the dynamics of the β-hairpin flap is still unknown. Understanding the structural and dynamical behaviour of the flap regions provides clear insight into the design of the effective inhibitors against Ldt_Mt2. Carbapenems, an specific class of β-lactam family, have been shown to inactivate this enzyme. Herein a comprehensive investigation of the flap dynamics of Ldt_Mt2 complex with substrate and three carbapenems namely, ertapenem, imipenem and meropenem is discussed and analyzed for the first account using 140 ns molecular dynamics simulations. The structural features (RMSD, RMSF and R_g) derived by MD trajectories were analyzed. Distance analysis, particularly tip-tip SER135-ASN167 index, identified conformational changes in terms of flap opening and closure within binding process. Principal component analysis (PCA) was employed to qualitatively understand the divergent effects of different inhibitors on the dominant motion of each residue. To probe different internal dynamics induced by ligand binding, dynamic cross-correlation marix (DCCM) analysis was used. The binding free energies of the selected complexes were assessed using MM-GBSA method and per residue free energy decomposition analysis were performed to characterize the contribution of the key residues to the total binding free energies.

Non-classical transpeptidases yield insight into new antibacterials

Bacterial survival requires an intact peptidoglycan layer, a three-dimensional exoskeleton that encapsulates the cytoplasmic membrane. Historically, the final steps of peptidoglycan synthesis are known to be carried out by D,D-transpeptidases, enzymes that are inhibited by the β-lactams, which constitute >50% of all antibacterials in clinical use. Here, we show that the carbapenem subclass of β-lactams are distinctly effective not only because they inhibit D,D-transpeptidases and are poor substrates for β-lactamases, but primarily because they also inhibit non-classical transpeptidases, namely the L,D-transpeptidases, which generate the majority of linkages in the peptidoglycan of mycobacteria. We have characterized the molecular mechanisms responsible for inhibition of L,D-transpeptidases of Mycobacterium tuberculosis and a range of bacteria including ESKAPE pathogens, and used this information to design, synthesize and test simplified carbapenems with potent antibacterial activity.

Structure and dynamics of the multi-domain resuscitation promoting factor RpfB from Mycobacterium tuberculosis

RpfB is multidomain protein that is crucial for Mycobacterium tuberculosis resuscitation from dormancy. This protein cleaves cell wall peptidoglycan, an essential bacterial cell wall polymer formed by glycan chains of β-(1-4)-linked-N-acetylglucosamine (GlcNAc) and N-acetylmuramic acid (MurNAc) cross-linked by short peptide stems. RpfB is structurally complex being composed of five distinct domains, namely a catalytic, a G5 and three DUF348 domains. Here, we have undertaken a combined experimental and computation structural investigations on the entire protein to gain insights into its structure-function relationships. CD spectroscopy and light scattering experiments have provided insights into the protein fold stability and into its oligomeric state. Using the available structure information, we modeled the entire protein structure, which includes the two DUF348 domains whose structure is experimentally unknown, and we analyzed the dynamic behavior of RpfB using molecular dynamics simulations. Present results highlight an intricate mutual influence of the dynamics of the different protein domains. These data provide interesting clues on the functional role of non-catalytic domains of RpfB and on the mechanism of peptidoglycan degradation necessary to resuscitation of M. tuberculosis.