Identification of lectin receptors for conserved SARS-CoV-2 glycosylation sites.

New SARS-CoV-2 variants are continuously emerging with critical implications for therapies or vaccinations. The 22 N-glycan sites of Spike remain highly conserved among SARS-CoV-2 variants, opening an avenue for robust therapeutic intervention. Here we used a comprehensive library of mammalian carbohydrate-binding proteins (lectins) to probe critical sugar residues on the full-length trimeric Spike and the receptor binding domain (RBD) of SARS-CoV-2. Two lectins, Clec4g and CD209c, were identified to strongly bind to Spike. Clec4g and CD209c binding to Spike was dissected and visualized in real time and at single-molecule resolution using atomic force microscopy. 3D modelling showed that both lectins can bind to a glycan within the RBD-ACE2 interface and thus interferes with Spike binding to cell surfaces. Importantly, Clec4g and CD209c significantly reduced SARS-CoV-2 infections. These data report the first extensive map and 3D structural modelling of lectin-Spike interactions and uncovers candidate receptors involved in Spike binding and SARS-CoV-2 infections. The capacity of CLEC4G and mCD209c lectins to block SARS-CoV-2 viral entry holds promise for pan-variant therapeutic interventions.

Molecular cloning, characterisation and molecular modelling of two novel T-synthases from mollusc origin.

The glycoprotein-N-acetylgalactosamine ß1,3-galactosyltransferase, known as T-synthase (EC 2.4.1.122), plays a crucial role in the synthesis of the T-antigen, which is the core 1 O-glycan structure. This enzyme transfers galactose from UDP-Gal to GalNAc-Ser/Thr. The T-antigen has significant functions in animal development, immune response, and recognition processes. Molluscs are a successful group of animals that inhabit various environments, such as freshwater, marine, and terrestrial habitats. They serve important roles in ecosystems as filter feeders and decomposers but can also be pests in agriculture and intermediate hosts for human and cattle parasites. The identification and characterization of novel carbohydrate active enzymes, such as T-synthase, can aid in the understanding of molluscan glycosylation abilities and their adaptation and survival abilities. Here, the T-synthase enzymes from the snail Pomacea canaliculata and the oyster Crassostrea gigas are identified, cloned, expressed, and characterized, with a focus on structural elucidation. The synthesized enzymes display core 1 ß1,3-galactosyltransferase activity using pNP-α-GalNAc as substrate and exhibit similar biochemical parameters as previously characterised T-synthases from other species. While the enzyme from C. gigas shares the same structural parameters with the other enzymes characterised so far, the T-synthase from P. canaliculata lacks the consensus sequence CCSD, which was previously considered indispensable.

Engineering A-type Dye-Decolorizing Peroxidases by Modification of a Conserved Glutamate Residue.

Dye-decolorizing peroxidases (DyPs) are recently identified microbial enzymes that have been used in several Biotechnology applications from wastewater treatment to lignin valorization. However, their properties and mechanism of action still have many open questions. Their heme-containing active site is buried by three conserved flexible loops with a putative role in modulating substrate access and enzyme catalysis. Here, we investigated the role of a conserved glutamate residue in stabilizing interactions in loop 2 of A-type DyPs. First, we did site saturation mutagenesis of this residue, replacing it with all possible amino acids in bacterial DyPs from Bacillus subtilis (BsDyP) and from Kitasatospora aureofaciens (KaDyP1), the latter being characterized here for the first time. We screened the resulting libraries of variants for activity towards ABTS and identified variants with increased catalytic efficiency. The selected variants were purified and characterized for activity and stability. We furthermore used Molecular Dynamics simulations to rationalize the increased catalytic efficiency and found that the main reason is the electron channeling becoming easier from surface-exposed tryptophans. Based on our findings, we also propose that this glutamate could work as a pH switch in the wild-type enzyme, preventing intracellular damage.

Understanding the riddle of amine oxidase flavoenzyme reactivity on the stereoisomers of N-methyl-dopa and N-methyl-tyrosine.

