Glutamate 139 of tropomyosin is critical for cardiac thin filament blocked-state stabilization

The cardiac thin filament proteins troponin and tropomyosin control actomyosin formation and thus cardiac contractility. Calcium binding to troponin changes tropomyosin position along the thin filament, allowing myosin head binding to actin required for heart muscle contraction. The thin filament regulatory proteins are hot spots for genetic mutations causing heart muscle dysfunction. While much of the thin filament structure has been characterized, critical regions of troponin and tropomyosin involved in triggering conformational changes remain unresolved. A poorly resolved region, helix-4 (H4) of troponin I, is thought to stabilize tropomyosin in a position on actin that blocks actomyosin interactions at low calcium concentrations during muscle relaxation. We have proposed that contact between glutamate 139 on tropomyosin and positively charged residues on H4 leads to blocking-state stabilization. In this study, we attempted to disrupt these interactions by replacing E139 with lysine (E139K) to define the importance of this residue in thin filament regulation. Comparison of mutant and wild-type tropomyosin was carried out using in-vitro motility assays, actin co-sedimentation, and molecular dynamics simulations to determine perturbations in troponin-tropomyosin function caused by the tropomyosin mutation. Motility assays revealed that mutant thin filaments moved at higher velocity at low calcium with increased calcium sensitivity demonstrating that tropomyosin residue 139 is vital for proper tropomyosin-mediated inhibition during relaxation. Similarly, molecular dynamic simulations revealed a mutation-induced decrease in interaction energy between tropomyosin-E139K and troponin I (R170 and K174). These results suggest that salt-bridge stabilization of tropomyosin position by troponin IH4 is essential to prevent actomyosin interactions during cardiac muscle relaxation.

Myosin loop-4 is critical for optimal tropomyosin repositioning on actin during muscle activation and relaxation

During force-generating steps of the muscle crossbridge cycle, the tip of the myosin motor, specifically loop-4, contacts the tropomyosin cable of actin filaments. In the current study, we determined the corresponding effect of myosin loop-4 on the regulatory positioning of tropomyosin on actin. To accomplish this, we compared high-resolution cryo-EM structures of myosin S1-decorated thin filaments containing either wild-type or a loop-4 mutant construct, where the seven-residue portion of myosin loop-4 that contacts tropomyosin was replaced by glycine residues, thus removing polar side chains from residues 366-372. Cryo-EM analysis of fully decorated actin-tropomyosin filaments with wild-type and mutant S1, yielded 3.4-3.6 Å resolution reconstructions, with even higher definition at the actin-myosin interface. Loop-4 densities both in wild-type and mutant S1 were clearly identified, and side chains were resolved in the wild-type structure. Aside from loop-4, actin and myosin structural domains were indistinguishable from each other when filaments were decorated with either mutant or wild-type S1. In marked contrast, the position of tropomyosin on actin in the two reconstructions differed by 3 to 4 Å. In maps of filaments containing the mutant, tropomyosin was located closer to the myosin-head and thus moved in the direction of the C-state conformation adopted by myosin-free thin filaments. Complementary interaction energy measurements showed that tropomyosin in the mutant thin filaments sits on actin in a local energy minimum, whereas tropomyosin is positioned by wild-type S1 in an energetically unfavorable location. We propose that the high potential energy associated with tropomyosin positioning in wild-type filaments favors an effective transition to B- and C-states following release of myosin from the thin filaments during relaxation.

Modulation of cardiac thin filament structure by phosphorylated troponin-I analyzed by protein-protein docking and molecular dynamics simulation

