PvP01-DB: computational structural and functional characterization of soluble proteome of PvP01 strain of Plasmodium vivax.

Despite Plasmodium vivax being the main offender in the majority of malarial infections, very little information is available about its adaptation and development in humans. Its capability for activating relapsing infections through its dormant liver stage and resistance to antimalarial drugs makes it as one of the major challenges in eradicating malaria. Noting the immediate necessity for the availability of a comprehensive and reliable structural and functional repository for P. vivax proteome, here we developed a web resource for the new reference genome, PvP01, furnishing information on sequence, structure, functions, active sites and metabolic pathways compiled and predicted using some of the state-of-the-art methods in respective fields. The PvP01 web resource comprises organized data on the soluble proteome consisting of 3664 proteins in blood and liver stages of malarial cycle. The current public resources represent only 163 proteins of soluble proteome of PvP01, with complete information about their molecular function, biological process and cellular components. Also, only 46 proteins of P. vivax have experimentally determined structures. In this milieu of extreme scarcity of structural and functional information, PvP01 web resource offers meticulously validated structures of 3664 soluble proteins. The sequence and structure-based functional characterization led to a quantum leap from 163 proteins available presently to whole soluble proteome offered through PvP01 web resource. We believe PvP01 web resource will serve the researchers in identifying novel protein drug targets and in accelerating the development of structure-based new drug candidates to combat malaria. Database Availability: http://www.scfbio-iitd.res.in/PvP01.

Discrete analyses of protein dynamics.

Protein structures are highly dynamic macromolecules. This dynamics is often analysed through experimental and/or computational methods only for an isolated or a limited number of proteins. Here, we explore large-scale protein dynamics simulation to observe dynamics of local protein conformations using different perspectives. We analysed molecular dynamics to investigate protein flexibility locally, using classical approaches such as RMSf, solvent accessibility, but also innovative approaches such as local entropy. First, we focussed on classical secondary structures and analysed specifically how ß-strand, ß-turns, and bends evolve during molecular simulations. We underlined interesting specific bias between ß-turns and bends, which are considered as the same category, while their dynamics show differences. Second, we used a structural alphabet that is able to approximate every part of the protein structures conformations, namely protein blocks (PBs) to analyse (i) how each initial local protein conformations evolve during dynamics and (ii) if some exchange can exist among these PBs. Interestingly, the results are largely complex than simple regular/rigid and coil/flexible exchange. AbbreviationsNeqnumber of equivalentPBProtein BlocksPDBProtein DataBankRMSfroot mean square fluctuationsCommunicated by Ramaswamy H. Sarma.

Computational Simulations Identified Two Candidate Inhibitors of Cdk5/p25 to Abrogate Tau-associated Neurological Disorders.

Deregulation of Cdk5 is a hallmark in neurodegenerative diseases and its complex with p25 forms Cdk5/p25, thereby causes severe neuropathological insults. Cdk5/p25 abnormally phosphorylates tau protein, and induces tau-associated neurofibrillary tangles in neurological disorders. Therefore, the pharmacological inhibition of Cdk5/p25 alleviates tau-associated neurological disorders. Herein, computational simulations probed two candidate inhibitors of Cdk5/p25. Structure-based pharmacophore investigated the essential complementary chemical features of ATP-binding site of Cdk5 in complex with roscovitine. Resultant pharmacophore harbored polar interactions with Cys83 and Asp86 residues and non-polar interactions with Ile10, Phe80, and Lys133 residues of Cdk5. The chemical space of selected pharmacophore was comprised of two hydrogen bond donors, one hydrogen bond acceptor, and three hydrophobic features. Decoy test validation of pharmacophore obtained highest Guner-Henry score (0.88) and enrichment factor score (7.23). The screening of natural product drug-like databases by validated pharmacophore retrieved 1126 compounds as candidate inhibitors of Cdk5/p25. The docking of candidate inhibitors filtered 10 molecules with docking score >80.00 and established polar and non-polar interactions with the ATP-binding site residues of Cdk5/p25. Finally, molecular dynamics simulation and binding free energy analyses identified two candidate inhibitors of Cdk5/p25. During 30 ns simulation, the candidate inhibitors established <3.0 Šroot mean square deviation and stable hydrogen bond interactions with the ATP-binding site residues of Cdk5/p25. The final candidate inhibitors obtained lowest binding free energies of -122.18 kJ/mol and - 117.26 kJ/mol with Cdk5/p25. Overall, we recommend two natural product candidate inhibitors to target the pharmacological inhibition of Cdk5/p25 in tau-associated neurological disorders.

