Three decades targeting falcipains to develop antiplasmodial agents: what have we learned and what can be done next?

Malaria is a devastating infectious disease that affects large swathes of human populations across the planet's tropical regions. It is caused by parasites of the genus Plasmodium, with Plasmodium falciparum being responsible for the most lethal form of the disease. During the intraerythrocytic stage in the human hosts, malaria parasites multiply and degrade hemoglobin (Hb) using a battery of proteases, which include two cysteine proteases, falcipains 2 and 3 (FP-2 and FP-3). Due to their role as major hemoglobinases, FP-2 and FP-3 have been targeted in studies aiming to discover new antimalarials and numerous inhibitors with activity against these enzymes, and parasites in culture have been identified. Nonetheless, cross-inhibition of human cysteine cathepsins remains a serious hurdle to overcome for these compounds to be used clinically. In this article, we have reviewed key functional and structural properties of FP-2/3 and described different compound series reported as inhibitors of these proteases during decades of active research in the field. Special attention is also paid to the wide range of computer-aided drug design (CADD) techniques successfully applied to discover new active compounds. Finally, we provide guidelines that, in our understanding, will help advance the rational discovery of new FP-2/3 inhibitors.

Genetic diversity of the melanocortin-1 receptor in an admixed population of Rio de Janeiro: Structural and functional impacts of Cys35Tyr variant.

The melanocortin-1 receptor (MC1R) is one of the key proteins involved in the regulation of melanin production and several polymorphisms have been associated with different phenotypes of skin and hair color in human and nonhuman species. Most of the knowledge is centered on more homogeneous populations and studies involving an admixed group of people should be encouraged due to the great importance of understanding the human color variation. This work evaluates the MC1R diversity and the possible impacts of MC1R variants in an admixed sample population of Rio de Janeiro, Brazil, which is a product of Native American, African, and European miscegenation. Sequencing of complete coding region and part of the 3´UTR of MC1R gene identified 31 variants including one insertion and three novel synonymous substitutions in sample population grouped according to skin, hair and eye pigmentation levels. In nonmetric multidimensional scaling analysis (NMDS), three main clusters were identified, in which the Brazilian dark skin group remained in the African cluster whereas the intermediate and the light skin color phenotype in the European one. None gathered with Asians since their immigration to Brazil was a recent event. In silico analyses demonstrated that Cys35Tyr, Ile155Thr and Pro256Ser, found in our population, have a negative effect on receptor function probably due to changes on the receptor structure. Notably, Cys35Tyr mutation could potentially impair agonist binding. Altogether, this work contributes to the understanding of the genetic background of color variation on an admixed population and gives insights into the damaging effects of MC1R variants.

A novel receptor for platelet-activating factor and lysophosphatidylcholine in Trypanosoma cruzi.

The lipid mediators, platelet-activating factor (PAF) and lysophosphatidylcholine (LPC), play relevant pathophysiological roles in Trypanosoma cruzi infection. Several species of LPC, including C18:1 LPC, which mimics the effects of PAF, are synthesized by T. cruzi. The present study identified a receptor in T. cruzi, which was predicted to bind to PAF, and found it to be homologous to members of the progestin and adiponectin family of receptors (PAQRs). We constructed a three-dimensional model of the T. cruzi PAQR (TcPAQR) and performed molecular docking to predict the interactions of the TcPAQR model with C16:0 PAF and C18:1 LPC. We knocked out T. cruzi PAQR (TcPAQR) gene and confirmed the identity of the expressed protein through immunoblotting and immunofluorescence assays using an anti-human PAQR antibody. Wild-type and knockout (KO) parasites were also used to investigate the in vitro cell differentiation and interactions with peritoneal mouse macrophages; TcPAQR KO parasites were unable to react to C16:0 PAF or C18:1 LPC. Our data are highly suggestive that PAF and LPC act through TcPAQR in T. cruzi, triggering its cellular differentiation and ability to infect macrophages.

Improving Blind Docking in DOCK6 through an Automated Preliminary Fragment Probing Strategy.

