Computational modelling of the open-state Kv 1.5 ion channel block by bupivacaine

Binding of R(+)-bupivacaine to open-state homology models of the mammalian K(v)1.5 membrane ion channel is studied using automated docking and molecular dynamics (MD) methods. Homology models of K(v)1.5 are built using the 3D structures of the KcsA and MthK channels as a template. The packing of transmembrane (TM) alpha-helices in the KcsA structure corresponds to a closed channel state. Opening of the channel may be reached by a conformational transition yielding a bent structure of the internal S6 helices. Our first model of the K(v) open state involves a PVP-type of bending hinge in the internal helices, while the second model corresponds to a Gly-type of bending hinge as found in the MthK channel. Ligand binding to these models is probed using the common local anaesthetic bupivacaine, where blocker binding from the intracellular side of the channel is considered. Conformational properties and partial atomic charges of bupivacaine are determined from quantum mechanical HF/6-31G* calculations with inclusion of solvent effects. The automated docking and MD calculations for the PVP-bend model predict that bupivacaine could bind either in the central cavity or in the PVP region of the channel pore. Linear interaction energy (LIE) estimates of the binding free energies for bupivacaine predict strongest binding to the PVP region. Surprisingly, no binding is predicted for the Gly-bend model. These results are discussed in light of electrophysiological data which show that the K(v)1.5 channel is unable to close when bupivacaine is bound.     

Molecular modeling study of the editing active site of Escherichia coli leucyl-tRNA synthetase: two amino acid binding sites in the editing domain

Aminoacyl-tRNA synthetases (aaRSs) strictly discriminate their cognate amino acids. Some aaRSs accomplish this via proofreading and editing mechanisms. Mursinna and coworkers recently reported that substituting a highly conserved threonine (T252) with an alanine within the editing domain of Escherichia coli leucyl-tRNA synthetase (LeuRS) caused LeuRS to cleave its cognate aminoacylated leucine from tRNA(Leu) (Mursinna et al., Biochemistry 2001;40:5376-5381). To achieve atomic level insight into the role of T252 in LeuRS and the editing reaction of aaRSs, a series of molecular modeling studies including homology modeling and automated docking simulations were carried out. A 3D structure of E. coli LeuRS was constructed via homology modeling using the X-ray structure of Thermus thermophilus LeuRS as a template because the E. coli LeuRS structure is not available from X-ray or NMR studies. However, both the X-ray T. thermophilus and homology-modeled E. coli structures were used in our studies. Amino acid binding sites in the proposed editing domain, which is also called the connective polypeptide 1 (CP1) domain, were investigated by automated docking studies. The root mean square deviation (RMSD) for backbone atoms between the X-ray and homology-modeled structures was 1.18 A overall and 0.60 A for the editing (CP1) domain. Automated docking studies of a leucine ligand into the editing domain were performed for both structures: homology structure of E. coli LeuRS and X-ray structure of T. thermophilus LeuRS for comparison. The results of the docking studies suggested that there are two possible amino acid binding sites in the CP1 domain for both proteins. The first site lies near a threonine-rich region that includes the highly conserved T252 residue, which is important for amino acid discrimination. The second site is located in a flexible loop region surrounded by residues E292, A293, M295, A296, and M298. The important T252 residue is at the bottom of the first binding pocket.     

Homology model and targeted mutagenesis identify critical residues for arachidonic acid inhibition of Kv4 channels

Polyunsaturated fatty acids such as arachidonic acid (AA) exhibit inhibitory modulation of Kv4 potassium channels. Molecular docking approaches using a Kv4.2 homology model predicted a membrane-embedded binding pocket for AA comprised of the S4-S5 linker on one subunit and several hydrophobic residues within S3, S5 and S6 from an adjacent subunit. The pocket is conserved among Kv4 channels. We tested the hypothesis that modulatory effects of AA on Kv4.2/KChIP channels require access to this site. Targeted mutation of a polar residue (K318) and a nonpolar residue (G314) within the S4-S5 linker as well as a nonpolar residue in S3 (V261) significantly impaired the effects of AA on K⁺ currents in Xenopus oocytes. These residues may be important in stabilizing (K318) or regulating access to (V261, G314) the negatively charged carboxylate moiety on the fatty acid. Structural specificity was supported by the lack of disruption of AA effects observed with mutations at residues located near, but not within the predicted binding pocket. Furthermore, we found that the crystal structure of the related Kv1.2/2.1 chimera lacks the structural features present in the proposed AA docking site of Kv4.2 and the Kv1.2/2.1 K⁺ currents were unaffected by AA. We simulated the mutagenic substitutions in our Kv4.2 model to demonstrate how specific mutations may disrupt the putative AA binding pocket. We conclude that AA inhibits Kv4 channel currents and facilitates current decay by binding within a hydrophobic pocket in the channel in which K318 within the S4-S5 linker is a critical residue for AA interaction.     

