Sequence, Structural Analysis and Metrics to Define the Unique Dynamic Features of the Flap Regions Among Aspartic Proteases

Extracting binding energies and binding modes from biomolecular simulations of fragment binding to endothiapepsin

Fragment-based drug discovery (FBDD) aims to discover a set of small binding fragments that may be subsequently linked together. Therefore, in-depth knowledge of the individual fragments' structural and energetic binding properties is essential. In addition to experimental techniques, the direct simulation of fragment binding by molecular dynamics (MD) simulations became popular to characterize fragment binding. However, former studies showed that long simulation times and high computational demands per fragment are needed, which limits applicability in FBDD. Here, we performed short, unbiased MD simulations of direct fragment binding to endothiapepsin, a well-characterized model system of pepsin-like aspartic proteases. To evaluate the strengths and limitations of short MD simulations for the structural and energetic characterization of fragment binding, we predicted the fragments' absolute free energies and binding poses based on the direct simulations of fragment binding and compared the predictions to experimental data. The predicted absolute free energies are in fair agreement with the experiment. Combining the MD data with binding mode predictions from molecular docking approaches helped to correctly identify the most promising fragments for further chemical optimization. Importantly, all computations and predictions were done within 5 days, suggesting that MD simulations may become a viable tool in FBDD projects.

pH-Dependent Structural Dynamics of Cathepsin D-Family Aspartic Peptidase of Clonorchis sinensis

Cathepsin D (CatD; EC 3.4.23.5) family peptidases of parasitic organisms are regarded as potential drug targets as they play critical roles in the physiology and pathobiology of parasites. Previously, we characterized the biochemical features of cathepsin D isozyme 2 (CatD2) in the carcinogenic liver fluke Clonorchis sinensis (CsCatD2). In this study, we performed all-atomic molecular dynamics simulations by applying different systems for the ligand-free/bound forms under neutral and acidic conditions to investigate the pH-dependent structural alterations and associated functional changes in CsCatD2. CsCatD2 showed several distinctive characteristics as follows: (1) acidic pH caused major conformational transitions from open to closed state in this enzyme; (2) during 30–36-ns simulations, acidic pH contributed significantly to the formation of rigid β-sheets around the catalytic residue Asp₂₁₉, higher occupancy (0% to 99%) of hydrogen bond than that of Asp₃₃, and enhanced stabilization of the CsCatD2-inhibtor complex; (3) neutral pH-induced displacement of the N-terminal part to hinder the accessibility of the active site and open allosteric site of this enzyme; and (4) the flap dynamics metrics, including distance (d₁), TriCα angles (θ₁ and θ₂), and dihedral angle (ϕ), account for the asymmetrical twisting motion of the active site of this enzyme. These findings provide an in-depth understanding of the pH-dependent structural dynamics of free and bound forms of CsCatD2 and basic information for the rational design of an inhibitor as a drug targeting parasitic CatD.

A target safety assessment of the potential toxicological risks of targeting plasmepsin IX/X for the treatment of malaria

The aspartic proteases plasmepsin IX/X are important antimalarial drug targets due to their specificity to the malaria parasite and their vital role as mediators of disease progression. Focusing on parasite-specific targets where no human homologue exists reduces the possibility of on-target drug toxicity. However, there is a risk of toxicity driven by inadequate selectivity for plasmepsins IX/X in Plasmodium over related mammalian aspartic proteases. Of these, CatD/E may be of most toxicological relevance as CatD is a ubiquitous lysosomal enzyme present in most cell types and CatE is found in the gut and in erythrocytes, the clinically significant site of malarial infection. Based on mammalian aspartic protease physiology and adverse drug reactions (ADRs) to FDA-approved human immunodeficiency virus (HIV) aspartic protease inhibitors, we predicted several potential toxicities including β-cell and congenital abnormalities, hypotension, hypopigmentation, hyperlipidaemia, increased infection risk and respiratory, renal, gastrointestinal, dermatological, and other epithelial tissue toxicities. These ADRs to the HIV treatments are likely to be a result of host aspartic protease inhibition due a lack of specificity for the HIV protease; plasmepsins are much more closely related to human CatD than to HIV proteinase. Plasmepsin IX/X inhibition presents an opportunity to specifically target Plasmodium as an effective antimalarial treatment, providing adequate selectivity can be obtained. Potential plasmepsin IX/X inhibitors should be assayed for inhibitory activity against the main human aspartic proteases and particularly CatD/E. An investigative rodent study conducted early in drug discovery would serve as an initial risk assessment of the potential hazards identified.