Probing protein surfaces to accurately predict the binding site and conformation of a small molecule is a challenge currently addressed through mainly two different approaches: blind docking and cavity detection-guided docking. Although cavity detection-guided blind docking has yielded high success rates, it is less practical when a large number of molecules must be screened against many detected binding sites. On the other hand, blind docking allows for simultaneous search of the whole protein surface, which however entails the loss of accuracy and speed. To bridge this gap, in this study, we developed and tested BLinDPyPr, an automated pipeline which uses FTMap and DOCK6 to perform a hybrid blind docking strategy. Through our algorithm, FTMap docked probe clusters are converted into DOCK6 spheres for determining binding regions. Because these spheres are solely derived from FTMap probes, their locations are contained in and specific to multiple potential binding pockets, which become the regions that are simultaneously probed and chosen by the search algorithm based on the properties of each candidate ligand. This method yields pose prediction results (45.2-54.3% success rates) comparable to those of site-specific docking with the classic DOCK6 workflow (49.7-54.3%) and is half as time-consuming as the conventional blind docking method with DOCK6.

Enzymatic and structural properties of human glutamine:fructose-6-phosphate amidotransferase 2 (hGFAT2).

Glycoconjugates play a central role in several cellular processes, and alteration in their composition is associated with numerous human pathologies. Substrates for cellular glycosylation are synthesized in the hexosamine biosynthetic pathway, which is controlled by the glutamine:fructose-6-phosphate amidotransfera-se (GFAT). Human isoform 2 GFAT (hGFAT2) has been implicated in diabetes and cancer; however, there is no information about structural and enzymatic properties of this enzyme. Here, we report a successful expression and purification of a catalytically active recombinant hGFAT2 (rhGFAT2) in Escherichia coli cells fused or not to a HisTag at the C-terminal end. Our enzyme kinetics data suggest that hGFAT2 does not follow the expected ordered bi-bi mechanism, and performs the glucosamine-6-phosphate synthesis much more slowly than previously reported for other GFATs. In addition, hGFAT2 is able to isomerize fructose-6-phosphate into glucose-6-phosphate even in the presence of equimolar amounts of glutamine, which results in unproductive glutamine hydrolysis. Structural analysis of a three-dimensional model of rhGFAT2, corroborated by circular dichroism data, indicated the presence of a partially structured loop in the glutaminase domain, whose sequence is present in eukaryotic enzymes but absent in the E. coli homolog. Molecular dynamics simulations suggest that this loop is the most flexible portion of the protein and plays a key role on conformational states of hGFAT2. Thus, our study provides the first comprehensive set of data on the structure, kinetics, and mechanics of hGFAT2, which will certainly contribute to further studies on the (patho)physiology of hGFAT2.

Theoretical studies and NMR assay of coumarins and neoflavanones derivatives as potential inhibitors of acetylcholinesterase.

Currently Alzheimer's disease (AD) is a devastating neurological disorder that mainly affects the elderly. The treatment of AD has as main objective to increase the levels of ACh in the synaptic cleft by inhibiting the cholinesterase enzymes, which are responsible for the degradation of ACh. Twenty one synthesized coumarins and neoflavanones (4-arylcoumarins) and theoretical studies were used to select the most promising ligands for in vitro experimental studies by Nuclear Magnetic Resonance. The eight compounds selected for the experimental study only 12b (effectiveness 68.54 ±â€¯3.22%) was promising AChE inhibitor. This compound (12b) presents substituents at positions 4, 5, 6, 7 and 8 in a coumarin nucleus, being the most significant characteristic in comparison to the other studied compounds. These results can be used for the design and synthesis of other possible derivatives with inhibitory potential of AChE.

Identification of (4-(9H-fluoren-9-yl) piperazin-1-yl) methanone derivatives as falcipain 2 inhibitors active against Plasmodium falciparum cultures.