How Structural Biologists and the Protein Data Bank Contributed to Recent FDA New Drug Approvals.

Discovery and development of 210 new molecular entities (NMEs; new drugs) approved by the US Food and Drug Administration 2010-2016 was facilitated by 3D structural information generated by structural biologists worldwide and distributed on an open-access basis by the PDB. The molecular targets for 94% of these NMEs are known. The PDB archive contains 5,914 structures containing one of the known targets and/or a new drug, providing structural coverage for 88% of the recently approved NMEs across all therapeutic areas. More than half of the 5,914 structures were published and made available by the PDB at no charge, with no restrictions on usage >10 years before drug approval. Citation analyses revealed that these 5,914 PDB structures significantly affected the very large body of publicly funded research reported in publications on the NME targets that motivated biopharmaceutical company investment in discovery and development programs that produced the NMEs.

Parsing the functional specificity of Siderocalin/Lipocalin 2/NGAL for siderophores and related small-molecule ligands.

Siderocalin/Lipocalin 2/Neutrophil Gelatinase Associated Lipocalin/24p3 is an innate immune system protein with bacteriostatic activity, acting by tightly binding and sequestering diverse catecholate and mixed-type ferric siderophores from enteric bacteria and mycobacteria. Bacterial virulence achieved through siderophore modifications, or utilization of alternate siderophores, can be explained by evasion of Siderocalin binding. Siderocalin has also been implicated in a wide variety of disease processes, though often in seemingly contradictory ways, and has been proposed to bind to a broader array of ligands beyond siderophores. Using structural, directed mutational, and binding studies, we have sought to rigorously test, and fully elucidate, the Siderocalin recognition mechanism. Several proposed ligands fail to meet rigorous binding criteria, including the bacterial siderophore pyochelin, the iron-chelating catecholamine hormone norepinephrine, and the bacterial second messenger cyclic diguanylate monophosphate. While possessing a remarkably rigid structure, in principle simplifying analyses of ligand recognition, understanding Scn recognition is complicated by the observed conformational and stoichiometric plasticity, and instability, of its bona fide siderophore ligands. Since the role of Siderocalin at the early host/pathogen interface is to compete for bacterial ferric siderophores, we also analyzed how bacterial siderophore binding proteins and enzymes alternately recognize siderophores that efficiently bind to, or evade, Siderocalin sequestration - including determining the crystal structure of Bacillus cereus YfiY bound to schizokinen. These studies combine to refine the potential physiological functions of Siderocalin by defining its multiplexed recognition mechanism.

Computational insights into the interaction of small molecule inhibitors with HRI kinase domain.

The Heme-Regulated Inhibitor (HRI) kinase regulates globin synthesis in a heme-dependent manner in reticulocytes and erythroid cells in bone marrow. Inhibitors of HRI have been proposed to lead to an increased amount of haemoglobin, benefitting anaemia patients. A series of indeno[1,2-c]pyrazoles were discovered to be the first known in vitro inhibitors of HRI. However, the structural mechanism of inhibition is yet to be understood. The aim of this study was to unravel the binding mechanism of these inhibitors using molecular dynamic simulations and docking. The docking scores were observed to correlate well with experimentally determined pIC50 values. The inhibitors were observed to bind in the ATP-binding site forming hydrogen bonds with the hinge region and van der Waals interactions with non-polar residues in the binding site. Further, quantitative structure-activity relationship (QSAR) studies were performed to correlate the structural features of the inhibitors with their biological activity. The developed QSAR models were found to be statistically significant in terms of internal and external predictabilities. The presence of chlorine atoms and the hydroxymethyl groups were found to correlate with higher activity. The identified binding modes and the descriptors can support future rational identification of more potent and selective small molecule inhibitors for this kinase which are of therapeutic importance in the context of various human pathological disorders.

Identification of anti-filarial leads against aspartate semialdehyde dehydrogenase of Wolbachia endosymbiont of Brugia malayi: combined molecular docking and molecular dynamics approaches.