Three consecutive arginines are important for the mycobacterial peptide deformylase enzyme activity

Genes encoding the peptide deformylase enzyme (def) are present in all eubacteria and are involved in the deformylation of the N-formyl group of newly synthesized polypeptides during protein synthesis. We compared the amino acid sequences of this enzyme in different mycobacterial species and found that they are highly conserved (76% homology with 62% identity); however, when this comparison was extended to other eubacterial homologs, it emerged that the mycobacterial proteins have an insertion region containing three consecutive arginine residues (residues 77-79 in Mycobacterium tuberculosis peptide deformylase (mPDF)). Here, we demonstrate that these three arginines are important for the activity of mPDF. Circular dichroism studies of wild-type mPDF and of mPDF containing individual conservative substitutions (R77K, R78K, or R79K) or combined substitutions incorporated into a triple mutant (R77K/R78K/R79K) indicate that such mutations cause mPDF to undergo structural alterations. Molecular modeling of mPDF suggests that the three arginines are distal to the active site. Molecular dynamics simulations of wild-type and mutant mPDF structures indicate that the arginines may be involved in the stabilization of substrate binding pocket residues for their proper interaction with peptide(s). Treatment with 5'-phosphothiorate-modified antisense oligodeoxyribonucleotides directed against different regions of def from M. tuberculosis inhibits growth of Mycobacterium smegmatis in culture. Taken together, these results hold out the possibility of future design of novel mycobacteria-specific PDF inhibitors.     

The NMR solution structure and function of RPA3313: a putative ribosomal transport protein from Rhodopseudomonas palustris

Protein function elucidation often relies heavily on amino acid sequence analysis and other bioinformatics approaches. The reliance is extended to structure homology modeling for ligand docking and protein–protein interaction mapping. However, sequence analysis of RPA3313 exposes a large, unannotated class of hypothetical proteins mostly from the Rhizobiales order. In the absence of sequence and structure information, further functional elucidation of this class of proteins has been significantly hindered. A high quality NMR structure of RPA3313 reveals that the protein forms a novel split ββαβ fold with a conserved ligand binding pocket between the first β‐strand and the N‐terminus of the α‐helix. Conserved residue analysis and protein–protein interaction prediction analyses reveal multiple protein binding sites and conserved functional residues. Results of a mass spectrometry proteomic analysis strongly point toward interaction with the ribosome and its subunits. The combined structural and proteomic analyses suggest that RPA3313 by itself or in a larger complex may assist in the transportation of substrates to or from the ribosome for further processing. Proteins 2016; 85:93–102. © 2016 Wiley Periodicals, Inc.     

Uncoupling the ATPase activity of the N-ethylmaleimide sensitive factor (NSF) from 20S complex disassembly.

The N-ethylmaleimide sensitive factor (NSF) plays a critical role in intracellular trafficking by disassembling soluble NSF attachment protein receptor (SNARE ) complexes. The NSF protomer consists of three domains (NSF-N, NSF-D1, and NSF-D2). Two domains (NSF-D1 and NSF-D2) contain a conserved approximately 230 amino acid cassette, which includes a distinctive motif termed the second region of homology (SRH) common to all ATPases associated with various cellular activities (AAA). In hexameric NSF, several SRH residues become trans elements of the ATP binding pocket. Mutation of two conserved arginine residues in the NSF-D1 SRH (R385A and R388A) did not effect basal or soluble NSF attachment protein (SNAP)-stimulated ATPase activity; however, neither mutant underwent ATP-dependent release from SNAP-SNARE complexes. A trans element of the NSF-D2 ATP binding site (K631) has been proposed to limit the ATPase activity of NSF-D2, but a K631D mutant retained wild-type activity. A mutation of the equivalent residue in NSF-D1 (D359K) also did not affect nucleotide hydrolysis activity but did limit NSF release from SNAP-SNARE complexes. These trans elements of the NSF-D1 ATP binding site (R385, R388, and D359) are not required for nucleotide hydrolysis but are important as nucleotide-state sensors. NSF-N mediates binding to the SNAP-SNARE complex. To identify the structural features required for binding, three conserved residues (R67, S73, and Q76) on the surface of NSF-N were mutated. R67E completely eliminated binding, while S73R and Q76E showed limited effect. This suggests that the surface important for SNAP binding site lies in the cleft between the NSF-N subdomains adjacent to a conserved, positively charged surface.     