Tropomyosin, controlled by troponin-linked Ca2+-binding, regulates muscle contraction by a macromolecular scale steric-mechanism that governs myosin-crossbridge-actin interactions. At low-Ca2+, C-terminal domains of troponin-I (TnI) trap tropomyosin in a position on thin filaments that interferes with myosin-binding, thus causing muscle relaxation. Steric inhibition is reversed at high-Ca2+ when TnI releases from F-actin-tropomyosin as Ca2+ and the TnI switch-peptide bind to the N-lobe of troponin-C (TnC). The opposite end of cardiac TnI contains a phosphorylation-sensitive ∼30 residue-long N-terminal peptide that is absent in skeletal muscle, and likely modifies these interactions in hearts. Here, PKA-dependent phosphorylation of serine 23 and 24 modulates Ca2+ and possibly switch-peptide binding to TnC, causing faster relaxation during the cardiac-cycle (lusitropy). The cardiac-specific N-terminal TnI domain is not captured in crystal structures of troponin or in cryo-EM reconstructions of thin filaments; thus, its global impact on thin filament structure and function is uncertain. Here, we used protein-protein docking and molecular dynamics simulation-based protocols to build a troponin model that was guided by and hence consistent with the recent seminal Yamada structure of Ca2+-activated thin filaments. We find that when present on thin filaments, phosphorylated Ser23/24 along with adjacent polar TnI residues interact closely with both tropomyosin and the N-lobe of TnC during our simulations. These interactions would likely bias tropomyosin to an off-state positioning on actin. In situ, such enhanced relaxation kinetics would promote cardiac lusitropy.

Immunoinformatic screening of Marburgvirus epitopes and computational investigations of epitope-allele complexes

Marburgvirus (MARV), a member of the Filovirus family, causes severe hemorrhagic fever in humans. Currently, there are no approved vaccines or post exposure treatment methods available against MARV. With the aim of identifying vaccine candidates against MARV, we employ different sequence-based computational methods to predict the MHC-I and MHC-II T-cell epitopes as well as B-cell epitopes for the complete MARV genome. We analyzed the variations in the predicted epitopes among four MARV variants, the Lake Victoria, Angola, Musoke, and Ravn. We used a consensus approach to identify several epitopes, including novel epitopes, and narrowed down the selection based on different parameters such as antigenicity and IC₅₀ values. The selected epitopes can be used in various vaccine constructs that give effective antibody responses. The MHC-I epitope-allele complexes for GP and NP with favorably low IC₅₀ values were investigated using molecular dynamics computations to determine the molecular details of the epitope-allele complexes. This study provides information for further experimental validation of the potential epitopes and the design and development of MARV vaccines.

Modeling Human Cardiac Thin Filament Structures

Striated muscle contraction is regulated in a calcium-dependent manner through dynamic motions of the tropomyosin/troponin polymer, a multicomponent complex wrapped around actin-containing thin filaments. Tropomyosin/troponin sterically blocks myosin-binding at low-calcium concentrations but moves to expose myosin-binding sites at high-calcium concentrations leading to force development. Understanding the key intermolecular interactions that define these dynamic motions will promote our understanding of mutation-induced contractile dysfunction that eventually leads to hypertrophic cardiomyopathy, dilated cardiomyopathy, and skeletal myopathies. Advancements in cryoelectron microscopy (cryoEM) have resulted in a partial elucidation of structures of the thin filament, revealing many atomic-level interactions between the component proteins and critical calcium-dependent conformational alterations. However, building models at the resolutions achieved can be challenging since landmarks in the maps are often missing or ambiguous. Therefore, current computational analyses including de novo structure prediction, protein-protein docking, molecular dynamics flexible fitting, and molecular dynamics simulations are needed to ensure good quality models. We review here our efforts to model the troponin T domain spanning the head-to-tail overlap domain of tropomyosin, improving previous models. Next, we refined the published cryoEM modeled structures, which had mistakenly compressed alpha helices, with a model that has expected helical parameters while matching densities in the cryoEM volume. Lastly, we used this model to reinterpret the interactions between tropomyosin and troponin I showing key features that hold the tropomyosin cable in its low-calcium, sterically blocking position. These revised thin filament models show improved intermolecular interactions in the key low- and high-calcium regulatory states, providing novel insights into function.

Ebola virus protein VP40 binding to Sec24c for transport to the plasma membrane