Enzymes are usually stereospecific against chiral substrates, which is commonly accepted for the amine oxidase family of enzymes as well. However, the FsqB (fumisoquin biosynthesis gene B) enzyme that belongs to the family of sarcosine oxidase and oxidizes L-N-methyl-amino acids, shows surprising activity for both enantiomers of N-methyl-dopa. The aim of this study is to understand the mechanism behind this behavior. Primary docking experiments showed that tyrosine and aspartate residues (121 and 315 respectively) are located on the ceiling of the active site of FsqB and may play a role in fixing the N-methyl-dopa via its catechol moiety and allowing both stereoisomers of this substrate to be in close proximity of the N5 atom of the isoalloxazine ring of the cofactor. Three experimental approaches were used to prove this hypothesis which are: (1) studying the oxidative ability of the variants Y121F and D315A on N-methyl-dopa substrates in comparison with N-methyl-tyrosine substrates; (2) studying the FsqB WT and variants catalyzed biotransformation via high-performance liquid chromatography (HPLC); (3) molecular dynamics simulations to characterize the underlying mechanisms of the molecular recognition. First, we found that the chemical characteristics of the catechol moiety of N-methyl-dopa are important to explain the differences between N-methyl-dopa and N-methyl-tyrosine. Furthermore, we found that Y121 and D315 are specific in FsqB and not found in the model enzyme sarcosine oxidase. The on-bench and theoretical mutagenesis studies show that Y121 residue has a major role in fixing the N-methyl-dopa substrates close to the N5 atom of the isoalloxazine ring of the cofactor. Simultaneously, D315 has a supportive role in this mechanism. Jointly, the experimental and theoretical approaches help to solve the riddle of FsqB amine oxidase substrate specificity.

The Effect of Microwaves on Protein Structure: Molecular Dynamics Approach.

The impact of microwave (MW) irradiation on protein folding, potentially inciting misfolding, was investigated by employing molecular dynamics (MD) simulations. Twenty-nine proteins were subjected to MD simulations under equilibrium (300 K) and MW conditions, where the rotational temperature was elevated to 700 K. The utilized replacement model captures the microwave effects of δ- and γ-relaxation processes (frequency range of ∼300 MHz to ∼20 GHz). The results disclosed that MW heating incited a shift toward more compact protein conformations, as indicated by decreased root-mean-square deviations, root-mean-square fluctuations, head-to-tail distances, and radii of gyration. This compaction was attributed to the intensification of intramolecular electrostatic interactions and hydrogen bonds within the protein caused by MW-destabilized hydrogen bonds between the protein and solvent. The solvent-accessible surface area (SASA), particularly that of polar amino-acid residues, shrank under MW conditions, corresponding to a reduced polarity of the water solvent. However, MW irradiation produced no significant alterations in protein secondary structures; hence, MW heating was observed to primarily affect the protein tertiary structures.

Caspase-Based Fusion Protein Technology: Substrate Cleavability Described by Computational Modeling and Simulation.

The Caspase-based fusion protein technology (CASPON) allows for universal cleavage of fusion tags from proteins of interest to reconstitute the native N-terminus. While the CASPON enzyme has been optimized to be promiscuous against a diversity of N-terminal peptides, the cleavage efficacy for larger proteins can be surprisingly low. We develop an efficient means to rationalize and predict the cleavage efficiency based on a structural representation of the intrinsically disordered N-terminal peptides and their putative interactions with the CASPON enzyme. The number of favorably interacting N-terminal conformations shows a very good agreement with the experimentally observed cleavage efficiency, in agreement with a conformational selection model. The method relies on computationally cheap molecular dynamics simulations to efficiently generate a diverse collection of N-terminal conformations, followed by a simple fitting procedure into the CASPON enzyme. It can be readily used to assess the CASPON cleavability a priori.

Amino Acid Residues Controlling Domain Interaction and Interdomain Electron Transfer in Cellobiose Dehydrogenase.

The function of cellobiose dehydrogenase (CDH) in biosensors, biofuel cells, and as a physiological redox partner of lytic polysaccharide monooxygenase (LPMO) is based on its role as an electron donor. Before donating electrons to LPMO or electrodes, an interdomain electron transfer from the catalytic FAD-containing dehydrogenase domain to the electron shuttling cytochrome domain of CDH is required. This study investigates the role of two crucial amino acids located at the dehydrogenase domain on domain interaction and interdomain electron transfer by structure-based engineering. The electron transfer kinetics of wild-type Myriococcum thermophilum CDH and its variants M309A, R698S, and M309A/R698S were analyzed by stopped-flow spectrophotometry and structural effects were studied by small-angle X-ray scattering. The data show that R698 is essential to pull the cytochrome domain close to the dehydrogenase domain and orient the heme propionate group towards the FAD, while M309 is an integral part of the electron transfer pathway - its mutation reducing the interdomain electron transfer 10-fold. Structural models and molecular dynamics simulations pinpoint the action of these two residues on the domain interaction and interdomain electron transfer.