Lymphatic filariasis is a debilitating vector borne parasitic disease that infects human lymphatic system by nematode Brugia malayi. Currently available anti-filarial drugs are effective only on the larval stages of parasite. So far, no effective drugs are available for humans to treat filarial infections. In this regard, aspartate semialdehyde dehydrogenase (ASDase) in lysine biosynthetic pathway from Wolbachia endosymbiont Brugia malayi represents an attractive therapeutic target for the development of novel anti-filarial agents. In this present study, molecular modeling combined with molecular dynamics simulations and structure-based virtual screening were performed to identify potent lead molecules against ASDase. Based on Glide score, toxicity profile, binding affinity and mode of interactions with the ASDase, five potent lead molecules were selected. The molecular docking and dynamics results revealed that the amino acid residues Arg103, Asn133, Cys134, Gln161, Ser164, Lys218, Arg239, His246, and Asn321 plays a crucial role in effective binding of Top leads into the active site of ASDase. The stability of the ASDase-lead complexes was confirmed by running the 30 ns molecular dynamics simulations. The pharmacokinetic properties of the identified lead molecules are in the acceptable range. Furthermore, density functional theory and binding free energy calculations were performed to rank the lead molecules. Thus, the identified lead molecules can be used for the development of anti-filarial agents to combat the pathogenecity of Brugia malayi.

Virtual analysis of structurally diverse synthetic analogs as inhibitors of snake venom secretory phospholipase A2.

Due to the toxic pathophysiological role of snake venom phospholipase A2 (PLA2 ), its compelling limitations to anti-venom therapy in humans and the need for alternative therapy foster considerable pharmacological interest towards search of PLA2 specific inhibitors. In this study, an integrated approach involving homology modeling, molecular dynamics and molecular docking studies on VRV-PL-V (Vipera russellii venom phospholipase A2 fraction-V) belonging to Group II-B secretory PLA2 from Daboia russelli pulchella is carried out in order to study the structure-based inhibitor design. The accuracy of the model was validated using multiple computational approaches. The molecular docking study of this protein was undertaken using different classes of experimentally proven, structurally diverse synthetic inhibitors of secretory PLA2 whose selection is based on IC50 value that ranges from 25 µM to 100 µM. Estimation of protein-ligand contacts by docking analysis sheds light on the importance of His 47 and Asp 48 within the VRV-PL-V binding pocket as key residue for hydrogen bond interaction with ligands. Our virtual analysis revealed that compounds with different scaffold binds to the same active site region. ADME analysis was also further performed to filter and identify the best potential specific inhibitor against VRV-PL-V. Additionally, the e-pharmacophore was generated for the best potential specific inhibitor against VRV-PL-V and reported here. The present study should therefore play a guiding role in the experimental design of VRV-PL-V inhibitors that may provide better therapeutic molecular models for PLA2 recognition and anti-ophidian activity.

Kinesin-5: cross-bridging mechanism to targeted clinical therapy.

Kinesin motor proteins comprise an ATPase superfamily that works hand in hand with microtubules in every eukaryote. The mitotic kinesins, by virtue of their potential therapeutic role in cancerous cells, have been a major focus of research for the past 28 years since the discovery of the canonical Kinesin-1 heavy chain. Perhaps the simplest player in mitotic spindle assembly, Kinesin-5 (also known as Kif11, Eg5, or kinesin spindle protein, KSP) is a plus-end-directed motor localized to interpolar spindle microtubules and to the spindle poles. Comprised of a homotetramer complex, its function primarily is to slide anti-parallel microtubules apart from one another. Based on multi-faceted analyses of this motor from numerous laboratories over the years, we have learned a great deal about the function of this motor at the atomic level for catalysis and as an integrated element of the cytoskeleton. These data have, in turn, informed the function of motile kinesins on the whole, as well as spearheaded integrative models of the mitotic apparatus in particular and regulation of the microtubule cytoskeleton in general. We review what is known about how this nanomotor works, its place inside the cytoskeleton of cells, and its small-molecule inhibitors that provide a toolbox for understanding motor function and for anticancer treatment in the clinic.

Computational studies of novel carbonyl-containing diazabicyclic ligands interacting with α4ß2 nicotinic acetylcholine receptor (nAChR) reveal alternative binding modes.

