Generation of AMBER force field parameters for zinc centres of M1 and M17 family aminopeptidases

Identification and Validation of Compounds Targeting Leishmania major Leucyl-Aminopeptidase M17

Leishmaniasis is a neglected tropical disease; there is currently no vaccine and treatment is reliant upon a handful of drugs suffering from multiple issues including toxicity and resistance. There is a critical need for development of new fit-for-purpose therapeutics, with reduced toxicity and targeting new mechanisms to overcome resistance. One enzyme meriting investigation as a potential drug target in Leishmania is M17 leucyl-aminopeptidase (LAP). Here, we aimed to chemically validate LAP as a drug target in L. major through identification of potent and selective inhibitors. Using RapidFire mass spectrometry, the compounds DDD00057570 and DDD00097924 were identified as selective inhibitors of recombinant Leishmania major LAP activity. Both compounds inhibited in vitro growth of L. major and L. donovani intracellular amastigotes, and overexpression of LmLAP in L. major led to reduced susceptibility to DDD00057570 and DDD00097924, suggesting that these compounds specifically target LmLAP. Thermal proteome profiling revealed that these inhibitors thermally stabilized two M17 LAPs, indicating that these compounds selectively bind to enzymes of this class. Additionally, the selectivity of the inhibitors to act on LmLAP and not against the human ortholog was demonstrated, despite the high sequence similarities LAPs of this family share. Collectively, these data confirm LmLAP as a promising therapeutic target for Leishmania spp. that can be selectively inhibited by drug-like small molecules.

Understanding the structure and function of Plasmodium aminopeptidases to facilitate drug discovery

Malaria continues to be the most widespread parasitic disease affecting humans globally. As parasites develop drug resistance at an alarming pace, it has become crucial to identify novel drug targets. Over the last decade, the metalloaminopeptidases have gained importance as potential targets for new antimalarials. These enzymes are responsible for removing the N-terminal amino acids from proteins and peptides, and their restricted specificities suggest that many perform unique and essential roles within the malaria parasite. This mini-review focuses on the recent progress in structure and functional data relating to the Plasmodium metalloaminopeptidases that have been validated or shown promise as new antimalarial drug targets.

Insights into the Inhibitory Mechanisms of the Covalent Drugs for DNMT3A

The perturbations of DNA methyltransferase 3 alpha (DNMT3A) may cause uncontrolled gene expression, resulting in cancers and tumors. The DNMT inhibitors Azacytidine (AZA) and Zebularine (ZEB) inhibit the DNMT family with no specificities, and consequently would bring side effects during the treatment. Therefore, it is vital to understand the inhibitory mechanisms in DNMT3A to inform the new inhibitor design for DNMTs. Herein, we carried out molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) simulations to investigate the inhibitory mechanisms of the AZA and ZEB. The results were compared to the methyl transfer of cytosine. We showed how the AZA might stop the methyl transfer process, whereas the ZEB might be stuck in a methyl-transferred intermediate (IM3). The IM3 state then fails the elimination due to the unique protein dynamics that result in missing the catalytic water chain. Our results brought atomic-level insights into the mechanisms of the two drugs in DNMT3A, which could benefit the new generation of drug design for the DNMTs.

Omicron Binding Mode: Contact Analysis and Dynamics of the Omicron Receptor-Binding Domain in Complex with ACE2