BACKGROUND: Falcipain 2 (FP-2) is the hemoglobin-degrading cysteine protease of Plasmodium falciparum most extensively targeted to develop novel antimalarials. However, no commercial antimalarial drugs based on FP-2 inhibition are available yet due to the low selectivity of most FP-2 inhibitors against the human cysteine proteases. METHODS: A structure-based virtual screening (SVBS) using Maybridge HitFinder™ compound database was conducted to identify potential FP-2 inhibitors. In vitro enzymatic and cell-growth inhibition assays were performed for the top-scoring compounds. Docking, molecular dynamics (MD) simulations and free energy calculations were employed to study the interaction of the best hits with FP-2 and other related enzymes. RESULTS AND CONCLUSIONS: Two hits based on 4-(9H-fluoren-9-yl) piperazin-1-yl) methanone scaffold, HTS07940 and HTS08262, were identified as inhibitors of FP-2 (half-maximal inhibitory concentration (IC50) = 64 µM and 14.7 µM, respectively) without a detectable inhibition against the human off-target cathepsin K (hCatK). HTS07940 and HTS08262 inhibited the growth of the multidrug-resistant P. falciparum strain FCR3 in culture (half-maximal inhibitory concentrations (IC50) = 2.91 µM and 34 µM, respectively) and exhibited only moderate cytotoxicity against HeLa cells (Half-maximal cytotoxic concentration (CC50) = 133 µM and 350 µM, respectively). Free energy calculations reproduced the experimental affinities of the hits for FP-2 and explained the selectivity with respect to hCatK. GENERAL SIGNIFICANCE: To the best of our knowledge, HTS07940 stands among the most selective FP-2 inhibitors identified by SBVS reported so far, displaying moderate antiplasmodial activity and low cytotoxicity against human cells. Hence, this compound constitutes a promising lead for the design of more potent and selective FP-2 inhibitors.

Dataset showing the impact of the protonation states on molecular dynamics of HIV protease.

