Determination of Ligand Pathways in Globins: APOLAR TUNNELS VERSUS POLAR GATES*

Although molecular dynamics simulations suggest multiple interior pathways for O₂ entry into and exit from globins, most experiments indicate well defined single pathways. In 2001, we highlighted the effects of large-to-small amino acid replacements on rates for ligand entry and exit onto the three-dimensional structure of sperm whale myoglobin. The resultant map argued strongly for ligand movement through a short channel from the heme iron to solvent that is gated by the distal histidine (His-64(E7)) near the solvent edge of the porphyrin ring. In this work, we have applied the same mutagenesis mapping strategy to the neuronal mini-hemoglobin from Cerebratulus lacteus (CerHb), which has a large internal tunnel from the heme iron to the C-terminal ends of the E and H helices, a direction that is 180° opposite to the E7 channel. Detailed comparisons of the new CerHb map with expanded results for Mb show unambiguously that the dominant (>90%) ligand pathway in CerHb is through the internal tunnel, and the major (>75%) ligand pathway in Mb is through the E7 gate. These results demonstrate that: 1) mutagenesis mapping can identify internal pathways when they exist; 2) molecular dynamics simulations need to be refined to address discrepancies with experimental observations; and 3) alternative pathways have evolved in globins to meet specific physiological demands.     

Escape of a Small Molecule from Inside T4 Lysozyme by Multiple Pathways

The T4 lysozyme L99A mutant is often used as a model system to study small-molecule binding to proteins, but pathways for ligand entry and exit from the buried binding site and the associated protein conformational changes have not been fully resolved. Here, molecular dynamics simulations were employed to model benzene exit from its binding cavity using the weighted ensemble (WE) approach to enhance sampling of low-probability unbinding trajectories. Independent WE simulations revealed four pathways for benzene exit, which correspond to transient tunnels spontaneously formed in previous simulations of apo T4 lysozyme. Thus, benzene unbinding occurs through multiple pathways partially created by intrinsic protein structural fluctuations. Motions of several α-helices and side chains were involved in ligand escape from metastable microstates. WE simulations also provided preliminary estimates of rate constants for each exit pathway. These results complement previous works and provide a semiquantitative characterization of pathway heterogeneity for binding of small molecules to proteins.     

NLDB: a database for 3D protein–ligand interactions in enzymatic reactions