Model of three-dimensional structure of vitamin D receptor and its binding mechanism with 1alpha,25-dihydroxyvitamin D(3).

Comparative modeling of the vitamin D receptor three-dimensional structure and computational docking of 1alpha,25-dihydroxyvitamin D(3) into the putative binding pocket of the two deletion mutant receptors: (207-423) and (120-422, Delta [164-207]) are reported and evaluated in the context of extensive mutagenic analysis and crystal structure of holo hVDR deletion protein published recently. The obtained molecular model agrees well with the experimentally determined structure. Six different conformers of 1alpha,25-dihydroxyvitamin D(3) were used to study flexible docking to the receptor. On the basis of values of conformational energy of various complexes and their consistency with functional activity, it appears that 1alpha,25-dihydroxyvitamin D(3) binds the receptor in its 6-s-trans form. The two lowest energy complexes obtained from docking the hormone into the deletion protein (207-423) differ in conformation of ring A and orientation of the ligand molecule in the VDR pocket. 1alpha,25-Dihydroxyvitamin D(3) possessing the A-ring conformation with axially oriented 1alpha-hydroxy group binds receptor with its 25-hydroxy substituent oriented toward the center of the receptor cavity, whereas ligand possessing equatorial conformation of 1alpha-hydroxy enters the pocket with A ring directed inward. The latter conformation and orientation of the ligand is consistent with the crystal structure of hVDR deletion mutant (118-425, Delta [165-215]). The lattice model of rVDR (120-422, Delta [164-207]) shows excellent agreement with the crystal structure of the hVDR mutant. The complex obtained from docking the hormone into the receptor has lower energy than complexes for which homology modeling was used. Thus, a simple model of vitamin D receptor with the first two helices deleted can be potentially useful for designing a general structure of ligand, whereas the advanced lattice model is suitable for examining binding sites in the pocket.     

Exploring the binding pocket for pyridopyrimidine ligands at the CCK1 receptor by molecular docking

Pyridopyrimidine-based analogues are among the most highly potent and selective antagonists of cholecystokinin receptor subtype-1 (CCK1R) described to date. To better understand the structural and chemical features responsible for the recognition mechanism, and to explore the binding pocket of these compounds, we performed automated molecular docking using GOLD2.2 software on some derivatives with structural diversity, and propose a putative binding conformation for each compound. The docking protocol was guided by the key role of the Asn333 residue, as revealed by site directed mutagenesis studies. The results suggest two putative binding modes located in the same pocket. Both are characterized by interaction with the main residues revealed by experiment, Asn333 and Arg336, and differ in the spatial position of the Boc-Trp moiety of these compounds. Hydrophobic contacts with residues Thr117, Phe107, Ile352 and Ile329 are also in agreement with experimental data. Despite the poor correlation obtained between the estimated binding energies and the experimental activity, the proposed models allow us to suggest a plausible explanation of the observed binding data in accordance with chemical characteristics of the compounds, and also to explain the observed diastereoselectivity of this family of antagonists towards CCK1R. The most reasonable selected binding conformations could be the starting point for future studies. Figure Superimposition of the two putative binding conformations revealed by molecular docking for pyridopyrimidine-based CCK1 antagonists     

A low molecular weight agonist signals by binding to the transmembrane domain of thyroid-stimulating hormone receptor (TSHR) and luteinizing hormone/chorionic gonadotropin receptor (LHCGR)