Accelerated Enveloping Distribution Sampling (AEDS) Allows for Efficient Sampling of Orthogonal Degrees of Freedom.

One of the most challenging aspects in the molecular simulation of proteins is the study of slowly relaxing processes, such as loop rearrangements or ligands that adopt different conformations in the binding site. State-of-the-art methods used to calculate binding free energies rely on performing several short simulations (lambda steps), in which the ligand is slowly transformed into the endstates of interest. This makes capturing the slowly relaxing processes even more difficult, as they would need to be observed in all of the lambda steps. One attractive alternative is the use of a reference state capable of sampling all of the endstates of interest in a single simulation. However, the energy barriers between the states can be too high to overcome, thus hindering the sampling of all of the relevant conformations. Accelerated enveloping distribution sampling (AEDS) is a recently developed reference state technique that circumvents the high-energy-barrier challenge by adding a boosting potential that flattens the energy landscape without distorting the energy minima. In the present work, we apply AEDS to the well-studied benchmark system T4 lysozyme L99A. The T4 lysozyme L99A mutant contains a hydrophobic pocket in which there is a valine (valine 111), whose conformation influences the binding efficiencies of the different ligands. Incorrectly sampling the dihedral angle can lead to errors in predicted binding free energies of up to 16 kJ mol-1. This protein constitutes an ideal scenario to showcase the advantages and challenges when using AEDS in the presence of slow relaxing processes. We show that AEDS is able to successfully sample the relevant degrees of freedom, providing accurate binding free energies, without the need of previous information of the system in the form of collective variables. Additionally, we showcase the capabilities of AEDS to efficiently screen a set of ligands. These results represent a promising first step toward the development of free-energy methods that can respond to more intricate molecular events.

Size Matters: Free-Energy Calculations of Amino Acid Adsorption over Pristine Graphene.

There is still growing interest in graphene interactions with proteins, both for its possible biological applications and due to concerns over detrimental effects at the cellular level. As with any process involving proteins, an understanding of amino acid composition is desirable. In this work, we systematically studied the adsorption process of amino acids onto pristine graphene via rigorous free-energy calculations. We characterized the free energy, potential energy, and entropy of the adsorption of all proteinogenic amino acids. The energetic components were further separated into pair interaction contributions. A linear correlation was found between the free energy and the solvent accessible surface area change during adsorption (ΔSASAads) over pristine graphene and uncharged amino acids. Free energies over pristine graphene were compared with adsorption onto graphene oxide, finding an almost complete loss of the favorability of amino acid adsorption onto graphene. Finally, the correlation with ΔSASAads was used to successfully predict the free energy of adsorption of several penta-l-peptides in different structural states and sequences. Due to the relative ease of calculating the ΔSASAads compared to free-energy calculations, it could prove to be a cost-effective predictor of the free energy of adsorption for proteins onto nonpolar surfaces.

Insights into the Binding Mode of Lipid A to the Anti-lipopolysaccharide Factor ALFPm3 from Penaeus monodon: An In Silico Study through MD Simulations.

The globally expanding threat of antibiotic resistance calls for the development of new strategies for abating Gram-negative bacterial infections. The use of extracorporeal blood cleansing devices with affinity sorbents to selectively capture bacterial lipopolysaccharide (LPS), which is the major constituent of Gram-negative bacterial outer membranes and the responsible agent for eliciting an exacerbated innate immune response in the host during infection, has received outstanding interest. For that purpose, molecules that bind tightly to LPS are required to functionalize the affinity sorbents. Particularly, anti-LPS factors (ALFs) are promising LPS-sequestrating molecules. Hence, in this work, molecular dynamics (MD) simulations are used to investigate the interaction mechanism and binding pose of the ALF isoform 3 from Penaeus monodon (ALFPm3), which is referred to as "AL3" for the sake of simplicity, and lipid A (LA, the component of LPS that represents its endotoxic principle). We concluded that hydrophobic interactions are responsible for AL3-LA binding and that LA binds to AL3 within the protein cavity, where it buries its aliphatic tails, whereas the negatively charged phosphate groups are exposed to the medium. AL3 residues that are key for its interaction with LA were identified, and their conservation in other ALFs (specifically Lys39 and Tyr49) was also analyzed. Additionally, based on the MD-derived results, we provide a picture of the possible AL3-LA interaction mechanism. Finally, an in vitro validation of the in silico predictions was performed. Overall, the insights gained from this work can guide the design of novel therapeutics for treating sepsis, since they may be significantly valuable for designing LPS-sequestrating molecules that could functionalize affinity sorbents to be used for extracorporeal blood detoxification.