NLDB (Natural Ligand DataBase; URL: http://nldb.hgc.jp) is a database of automatically collected and predicted 3D protein–ligand interactions for the enzymatic reactions of metabolic pathways registered in KEGG. Structural information about these reactions is important for studying the molecular functions of enzymes, however a large number of the 3D interactions are still unknown. Therefore, in order to complement such missing information, we predicted protein–ligand complex structures, and constructed a database of the 3D interactions in reactions. NLDB provides three different types of data resources; the natural complexes are experimentally determined protein–ligand complex structures in PDB, the analog complexes are predicted based on known protein structures in a complex with a similar ligand, and the ab initio complexes are predicted by docking simulations. In addition, NLDB shows the known polymorphisms found in human genome on protein structures. The database has a flexible search function based on various types of keywords, and an enrichment analysis function based on a set of KEGG compound IDs. NLDB will be a valuable resource for experimental biologists studying protein–ligand interactions in specific reactions, and for theoretical researchers wishing to undertake more precise simulations of interactions.     

Molecular dynamics simulation of water in cytochrome c oxidase reveals two water exit pathways and the mechanism of transport

We have examined the network of connected internal cavities in cytochrome c oxidase along which water produced at the catalytic center is removed from the enzyme. Using combination of structural analysis, molecular dynamics simulations, and free energy calculations we have identified two exit pathways that connect the Mg²⁺ ion cavity to the outside of the enzyme. Each pathway has a well-defined bottleneck, which determines the overall rate of water traffic along the exit pathway, and a specific cooperative mechanism of passing it. One of the pathways is going via Arg438/439 (in bovine numbering) toward the Cu_A center, approaching closely its His204_B ligand and Lys171_B residue; and the other is going toward Asp364 and Thr294. Comparison of the pathways among different aa₃-type enzymes shows that they are well conserved. Possible connections of the finding to redox-coupled proton pumping mechanism are discussed. We propose specific mutations near the bottlenecks of the exit pathways that can test some of our hypotheses.     

Ligand-release pathways in the pheromone-binding protein of Bombyx mori

Pheromone-binding proteins (PBP) supply olfactory neuron cells with pheromones by binding the ligands they are tailored for and carrying them to their receptor. The function of a PBP as an efficient carrier requires fast ligand uptake and release. The molecular basis of the ligand-binding mechanism was addressed here for the intriguing case of the PBP of the silk moth Bombyx mori. This PBP completely encapsulates its ligand bombykol without displaying any obvious ligand entrance/exit sites. Here, two opposite dissociation routes were identified as the most likely entrance/exit paths by replica-exchange molecular dynamics, essential dynamics, and force-probe molecular dynamics simulations. One of the paths runs along a flexible front lid; the other along the termini at the back. Calculated forces and energies suggest that both routes are physiologically relevant. The multiplicity of pathways may reduce or tune the entropic barrier for ligand binding.     

Free-energy relationships in ion channels activated by voltage and ligand

Many ion channels are modulated by multiple stimuli, which allow them to integrate a variety of cellular signals and precisely respond to physiological needs. Understanding how these different signaling pathways interact has been a challenge in part because of the complexity of underlying models. In this study, we analyzed the energetic relationships in polymodal ion channels using linkage principles. We first show that in proteins dually modulated by voltage and ligand, the net free-energy change can be obtained by measuring the charge-voltage (Q-V) relationship in zero ligand condition and the ligand binding curve at highly depolarizing membrane voltages. Next, we show that the voltage-dependent changes in ligand occupancy of the protein can be directly obtained by measuring the Q-V curves at multiple ligand concentrations. When a single reference ligand binding curve is available, this relationship allows us to reconstruct ligand binding curves at different voltages. More significantly, we establish that the shift of the Q-V curve between zero and saturating ligand concentration is a direct estimate of the interaction energy between the ligand- and voltage-dependent pathway. These free-energy relationships were tested by numerical simulations of a detailed gating model of the BK channel. Furthermore, as a proof of principle, we estimate the interaction energy between the ligand binding and voltage-dependent pathways for HCN2 channels whose ligand binding curves at various voltages are available. These emerging principles will be useful for high-throughput mutagenesis studies aimed at identifying interaction pathways between various regulatory domains in a polymodal ion channel.     

Long-Timescale Molecular-Dynamics Simulations of the Major Urinary Protein Provide Atomistic Interpretations of the Unusual Thermodynamics of Ligand Binding

The mouse major urinary protein (MUP) has proved to be an intriguing test bed for detailed studies on protein-ligand recognition. NMR, calorimetric, and modeling investigations have revealed that the thermodynamics of ligand binding involve a complex interplay between competing enthalpic and entropic terms. We performed six independent, 1.2 μs molecular-dynamics simulations on MUP—three replicates on the apo-protein, and three on the complex with the pheromone isobutylmethoxypyrazine. Our findings provide the most comprehensive picture to date of the structure and dynamics of MUP, and how they are modulated by ligand binding. The mechanical pathways by which amino acid side chains can transmit information regarding ligand binding to surface loops and either increase or decrease their flexibility (entropy-entropy compensation) are identified. Dewetting of the highly hydrophobic binding cavity is confirmed, and the results reveal an aspect of ligand binding that was not observed in earlier, shorter simulations: bound ligand retains extensive rotational freedom. Both of these features have significant implications for interpretations of the entropic component of binding. More generally, these simulations test the ability of current molecular simulation methods to produce a reliable and reproducible picture of protein dynamics on the microsecond timescale.     

The K-path entrance in cytochrome c oxidase is defined by mutation of E101 and controlled by an adjacent ligand binding domain

Three mutant forms of Rhodobacter sphaeroides cytochrome c oxidase (RsCcO) were created to test for multiple K-path entry sites (E101W), the existence of an “upper ligand site” (M350?W), and the nature and binding specificity of the “lower ligand site” (P315W/E101A) in the region of a crystallographically-defined deoxycholate at the K-path entrance. The effects of inhibitory and stimulatory detergents (dodecyl maltoside and Tween20) on these mutants are presented, as well as competition with other ligands, including the potentially physiologically relevant ligands cholesterol and retinoic acid. Ligands are shown to be able to compete with natural lipids to affect the activity of membrane-bound RsCcO. Results point to a single K-path entrance site at E101, with a single ligand binding pocket proximal to the entrance. The affinity of this pocket for amphipathic ligands is enhanced by removal of the E101 carboxyl and blocked by substituting a tryptophan in this area. A new crystal structure of the E101A mutant of RsCcO is presented that illustrates the structural basis of these results, showing that the loss of the E101 carboxyl creates a more hydrophobic groove consistent with altered ligand affinities.     

Unveiling Prolyl Oligopeptidase Ligand Migration by Comprehensive Computational Techniques

Prolyl oligopeptidase (POP) is a large 80 kDa protease, which cleaves oligopeptides at the C-terminal side of proline residues and constitutes an important pharmaceutical target. Despite the existence of several crystallographic structures, there is an open debate about migration (entrance and exit) pathways for ligands, and their coupling with protein dynamics. Recent studies have shown the capabilities of molecular dynamics and classical force fields in describing spontaneous binding events and nonbiased ligand migration pathways. Due to POP’s size and to the buried nature of its active site, an exhaustive sampling by means of conventional long enough molecular dynamics trajectories is still a nearly impossible task. Such a level of sampling, however, is possible with the breakthrough protein energy landscape exploration technique. Here, we present an exhaustive sampling of POP with a known inhibitor, Z-pro-prolinal. In >3000 trajectories Z-pro-prolinal explores all the accessible surface area, showing multiple entrance events into the large internal cavity through the pore in the β-propeller domain. Moreover, we modeled a natural substrate binding and product release by predicting the entrance of an undecapeptide substrate, followed by manual active site cleavage and nonbiased exit of one of the products (a dipeptide). The product exit shows preference from a flexible 18-amino acid residues loop, pointing to an overall mechanism where entrance and exit occur in different sites.     

Theoretical Investigations of Nitric Oxide Channeling in Mycobacterium tuberculosis Truncated Hemoglobin N

Mycobacterium tuberculosis group I truncated hemoglobin trHbN catalyzes the oxidation of nitric oxide (•NO) to nitrate with a second-order rate constant k ≈ 745 μM⁻¹ s⁻¹ at 23°C (nitric oxide dioxygenase reaction). It was proposed that this high efficiency is associated with the presence of hydrophobic tunnels inside trHbN structure that allow substrate diffusion to the distal heme pocket. In this work, we investigated the mechanisms of •NO diffusion within trHbN tunnels in the context of the nitric oxide dioxygenase reaction using two independent approaches. Molecular dynamics simulations of trHbN were performed in the presence of explicit •NO molecules. Successful •NO diffusion from the bulk solvent to the distal heme pocket was observed in all simulations performed. The simulations revealed that •NO interacts with trHbN at specific surface sites, composed of hydrophobic residues located at tunnel entrances. The entry and the internal diffusion of •NO inside trHbN were performed using the Long, Short, and EH tunnels reported earlier. The Short tunnel was preferentially used by •NO to reach the distal heme pocket. This preference is ascribed to its hydrophobic funnel-shape entrance, covering a large area extending far from the tunnel entrance. This funnel-shape entrance triggers the frequent formation of solvent-excluded cavities capable of hosting up to three •NO molecules, thereby accelerating •NO capture and entry. The importance of hydrophobicity of entrances for •NO capture is highlighted by a comparison with a polar mutant for which residues at entrances were mutated with polar residues. A complete map of •NO diffusion pathways inside trHbN matrix was calculated, and •NO molecules were found to diffuse from Xe cavity to Xe cavity. This scheme was in perfect agreement with the three-dimensional free-energy distribution calculated using implicit ligand sampling. The trajectories showed that •NO significantly alters the dynamics of the key amino acids of Phe⁶²(E15), a residue proposed to act as a gate controlling ligand traffic inside the Long tunnel, and also of Ile¹¹⁹(H11), at the entrance of the Short tunnel. It is noteworthy that •NO diffusion inside trHbN tunnels is much faster than that reported previously for myoglobin. The results presented in this work shed light on the diffusion mechanism of apolar gaseous substrates inside protein matrix.     

Nitric Oxide Dynamics in Truncated Hemoglobin: Docking Sites, Migration Pathways, and Vibrational Spectroscopy from Molecular Dynamics Simulations

Atomistic simulations of nitric oxide (NO) dynamics and migration in the trHbN of Mycobacterium tuberculosis are reported. From extensive molecular dynamics simulations (48 ns in total), the structural and energetic properties of the ligand docking sites in the protein have been characterized and a connectivity network between the ligand docking sites has been built. Several novel migration and exit pathways are found and are analyzed in detail. The interplay between a hydrogen-bonding network involving residues Tyr³³ and Gln⁵⁸ and the bound O₂ ligand is discussed and the role of Phe⁶² residue in ligand migration is examined. It is found that Phe⁶² is directly involved in controlling ligand migration. This is reminiscent of His⁶⁴ in myoglobin, which also plays a central role in CO migration pathways. Finally, infrared spectra of the NO molecule in different ligand docking sites of the protein are calculated. The pocket-specific spectra are typically blue-shifted by 5-10 cm⁻¹, which should be detectable in future spectroscopic experiments.     

Conformational plasticity and dynamic interactions of the N-terminal domain of the chemokine receptor CXCR1

The dynamic interactions between G protein-coupled receptors (GPCRs) and their cognate protein partners are central to several cell signaling pathways. For example, the association of CXC chemokine receptor 1 (CXCR1) with its cognate chemokine, interleukin-8 (IL8 or CXCL8) initiates pathways leading to neutrophil-mediated immune responses. The N-terminal domain of chemokine receptors confers ligand selectivity, but unfortunately the conformational dynamics of this intrinsically disordered region remains unresolved. In this work, we have explored the interaction of CXCR1 with IL8 by microsecond time scale coarse-grain simulations, complemented by atomistic models and NMR chemical shift predictions. We show that the conformational plasticity of the apo-receptor N-terminal domain is restricted upon ligand binding, driving it to an open C-shaped conformation. Importantly, we corroborated the dynamic complex sampled in our simulations against chemical shift perturbations reported by previous NMR studies and show that the trends are similar. Our results indicate that chemical shift perturbation is often not a reporter of residue contacts in such dynamic associations. We believe our results represent a step forward in devising a strategy to understand intrinsically disordered regions in GPCRs and how they acquire functionally important conformational ensembles in dynamic protein-protein interfaces.     

Homology modeling of phosphoryl thymidine kinase of enterohemorrhagic Escherichia coli OH: 157

Enterohemorrhagic Escherichia coli (EHEC) are source of emerging infectious disease in India. Escherichia coli O157:H7 is an EHEC strain showing multiple antibiotic resistances and the cause of infantile diarrhea and hemolytic uremic syndrome worldwide. A novel strategy to counteract multiple antibiotic resistant organisms is to design drugs which specifically target metabolic pathways such as thiamine biosynthetic pathways found exclusively in prokaryotes. Homology modeling was used for model building of a terminal thiamine biosynthesis enzyme phosphoryl thymidine kinase (Thi E) using Geno3D, Swiss Model and Modeller. The best model was selected based on overall stereochemical quality. The potential ligand binding sites in the model were identified by CASTp server. The validated theoretical model of the 3D structure of the thiE protein of E. coli O157:H7 was predicted using a thiamine phosphate pyrophosphatase from Pyrococcus furiosus (PDB ID: 1X13_A) as template. The active pockets of ligand binding sites in the enzyme were identified. In this study, phosphoryl thymidine kinase (thi E), a terminal enzyme in the thiamine biosynthesis pathway in the pathogen has been modeled to be used in future as a potential drug target by the design of suitable inhibitors.     

A combined computational-experimental analyses of selected metabolic enzymes in Pseudomonas species

Comparative genomic analysis has revolutionized our ability to predict the metabolic subsystems that occur in newly sequenced genomes, and to explore the functional roles of the set of genes within each subsystem. These computational predictions can considerably reduce the volume of experimental studies required to assess basic metabolic properties of multiple bacterial species. However, experimental validations are still required to resolve the apparent inconsistencies in the predictions by multiple resources. Here, we present combined computational-experimental analyses on eight completely sequenced Pseudomonas species. Comparative pathway analyses reveal that several pathways within the Pseudomonas species show high plasticity and versatility. Potential bypasses in 11 metabolic pathways were identified. We further confirmed the presence of the enzyme O-acetyl homoserine (thiol) lyase (EC: 2.5.1.49) in P. syringae pv. tomato that revealed inconsistent annotations in KEGG and in the recently published SYSTOMONAS database. These analyses connect and integrate systematic data generation, computational data interpretation, and experimental validation and represent a synergistic and powerful means for conducting biological research.     

Role of substrate recognition in modulating strigolactone receptor selectivity in witchweed

Witchweed, or Striga hermonthica, is a parasitic weed that destroys billions of dollars'' worth of crops globally every year. Its germination is stimulated by strigolactones exuded by its host plants. Despite high sequence, structure, and ligand-binding site conservation across different plant species, one strigolactone receptor in witchweed, ShHTL7, uniquely exhibits a picomolar EC50 for downstream signaling. Previous biochemical and structural analyses have hypothesized that this unique ligand sensitivity can be attributed to a large binding pocket volume in ShHTL7 resulting in enhanced ability to bind substrates, but additional structural details of the substrate-binding process would help explain its role in modulating the ligand selectivity. Using long-timescale molecular dynamics simulations, we demonstrate that mutations at the entrance of the binding pocket facilitate a more direct ligand-binding pathway to ShHTL7, whereas hydrophobicity at the binding pocket entrance results in a stable “anchored” state. We also demonstrate that several residues on the D-loop of AtD14 stabilize catalytically inactive conformations. Finally, we show that strigolactone selectivity is not modulated by binding pocket volume. Our results indicate that while ligand binding is not the sole modulator of strigolactone receptor selectivity, it is a significant contributing factor. These results can be used to inform the design of selective antagonists for strigolactone receptors in witchweed.     

Signaling pathways of PDZ2 domain: a molecular dynamics interaction correlation analysis

PDZ domains are found in many signaling proteins. One of their functions is to provide scaffolds for forming membrane-associated protein complexes by binding to the carboxyl termini of their partners. PDZ domains are thought also to play a signal transduction role by propagating the information that binding has occurred to remote sites. In this study, a molecular dynamics (MD) simulation-based approach, referred to as an interaction correlation analysis, is applied to the PDZ2 domain to identify the possible signal transduction pathways. A residue correlation matrix is constructed from the interaction energy correlations between all residue pairs obtained from the MD simulations. Two continuous interaction pathways, starting at the ligand binding pocket, are identified by a hierarchical clustering analysis of the residue correlation matrix. One pathway is mainly localized at the N-terminal side of helix alpha1 and the adjacent C-terminus of loop beta1-beta2. The other pathway is perpendicular to the central beta-sheet and extends toward the side of PDZ2 domain opposite to the ligand binding pocket. The results complement previous studies based on multiple sequence analysis, NMR, and MD simulations. Importantly, they reveal the energetic origin of the long-range coupling. The PDZ2 results, as well as the earlier rhodopsin analysis, show that the interaction correlation analysis is a robust approach for determining pathways of intramolecular signal transduction.     

Biochemical characterization of recombinant nucleoside hydrolase from Mycobacterium tuberculosis H37Rv

Tuberculosis (TB) is a major global health threat. There is a need for the development of more efficient drugs for the sterilization of the disease’s causative agent, Mycobacterium tuberculosis (MTB). A more comprehensive understanding of the bacilli’s nucleotide metabolic pathways could aid in the development of new anti-mycobacterial drugs. Here we describe expression and purification of recombinant iunH-encoded nucleoside hydrolase from MTB (MtIAGU-NH). Glutaraldehyde cross-linking results indicate that MtIAGU-NH predominates as a monomer, presenting varied oligomeric states depending upon binding of ligands. Steady-state kinetics results show that MtIAGU-NH has broad substrate specificity, accepting inosine, adenosine, guanosine, and uridine as substrates. Inosine and adenosine displayed positive homotropic cooperativity kinetics, whereas guanosine and uridine displayed hyperbolic saturation curves. Measurements of kinetics of ribose binding to MtIAGU-NH by fluorescence spectroscopy suggest two pre-existing forms of enzyme prior to ligand association. The intracellular concentrations of inosine, uridine, hypoxanthine, and uracil were determined and thermodynamic parameters estimated. Thermodynamic activation parameters (E_a, ΔG^#, ΔS^#, ΔH^#) for MtIAGU-NH-catalyzed chemical reaction are presented. Results from mass spectrometry, isothermal titration calorimetry (ITC), pH-rate profile experiment, multiple sequence alignment, and molecular docking experiments are also presented. These data should contribute to our understanding of the biological role played by MtIAGU-NH.     

Differential activities of the Drosophila JAK/STAT pathway ligands Upd, Upd2 and Upd3

JAK/STAT signalling in vertebrates is activated by multiple cytokines and growth factors. By contrast, the Drosophila genome encodes for only three related JAK/STAT ligands, Upd, Upd2 and Upd3. Identifying the differences between these three ligands will ultimately lead to a greater understanding of this disease-related signalling pathway and its roles in development. Here, we describe the analysis of the least well characterised of the Upd-like ligands, Upd3. We show that in tissue culture-based assays Upd3-GFP is secreted from cells and appears to interact with the extracellular matrix (ECM) in a similar manner to Upd, while still non-autonomously activating JAK/STAT signalling. Quantification of each of the Upd-like ligands in conditioned media has allowed us to determine the activity of equal amounts of each ligand on JAK/STAT ex vivo and reveals that Upd is the most potent ligand in this system. Finally, investigations into the effects of ectopic expression of Upd3 in vivo have confirmed its ability to activate pathway signalling at long-distance.