Outbreaks of the Ebola virus (EBOV) continue to occur and while a vaccine and treatment are now available, there remains a dearth of options for those who become sick with EBOV disease. An understanding at the atomic and molecular level of the various steps in the EBOV replication cycle can provide molecular targets for disrupting the virus. An important step in the EBOV replication cycle is the transport of EBOV structural matrix VP40 protein molecules to the plasma membrane inner leaflet, which involves VP40 binding to the host cell's Sec24c protein. Though some VP40 residues involved in the binding are known, the molecular details of VP40-Sec24c binding are not known. We use various molecular computational techniques to investigate the molecular details of how EBOV VP40 binds with the Sec24c complex of the ESCRT-I pathway. We employed different docking programs to identify the VP40-binding site on Sec24c and then performed molecular dynamics simulations to determine the atomic details and binding interactions of the complex. We also investigated how the inter-protein interactions of the complex are affected upon mutations of VP40 amino acids in the Sec24c-binding region. Our results provide a molecular basis for understanding previous coimmunoprecipitation experimental studies. In addition, we found that VP40 can bind to a site on Sec24c that can also bind Sec23 and suggests that VP40 may use the COPII transport mechanism in a manner that may not need the Sec23 protein in order for VP40 to be transported to the plasma membrane.

C-terminal troponin-I residues trap tropomyosin in the muscle thin filament blocked-state

Tropomyosin and troponin regulate muscle contraction by participating in a macromolecular scale steric-mechanism to control myosin-crossbridge - actin interactions and consequently contraction. At low-Ca²⁺, the C-terminal 30% of troponin subunit-I (TnI) is proposed to trap tropomyosin in a position on thin filaments that sterically interferes with myosin-binding, thus causing muscle relaxation. In contrast, at high-Ca²⁺, inhibition is released after the C-terminal domains dissociate from F-actin-tropomyosin as its component switch-peptide domain binds to the N-lobe of troponin-C (TnC). Recent, paradigm-shifting, cryo-EM reconstructions by the Namba group have revealed density attributed to TnI along cardiac muscle thin filaments at both low- and high-Ca²⁺ concentration. Modeling the reconstructions showed expected high-Ca²⁺ hydrophobic interactions of the TnI switch-peptide and TnC. However, under low-Ca²⁺ conditions, sparse interactions of TnI and tropomyosin, and in particular juxtaposition of non-polar switch-peptide residues and charged tropomyosin amino acids in the published model seem difficult to reconcile with an expected steric-blocking conformation. This anomaly is likely due to inaccurate fitting of tropomyosin into the cryo-EM volume. In the current study, the low-Ca²⁺ cryo-EM volume was fitted with a more accurate tropomyosin model and representation of cardiac TnI. Our results show that at low-Ca²⁺ a cluster of hydrophobic residues at the TnI switch-peptide and adjacent H₄ helix (Ala149, Ala151, Met 154, Leu159, Gly160, Ala161, Ala163, Leu167, Leu169, Ala171, Leu173) draw-in tropomyosin surface residues (Ile143, Ile146, Ala151, Ile154), presumably attracting the entire tropomyosin cable to its myosin-blocking position on actin. The modeling confirms that neighboring TnI "inhibitory domain" residues (Arg145, Arg148) bind to thin filaments at actin residue Asp25, as previously suggested. ClusPro docking of TnI residues 137-184 to actin-tropomyosin, including the TnI inhibitory-domain, switch-peptide and Helix H₄, verified the modeled configuration. Our residue-to-residue contact-mapping of the TnI-tropomyosin association lends itself to experimental validation and functional localization of disease-bearing mutations.

Cryo-EM and Molecular Docking Shows Myosin Loop 4 Contacts Actin and Tropomyosin on Thin Filaments