Identification of Activating Mutations in the Transmembrane and Extracellular Domains of EGFR.

The epidermal growth factor receptor (EGFR) is frequently mutated in human cancer, most notably non-small-cell lung cancer and glioblastoma. While many frequently occurring EGFR mutations are known to confer constitutive EGFR activation, the situation is less clear for rarely detected variants. In fact, more than 1000 distinct EGFR mutations are listed in the Catalogue of Somatic Mutations in Cancer (COSMIC), but for most of them, the functional consequence is unknown. To identify additional, previously unknown activating mutations in EGFR, we screened a randomly mutated EGFR library for constitutive EGFR phosphorylation using a recently developed high-throughput approach termed PhosphoFlowSeq. Enrichment of the well-known activating mutations S768I, T790M, and L858R validated the experimental approach. Importantly, we also identified the activating mutations S442I and L658Q located in the extracellular and transmembrane domains of EGFR, respectively. To the best of our knowledge, neither S442I nor L658Q has been associated with an activating phenotype before. However, both have been detected in cancer samples. Interestingly, molecular dynamics (MD) simulations suggest that the L658Q mutation located in the hydrophobic transmembrane region forms intermolecular hydrogen bonds, thereby promoting EGFR dimerization and activation. Based on these findings, we screened the COSMIC database for additional hydrophilic mutations in the EGFR transmembrane region and indeed detected moderate constitutive activation of EGFR-G652R. Together, this study demonstrates that unbiased screening for activating mutations in EGFR not only yields well-established substitutions located in the kinase domain but also activating mutations in other regions of EGFR, including the extracellular and transmembrane domains.

PROFICS: A bacterial selection system for directed evolution of proteases.

Proteases serve as important tools in biotechnology and as valuable drugs or drug targets. Efficient protein engineering methods to study and modulate protease properties are thus of great interest for a plethora of applications. We established PROFICS (PRotease Optimization via Fusion-Inhibited Carbamoyltransferase-based Selection), a bacterial selection system, which enables the optimization of proteases for biotechnology, therapeutics or diagnosis in a simple overnight process. During the PROFICS process, proteases are selected for their ability to specifically cut a tag from a reporter enzyme and leave a native N-terminus. Precise and efficient cleavage after the recognition sequence reverses the phenotype of an Escherichia coli knockout strain deficient in an essential enzyme of pyrimidine synthesis. A toolbox was generated to select for proteases with different preferences for P1' residues (the residue immediately following the cleavage site). The functionality of PROFICS is demonstrated with viral proteases and human caspase-2. PROFICS improved caspase-2 activity up to 25-fold after only one round of mutation and selection. Additionally, we found a significantly improved tolerance for all P1' residues caused by a mutation in a substrate interaction site. We showed that this improved activity enables cells containing the new variant to outgrow cells containing all other mutants, facilitating its straightforward selection. Apart from optimizing enzymatic activity and P1' tolerance, PROFICS can be used to reprogram specificities, erase off-target activity, optimize expression via tags/codon usage, or even to screen for potential drug-resistance-conferring mutations in therapeutic targets such as viral proteases in an unbiased manner.

On the effects of induced polarizability at the water-graphene interface via classical charge-on-spring models.

Molecular models of the water-graphene interaction are essential to describe graphene in condensed phases. Different challenges are associated with the generation of these models, in particular π-π and dispersion interactions; thus quantum and classical models have been developed and due to the numerical efficiency of the latter, they have been extensively employed. In this work, we have systematically studied, via molecular dynamics, two polarizable graphene models, denominated CCCP and CCCPD, employing the charge-on-spring model of the GROMOS forcefield, both being compatible with the polarizable water models COS/G2 and COS/D2, respectively. These models were compared with non-polarizable graphene and SPC water models. We focused the study on the water-graphene interface in two distinct systems and under the influence of an electric field: one composed of graphene immersed in water and the other composed of graphene with a water droplet above it. In the former, the orientation of water close to the graphene layer is affected by polarizable graphene in comparison to non-polarizable graphene. This effect is emphasised when an electric field is applied. In the latter, carbon polarizability reduced water contact angles, but graphene retained its hydrophobicity and the computed angles are within the experimental data. Given the significant extra computational cost, the use of polarizable models instead of the traditional fixed-charged approach for the graphene-water interaction may be justified when polarizability effects are relevant, for example, when applying relatively strong fields or in very anisotropic systems, such as the vacuum-bulk interface, as these models are more responsive to such conditions.

