Exploring the Flap Pocket of the Antimalarial Target Plasmepsin II: The “55 % Rule” Applied to Enzymes

CageCavityCalc (C3): A Computational Tool for Calculating and Visualizing Cavities in Molecular Cages

Organic(porous) and metal-organic cages are promising biomimetic platforms with diverse applications spanning recognition, sensing, and catalysis. The key to the emergence of these functions is the presence of well-defined inner cavities capable of binding a wide range of guest molecules and modulating their properties. However, despite the myriad cage architectures currently available, the rational design of structurally diverse and functional cages with specific host-guest properties remains challenging. Efficiently predicting such properties is critical for accelerating the discovery of novel functional cages. Herein, we introduce CageCavityCalc (C3), a Python-based tool for calculating the cavity size of molecular cages. The code is available on GitHub at https://github.com/VicenteMartiCentelles/CageCavityCalc. C3 utilizes a novel algorithm that enables the rapid calculation of cavity sizes for a wide range of molecular structures and porous systems. Moreover, C3 facilitates easy visualization of the computed cavity size alongside hydrophobic and electrostatic potentials, providing insights into host-guest interactions within the cage. Furthermore, the calculated cavity can be visualized using widely available visualization software, such as PyMol, VMD, or ChimeraX. To enhance user accessibility, a PyMol plugin has been created, allowing nonspecialists to use this tool without requiring computer programming expertise. We anticipate that the deployment of this computational tool will significantly streamline cage cavity calculations, thereby accelerating the discovery of functional cages.

Exploring Aspartic Protease Inhibitor Binding to Design-Selective Antimalarials

Selectivity is a major issue in the development of drugs targeting pathogen aspartic proteases. Here, we explore the selectivity-determining factors by studying specifically designed malaria aspartic protease (plasmepsin) open-flap inhibitors. Metadynamics simulations are used to uncover the complex binding/unbinding pathways of these inhibitors and describe the critical transition states in atomistic resolution. The simulation results are compared with experimentally determined enzymatic activities. Our findings demonstrate that plasmepsin inhibitor selectivity can be achieved by targeting the flap loop with hydrophobic substituents that enable ligand binding under the flap loop, as such a behavior is not observed for several other aspartic proteases. The ability to estimate the selectivity of compounds before they are synthesized is of considerable importance in drug design; therefore, we expect that our approach will be useful in selective inhibitor designs against not only aspartic proteases but also other enzyme classes.

Molecular Recognition of Nerve Agents and Their Organophosphorus Surrogates: Toward Supramolecular Scavengers and Catalysts

Nerve agents are tetrahedral organophosphorus compounds (OPs) that were developed in the last century to irreversibly inhibit acetylcholinesterase (AChE) and therefore impede neurological signaling in living organisms. Exposure to OPs leads to a rapid development of symptoms from excessive salivation, nasal congestion and chest pain to convulsion and asphyxiation which if left untreated may lead to death. These potent toxins are prepared on a large scale from inexpensive staring materials, making it feasible for terrorist groups or states to use them against military and civilians. The existing antidotes provide limited protection and are difficult to apply to a large number of affected individuals. While new prophylactics are currently being developed, there is still need for therapeutics capable of both preventing and reversing the effects of OP poisoning. In this review, we describe how the science of molecular recognition can expand the pallet of tools for rapid and safe sequestration of nerve agents.

Expanding Coefficient: A Parameter To Assess the Stability of Induced-Fit Complexes

Here we propose a new parameter, the Expanding Coefficient (EC), that can be correlated with the thermodynamic stability of supramolecular complexes governed by weak secondary interactions and obeying the induced-fit model. The EC values show a good linear relationship with the log Kapp of the respective pseudorotaxane complexes investigated. According to Cram's Principle of Preorganization, the EC can be considered an approximate mechanical measure of the host's reorganization energy cost upon adopting the final bound geometry.

Fentanyl Family at the Mu-Opioid Receptor: Uniform Assessment of Binding and Computational Analysis