The data described here supports the research article "Unraveling HIV Protease Flaps Dynamics by Constant pH Molecular Dynamics Simulations" (Soares et al., 2016) [1]. The data involves both standard Molecular Dynamics (MD) and Constant pH Molecular Dynamics (CpHMD) to elucidate the effect of protonation states of catalytic dyad on the HIV-PR conformation. The data obtained from MD simulation demonstrate that the protonation state of the two aspartic acids (Asp25/Asp25') has a strong influence on the dynamics of the HIV-PR. Regarding the CpHMD simulation, we performed pka calculations for HIV-PR and the data indicate that only one catalytic aspartate should be protonated.

Unraveling HIV protease flaps dynamics by Constant pH Molecular Dynamics simulations.

The active site of HIV protease (HIV-PR) is covered by two flaps. These flaps are known to be essential for the catalytic activity of the HIV-PR, but their exact conformations at the different stages of the enzymatic pathway remain subject to debate. Understanding the correct functional dynamics of the flaps might aid the development of new HIV-PR inhibitors. It is known that, the HIV-PR catalytic efficiency is pH-dependent, likely due to the influence of processes such as charge transfer and protonation/deprotonation of ionizable residues. Several Molecular Dynamics (MD) simulations have reported information about the HIV-PR flaps. However, in MD simulations the protonation of a residue is fixed and thus it is not possible to study the correlation between conformation and protonation state. To address this shortcoming, this work attempts to capture, through Constant pH Molecular Dynamics (CpHMD), the conformations of the apo, substrate-bound and inhibitor-bound HIV-PR, which differ drastically in their flap arrangements. The results show that the HIV-PR flaps conformations are defined by the protonation of the catalytic residues Asp25/Asp25' and that these residues are sensitive to pH changes. This study suggests that the catalytic aspartates can modulate the opening of the active site and substrate binding.

Allosteric regulation of the Plasmodium falciparum cysteine protease falcipain-2 by heme.

During the erythrocytic cycle of Plasmodium falciparum malaria parasites break down host hemoglobin, resulting in the release of free heme (ferriprotoporphyrin IX). Heme is a generator of free radicals that cause oxidative stress, but it is detoxified by crystallization into hemozoin inside the food vacuole. We evaluated the interaction of heme and heme analogues with falcipain-2, a P. falciparum food vacuole cysteine protease that plays a key role in hemoglobin digestion. Heme bound to falcipain-2 with a 1:1 stoichiometry, and heme inhibited falcipain-2 activity against both human hemoglobin and chromogenic peptide substrates through a noncompetitive-like mechanism. A series of porphyrin analogues was screened for inhibition of falcipain-2, demonstrating a minor contribution of iron to heme-falcipain-2 interaction, and revealing dependence on both propionic and vinyl groups for inhibition of falcipain-2 by heme. Docking and molecular dynamics simulation unveiled a novel, inducible heme-binding moiety in falcipain-2 adjacent to the catalytic site. Kinetic data suggested that the noncompetitive-like inhibition was substrate inhibition induced by heme. Collectively these data suggest that binding of heme to falcipain-2 may limit the accumulation of free heme in the parasite food vacuole, providing a means of heme detoxification in addition to crystallization into hemozoin.

Structural and functional analysis of a platelet-activating lysophosphatidylcholine of Trypanosoma cruzi.

BACKGROUND: Trypanosoma cruzi is the causative agent of the life-threatening Chagas disease, in which increased platelet aggregation related to myocarditis is observed. Platelet-activating factor (PAF) is a potent intercellular lipid mediator and second messenger that exerts its activity through a PAF-specific receptor (PAFR). Previous data from our group suggested that T. cruzi synthesizes a phospholipid with PAF-like activity. The structure of T. cruzi PAF-like molecule, however, remains elusive. METHODOLOGY/PRINCIPAL FINDINGS: Here, we have purified and structurally characterized the putative T. cruzi PAF-like molecule by electrospray ionization-tandem mass spectrometry (ESI-MS/MS). Our ESI-MS/MS data demonstrated that the T. cruzi PAF-like molecule is actually a lysophosphatidylcholine (LPC), namely sn-1 C18:1(delta 9)-LPC. Similar to PAF, the platelet-aggregating activity of C18:1-LPC was abrogated by the PAFR antagonist, WEB 2086. Other major LPC species, i.e., C16:0-, C18:0-, and C18:2-LPC, were also characterized in all T. cruzi stages. These LPC species, however, failed to induce platelet aggregation. Quantification of T. cruzi LPC species by ESI-MS revealed that intracellular amastigote and trypomastigote forms have much higher levels of C18:1-LPC than epimastigote and metacyclic trypomastigote forms. C18:1-LPC was also found to be secreted by the parasite in extracellular vesicles (EV) and an EV-free fraction. A three-dimensional model of PAFR was constructed and a molecular docking study was performed to predict the interactions between the PAFR model and PAF, and each LPC species. Molecular docking data suggested that, contrary to other LPC species analyzed, C18:1-LPC is predicted to interact with the PAFR model in a fashion similar to PAF. CONCLUSIONS/SIGNIFICANCE: Taken together, our data indicate that T. cruzi synthesizes a bioactive C18:1-LPC, which aggregates platelets via PAFR. We propose that C18:1-LPC might be an important lipid mediator in the progression of Chagas disease and its biosynthesis could eventually be exploited as a potential target for new therapeutic interventions.

Conformational changes in human Hsp70 induced by high hydrostatic pressure produce oligomers with ATPase activity but without chaperone activity.

We investigated the folding of the 70 kDa human cytosolic inducible protein (Hsp70) in vitro using high hydrostatic pressure as a denaturing agent. We followed the structural changes in Hsp70 induced by high hydrostatic pressure using tryptophan fluorescence, molecular dynamics, circular dichroism, high-performance liquid chromatography gel filtration, dynamic light scattering, ATPase activity, and chaperone activity. Although monomeric, Hsp70 is very sensitive to hydrostatic pressure; after pressure had been removed, the protein did not return to its native sate but instead formed oligomeric species that lost chaperone activity but retained ATPase activity.

Anti-prion activity of a panel of aromatic chemical compounds: in vitro and in silico approaches.

The prion protein (PrP) is implicated in the Transmissible Spongiform Encephalopathies (TSEs), which comprise a group of fatal neurodegenerative diseases affecting humans and other mammals. Conversion of cellular PrP (PrP(C)) into the scrapie form (PrP(Sc)) is the hallmark of TSEs. Once formed, PrP(Sc) aggregates and catalyzes PrP(C) misfolding into new PrP(Sc) molecules. Although many compounds have been shown to inhibit the conversion process, so far there is no effective therapy for TSEs. Besides, most of the previously evaluated compounds failed in vivo due to poor pharmacokinetic profiles. In this work we propose a combined in vitro/in silico approach to screen for active anti-prion compounds presenting acceptable drugability and pharmacokinetic parameters. A diverse panel of aromatic compounds was screened in neuroblastoma cells persistently infected with PrP(Sc) (ScN2a) for their ability to inhibit PK-resistant PrP (PrP(Res)) accumulation. From ∼200 compounds, 47 were effective in decreasing the accumulation of PrP(Res) in ScN2a cells. Pharmacokinetic and physicochemical properties were predicted in silico, allowing us to obtain estimates of relative blood brain barrier permeation and mutagenicity. MTT reduction assays showed that most of the active compounds were non cytotoxic. Compounds that cleared PrP(Res) from ScN2a cells, were non-toxic in the MTT assay, and presented a good pharmacokinetic profile were investigated for their ability to inhibit aggregation of an amyloidogenic PrP peptide fragment (PrP(109-149)). Molecular docking results provided structural models and binding affinities for the interaction between PrP and the most promising compounds. In summary, using this combined in vitro/in silico approach we have identified new small organic anti-scrapie compounds that decrease the accumulation of PrP(Res) in ScN2a cells, inhibit the aggregation of a PrP peptide, and possess pharmacokinetic characteristics that support their drugability. These compounds are attractive candidates for prion disease therapy.

Impact of M36I polymorphism on the interaction of HIV-1 protease with its substrates: insights from molecular dynamics.

BACKGROUND: Over the last decades, a vast structural knowledge has been gathered on the HIV-1 protease (PR). Noticeably, most of the studies focused the B-subtype, which has the highest prevalence in developed countries. Accordingly, currently available anti-HIV drugs target this subtype, with considerable benefits for the corresponding patients. RESULTS: Herein, we used molecular dynamics simulations to investigate the role of this polymorphism on the interaction of PR with six of its natural cleavage-sites substrates. CONCLUSIONS: With multiple approaches and analyses we identified structural and dynamical determinants associated with the changes found in the binding affinity of the M36I variant. This mutation influences the flexibility of both PR and its complexed substrate. The observed impact of M36I, suggest that combination with other non-B subtype polymorphisms, could lead to major effects on the interaction with the 12 known cleavage sites, which should impact the virion maturation.

The structural dynamics of the flavivirus fusion peptide-membrane interaction.

Membrane fusion is a crucial step in flavivirus infections and a potential target for antiviral strategies. Lipids and proteins play cooperative roles in the fusion process, which is triggered by the acidic pH inside the endosome. This acidic environment induces many changes in glycoprotein conformation and allows the action of a highly conserved hydrophobic sequence, the fusion peptide (FP). Despite the large volume of information available on the virus-triggered fusion process, little is known regarding the mechanisms behind flavivirus-cell membrane fusion. Here, we evaluated the contribution of a natural single amino acid difference on two flavivirus FPs, FLA(G) ((98)DRGWGNGCGLFGK(110)) and FLA(H) ((98)DRGWGNHCGLFGK(110)), and investigated the role of the charge of the target membrane on the fusion process. We used an in silico approach to simulate the interaction of the FPs with a lipid bilayer in a complementary way and used spectroscopic approaches to collect conformation information. We found that both peptides interact with neutral and anionic micelles, and molecular dynamics (MD) simulations showed the interaction of the FPs with the lipid bilayer. The participation of the indole ring of Trp appeared to be important for the anchoring of both peptides in the membrane model, as indicated by MD simulations and spectroscopic analyses. Mild differences between FLA(G) and FLA(H) were observed according to the pH and the charge of the target membrane model. The MD simulations of the membrane showed that both peptides adopted a bend structure, and an interaction between the aromatic residues was strongly suggested, which was also observed by circular dichroism in the presence of micelles. As the FPs of viral fusion proteins play a key role in the mechanism of viral fusion, understanding the interactions between peptides and membranes is crucial for medical science and biology and may contribute to the design of new antiviral drugs.

GSAFold: a new application of GSA to protein structure prediction.

The folding process defines three-dimensional protein structures from their amino acid chains. A protein's structure determines its activity and properties; thus knowing such conformation on an atomic level is essential for both basic and applied studies of protein function and dynamics. However, the acquisition of such structures by experimental methods is slow and expensive, and current computational methods mostly depend on previously known structures to determine new ones. Here we present a new software called GSAFold that applies the generalized simulated annealing (GSA) algorithm on ab initio protein structure prediction. The GSA is a stochastic search algorithm employed in energy minimization and used in global optimization problems, especially those that depend on long-range interactions, such as gravity models and conformation optimization of small molecules. This new implementation applies, for the first time in ab initio protein structure prediction, an analytical inverse for the Visitation function of GSA. It also employs the broadly used NAMD Molecular Dynamics package to carry out energy calculations, allowing the user to select different force fields and parameterizations. Moreover, the software also allows the execution of several simulations simultaneously. Applications that depend on protein structures include rational drug design and structure-based protein function prediction. Applying GSAFold in a test peptide, it was possible to predict the structure of mastoparan-X to a root mean square deviation of 3.00 Å.

Hybrid QM/MM Molecular Dynamics Study of Benzocaine in a Membrane Environment: How Does a Quantum Mechanical Treatment of Both Anesthetic and Lipids Affect Their Interaction.

Biomolecular dynamics studies using a QM/MM approach have been largely used especially to study enzymatic reactions. However, to the best of our knowledge, the very same approach has not been used to study the water/membrane interface using a quantum mechanical treatment for the lipids. Since a plethora of biochemical processes take place in this region, we believe that it is of primary importance to understand, at the level of molecular orbitals, the behavior of a drug in such an odd environment. In this work, we take advantage of an integration of the CPMD and the GROMACS code, using the Car-Parrinello method, to treat the benzocaine local anesthetic as well as two of the membrane lipids and the GROMOS force field to treat the remaining lipids and the water molecules.

Conformational variability of organophosphorus hydrolase upon soman and paraoxon binding.

The bacterial enzyme organophosphorus hydrolase (OPH) exhibits both catalytic and substrate promiscuity. It hydrolyzes bonds in a variety of phosphotriester (P-O), phosphonothioate (P-S), phosphofluoridate (P-F), and phosphonocyanate (F-CN) compounds. However, its catalytic efficiency varies markedly for different substrates, limiting the broad-range application of OPH as catalyst in the bioremediation of pesticides and chemical war agents. In the present study, pK(a) calculations and multiple explicit-solvent molecular dynamics (MD) simulations were performed to characterize and contrast the structural dynamics of OPH bound to two substrates hydrolyzed with very distinct catalytic efficiencies: the nerve agent soman (O-pinacolylmethylphosphonofluoridate) and the pesticide paraoxon (diethyl p-nitrophenyl phosphate). pK(a) calculations for the substrate-bound and unbound enzyme showed a significant pK(a) shift from standard values (ΔpK(a) = ±3 units) for residues His254 and Arg275. MD simulations of protonated His254 revealed a dynamic hydrogen bond network connecting the catalytic residue Asp301 via His254 to Asp232, Asp233, Arg275, and Asp235, and is consistent with a previously postulated proton relay mechanism to ferry protons away from the active site with substrates that do not require activation of the leaving group. Hydrogen bonds between Asp301 and His254 were persistent in the OPH-paraoxon complex but not in the OPH-soman one, suggesting a potential role for such interaction in the more efficient hydrolysis of paraoxon over soman by OPH. These results are in line with previous mutational studies of residue His254, which led to an increase of the catalytic efficiency of OPH over soman yet decreased its efficiency for paraoxon. In addition, comparative analysis of the molecular trajectories for OPH bound to soman and paraoxon suggests that binding of the latter facilitates the conformational transition of OPH from the open to the closed substate promoting a tighter binding of paraoxon.

Computational perspectives into plasmepsins structure-function relationship: implications to inhibitors design.

The development of efficient and selective antimalariais remains a challenge for the pharmaceutical industry. The aspartic proteases plasmepsins, whose inhibition leads to parasite death, are classified as targets for the design of potent drugs. Combinatorial synthesis is currently being used to generate inhibitor libraries for these enzymes, and together with computational methodologies have been demonstrated capable for the selection of lead compounds. The high structural flexibility of plasmepsins, revealed by their X-ray structures and molecular dynamics simulations, made even more complicated the prediction of putative binding modes, and therefore, the use of common computational tools, like docking and free-energy calculations. In this review, we revised the computational strategies utilized so far, for the structure-function relationship studies concerning the plasmepsin family, with special focus on the recent advances in the improvement of the linear interaction estimation (LIE) method, which is one of the most successful methodologies in the evaluation of plasmepsin-inhibitor binding affinity.

Free Energy Profiles along Consensus Normal Modes Provide Insight into HIV-1 Protease Flap Opening.

Describing biological macromolecular energetics from computer simulations can pose major challenges, and often necessitates enhanced conformational sampling. We describe the calculation of conformational free-energy profiles along carefully chosen collective coordinates: "consensus" normal modes, developed recently as robust alternatives to conventional normal modes. In an application to the HIV-1 protease, we obtain efficient sampling of significant flap opening movements governing inhibitor binding from relatively short simulations, in close correspondence with experimental results.