On 26 November 2021, the WHO classified the Omicron variant of the SARS-CoV-2 virus (B.1.1.529 lineage) as a variant of concern (VOC) (COVID-19 Variant Data, Department of Health, 2022). The Omicron variant contains as many as 26 unique mutations of effects not yet determined (Venkatakrishnan, A., Open Science Framework, 2021). Out of its total of 34 Spike protein mutations, 15 are located on the receptor-binding domain (S-RBD) (Stanford Coronavirus Antiviral & Resistance Database, 2022) that directly contacts the angiotensin-converting enzyme 2 (ACE2) host receptor and is also a primary target for antibodies. Here, we studied the binding mode of the S-RBD domain of the Spike protein carrying the Omicron mutations and the globular domain of human ACE2 using molecular dynamics (MD) simulations. We identified new and key Omicron-specific interactions such as R493 (of mutation Q493R), which forms salt bridges both with E35 and D38 of ACE2, Y501 (N501Y), which forms an edge-to-face aromatic interaction with Y41, and Y505 (Y505H), which makes an H-bond with E37 and K353. The glycan chains of ACE2 also bind differently in the WT and Omicron variants in response to different charge distributions on the surface of Spike proteins. However, while the Omicron mutations considerably improve the overall electrostatic fit of the two interfaces, the total number of specific and favorable interactions between the two does not increase. The dynamics of the complexes are highly affected too, making the Omicron S-RBD:ACE2 complex more rigid; the two main interaction sites, Patches I and II, isolated in the WT complex, become connected in the Omicron complex through the alternating interaction of R493 and R498 with E35 and D38.

A metal ion-dependent conformational switch modulates activity of the Plasmodium M17 aminopeptidase

The metal-dependent M17 aminopeptidases are conserved throughout all kingdoms of life. This large enzyme family is characterized by a conserved binuclear metal center and a distinctive homohexameric arrangement. Recently, we showed that hexamer formation in Plasmodium M17 aminopeptidases was controlled by the metal ion environment, although the functional necessity for hexamer formation is still unclear. To further understand the mechanistic role of the hexameric assembly, here we undertook an investigation of the structure and dynamics of the M17 aminopeptidase from Plasmodium falciparum, PfA-M17. We describe a novel structure of PfA-M17, which shows that the active sites of each trimer are linked by a dynamic loop, and loop movement is coupled with a drastic rearrangement of the binuclear metal center and substrate-binding pocket, rendering the protein inactive. Molecular dynamics simulations and biochemical analyses of PfA-M17 variants demonstrated that this rearrangement is inherent to PfA-M17, and that the transition between the active and inactive states is metal dependent and part of a dynamic regulatory mechanism. Key to the mechanism is a remodeling of the binuclear metal center, which occurs in response to a signal from the neighboring active site and serves to moderate the rate of proteolysis under different environmental conditions. In conclusion, this work identifies a precise mechanism by which oligomerization contributes to PfA-M17 function. Furthermore, it describes a novel role for metal cofactors in the regulation of enzymes, with implications for the wide range of metalloenzymes that operate via a two-metal ion catalytic center, including DNA processing enzymes and metalloproteases.

8-Alkylmercaptocaffeine derivatives: antioxidant, molecular docking, and in-vitro cytotoxicity studies

Due to their unique pharmacological characteristics, methylxanthines are known as therapeutic agents in a fascinating range of medicinal scopes. In this report, we aimed to examine some biological effects of previously synthesized 8-alkylmercaptocaffeine derivatives. Cytotoxic and antioxidative activity of 8-alkylmercaptocaffeine derivatives were measured in malignant A549, MCF7, and C152 cell lines. Assessment of cGMP levels and caspase-3 activity were carried out using a colorimetric competitive ELISA kit. Computational approaches were employed to discover the inhibitory mechanism of synthesized compounds. Among the twelve synthesized derivatives, three compounds (C1, C5, and C7) bearing propyl, heptyl, and 3-methyl-butyl moieties showed higher and more desirable cytotoxic activity against all the studied cell lines (IC₅₀ < 100 µM). Furthermore, C5 synergistically enhanced cisplatin-induced cytotoxicity in MCF-7 cells (CI < 1). Both C5 and C7 significantly increased caspase-3 activity and intracellular cGMP levels at specific time intervals in all studied cell lines (P < 0.05). However, these derivatives did not elevate LDH leakage (P > 0.05) and exhibited no marked ameliorating effects on oxidative damage (P > 0.05). Computational studies showed that H-bond formation between the nitrogen atom in pyrazolo[4,3-D] pyrimidine moiety with Gln817 and creating a hydrophobic cavity result in the stability of the alkyl group in the PDE5A active site. We found that synthesized 8-alkylmercaptocaffeine derivatives induced cell death in different cancer cells through the cGMP pathway. These findings will help us to get a deeper insight into the role of methylxanthines as useful alternatives to conventional cancer therapeutics.