The motor protein myosin drives muscle and nonmuscle motility by binding to and moving along actin of thin filaments. Myosin binding to actin also modulates interactions of the regulatory protein, tropomyosin, on thin filaments, and conversely tropomyosin affects myosin binding to actin. Insight into this reciprocity will facilitate a molecular level elucidation of tropomyosin regulation of myosin interaction with actin in muscle contraction, and in turn, promote better understanding of nonmuscle cell motility. Indeed, experimental approaches such as fiber diffraction, cryoelectron microscopy, and three-dimensional reconstruction have long been used to define regulatory interaction of tropomyosin and myosin on actin at a structural level. However, their limited resolution has not proven sufficient to determine tropomyosin and myosin contacts at an atomic-level and thus to fully substantiate possible functional contributions. To overcome this deficiency, we have followed a hybrid approach by performing new cryogenic electron microscopy reconstruction of myosin-S1-decorated F-actin-tropomyosin together with atomic scale protein-protein docking of tropomyosin to the EM models. Here, cryo-EM data were derived from filaments reconstituted with α1-actin, cardiac αα-tropomyosin, and masseter muscle β-myosin complexes; masseter myosin, which shares sequence identity with β-cardiac myosin-heavy chain, was used because of its stability in vitro. The data were used to build an atomic model of the tropomyosin cable that fits onto the actin filament between the tip of the myosin head and a cleft on the innermost edge of actin subunits. The docking and atomic scale fitting showed multiple discrete interactions of myosin loop 4 and acidic residues on successive 39-42 residue-long tropomyosin pseudorepeats. The contacts between S1 and tropomyosin on actin appear to compete with and displace ones normally found between actin and tropomyosin on myosin-free thin filaments in relaxed muscle, thus restructuring the filament during myosin-induced activation.

Protein-Protein Docking Reveals Dynamic Interactions of Tropomyosin on Actin Filaments

Experimental approaches such as fiber diffraction and cryo-electron microscopy reconstruction have defined regulatory positions of tropomyosin on actin but have not, as yet, succeeded at determining key atomic-level contacts between these proteins or fully substantiated the dynamics of their interactions at a structural level. To overcome this deficiency, we have previously employed computational approaches to deduce global dynamics of thin filament components by energy landscape determination and molecular dynamics simulations. Still, these approaches remain computationally challenging for any complex and large macromolecular assembly like the thin filament. For example, tropomyosin cable wrapping around actin of thin filaments features both head-to-tail polymeric interactions and local twisting, both of which depart from strict superhelical symmetry. This produces a complex energy surface that is difficult to model and thus to evaluate globally. Therefore, at this stage of our understanding, assessing global molecular dynamics can prove to be inherently impractical. As an alternative, we adopted a "divide and conquer" protocol to investigate actin-tropomyosin interactions at an atomistic level. Here, we first employed unbiased protein-protein docking tools to identify binding specificity of individual tropomyosin pseudorepeat segments over the actin surface. Accordingly, tropomyosin "ligand" segments were rotated and translated over potential "target" binding sites on F-actin where the corresponding interaction energetics of billions of conformational poses were ranked by the programs PIPER and ClusPro. These data were used to assess favorable interactions and then to rebuild models of seamless and continuous tropomyosin cables over the F-actin substrate, which were optimized further by flexible fitting routines and molecular dynamics. The models generated azimuthally distinct regulatory positions for tropomyosin cables along thin filaments on actin dominated by stereo-specific head-to-tail overlap linkage. The outcomes are in good agreement with current cryo-electron microscopy topology and consistent with long-thought residue-to-residue interactions between actin and tropomyosin.

Docking Troponin T onto the Tropomyosin Overlapping Domain of Thin Filaments

Complete description of thin filament conformational transitions accompanying muscle regulation requires ready access to atomic structures of actin-bound tropomyosin-troponin. To date, several molecular-docking protocols have been employed to identify troponin interactions on actin-tropomyosin because high-resolution experimentally determined structures of filament-associated troponin are not available. However, previously published all-atom models of the thin filament show chain separation and corruption of components during our molecular dynamics simulations of the models, implying artifactual subunit organization, possibly due to incorporation of unorthodox tropomyosin-TnT crystal structures and complex FRET measurements during model construction. For example, the recent Williams et al. (2016) atomistic model of the thin filament displays a paucity of salt bridges and hydrophobic complementarity between the TnT tail (TnT1) and tropomyosin, which is difficult to reconcile with the high, 20 nM K_d binding of TnT onto tropomyosin. Indeed, our molecular dynamics simulations show the TnT1 component in their model partially dissociates from tropomyosin in under 100 ns, whereas actin-tropomyosin and TnT1 models themselves remain intact. We therefore revisited computational work aiming to improve TnT1-thin filament models by employing unbiased docking methodologies, which test billions of trial rotations and translations of TnT1 over three-dimensional grids covering end-to-end bonded tropomyosin alone or tropomyosin on F-actin. We limited conformational searches to the association of well-characterized TnT1 helical domains and either isolated tropomyosin or actin-tropomyosin yet avoided docking TnT domains that lack known or predicted structure. The docking programs PIPER and ClusPro were used, followed by interaction energy optimization and extensive molecular dynamics. TnT1 docked to either side of isolated tropomyosin but uniquely onto one location of actin-bound tropomyosin. The antiparallel interaction with tropomyosin contained abundant salt bridges and intimately integrated hydrophobic networks joining TnT1 and the tropomyosin N-/C-terminal overlapping domain. The TnT1-tropomyosin linkage yields well-defined molecular crevices. Interaction energy measurements strongly favor this TnT1-tropomyosin design over previously proposed models.