Reaction intermediate rotation during the decarboxylation of coproheme to heme b in C. diphtheriae.

Monoderm bacteria utilize coproheme decarboxylases (ChdCs) to generate heme b by a stepwise decarboxylation of two propionate groups of iron coproporphyrin III (coproheme), forming two vinyl groups. This work focuses on actinobacterial ChdC from Corynebacterium diphtheriae (CdChdC) to elucidate the hydrogen peroxide-mediated decarboxylation of coproheme via monovinyl monopropionyl deuteroheme (MMD) to heme b, with the principal aim being to understand the reorientation mechanism of MMD during turnover. Wild-type CdChdC and variants, namely H118A, H118F, and A207E, were studied by resonance Raman and ultraviolet-visible spectroscopy, mass spectrometry, and molecular dynamics simulations. As actinobacterial ChdCs use a histidine (H118) as a distal base, we studied the H118A and H118F variants to elucidate the effect of 1) the elimination of the proton acceptor and 2) steric constraints within the active site. The A207E variant mimics the proximal H-bonding network found in chlorite dismutases. This mutation potentially increases the rigidity of the proximal site and might impair the rotation of the reaction intermediate MMD. We found that both wild-type CdChdC and the variant H118A convert coproheme mainly to heme b upon titration with H2O2. Interestingly, the variant A207E mostly accumulates MMD along with small amounts of heme b, whereas H118F is unable to produce heme b and accumulates only MMD. Together with molecular dynamics simulations, the spectroscopic data provide insight into the reaction mechanism and the mode of reorientation of MMD, i.e., a rotation in the active site versus a release and rebinding.

Exploring the structure and dynamics of proteins in soil organic matter.

Alongside inorganic materials, water, and air, soil organic matter (SOM) is one of the major components of soil and has tremendous influence on the environment given its vital role in the carbon cycle. Many soil dwelling organisms like plants, fungi and bacteria excrete proteins, whose interaction with SOM is poorly understood on an atomistic level. In this study, molecular dynamics simulations were used to investigate selected proteins in soil models of different complexity from simple co-solvent molecules to Leonardite humic acids (LHA). We analyzed the proteins in terms of their structural stability, the nature and strength of the interactions with their surroundings, as well as their aggregation behavior. Upon insertion of proteins in complex SOM models, their structural stability decreased, although no unfolding or disruption of secondary structure was observed. The interactions of proteins and SOM were primarily governed by electrostatic forces, often in form of hydrogen bonds. However, also weaker van der Waals forces made a significant contribution to the total interaction energies. Moreover, we showed that even though the molecular structure and size of SOM molecules varied, the functional groups of SOM ordered around the protein in a similar pattern. Finally, the number of aggregates formed by proteins and SOM molecules was shown to be primarily proportional to the size of the latter. Strikingly, for varying protein net charges no changes in the formation of aggregates with the strongly negatively charged LHA were observed.

On the use of multiple-time-step algorithms to save computing effort in molecular dynamics simulations of proteins.

Computer simulation of proteins in aqueous solution at the atomic level of resolution is still limited in time span and system size due to limited computing power available and thus employs a variety of time-saving techniques that trade some accuracy against computational effort. Examples of such time-saving techniques are the application of constraints to particular degrees of freedom or the use of a multiple-time-step (MTS) algorithm distinguishing between particular forces when integrating Newton's equations of motion. The application of two types of MTS algorithms to bond-stretching forces versus the remaining forces in molecular dynamics (MD) simulations of a protein in aqueous solution or of liquid water is investigated and the results in terms of total energy conservation and the influence on various other properties are compared to those of MD simulations of the same systems using bond-length, and for water bond-angle, constraints. At comparable computational effort, the use of bond-length constraints in proteins leads to better energy conservation and less distorted properties than the two MTS algorithms investigated.

On the Effect of the Various Assumptions and Approximations used in Molecular Simulations on the Properties of Bio-Molecular Systems: Overview and Perspective on Issues.