Potential effective inhibitory compounds against Prostate Specific Membrane Antigen (PSMA): A molecular docking and molecular dynamics study

One of the most prevalent cancers in men is prostate cancer and could be managed with immunotoxins or antibody treatment. Because of the substantial rise of the Prostate-Specific Antigen and the Prostate-Specific Membrane Antigen (PSMA), cancer vaccination should be rendered with these antigens. Through pharmacodynamic experiments in a library of natural compounds from ZINC database, the current research sought to identify compounds that could suppress PSMA protein. To test the most productive compounds for further research, the Library has been scanned with Pharmacophore and ADMET analysis followed by molecular docking methods in the first phase. After selecting 15 ligands with the best pose related to docking results, to evaluate the stability of the ligand-protein bounds of the compounds, a molecular dynamics simulation considering the effect of the presence of zinc ions on the protein structure was performed. The measurement of ligand binding modes and free energy has shown that four compounds, including Z10, Z06, Z01, and Z03, have formed critical interactions with the active site's residues. Besides, multiple approaches were employed to determine their inhibition rating and describe the variables that facilitate the attachment of ligands to the protein active site. The results are obtained from the MMPBSA/GBSA analysis of four selected small molecules (Z10, Z06, Z01, and Z03), which are very close to the IC50 value of reference ligand (DCIBzl); they are -13.85 kcal/mol, -12.58 kcal/mol, -10.71 kcal/mol and -9.39 kcal/mol respectively. Finally, we evaluate the results obtained from selected ligands using hydrogen bond and decomposition analyzes. We have examined the effective interactions between ligands and S1/S1'pockets in protein. Our computational results illustrate the design of more efficient inhibitors of PSMA.

A dual-functional aminopeptidase from Streptomyces canus T20 and its application in the preparation of small rice peptides

An aminopeptidase that derived from Streptomyces canus T20 (ScAP) was successfully expressed and characterized in Escherichia coli BL21 (DE3) for first time. The specific activity was 6000 U/mg, which is highest in Streptomyces aminopeptidases. Its optimal conditions were 60 °C and pH 8.0, respectively. ScAP exhibited excellent thermal and alkaline pH stability, retained 80.0% maximal activity at 50 °C for 200 h or at pH 9.0 for 24 h. Its activity observed to be complete inhibited by 0.1 mM EDTA and enhanced by Ca²⁺ and Co²⁺ to 115.4% and 104.0% respectively. ScAP also has exhibited high specificity towards rice protein on preparation of small peptides. The yield of small rice peptides achieved 66.5%, which is highest by far. Besides, ScAP have significant debittering effect on rice peptides. Results showed that bitter intensity score decreased by 49.0% with optimum condition (0.048% ScAP at 50 °C for 6 h). Therefore, ScAP as dual functional aminopeptidase of hydrolytic and debittering might have a potential application in the production of high yield and low bitterness of small rice peptides.

Development of parameters compatible with the CHARMM36 force field for [Fe4S4]2+ clusters and molecular dynamics simulations of adenosine-5'-phosphosulfate reductase in GROMACS 2019

DFT calculations were used to obtain parameters compatible with the CHARMM36 force field for iron-sulfur clusters (Fe-S) of the type [Fe₄S₄]²⁺ that are coordinated to dissimilatory adenosine-5'-phosphosulfate reductase (APSrAB). Classical molecular dynamics (MD) simulations were performed on two APSrAB systems to validate the parameters and verify the stability of the studied systems. The time analysis of the parameters inserted into the force field was in reasonable agreement with the experimental X-ray diffraction data. The analysis of the time evolution of the studied systems indicated that these systems and, in particular, the clusters in their respective cavities had a good stability and were in agreement with what was observed in previous works. The parameters obtained provide the basis for the study of APSrAB as well as other systems that contain [Fe₄S₄]²⁺ through the CHARMM36 force field.