Conformational Flexibility of the Protein-Protein Interfaces of the Ebola Virus VP40 Structural Matrix Filament

The Ebola virus (EBOV) is a virulent pathogen that causes severe hemorrhagic fever with a high fatality rate in humans. The EBOV transformer protein VP40 plays crucial roles in viral assembly and budding at the plasma membrane of infected cells. One of VP40's roles is to form the long, flexible, pleomorphic filamentous structural matrix for the virus. Each filament contains three unique interfaces: monomer NTD-NTD to form a dimer, dimer-to-dimer NTD-NTD oligomerization to form a hexamer, and end-to-end hexamer CTD-CTD to build the filament. However, the atomic-level details of conformational flexibility of the VP40 filament are still elusive. In this study, we have performed explicit-solvent, all-atom molecular dynamic simulations to explore the conformational flexibility of the three different interface structures of the filament. Using dynamic network analysis and other calculational methods, we find that the CTD-CTD hexamer interface with weak interdomain amino acid communities is the most flexible, and the NTD-NTD oligomer interface with strong interdomain communities is the least flexible. Our study suggests that the high flexibility of the CTD-CTD interface may be essential for the supple bending of the Ebola filovirus, and such flexibility may present a target for molecular interventions to disrupt the Ebola virus functioning.

The adipokinetic hormones and their cognate receptor from the desert locust, Schistocerca gregaria: solution structure of endogenous peptides and models of their binding to the receptor

Background

Neuropeptides exert their activity through binding to G protein-coupled receptors (GPCRs). GPCRs are well-known drug targets in the pharmaceutical industry and are currently discussed as targets to control pest insects. Here, we investigate the neuropeptide adipokinetic hormone (AKH) system of the desert locust Schistocerca gregaria. The desert locust is known for its high reproduction, and for forming devastating swarms consisting of billions of individual insects. It is also known that S. gregaria produces three different AKHs as ligands but has only one AKH receptor (AKHR). The AKH system is known to be essential for metabolic regulation, which is necessary for reproduction and flight activity.

Methods

Nuclear magnetic resonance techniques (NMR) in a dodecylphosphocholin (DPC) micelle solution were used to determine the structure of the three AKHs. The primary sequence of the S. gregaria AKHR was used to construct a 3D molecular model. Next, the three AKHs were individually docked to the receptor, and dynamic simulation of the whole ligand–receptor complex in a model membrane was performed.

Results

Although the three endogenous AKHs of S. gregaria have quite different amino acids sequences and chain length (two octa- and one decapeptide), NMR experiments assigned a turn structure in DPC micelle solution for all. The GPCR-ModSim program identified human kappa opioid receptor to be the best template after which the S. gregaria AKHR was modeled. All three AKHs were found to have the same binding site on this receptor, interact with similar residues of the receptor and have comparable binding constants. Molecular switches were also identified; the movement of the receptor could be visually shown when ligands (AKHs) were docked and the receptor was activated.

Conclusions

The study proposes a model of binding of the three endogenous ligands to the one existing AKHR in the desert locust and paves the way to use such a model for the design of peptide analogs and finally, peptide mimetics, in the search for novel species-specific insecticides based on receptor–ligand interaction.