Weak spots inhibition in the Mycobacterium tuberculosis antigen 85C target for antitubercular drug design through selective irreversible covalent inhibitor-SER124

Mycobacterium tuberculosis (Mtb) encoded secreted antigen 85 enzymes (Ag85A/Ag85B/Ag85C) play that critical roles in the virulence, survival and drug-resistant TB of the pathogen. Ag85 proteins are potential antitubercular drug targets because they are essential in the catalytic synthesis of trehalose moieties and mycolic acid attachment to the Mtb cell wall. Recently, experimental protocols led to the discovery of a selective covalent Ag85 inhibitor, β-isomer monocyclic enolphosphorus Cycliphostin (CyC_8β) compound, which targets the Ag85 serine 124 to exhibit a promising therapeutic activity. For the first time, our study unravelled the structural features among Mtb Ag85C homologs and motions and dynamics of Ag85C when the CyC_8β bound covalently and in open model conformations to the protein using bioinformatics tools and integrated Molecular dynamics simulations. Comparative Ag85C sequence analysis revealed conserved regions; 70% active site, 90% Adeniyi loop L1 and 50% loop L2, which acts as a switch between open and closed conformations. The average C-α atoms RMSD (2.05 Å) and RMSF (0.9 Å) revealed instability and high induced flexibility in the CyC_8β covalent-bound compared to the apo and open model systems, which displayed more stability and lower fluctuations. DSSP showed structural transitions of α-helices to bend and loops to 3₁₀-helices in the bound systems. SASA of CyC_8β covalent bound showed active site hydrophobic residues exposure to huge solvent. Therefore, these findings present the potential opportunity hotspots in Ag85C protein that would aid the structure-based design of novel chemical entities capable of resulting in potent antitubercular drugs.Communicated by Ramaswamy H. Sarma.

Molecular dynamic mechanism(s) of inhibition of bioactive antiviral phytochemical compounds targeting cytochrome P450 3A4 and P-glycoprotein

P-glycoprotein (ABCB1) and cytochrome P450 3A4 (CYP3A4) metabolize almost all known human immunodeficiency virus' protease inhibitor drugs (PIs). Over induction of these proteins' activities has been linked to rapid metabolism of PIs which are then pumped out of the circulatory system, eventually leading to drug-resistance in HIV-positive patients. This study aims to determine, with the use of computational tools, the inhibitory potential of four phytochemical compounds (PCs) (epigallocatechin gallate (EGCG), kaempferol-7-glucoside (K7G), luteolin (LUT) and ellagic acid (EGA)) in inhibiting the activities of these drug-metabolizing proteins. The comparative analysis of the MM/GBSA results revealed that the binding affinity (ΔG_bind) of EGCG and K7G for CYP3A4 and ABCB1 are higher than LUT and EGA and fall between the ΔG_bind of the inhibitors of CYP3A4 and ABCB1 (Ritonavir (strong inhibitor) and Lopinavir (moderate inhibitor)). The structural analysis (RMSD, RMSF, RoG and protein-ligand interaction plots) also confirmed that EGCG and K7G showed similar inhibitory activities with the inhibitors. The study has shown that EGCG and K7G have inhibitory activities against the two proteins and assumes they could decrease intracellular efflux of PIs, consequently increasing the optimal concentration of PIs in the systemic circulation.Communicated by Ramaswamy H. Sarma.

Delving into the Characteristic Features of "Menace" Mycobacterium tuberculosis Homologs: A Structural Dynamics and Proteomics Perspectives