Interactions of 21 fentanyl derivatives with μ-opioid receptor (μOR) were studied using experimental and theoretical methods. Their binding to μOR was assessed with radioligand competitive binding assay. A uniform set of binding affinity data contains values for two novel and one previously uncharacterized derivative. The data confirms trends known so far and thanks to their uniformity, they facilitate further comparisons. In order to provide structural hypotheses explaining the experimental affinities, the complexes of the studied derivatives with μOR were modeled and subject to molecular dynamics simulations. Five common General Features (GFs) of fentanyls' binding modes stemmed from these simulations. They include: GF1) the ionic interaction between D147 and the ligands' piperidine NH⁺ moiety; GF2) the N-chain orientation towards the μOR interior; GF3) the other pole of ligands is directed towards the receptor outlet; GF4) the aromatic anilide ring penetrates the subpocket formed by TM3, TM4, ECL1 and ECL2; GF5) the 4-axial substituent (if present) is directed towards W318. Except for the ionic interaction with D147, the majority of fentanyl-μOR contacts is hydrophobic. Interestingly, it was possible to find nonlinear relationships between the binding affinity and the volume of the N-chain and/or anilide's aromatic ring. This kind of relationships is consistent with the apolar character of interactions involved in ligand⁻receptor binding. The affinity reaches the optimum for medium size while it decreases for both large and small substituents. Additionally, a linear correlation between the volumes and the average dihedral angles of W293 and W133 was revealed by the molecular dynamics study. This seems particularly important, as the W293 residue is involved in the activation processes. Further, the Y326 (OH) and D147 (Cγ) distance found in the simulations also depends on the ligands' size. In contrast, neither RMSF measures nor D114/Y336 hydrations show significant structure-based correlations. They also do not differentiate studied fentanyl derivatives. Eventually, none of 14 popular scoring functions yielded a significant correlation between the predicted and observed affinity data (R < 0.30, n = 28).

2-Aminoquinazolin-4(3H)-one based plasmepsin inhibitors with improved hydrophilicity and selectivity

2-Aminoquinazolin-4(3H)-ones were previously discovered as perspective leads for antimalarial drug development targeting the plasmepsins. Here we report the lead optimization studies with the aim to reduce inhibitor lipophilicity and increase selectivity versus the human aspartic protease Cathepsin D. Exploiting the solvent exposed area of the enzyme provides an option to install polar groups (R¹) the 5-position of 2-aminoquinazolin-4(3H)-one to inhibitors such as carboxylic acid without scarifying enzymatic potency. Moreover, introduction of R¹ substituents increased selectivity factors of compounds in this series up to 100-fold for Plm II, IV vs CatD inhibition. The introduction of flap pocket substituent (R²) at 7-postion of 2-aminoquinazolin-4(3H)-one allows to remove Ph group from THF ring without notably impairing Plm inhibitory potency. Based on these findings, inhibitors were developed, which show Plm II and IV inhibitory potency in low nanomolar range and remarkable selectivity against Cathepsin D along with decreased lipophilicity and increased solubility.

Rational Tuning of Visual Cycle Modulator Pharmacodynamics

Modulators of the visual cycle have been developed for treatment of various retinal disorders. These agents were designed to inhibit retinoid isomerase [retinal pigment epithelium-specific 65 kDa protein (RPE65)], the rate-limiting enzyme of the visual cycle, based on the idea that attenuation of visual pigment regeneration could reduce formation of toxic retinal conjugates. Of these agents, certain ones that contain primary amine groups can also reversibly form retinaldehyde Schiff base adducts, which contributes to their retinal protective activity. Direct inhibition of RPE65 as a therapeutic strategy is complicated by adverse effects resulting from slowed chromophore regeneration, whereas effective retinal sequestration can require high drug doses with potential off-target effects. We hypothesized that the RPE65-emixustat crystal structure could help guide the design of retinaldehyde-sequestering agents with varying degrees of RPE65 inhibitory activity. We found that addition of an isopropyl group to the central phenyl ring of emixustat and related compounds resulted in agents effectively lacking in vitro retinoid isomerase inhibitory activity, whereas substitution of the terminal 6-membered ring with branched moieties capable of stronger RPE65 interaction potentiated inhibition. The isopropyl derivative series produced discernible visual cycle suppression in vivo, albeit much less potently than compounds with a high affinity for the RPE65 active site. These agents were distributed into the retina and formed Schiff base adducts with retinaldehyde. Except for one compound [3-amino-1-(3-isopropyl-5-((2,6,6-trimethylcyclohex-1-en-1-yl)methoxy)phenyl)propan-1-ol (MB-007)], these agents conferred protection against retinal phototoxicity, suggesting that both direct RPE65 inhibition and retinal sequestration are mechanisms of potential therapeutic relevance.

IDD388 Polyhalogenated Derivatives as Probes for an Improved Structure-Based Selectivity of AKR1B10 Inhibitors