Many cognate low molecular weight (LMW) agonists bind to seven transmembrane-spanning receptors within their transmembrane helices (TMHs). The thienopyrimidine org41841 was identified previously as an agonist for the luteinizing hormone/chorionic gonadotropin receptor (LHCGR) and suggested to bind within its TMHs because it did not compete for LH binding to the LHCGR ectodomain. Because of its high homology with LHCGR, we predicted that thyroid-stimulating hormone receptor (TSHR) might be activated by org41841 also. We show that org41841 is a partial agonist for TSHR but with lower potency than for LHCGR. Analysis of three-dimensional molecular models of TSHR and LHCGR predicted a binding pocket for org41841 in common clefts between TMHs 3, 4, 5, 6, and 7 and extracellular loop 2 in both receptors. Evidence for this binding pocket was obtained in signaling studies with chimeric receptors that exhibited improved responses to org41841. Furthermore, a key receptor-ligand interaction between the highly conserved negatively charged E3.37 and the amino group of org41841 predicted by docking of the ligand into the three-dimensional TSHR model was experimentally confirmed. These findings provide the first evidence that, in contrast to the ectodomain binding of cognate ligands, a LMW agonist can bind to and activate glycoprotein hormone receptors via interaction with their transmembrane domain.     

The structure of a resuscitation-promoting factor domain from Mycobacterium tuberculosis shows homology to lysozymes

Resuscitation-promoting factor (RPF) proteins reactivate stationary-phase cultures of (G+C)-rich Gram-positive bacteria including the causative agent of tuberculosis, Mycobacterium tuberculosis. We report the solution structure of the RPF domain from M. tuberculosis Rv1009 (RpfB) solved by heteronuclear multidimensional NMR. Structural homology with various glycoside hydrolases suggested that RpfB cleaved oligosaccharides. Biochemical studies indicate that a conserved active site glutamate is important for resuscitation activity. These data, as well as the presence of a clear binding pocket for a large molecule, indicate that oligosaccharide cleavage is probably the signal for revival from dormancy.     

Identification of hotspot regions of MurB oxidoreductase enzyme using homology modeling,molecular dynamics and molecular docking techniques

Despite the availability of effective chemotherapy and a moderately protective vaccine, new anti-tuberculosis agents are urgently needed to decrease the global incidence of tuberculosis (TB) disease. The MurB gene belongs to the bacterial cell wall biosynthesis pathway and is an essential drug target in Mycobacterium tuberculosis (Mtb) that has no mammalian counterparts. Here, we present an integrated approach involving homology modeling, molecular dynamics and molecular docking studies on Mtb-MurB oxidoreductase enzyme. A homology model of Mtb-MurB enzyme was built for the first time in order to carry out structure-based inhibitor design. The accuracy of the model was validated using different techniques. The molecular docking study on this enzyme was undertaken using different classes of well known MurB inhibitors. Estimation of binding free energy by docking analysis indicated the importance of Tyr155, Arg156, Ser237, Asn241 and His304 residues within the Mtb-MurB binding pocket. Our computational analysis is in good agreement with experimental results of site-directed mutagenesis. The present study should therefore play a guiding role in the experimental design of Mtb-MurB inhibitors for in vitro/in vivo analysis.     

Design,synthesis, biological evaluation and molecular modelling studies of novel quinoline derivatives against Mycobacterium tuberculosis

We herein describe the synthesis and antimycobacterial activity of a series of 27 different derivatives of 3-benzyl-6-bromo-2-methoxy-quinolines and amides of 2-[(6-bromo-2-methoxy-quinolin-3-yl)-phenyl-methyl]-malonic acid monomethyl ester. The antimycobacterial activity of these compounds was evaluated in vitro against Mycobacterium tuberculosis H37Rv for nine consecutive days upon a fixed concentration (6.25 μg/mL) at day one in Bactec assay and compared to untreated TB cell culture as well as one with isoniazide treated counterpart, under identical experimental conditions. The compounds 3, 8, 17 and 18 have shown 92–100% growth inhibition of mycobacterial activity, with minimum inhibitory concentration (MIC) of 6.25 μg/mL. Based on our molecular modelling and docking studies on well-known diarylquinoline antitubercular drug R207910, the presence of phenyl, naphthyl and halogen moieties seem critical. Comparison of docking studies on different stereoisomers of R207910 as well as compounds from our data set, suggests importance of electrostatic interactions. Further structural analysis of docking studies on our compounds suggests attractive starting point to find new lead compounds with potential improvements.     