The global increase in the morbidity/mortality rate of Mycobacterial infections, predominantly renascent tuberculosis, leprosy, and Buruli ulcers have become worrisome over the years. More challenging is the incidence of resistance mediated by mutant Mycobacterium strains against front-line antitubercular drugs. Homologous to all Mycobacteria species is the GlcNAc-6-phosphate deacetylase (NagA) which catalyzes essential amino sugars synthesis required for cell wall architecture, hence, metamorphosing into an important pharmacological target for curtailing virulence and drug-resistance. This study used integrated bioinformatics methods, MD simulations, and DynaMut and PolyPhen2 to; explore unique features, monitor dynamics, and analyze the functional impact of non-synonymous single-nucleotide polymorphisms of the six NagA of most ruinous Mycobacterium species; tuberculosis (Mtb), smegmatis (MS), marinum (MM), ulcerans, africanum, and microti respectively. This approach is essential for multi-targeting and could result in the identification of potential polypharmacological antitubercular compounds. Comparative sequential analyses revealed ≤ 50% of the overall structure, including the catalytic Asp267 and reactive Cys131, remained conserved. Interestingly, MS-NagA and MM-NagA possess unique hydrophobic isoleucine (Ile) residues at their active sites in contrast to leucine (Leu) found in other variants. More so, unique to the active sites of the NagA is a 'subunit loop' that covers the active site; probably crucial in binding (entry and exit) mechanisms of targeted NagA inhibitors. Relatively, nsSNP mutations exerted a destabilizing effect on the native NagA conformation. Structural and dynamical insights provided, basically pin-pointed the "Achilles' heel" explorable for the rational drug design of target-specific 'NagA' inhibitors potent against a wide range of mycobacterial diseases.

A flap motif in human serine hydroxymethyltransferase is important for structural stabilization, ligand binding, and control of product release

Human cytosolic serine hydroxymethyltransferase (hcSHMT) is a promising target for anticancer chemotherapy and contains a flexible "flap motif" whose function is yet unknown. Here, using size-exclusion chromatography, analytical ultracentrifugation, small-angle X-ray scattering (SAXS), molecular dynamics (MD) simulations, and ligand-binding and enzyme-kinetic analyses, we studied the functional roles of the flap motif by comparing WT hcSHMT with a flap-deleted variant (hcSHMT/Δflap). We found that deletion of the flap results in a mixture of apo-dimers and holo-tetramers, whereas the WT was mostly in the tetrameric form. MD simulations indicated that the flap stabilizes structural compactness and thereby enhances oligomerization. The hcSHMT/Δflap variant exhibited different catalytic properties in (6S)-tetrahydrofolate (THF)-dependent reactions compared with the WT but had similar activity in THF-independent aldol cleavage of β-hydroxyamino acid. hcSHMT/Δflap was less sensitive to THF inhibition than the WT (K_i of 0.65 and 0.27 mm THF at pH 7.5, respectively), and the THF dissociation constant of the WT was also 3-fold lower than that of hcSHMT/Δflap, indicating that the flap is important for THF binding. hcSHMT/Δflap did not display the burst kinetics observed in the WT. These results indicate that, upon removal of the flap, product release is no longer the rate-limiting step, implying that the flap is important for controlling product release. The findings reported here improve our understanding of the functional roles of the flap motif in hcSHMT and provide fundamental insight into how a flexible loop can be involved in controlling the enzymatic reactions of hcSHMT and other enzymes.

A dual target of Plasmepsin IX and X: Unveiling the atomistic superiority of a core chemical scaffold in malaria therapy

Plasmepsin IX and X, members of the prominent aspartic family of proteases whose function were hitherto unknown have only recently been established as key mediators of erythrocyte invasion and egress of the virulent malarial parasite. Inhibitor 49c, a potent antimalarial peptidomimetic inhibitor initially developed to target Plasmepsin II has lately been proven to exhibit potent inhibitory activity against Plasmepsin IX and X. However, the molecular and structural dynamics supporting its inhibitory activity remain inconclusive. Hindering the motion of the flap and hinge region of an aspartic protease remains essential for disabling the catalytic activity of the enzyme. Integrating molecular dynamic simulations coupled with other advanced biocomputational tools, we reveal the enhanced structural mechanistic competence of 49c in complex with Plasmepsin IX and X relative to Pepstatin. Pepstatin, a known aspartic protease inhibitor which actively hinders the opening and closing of the flap tip and flexible loop and consequently limits access to the catalytic aspartic residues, however, its administration has been related to elevated levels of toxicity. Thermodynamic calculations reveal a higher relative binding free energy associated with Plasmepsin IX and X in complex with 49c as opposed to Pepstatin. A relatively compact and structurally rigid 49c bound complexes sequel into the restriction of the flap and hinge residues by restraining cohesive movement, consequently hindering their "twisting motion" from transpiring. Findings unveil an atomistic perspective into the structural superiority of 49c in complex with Plasmepsin IX and X.