Human enzyme aldo-keto reductase family member 1B10 (AKR1B10) has evolved as a tumor marker and promising antineoplastic target. It shares high structural similarity with the diabetes target enzyme aldose reductase (AR). Starting from the potent AR inhibitor IDD388, we have synthesized a series of derivatives bearing the same halophenoxyacetic acid moiety with an increasing number of bromine (Br) atoms on its aryl moiety. Next, by means of IC₅₀ measurements, X-ray crystallography, WaterMap analysis, and advanced binding free energy calculations with a quantum-mechanical (QM) approach, we have studied their structure-activity relationship (SAR) against both enzymes. The introduction of Br substituents decreases AR inhibition potency but improves it in the case of AKR1B10. Indeed, the Br atoms in ortho position may impede these drugs to fit into the AR prototypical specificity pocket. For AKR1B10, the smaller aryl moieties of MK181 and IDD388 can bind into the external loop A subpocket. Instead, the bulkier MK184, MK319, and MK204 open an inner specificity pocket in AKR1B10 characterized by a π-π stacking interaction of their aryl moieties and Trp112 side chain in the native conformation (not possible in AR). Among the three compounds, only MK204 can make a strong halogen bond with the protein (-4.4 kcal/mol, using QM calculations), while presenting the lowest desolvation cost among all the series, translated into the most selective and inhibitory potency AKR1B10 (IC₅₀ = 80 nM). Overall, SAR of these IDD388 polyhalogenated derivatives have unveiled several distinctive AKR1B10 features (shape, flexibility, hydration) that can be exploited to design novel types of AKR1B10 selective drugs.

Fragment-Based Discovery of 2-Aminoquinazolin-4(3H)-ones As Novel Class Nonpeptidomimetic Inhibitors of the Plasmepsins I, II, and IV

2-Aminoquinazolin-4(3H)-ones were identified as a novel class of malaria digestive vacuole plasmepsin inhibitors by using NMR-based fragment screening against Plm II. Initial fragment hit optimization led to a submicromolar inhibitor, which was cocrystallized with Plm II to produce an X-ray structure of the complex. The structure showed that 2-aminoquinazolin-4(3H)-ones bind to the open flap conformation of the enzyme and provided clues to target the flap pocket. Further improvement in potency was achieved via introduction of hydrophobic substituents occupying the flap pocket. Most of the 2-aminoquinazolin-4(3H)-one based inhibitors show a similar activity against digestive Plms I, II, and IV and >10-fold selectivity versus CatD, although varying the flap pocket substituent led to one Plm IV selective inhibitor. In cell-based assays, the compounds show growth inhibition of Plasmodium falciparum 3D7 with IC50 ∼ 1 μM. Together, these results suggest 2-aminoquinazolin-4(3H)-ones as perspective leads for future development of an antimalarial agent.

Limitations and extensions of the lock-and-key principle: differences between gas state, solution and solid state structures

The lock-and-key concept is discussed with respect to necessary extensions. Formation of supramolecular complexes depends not only, and often not even primarily on an optimal geometric fit between host and guest. Induced fit and allosteric interactions have long been known as important modifications. Different binding mechanisms, the medium used and pH effects can exert a major influence on the affinity. Stereoelectronic effects due to lone pair orientation can lead to variation of binding constants by orders of magnitude. Hydrophobic interactions due to high-energy water inside cavities modify the mechanical lock-and-key picture. That optimal affinities are observed if the cavity is only partially filled by the ligand can be in conflict with the lock-and-key principle. In crystals other forces than those between host and guest often dominate, leading to differences between solid state and solution structures. This is exemplified in particular with calixarene complexes, which by X-ray analysis more often than other hosts show guest molecules outside their cavity. In view of this the particular problems with the identification of weak interactions in crystals is discussed.

Molecular recognition in chemical and biological systems

Structure-based ligand design in medicinal chemistry and crop protection relies on the identification and quantification of weak noncovalent interactions and understanding the role of water. Small-molecule and protein structural database searches are important tools to retrieve existing knowledge. Thermodynamic profiling, combined with X-ray structural and computational studies, is the key to elucidate the energetics of the replacement of water by ligands. Biological receptor sites vary greatly in shape, conformational dynamics, and polarity, and require different ligand-design strategies, as shown for various case studies. Interactions between dipoles have become a central theme of molecular recognition. Orthogonal interactions, halogen bonding, and amide⋅⋅⋅π stacking provide new tools for innovative lead optimization. The combination of synthetic models and biological complexation studies is required to gather reliable information on weak noncovalent interactions and the role of water.