Membrane pore formation and ion selectivity of the Ebola virus delta peptide

The Ebola virus delta peptide homo-oligomerizes in the host cell membrane to form amphipathic pores that alter the membrane properties.

A cylindrical assembly model and dynamics of the Ebola virus VP40 structural matrix

The Ebola filovirus causes severe hemorrhagic fever with a high fatality rate in humans. The primary structural matrix protein VP40 displays transformer-protein characteristics and exists in different conformational and oligomeric states. VP40 plays crucial roles in viral assembly and budding at the plasma membrane of the infected cells and is capable of forming virus-like particles without the need for other Ebola proteins. However, no experimental three-dimensional structure for any filovirus VP40 cylindrical assembly matrix is currently available. Here, we use a protein-protein docking approach to develop cylindrical assembly models for an Ebola virion and also for a smaller structural matrix that does not contain genetic material. These models match well with the 2D averages of cryo-electron tomograms of the authentic virion. We also used all-atom molecular dynamics simulations to investigate the stability and dynamics of the cylindrical models and the interactions between the side-by-side hexamers to determine the amino acid residues that are especially important for stabilizing the hexamers in the cylindrical ring configuration matrix assembly. Our models provide helpful information to better understand the assembly processes of filoviruses and such structural studies may also lead to the design and development of antiviral drugs.

Data for the homology modelling of the red pigment-concentrating hormone receptor (Dappu-RPCHR) of the crustacean Daphnia pulex, and docking of its cognate agonist (Dappu-RPCH)

The data presented in this article are related to the publication “Interaction of the red pigment-concentrating hormone of the crustacean Daphnia pulex, with its cognate receptor, Dappu-RPCHR: A nuclear magnetic resonance and modeling study” (Jackson et al., 2017) [1]. This article contains the data for homology modeling of the red pigment-concentrating hormone (RPCH) receptor of the water flea, Daphnia pulex (Dappu-RPCHR), which was constructed from its primary sequence. This is the first 3D model of a crustacean G-protein coupled receptor. Docking of the agonist, pGlu-Val-Asn-Phe-Ser-Thr-Ser-Trp amide (Dappu-RPCH), was used to find a binding pocket on the receptor and compared to the binding pocket of the adipokinetic hormone (AKH) receptor from the malaria mosquito. Data for the receptor, with and without loop refinement, together with the docked agonist, are presented.

Interaction of the red pigment-concentrating hormone of the crustacean Daphnia pulex, with its cognate receptor, Dappu-RPCHR: A nuclear magnetic resonance and modeling study

The primary sequence of the red pigment-concentrating hormone (RPCH) receptor of the water flea, Daphnia pulex, was used in homology modeling to construct the first 3D model of a crustacean G-protein coupled receptor, Dappu-RPCHR. This receptor was found to belong to the class A subfamily of GPCRs with a disulfide bridge between Cys⁷² and Cys¹⁵⁰ and an ionic lock between Arg⁹⁷ and Thr²²⁴ and Thr²²⁰. NMR restrained molecular dynamics was used to determine the structure of an agonist, Dappu-RPCH, in a membrane-mimicking environment. The agonist was found to be flexible but has two main conformations in solution, both having β-turns. Docking of the predominant structure was used to find a binding pocket on the receptor. The pocket's spatial location was similar to that of the AKH receptor of Anopheles gambiae. The binding affinity was -69kcalmol^-1 with the N-terminus of Dappu-RPCH inserted between helices 4 and 6, and the C-terminus interacting with extra-cellular loop, ECL2. Upon binding, H-bonding to the peptide may activate the receptor. This development of the first Dappu-RPCH/Dappu-RPCHR model could be useful for understanding ligand-receptor interactions in crustaceans.

Identification of steroid-like natural products as antiplasmodial agents by 2D and 3D similarity-based virtual screening