We have carried out computational studies on interactions of diazabicyclic amide analogs with α4ß2 nAChR using homology modeling, docking and pharmacophore elucidation techniques. We have found alternative ligand binding modes in most cases. All these diverse poses exhibit the quintessential hydrogen-bonding interaction between the ligand basic nitrogen and the backbone carbonyl oxygen atom of the highly conserved Trp-149. This hydrogen bond was always found to be shorter than the one contracted by the ligand carbonyl group and a second hydrogen-bond made by the cationic center with Tyr-93 of the principal face of the protein. In most of the poses observed, cation-π interactions involved three aromatic residues located in the principal face of the protein: Trp-149, Tyr-190 and Tyr-197. The latter amino acid residue appears to often donate a hydrogen-bond to the ligand carbonyl oxygen atom. We also describe two rings of alternative receptor-based hydrogen-bond donor features equidistantly separated from the carbonyl oxygen of the highly conserved Trp-149 approximately by 5 and 8Å, respectively. These findings could be exploited to design diverse and selective novel chemical libraries for the treatment of diseases and conditions where the α4ß2 nAChR is disrupted, such as Alzheimer disease, Parkinson's disease and l-dopa-induced dyskinesia (LID).

QM-polarized ligand docking accurately predicts the trend in binding affinity of a series of arylmethylene quinuclidine-like derivatives at the α4ß2 and α3ß4 nicotinic acetylcholine receptors (nAChRs).

Compounds containing a quinuclidine scaffold are promising drug candidates for pharmacological management of the central nervous system (CNS) pathologies implicating nAChRs. We have carried out binding affinity and in-silico docking studies of arylmethylene quinuclidine-like derivatives at the α4ß2 receptor using in-vitro receptor binding assay and comparative modeling, respectively. We found that introducing a hydrogen-bond acceptor into the 3-benzylidene quinuclidine derivative resulted in a 266-fold increase in binding affinity and confers agonism properties. By contrast, addition of a phenyl group to 3-benzylidene quinuclidine derivative only results in an 18-fold increase in binding affinity, without conferring agonism. We also found that docking into the orthosteric binding site of the α4ß2 nAChR is consistent with the fact that the basic nitrogen atom donates a hydrogen-bond to the carbonyl group of the highly conserved Trp-149, as initially observed by Dougherty and co-workers.(1) The experimentally-observed trend in binding affinity at both α4ß2 and α3ß4 nAChRs was accurately and independently confirmed by quantum mechanics (QM)-polarized docking. The reduction in binding affinity to the α3ß4 subtype primarily results from a dampening of both coulombic and cation-π interactions.

Discovery of the first inhibitors of bacterial enzyme D-aspartate ligase from Enterococcus faecium (Aslfm).

The D-aspartate ligase of Enterococcus faecium (Aslfm) is an attractive target for the development of narrow-spectrum antibacterial agents that are active against multidrug-resistant E. faecium. Although there is currently little available information regarding the structural characteristics of Aslfm, we exploited the knowledge that this enzyme belongs to the ATP-grasp superfamily to target its ATP binding site. In the first design stage, we synthesized and screened a small library of known ATP-competitive inhibitors of ATP-grasp enzymes. A series of amino-oxazoles derived from bacterial biotin carboxylase inhibitors showed low micromolar activity. The most potent inhibitor compound 12, inhibits Aslfm with a Ki value of 2.9 µM. In the second design stage, a validated ligand-based pharmacophore modeling approach was used, taking the newly available inhibition data of an initial series of compounds into account. Experimental evaluation of the virtual screening hits identified two novel structural types of Aslfm inhibitors with 7-amino-9H-purine (18) and 7-amino-1H-pyrazolo[3,4-d]pyrimidine (30 and 34) scaffolds, and also with Ki values in the low micromolar range. Investigation the inhibitors modes of action confirmed that these compounds are competitive with respect to the ATP molecule. The binding of inhibitors to the target enzyme was also studied using isothermal titration calorimetry (ITC). Compounds 6, 12, 18, 30 and 34 represent the first inhibitors of Aslfm reported to date, and are an important step forward in combating infections due to E. faecium.

Unanswered questions about the structure of cytochrome bc1 complexes.

X-ray crystal structures of bc1 complexes obtained over the last 15 years have provided a firm structural basis for our understanding of the complex. For the most part there is good agreement between structures from different species, different crystal forms, and with different inhibitors bound. In this review we focus on some of the remaining unexplained differences, either between the structures themselves or the interpretations of the structural observations. These include the structural basis for the motion of the Rieske iron-sulfur protein in response to inhibitors, a possible conformational change involving tyrosine132 of cytochrome (cyt) b, the presence of cis-peptides at the beginnings of transmembrane helices C, E, and H, the structural insight into the function of the so-called "Core proteins", different modelings of the retained signal peptide, orientation of the low-potential heme b, and chirality of the Met ligand to heme c1. This article is part of a Special Issue entitled: Respiratory complex III and related bc complexes.