Computer simulations of molecular systems enable structure-energy-function relationships of molecular processes to be described at the sub-atomic, atomic, supra-atomic or supra-molecular level and plays an increasingly important role in chemistry, biology and physics. To interpret the results of such simulations appropriately, the degree of uncertainty and potential errors affecting the calculated properties must be considered. Uncertainty and errors arise from (1) assumptions underlying the molecular model, force field and simulation algorithms, (2) approximations implicit in the interatomic interaction function (force field), or when integrating the equations of motion, (3) the chosen values of the parameters that determine the accuracy of the approximations used, and (4) the nature of the system and the property of interest. In this overview, advantages and shortcomings of assumptions and approximations commonly used when simulating bio-molecular systems are considered. What the developers of bio-molecular force fields and simulation software can do to facilitate and broaden research involving bio-molecular simulations is also discussed.

Fighting Against Bacterial Lipopolysaccharide-Caused Infections through Molecular Dynamics Simulations: A Review.

Lipopolysaccharide (LPS) is the primary component of the outer leaflet of Gram-negative bacterial outer membranes. LPS elicits an overwhelming immune response during infection, which can lead to life-threatening sepsis or septic shock for which no suitable treatment is available so far. As a result of the worldwide expanding multidrug-resistant bacteria, the occurrence and frequency of sepsis are expected to increase; thus, there is an urge to develop novel strategies for treating bacterial infections. In this regard, gaining an in-depth understanding about the ability of LPS to both stimulate the host immune system and interact with several molecules is crucial for fighting against LPS-caused infections and allowing for the rational design of novel antisepsis drugs, vaccines and LPS sequestration and detection methods. Molecular dynamics (MD) simulations, which are understood as being a computational microscope, have proven to be of significant value to understand LPS-related phenomena, driving and optimizing experimental research studies. In this work, a comprehensive review on the methods that can be combined with MD simulations, recently applied in LPS research, is provided. We focus especially on both enhanced sampling methods, which enable the exploration of more complex systems and access to larger time scales, and free energy calculation approaches. Thereby, apart from outlining several strategies for surmounting LPS-caused infections, this work reports the current state-of-the-art of the methods applied with MD simulations for moving a step forward in the development of such strategies.

Efficient In Silico Saturation Mutagenesis of a Member of the Caspase Protease Family.

Rational-design methods have proven to be a valuable toolkit in the field of protein design. Numerical approaches such as free-energy calculations or QM/MM methods are fit to widen the understanding of a protein-sequence space but require large amounts of computational time and power. Here, we apply an efficient method for free-energy calculations that combines the one-step perturbation (OSP) with the third-power-fitting (TPF) approach. It is fit to calculate full free energies of binding from three different end states only. The nonpolar contribution to the free energies are calculated for a set of chosen amino acids from a single simulation of a judiciously chosen reference state. The electrostatic contributions, on the other hand, are predicted from simulations of the neutral and charged end states of the individual amino acids. We used this method to perform in silico saturation mutagenesis of two sites in human Caspase-2. We calculated relative binding free energies toward two different substrates that differ in their P1' site and in their affinity toward the unmutated protease. Although being approximate, our approach showed very good agreement upon validation against experimental data. 76% of the predicted relative free energies of amino acid mutations was found to be true positives or true negatives. We observed that this method is fit to discriminate amino acid mutations because the rate of false negatives is very low (<1.5%). The approach works better for a substrate with medium/low affinity with a Matthews correlation coefficient (MCC) of 0.63, whereas for a substrate with very low affinity, the MCC was 0.38. In all cases, the combined TPF + OSP approach outperformed the linear interaction energy method.

In silico identification of noncompetitive inhibitors targeting an uncharacterized allosteric site of falcipain-2.

Falcipain-2 (FP-2) is a Plasmodium falciparum hemoglobinase widely targeted in the search for antimalarials. FP-2 can be allosterically modulated by various noncompetitive inhibitors that have been serendipitously identified. Moreover, the crystal structures of two inhibitors bound to an allosteric site, termed site 6, of the homolog enzyme human cathepsin K (hCatK) suggest that the equivalent region in FP-2 might play a similar role. Here, we conduct the rational identification of FP-2 inhibitors through virtual screenings (VS) of compounds into several pocket-like conformations of site 6, sampled during molecular dynamics (MD) simulations of the free enzyme. Two noncompetitive inhibitors, ZINC03225317 and ZINC72290660, were confirmed using in vitro enzymatic assays and their poses into site 6 led to calculated binding free energies matching the experimental ones. Our results provide strong evidence about the allosteric inhibition of FP-2 through binding of small molecules to site 6, thus opening the way toward the discovery of new inhibitors against this enzyme.