Understanding structure-based dynamic interactions of antihypertensive peptides extracted from food sources

Functional foods are emerging as essential healthy nutritional component due to their abundant wellbeing benefits. Especially the food-derived peptides are considered as key components for playing their biologically active roles. One such robust therapeutics that already exploited with food peptides that help treating high blood pressure via targeting Angiotensin-Converting Enzyme (ACE). This in silico study demonstrated the inhibitory potential of antihypertensive peptides derived from food sources. This study involves an intensive structure-based analysis of enzyme-peptide interactions using Molecular Dynamics (MD) simulations. Interestingly, this study will help us to get deeper understanding on how food peptides achieve successful inhibition of ACE. In this study, the peptide-enzyme complexes revealed two binding pockets, A and B, on either side of the active site Zn atom. Pocket B has a smaller binding site volume than pocket A, comprised of β-sheets and the active site opening cleft. The interface of the binding sites showed that the enzyme structure was negative to neutral charge, and the peptide structure was positive to neutral charge. The dynamics of complex structures of seven highly potential peptides were performed for 20 ns each at 300 K. Comparative analysis of RMSD, RMSF and binding energies show the enzyme-peptide complexes and the overall stability of apo-enzyme. Importantly, two peptides AFKAWAVAR and IWHHTF showed the highest variation in their RMSD as compared to the apo-enzyme. This study will further be useful for the assessment of the characteristics to predict novel inhibitory peptides that can be generated from food proteins.Communicated by Ramaswamy H. Sarma.

Novel Human Aminopeptidase N Inhibitors: Discovery and Optimization of Subsite Binding Interactions

Aminopeptidase N (APN/CD13) is a zinc-dependent M1 aminopeptidase that contributes to cancer progression by promoting angiogenesis, metastasis, and tumor invasion. We have previously identified hydroxamic acid-containing analogues that are potent inhibitors of the APN homologue from the malarial parasite Plasmodium falciparum M1 aminopeptidase (PfA-M1). Herein, we describe the rationale that underpins the repurposing of PfA-M1 inhibitors as novel APN inhibitors. A series of novel hydroxamic acid analogues were developed using a structure-based design approach and evaluated their inhibition activities against APN. N-(2-(Hydroxyamino)-2-oxo-1-(3',4',5'-trifluoro-[1,1'-biphenyl]-4-yl)ethyl)-4-(methylsulfonamido)benzamide (6ad) proved to be an extremely potent inhibitor of APN activity in vitro, selective against other zinc-dependent enzymes such as matrix metalloproteases, and possessed limited cytotoxicity against Ad293 cells and favorable physicochemical and metabolic stability properties. The combined results indicate that compound 6ad may be a useful lead for the development of anticancer agents.

New Zinc Ion Parameters Suitable for Classical MD Simulations of Zinc Metallopeptidases

The main aim of this work was to find parameters for the zinc ion in human dipeptidyl peptidase III (DPP III) active site that would enable its reliable modeling. Since the parameters publicly available failed to reproduce the zinc ion coordination in the enzyme, we developed a new set of the hybrid bonded/nonbonded parameters for the zinc ion suitable for molecular modeling of the human DPP III, dynamics, and ligand binding. The parameters allowed exchange of the water molecules coordinating the zinc ion and proved to be robust enough to enable reliable modeling not only of human DPP III and its orthologues but also of the other zinc-dependent peptidases with the zinc ion coordination similar to that in dipeptidyl peptidases III, i.e., peptidases with the zinc ion coordinated with two histidines and one glutamate. The new parameters were tested on a set of 21 different systems comprising 8 different peptidases, 5 DPP III orthologues, thermolysin, neprilysin, and aminopeptidase N, and the results are summarized in the second part of the article.