Hydrocarbon binding by proteins: structures of protein binding sites for ≥C10 linear alkanes or long-chain alkyl and alkenyl groups

In order to identify potential de novo enzyme templates for the cleavage of C–C single bonds in long-chain hydrocarbons, we analyzed protein structures that bind substrates containing alkyl and alkenyl functional groups. A survey of ligand-containing protein structures deposited in the Protein Data Bank resulted in 874 entries, consisting of 194 unique ligands that have ≥10 carbons in a linear chain. Fatty acids and phospholipids are the most abundant types of ligands. Hydrophobic amino acids forming α-helical structures frequently line the binding pockets. Occupation of these binding sites was evaluated by calculating both the buried surface area and volume employed by the ligands; these quantities are similar to those computed for drug–protein complexes. Surface complementarity is relatively low due to the nonspecific nature of the interaction between the long-chain hydrocarbons and the hydrophobic amino acids. The selected PDB structures were annotated on the basis of their SCOP and EC identification numbers, which will facilitate design template searches based on structural and functional homologies. Relatively low surface complementarity and ∼55% volume occupancy, also observed in synthetic-host, alkane-guest systems, suggest general principles for the recognition of long-chain linear hydrocarbons.

Flap dynamics of plasmepsin proteases: insight into proposed parameters and molecular dynamics

Herein, for the first time, we report the flap opening and closing in Plasmepsin proteases – plasmepsin II (PlmII) was used as a prototype model.

Potent inhibitors of malarial aspartic proteases, the plasmepsins, by hydroformylation of substituted 7-azanorbornenes

The increasing prevalence of multidrug-resistant strains of the malarial parasite Plasmodium falciparum requires the urgent development of new therapeutic agents with novel modes of action. The vacuolar malarial aspartic proteases plasmepsin (PM) I, II, and IV are involved in hemoglobin degradation and play a central role in the growth and maturation of the parasite in the human host. We report the structure-based design, synthesis, and in vitro evaluation of a new generation of PM inhibitors featuring a highly decorated 7-azabicyclo[2.2.1]heptane core. While this protonated central core addresses the catalytic Asp dyad, three substituents bind to the flap, the S1/S3, and the S1' pockets of the enzymes. A hydroformylation reaction is the key synthetic step for the introduction of the new vector reaching into the S1' pocket. The configuration of the racemic ligands was confirmed by extensive NMR and X-ray crystallographic analysis. In vitro biological assays revealed high potency of the new inhibitors against the three plasmepsins (IC(50) values down to 6 nM) and good selectivity towards the closely related human cathepsins D and E. The occupancy of the S1' pocket makes an essential contribution to the gain in binding affinity and selectivity, which is particularly large in the case of the PM IV enzyme. Designing non-peptidic ligands for PM II is a valid route to generate compounds that inhibit the entire family of vacuolar plasmepsins.

Structural studies of vacuolar plasmepsins

Plasmepsins (PMs) are pepsin-like aspartic proteases present in different species of parasite Plasmodium. Four Plasmodium spp. (P. vivax, P. ovale, P. malariae, and the most lethal P. falciparum) are mainly responsible for causing human malaria that affects millions worldwide. Due to the complexity and rate of parasite mutation coupled with regional variations, and the emergence of P. falciparum strains which are resistant to antimalarial agents such as chloroquine and sulfadoxine/pyrimethamine, there is constant pressure to find new and lasting chemotherapeutic drug therapies. Since many proteases represent therapeutic targets and PMs have been shown to play an important role in the survival of parasite, these enzymes have recently been identified as promising targets for the development of novel antimalarial drugs. The genome of P. falciparum encodes 10 PMs (PMI, PMII, PMIV-X and histo-aspartic protease (HAP)), 4 of which (PMI, PMII, PMIV and HAP) reside within the food vacuole, are directly involved in degradation of human hemoglobin, and share 50-79% amino acid sequence identity. This review focuses on structural studies of only these four enzymes, including their orthologs in other Plasmodium spp.. Almost all original crystallographic studies were performed with PMII, but more recent work on PMIV, PMI, and HAP resulted in a more complete picture of the structure-function relationship of vacuolar PMs. Many structures of inhibitor complexes of vacuolar plasmepsins, as well as their zymogens, have been reported in the last 15 years. Information gained by such studies will be helpful for the development of better inhibitors that could become a new class of potent antimalarial drugs. This article is part of a Special Issue entitled: Proteolysis 50 years after the discovery of lysosome.