Computational modelling of the open-state Kv1.5 ion channel block by bupivacaine

Binding of R(+)-bupivacaine to open-state homology models of the mammalian K_v1.5 membrane ion channel is studied using automated docking and molecular dynamics (MD) methods. Homology models of K_v1.5 are built using the 3D structures of the KcsA and MthK channels as a template. The packing of transmembrane (TM) α-helices in the KcsA structure corresponds to a closed channel state. Opening of the channel may be reached by a conformational transition yielding a bent structure of the internal S6 helices. Our first model of the K_v open state involves a PVP-type of bending hinge in the internal helices, while the second model corresponds to a Gly-type of bending hinge as found in the MthK channel. Ligand binding to these models is probed using the common local anaesthetic bupivacaine, where blocker binding from the intracellular side of the channel is considered. Conformational properties and partial atomic charges of bupivacaine are determined from quantum mechanical HF/6-31G* calculations with inclusion of solvent effects. The automated docking and MD calculations for the PVP-bend model predict that bupivacaine could bind either in the central cavity or in the PVP region of the channel pore. Linear interaction energy (LIE) estimates of the binding free energies for bupivacaine predict strongest binding to the PVP region. Surprisingly, no binding is predicted for the Gly-bend model. These results are discussed in light of electrophysiological data which show that the K_v1.5 channel is unable to close when bupivacaine is bound.     

Role of alphaPhe-291 residue in the phosphate-binding subdomain of catalytic sites of Escherichia coli ATP synthase

The role of αPhe-291 residue in phosphate binding by Escherichia coli F₁F₀-ATP synthase was examined. X-ray structures of bovine mitochondrial enzyme suggest that this residue resides in close proximity to the conserved βR246 residue. Herein, we show that mutations αF291D and αF291E in E. coli reduce the ATPase activity of F₁F₀ membranes by 350-fold. Yet, significant oxidative phosphorylation activity is retained. In contrast to wild-type, ATPase activities of mutants were not inhibited by MgADP-azide, MgADP-fluoroaluminate, or MgADP-fluoroscandium. Whereas, 7-chloro-4-nitrobenzo-2-oxa-1,3-diazole (NBD-Cl) inhibited wild-type ATPase essentially completely, ATPase in mutants was inhibited maximally by ∼75%, although reaction still occurred at residue βTyr-297, proximal to αPhe-291 in the phosphate-binding pocket. Inhibition characteristics supported the conclusion that NBD-Cl reacts in βE (empty) catalytic sites, as shown previously by X-ray structure analysis. Phosphate protected against NBD-Cl inhibition in wild-type but not in mutants. In addition, our data suggest that the interaction of αPhe-291 with phosphate during ATP hydrolysis or synthesis may be distinct.     

Homology modeling and atomic level binding study of Leishmania MAPK with inhibitors

The current therapy for leishmaniasis is not sufficient and it has two severe drawbacks, host-toxicity and drug resistance. The substantial knowledge of parasite biology is not yet translating into novel drugs for leishmaniasis. Based on this observation, a 3D structural model of Leishmania mitogen-activated protein kinase (MAPK) homologue has been developed, for the first time, by homology modeling and molecular dynamics simulation techniques. The model provided clear insight in its structure features, i.e. ATP binding pocket, phosphorylation lip, and common docking site. Sequence-structure homology recognition identified Leishmania CRK3 (LCRK3) as a distant member of the MAPK superfamily. Multiple sequence alignment and 3D structure model provided the putative ATP binding pocket of Leishmania with respect to human ERK2 and LCRK3. This analysis was helpful in identifying the binding sites and molecular function of the Leishmania specific MAPK homologue. Molecular docking study was performed on this 3D structural model, using different classes of competitive ATP inhibitors of LCRK3, to check whether they exhibit affinity and could be identified as Leishmania MAPK specific inhibitors. It is well known that MAP kinases are extracellular signal regulated kinases ERK1 and ERK2, which are components of the Ras-MAPK signal transduction pathway which is complexed with HDAC4 protein, and their inhibition is of significant therapeutic interest in cancer biology. In order to understand the mechanism of action, docking of indirubin class of molecules to the active site of histone deacetylase 4 (HDAC4) protein is performed, and the binding affinity of the protein-ligand interaction was computed. The new structural insights obtained from this study are all consistent with the available experimental data, suggesting that the homology model of the Leishmania MAPK and its ligand interaction modes are reasonable. Further the comparative molecular electrostatic potential and cavity depth analysis of Leishmania MAPK and human ERK2 suggested several important differences in its ATP binding pocket. Such differences could be exploited in the future for designing Leishmania specific MAPK inhibitors.     