The emergence of drug resistance in Plasmodium falciparum to available antimalarial drugs has challenged current antimalarial treatments. New antimalarials, particularly those with novel mechanisms of action and no cross resistance to current drugs, are therefore urgently needed. To identify new growth inhibitors of Plasmodium falciparum, 2D and 3D similarity-based virtual screening methods were employed in parallel with an in-house database of steroid-type natural products using fusidic acid as a search query. The resulting hit compounds were further filtered based on the predicted partition coefficient, log P. The virtual screening strategy resulted in the identification of nine new compounds that inhibited parasite growth with IC₅₀ values of <20 μM. Four compounds exhibited IC₅₀ values in the range of 1.39-3.45 μM and three of which showed a promising selectivity index. Further, the predicted ADME properties of the four most active compounds were found to be comparable to fusidic acid. These compounds can be further explored using structural modifications in the identification and development of more potent parasite growth inhibitors with improved selectivity.

Identification of New Human Malaria Parasite Plasmodium falciparum Dihydroorotate Dehydrogenase Inhibitors by Pharmacophore and Structure-Based Virtual Screening

Plasmodium falciparum dihydroorotate dehydrogenase (PfDHODH), a key enzyme in the de novo pyrimidine biosynthesis pathway, which the Plasmodium falciparum relies on exclusively for survival, has emerged as a promising target for antimalarial drugs. In an effort to discover new and potent PfDHODH inhibitors, 3D-QSAR pharmacophore models were developed based on the structures of known PfDHODH inhibitors and the validated Hypo1 model was used as a 3D search query for virtual screening of the National Cancer Institute database. The virtual hit compounds were further filtered based on molecular docking and Molecular Mechanics/Generalized Born Surface Area binding energy calculations. The combination of the pharmacophore and structure-based virtual screening resulted in the identification of nine new compounds that showed >25% inhibition of PfDHODH at a concentration of 10 μM, three of which exhibited IC50 values in the range of 0.38-20 μM. The most active compound, NSC336047, displayed species-selectivity for PfDHODH over human DHODH and inhibited parasite growth with an IC50 of 26 μM. In addition to this, 13 compounds inhibited parasite growth with IC50 values of ≤ 50 μM, 4 of which showed IC50 values in the range of 5-12 μM. These compounds could be further explored in the identification and development of more potent PfDHODH and parasite growth inhibitors.

Computer-Aided Drug Discovery Approaches against the Tropical Infectious Diseases Malaria, Tuberculosis, Trypanosomiasis, and Leishmaniasis

Despite the tremendous improvement in overall global health heralded by the adoption of the Millennium Declaration in the year 2000, tropical infections remain a major health problem in the developing world. Recent estimates indicate that the major tropical infectious diseases, namely, malaria, tuberculosis, trypanosomiasis, and leishmaniasis, account for more than 2.2 million deaths and a loss of approximately 85 million disability-adjusted life years annually. The crucial role of chemotherapy in curtailing the deleterious health and economic impacts of these infections has invigorated the search for new drugs against tropical infectious diseases. The research efforts have involved increased application of computational technologies in mainstream drug discovery programs at the hit identification, hit-to-lead, and lead optimization stages. This review highlights various computer-aided drug discovery approaches that have been utilized in efforts to identify novel antimalarial, antitubercular, antitrypanosomal, and antileishmanial agents. The focus is largely on developments over the past 5 years (2010-2014).

Synthesis of fusidic acid bioisosteres as antiplasmodial agents and molecular docking studies in the binding site of elongation factor-G

Structural modifications through bioisosteric approach yielded fusidic acid analogues with 2–35 folds increase in antiplasmodial activity as compared to fusidic acid.

Synthesis of new verapamil analogues and their evaluation in combination with rifampicin against Mycobacterium tuberculosis and molecular docking studies in the binding site of efflux protein Rv1258c

New verapamil analogues were synthesized and their inhibitory activities against Mycobacterium tuberculosis H37Rv determined in vitro alone and in combination with rifampicin (RIF). Some analogues showed comparable activity to verapamil and exhibited better synergies with RIF. Molecular docking studies of the binding sites of Rv1258c, a M. tuberculosis efflux protein previously implicated in intrinsic resistance to RIF, suggested a potential rationale for the superior synergistic interactions observed with some analogues.