Three-dimensional structure and function study on the active region in the extracellular ligand-binding domain of human IL-6 receptor

In this study the three-dimensional (3-D) model of the ligand-binding domain (V106-P322) of human interleukin-6 receptor (hIL-6 R) was constructed by computer-guided homology modeling technique using the crystal structure of the ligand-binding domain (K52-L251) of human growth hormone receptor (hGHR) as templet. Furthermore, the active binding region of the 3-D model of hIL-6R with the ligand (hIL-6) was predicted. In light of the structural characteristics of the active region, a hydrophobic pocket shielded by two hydrophilic residues (E115 and E505) of the region was identified by a combination of molecular modelling and the site-directed or double-site mutation of the twelve crucial residues in the ligand-binding domain of hIL-6R (V106-P322). We observed and analyzed the effects of these mutants on the spatial conformation of the pocket-like region of hIL-6 R. The results indicated that any site-directed mutation of the five Cys residues (four conservative Cys residues: Cys121, Cys132, Cys165, Cys176; near membrane Cys residue: Cys193) or each double-site mutation of the five residues in WSEWS motif of hIL-6R (V106-P322) makes the corresponding spatial conformation of the pocket region block the linkage between hIL-6 R and hIL-6. However, the influence of the site-directed mutation of Cys211 and Cys277 individually on the conformation of the pocket region benefits the interaction between hIL-6R and hIL-6. Our study suggests that the predicted hydrophobic pocket in the 3-D model of hIL-6R (V106-P322) is the critical molecular basis for the binding of hIL-6R with its ligand, and the active pocket may be used as a target for designing small hIL-6R-inhibiting molecules in our further study.     

Molecular modeling aided design of nicotinic acid receptor GPR109A agonists

A homology model of the nicotinic acid receptor GPR109A was constructed based on the X-ray crystal structure of bovine rhodopsin. An HTS hit was docked into the homology model. Characterization of the binding pocket by a grid-based surface calculation of the docking model suggested that a larger hydrophobic body plus a polar tail would improve interaction between the ligand and the receptor. The designed compounds were synthesized, and showed significantly improved binding affinity and activation of GPR109A.     

Proteomic identification of M. tuberculosis protein kinase substrates: PknB recruits GarA, a FHA domain-containing protein, through activation loop-mediated interactions

Genes for functional Ser/Thr protein kinases (STPKs) are ubiquitous in prokaryotic genomes, but little is known about their physiological substrates and their actual involvement in bacterial signal transduction pathways. We report here the identification of GarA (Rv1827), a Forkhead-associated (FHA) domain-containing protein, as a putative physiological substrate of PknB, an essential Ser/Thr protein kinase from Mycobacterium tuberculosis. Using a global proteomic approach, GarA was found to be the best detectable substrate of the PknB catalytic domain in non-denatured whole-cell protein extracts from M. tuberculosis and the saprophyte Mycobacterium smegmatis. Enzymological and binding studies of the recombinant proteins demonstrate that docking interactions between the activation loop of PknB and the C-terminal FHA domain of GarA are required to enable efficient phosphorylation at a single N-terminal threonine residue, Thr22, of the substrate. The predicted amino acid sequence of the garA gene, including both the N-terminal phosphorylation motif and the FHA domain, is strongly conserved in mycobacteria and other related actinomycetes, suggesting a functional role of GarA in putative STPK-mediated signal transduction pathways. The ensuing model of PknB-GarA interactions suggests a substrate recruitment mechanism that might apply to other mycobacterial kinases bearing multiple phosphorylation sites in their activation loops.