Contribution of aromatic-aromatic interactions to the anomalous pK(a) of tyrosine-9 and the C-terminal dynamics of glutathione S-transferase A1-1

Most cytosolic glutathione S-transferases (GSTs) exploit a hydrogen bond between an active site Tyr and the bound glutathione (GSH) cofactor to lower the pK(a) of the GSH and generate the nucleophilic thiolate anion, GS(-). In human (hGSTA1-1) and rat (rGSTA1-1) homologues, the active site Tyr-9 has a low pK(a) of 8.1-8.3, for which the functional significance is unknown. Crystal structures of GSTA1-1 suggest that weakly polar interactions between the electropositive ring edge of Phe-10 and the pi-cloud of Tyr-9, in the apoenzyme, could stabilize the tyrosinate anion and also modulate the pK(a) of GSH. Upon binding a product GSH conjugate, Phe-10 moves away from Tyr-9, allowing the highly dynamic C-terminus to "close" over the active site. To explore the role of Phe-10 in modulating the Tyr-9 pK(a) and in ligand binding, rGSTA1-1 mutants F10Y, F10L, and F10A were characterized. The pK(a)s of Tyr-9 in the apoenzymes were 8.2 +/- 0.2, 8.7 +/- 0.2, and 9.3 +/- 0.1, respectively, for F10Y, F10L, and F10A, compared to 8.3 +/- 0.2 for the "wild type". The experimentally determined pK(a)s qualitatively paralleled the energies required to remove a proton predicted by ab initio calculations using model compounds constrained to the coordinates of rGSTA1-1. The pK(a) of GSH in the binary complex was significantly less affected by these substitutions. In contrast, F220I and F220Y C-terminal mutations caused the pK(a) of Tyr-9 to decrease modestly. For the binary complex with S-hexyl-GSH, which induces the "closed" conformation, Tyr-9 retains a low pK(a) and the Phe-10 substitutions have significant effects. Presumably, Phe-10 plays a critical structural role in stabilizing the closed conformation. The mutations F10L and F10A also slowed the rate of GSH conjugate binding by 10-20-fold, as measured by stopped-flow fluorescence. The effects of Phe-10 substitution were large for both steps of the biphasic binding reaction, suggesting the importance of aromatic interactions throughout the reaction coordinate. A unified view of the C-terminal dynamics of GSTA1-1 is discussed, which emphasizes the coupling between Tyr-9 ionization, active site solvation, and C-terminal dynamics.     

Vaccinia topoisomerase mutants illuminate conformational changes during closure of the protein clamp and assembly of a functional active site.

We present a mutational analysis of vaccinia topoisomerase that highlights the contributions of five residues in the catalytic domain (Phe-88 and Phe-101 in helix alpha1, Ser-204 in alpha5, and Lys-220 and Asn-228 in alpha6) to the DNA binding and transesterification steps. When augmented by structural information from exemplary type IB topoisomerases and tyrosine recombinases in different functional states, the results suggest how closure of the protein clamp around duplex DNA and assembly of a functional active site might be orchestrated by internal conformational changes in the catalytic domain. Lys-220 is a constituent of the active site, and a positive charge at this position is required for optimal DNA cleavage. Ser-204 and Asn-228 appear not to be directly involved in reaction chemistry at the scissile phosphodiester. We propose that (i) Asn-228 recruits the Tyr-274 nucleophile to the active site by forming a hydrogen bond to the main chain of the tyrosine-containing alpha8 helix and that (ii) contacts between Ser-204 and the DNA backbone upstream of the cleavage site trigger a separate conformational change required for active site assembly. Mutations of Phe-88 and Phe-101 affect DNA binding, most likely at the clamp closure step, which we posit to entail a distortion of helix alpha1.     

Site-directed mutagenesis in haemoglobin. Functional role of tyrosine-42(C7) alpha at the alpha 1-beta 2 interface

To clarify the functional role of Tyr-42(C7) alpha, which forms a hydrogen bond with Asp-99(G1) beta at the alpha 1-beta 2 interface of human deoxyhaemoglobin, we engineered two artificial mutant haemoglobins (Hb), in which Tyr-42 alpha was replaced by Phe (Hb Phe-42 alpha) or His (Hb His-42 alpha), and investigated their oxygen binding properties together with structural consequences of the mutations by using various spectroscopic probes. Like most of the natural Asp-99 beta mutants, Hb Phe-42 alpha showed a markedly increased oxygen affinity, a reduced Bohr effect and diminished co-operativity. Structural probes such as ultraviolet-region derivative and oxy-minus-deoxy difference spectra, resonance Raman scattering and proton nuclear magnetic resonance spectra indicate that, in Hb Phe-42 alpha, the deoxy T quaternary structure is highly destabilized and the strain imposed on the Fe-N epsilon (proximal His) bond is released, stabilizing the oxy tertiary structure. In contrast with Hb Phe-42 alpha, Hb His-42 alpha showed an intermediately impaired function and only moderate destabilization of the T-state, which can be explained by the formation of a new, weak hydrogen bond between His-42 alpha and Asp-99 beta. Spectroscopic data were consistent with this assumption. The present study proves that the hydrogen bond between Tyr-42 alpha and Asp-99 beta plays a key role in stabilizing the deoxy T structure and consequently in co-operative oxygen binding.     

Comparison of the binding affinity of chlorogenic acid with two serum albumins

Chlorogenic acid (CA) is a well-known ester of caffeic acid present in some food. It is also an active component in traditional Chinese medicines which are used to treat various diseases, but the molecular basis of CA is not clear. In the present work, the proton selective relaxation rate and the affinity index were used to investigate the interaction of CA with human serum albumin and bovine serum albumin under the same buffer conditions. The results indicated that the binding affinity of chlorogenic acid to BSA was stronger than that to HSA. The binding site of the ligand-protein complex was elucidated by molecular docking, and the specific interaction was observed from those hydrogen bonds formed by the ligand and active residues. Using a combination of TR-NOE detection, the optimal ligand conformation was illustrated. Further conformational analysis of the complex revealed that the ability of hydrogen bond formation by polar side chain residues in the binding site of BSA might contribute to the greater binding affinity. The results provide a better understanding of CA binding and should contribute towards the design of modifications of CA for therapeutic purposes.     

Computational identification of novel entry inhibitor scaffolds mimicking primary receptor CD4 of HIV-1 gp120

Virtual screening of novel entry inhibitor scaffolds mimicking primary receptor CD4 of HIV-1 gp120 was carried out in conjunction with evaluation of their potential inhibitory activity by molecular modeling. To do this, pharmacophore models presenting different sets of the hotspots of cellular receptor CD4 for its interaction with gp120 were generated. These models were used as the templates for identification of CD4-mimetic candidates by the pepMMsMIMIC screening platform. Complexes of these candidates with gp120 were built by high-throughput ligand docking and their stability was estimated by molecular dynamics simulations and binding free energy calculations. As a result, five top hits that exhibited strong attachment to the two well-conserved hotspots of the gp120 CD4-binding site were selected for the final analysis. In analogy to CD4, the identified compounds make hydrogen bonds with Asp-368_gp120 and multiple van der Waals contacts with the gp120 residues that bind to Phe-43_CD4, resulting in destruction of the critical interactions of gp120 with Phe-43_CD4 and Arg-59_CD4. The complexes of the CD4-mimetic candidates with gp120 show relative conformational stability within the molecular dynamics simulations and expose the high percentage occupancies of intermolecular hydrogen bonds, in line with the data on the binding free energy calculations. In light of these findings, the identified compounds are considered as good scaffolds for the development of new functional antagonists of viral entry with broad HIV-1 neutralization.     

Roles of Arg-97 and Phe-113 in regulation of distal ligand binding to heme in the sensor domain of Ec DOS protein. Resonance Raman and mutation study

The direct oxygen sensor protein isolated from Escherichia coli (Ec DOS) is a heme-based signal transducer protein responsible for phosphodiesterase (PDE) activity. Binding of O(2), CO, or NO to a reduced heme significantly enhances the PDE activity toward 3',5'-cyclic diguanylic acid. We report stationary and time-resolved resonance Raman spectra of the wild-type and several mutants (Glu-93 --> Ile, Met-95 --> Ala, Arg-97 --> Ile, Arg-97 --> Ala, Arg-97 --> Glu, Phe-113 --> Leu, and Phe-113 --> Thr) of the heme-containing PAS domain of Ec DOS. For the CO- and NO-bound forms, both the hydrogen-bonded and non-hydrogen-bonded conformations were found, and in the former Arg-97 forms a hydrogen bond with the heme-bound external ligand. The resonance Raman results revealed significant interactions of Arg-97 and Phe-113 with a ligand bound to the sixth coordination site of the heme and profound structural changes in the heme propionates upon dissociation of CO. Mutation of Phe-113 perturbed the PDE activities, and the mutation of Arg-97 and Phe-113 significantly influenced the transient binding of Met-95 to the heme upon photodissociation of CO. This suggests that the electrostatic interaction of Arg-97 and steric interaction of Phe-113 are crucial for regulating the competitive recombination of Met-95 and CO to the heme. On the basis of these results, we propose a model for the role of the heme propionates in communicating the heme structural changes to the protein moiety.     

A common site of the Fc receptor gamma subunit interacts with the unrelated immunoreceptors FcalphaRI and FcepsilonRI

The transmembrane (TM) region of the Fc receptor-gamma (FcRgamma) chain is responsible for the association of this ubiquitous signal transduction subunit with many immunoreceptor ligand binding chains, making FcRgamma key to a number of leukocyte activities in immunity and disease. Some receptors contain a TM arginine residue that interacts with Asp-11 of the FcRgamma subunit, but otherwise the molecular basis for the FcRgamma subunit interactions is largely unknown. This study reports residues in the TM region of the FcRgamma subunit are important for association with the high affinity IgE receptor FcepsilonRI and a leukocyte receptor cluster member, the IgA receptor FcalphaRI. FcRgamma residue Leu-21 was essential for surface expression of FcepsilonRIalpha/gamma2 and Tyr-8, Leu-14, and Phe-15 contributed to expression. Likewise, detergent-stable FcRgamma association with FcalphaRI was also dependent on Leu-14 and Leu-21 and in addition required residues Tyr-17, Tyr-25, and Cys-26. Modeling the TM regions of the FcRgamma dimer indicated these residues interacting with both FcalphaRI and FcepsilonRI are near the interface between the two FcRgamma TM helices. Furthermore, the FcRgamma residues interacting with FcalphaRI form a leucine zipper-like interface with mutagenesis confirming a complementary interface comprising FcalphaRI residues Leu-217, Leu-220, and Leu-224. The dependence of these two nonhomologous receptor interactions on FcRgamma Leu-14 and Leu-21 suggests that all the associated Fc receptors and the activating leukocyte receptor cluster members interact with this one site. Taken together these data provide a molecular basis for understanding how disparate receptor families assemble with the FcRgamma subunit.     

Thermodynamic consequences of disrupting a water-mediated hydrogen bond network in a protein:pheromone complex

The mouse pheromones (+/-)-2-sec-butyl-4,5-dihydrothiazole (SBT) and 6-hydroxy-6-methyl-3-heptanone (HMH) bind into an occluded hydrophobic cavity in the mouse major urinary protein (MUP-1). Although the ligands are structurally unrelated, in both cases binding is accompanied by formation of a similar buried, water-mediated hydrogen bond network between the ligand and several backbone and side chain groups on the protein. To investigate the energetic contribution of this hydrogen bond network to ligand binding, we have applied isothermal titration calorimetry to measure the binding thermodynamics using several MUP mutants and ligand analogs. Mutation of Tyr-120 to Phe, which disrupts a hydrogen bond from the phenolic hydroxyl group of Tyr-120 to one of the bound water molecules, results in a substantial loss of favorable binding enthalpy, which is partially compensated by a favorable change in binding entropy. A similar thermodynamic effect was observed when the hydrogen bonded nitrogen atom of the heterocyclic ligand was replaced by a methyne group. Several other modifications of the protein or ligand had smaller effects on the binding thermodynamics. The data provide supporting evidence for the role of the hydrogen bond network in stabilizing the complex.     

Molecular insights into the binding selectivity of a synthetic ligand DAAG to Fc fragment of IgG

Affinity chromatography with synthetic ligands has been focused as the potential alternative to protein A‐based chromatography for antibody capture because of its comparable selectivity and efficiency. Better understanding on the molecular interactions between synthetic ligand and antibody is crucial for improving and designing novel ligands. In this work, the molecular interaction mechanism between Fc fragment of IgG and a synthetic ligand (DAAG) was studied with molecular docking and dynamics simulation. The docking results on the consensus binding site (CBS) indicated that DAAG could bind to the CBS with the favorable orientation like a tripod for the top‐ranked binding complexes. The ligand‐Fc fragment complexes were then tested by molecular dynamics simulation at neutral condition (pH 7.0) for 10 ns. The results indicated that the binding of DAAG on the CBS of Fc fragment was achieved by the multimodal interactions, combining the hydrophobic interaction, electrostatic interaction, hydrogen bond, and so on. It was also found that multiple secondary interactions endowed DAAG with an excellent selectivity to Fc fragment. In addition, molecular dynamics simulation conducted at acidic condition (pH 3.0) showed that the departure of DAAG ligand from the surface of Fc fragment was the result of reduced interaction energies. The binding modes between DAAG and CBS not only shed light on the molecular mechanisms of DAAG for antibody purification but also provide useful information for the improvement of ligand design. Copyright © 2014 John Wiley & Sons, Ltd.     

Proton nuclear magnetic resonance analyses of the molecular conformations of unique long neurotoxins bearing Phe-25: Astrotia stokesii b, Astrotia stokesii c, and Acanthophis antarcticus b

The 270-MHz proton NMR spectra of the unique long neurotoxins bearing Phe-25, Astrotia stokesii b (As b) and Astrotia stokesii c (As c) from Astrotia stokesii, and Acanthophis antarcticus b (Aa b) from Acanthophis antarcticus, have been analyzed. The aromatic proton resonances of Phe-25 in As b and Aa b were assigned on the basis of the nuclear Overhauser effects observed on irradiation of slowly exchanging amide protons. Phe-25 was found to be involved in hydrophobic interactions with Ile/Val-42, Ala-46 and Ile-58 in As b and As c, and with Ala-46 and Val-58 in Aa b. These hydrophobic interactions, instead of the hydrogen bond between Tyr-25 and Glu-42 found in other neurotoxins, appear to be important for maintenance of the biologically active tertiary structure. The pH dependency of the chemical shift and intensity of the Trp-72 N-1 proton resonance of As b indicates that the indole ring is not fully exposed to the solvent and that the extra tail segment of this long neurotoxin interacts with the main part of the molecule.     

Evidence for the involvement of tyrosine-69 in the control of stereospecificity of porcine pancreatic phospholipase A2

We have studied the role of Tyr-69 of porcine pancreatic phospholipase A2 in catalysis and substrate binding, using site-directed mutagenesis. A mutant was constructed containing Phe at position 69. Kinetic characterization revealed that the Phe-69 mutant has retained enzymatic activity on monomeric and micellar substrates, and that the mutation has only minor effects on kcat and Km. This shows that Tyr-69 plays no role in the true catalytic events during substrate hydrolysis. In contrast, the mutation has a profound influence on the stereospecificity of the enzyme. Whereas the wild-type phospholipase A2 is only able to catalyse the degradation of sn-3 phospholipids, the Phe-69 mutant hydrolyses both the sn-3 isomers and, at a low (1-2%) rate, the sn-1 isomers. Despite the fact that the stereospecificity of the mutant phospholipase has been altered, Phe-69 phospholipase still requires Ca2+ ions as a cofactor and also retains its specificity for the sn-2 ester bond. Our data suggest that in porcine pancreatic phospholipase A2 the hydroxyl group of Tyr-69 serves to fix and orient the phosphate group of phospholipid monomers by hydrogen bonding. Because no such interaction can occur between the Phe-69 side-chain and the phosphate moiety of the substrate monomer, the mutant enzyme loses part of its stereospecificity but not its positional specificity.     

Stability of the beta-sheet of the WW domain: A molecular dynamics simulation study.

The WW domain consists of approximately 40 residues, has no disulfide bridges, and forms a three-stranded antiparallel beta-sheet that is monomeric in solution. It thus provides a model system for studying beta-sheet stability in native proteins. We performed molecular dynamics simulations of two WW domains, YAP65 and FBP28, with very different stability characteristics, in order to explore the initial unfolding of the beta-sheet. The less stable YAP domain is much more sensitive to simulation conditions than the FBP domain. Under standard simulation conditions in water (with or without charge-balancing counterions) at 300 K, the beta-sheet of the YAP WW domain disintegrated at early stages of the simulations. Disintegration commenced with the breakage of a hydrogen bond between the second and third strands of the beta-sheet due to an anticorrelated transition of the Tyr-28 psi and Phe-29 phi angles. Electrostatic interactions play a role in this event, and the YAP WW domain structure is more stable when simulated with a complete explicit model of the surrounding ionic strength. Other factors affecting stability of the beta-sheet are side-chain packing, the conformational entropy of the flexible chain termini, and the binding of cognate peptide.     

Structural studies of the streptavidin binding loop.

The streptavidin-biotin complex provides the basis for many important biotechnological applications and is an interesting model system for studying high-affinity protein-ligand interactions. We report here crystallographic studies elucidating the conformation of the flexible binding loop of streptavidin (residues 45 to 52) in the unbound and bound forms. The crystal structures of unbound streptavidin have been determined in two monoclinic crystal forms. The binding loop generally adopts an open conformation in the unbound species. In one subunit of one crystal form, the flexible loop adopts the closed conformation and an analysis of packing interactions suggests that protein-protein contacts stabilize the closed loop conformation. In the other crystal form all loops adopt an open conformation. Co-crystallization of streptavidin and biotin resulted in two additional, different crystal forms, with ligand bound in all four binding sites of the first crystal form and biotin bound in only two subunits in a second. The major change associated with binding of biotin is the closure of the surface loop incorporating residues 45 to 52. Residues 49 to 52 display a 3(10) helical conformation in unbound subunits of our structures as opposed to the disordered loops observed in other structure determinations of streptavidin. In addition, the open conformation is stabilized by a beta-sheet hydrogen bond between residues 45 and 52, which cannot occur in the closed conformation. The 3(10) helix is observed in nearly all unbound subunits of both the co-crystallized and ligand-free structures. An analysis of the temperature factors of the binding loop regions suggests that the mobility of the closed loops in the complexed structures is lower than in the open loops of the ligand-free structures. The two biotin bound subunits in the tetramer found in the MONO-b1 crystal form are those that contribute Trp 120 across their respective binding pockets, suggesting a structural link between these binding sites in the tetramer. However, there are no obvious signatures of binding site communication observed upon ligand binding, such as quaternary structure changes or shifts in the region of Trp 120. These studies demonstrate that while crystallographic packing interactions can stabilize both the open and closed forms of the flexible loop, in their absence the loop is open in the unbound state and closed in the presence of biotin. If present in solution, the helical structure in the open loop conformation could moderate the entropic penalty associated with biotin binding by contributing an order-to-disorder component to the loop closure.     

Energetics of hydrogen bond switch, residue burial and cavity analysis reveals molecular basis of improved heparin binding to antithrombin

Antithrombin III (ATIII) is the main inhibitor of the coagulation proteases like factor Xa and thrombin. Anticoagulant activity of ATIII is increased by several thousand folds when activated by vascular wall heparan sulfate proteoglycans (HSPGs) and pharmaceutical heparins. ATIII isoforms in human plasma, alpha-ATIII and beta-ATIII differ in the amount of glycosylation which is the basis of differences in their heparin binding affinity and function. Crystal structures and site directed mutagenesis studies have mapped the heparin binding site in ATIII, however the hydrogen bond switch and energetics of interaction during the course of heparin dependent conformational change remains largely unclear. An analysis of heparin bound conformational states of ATIII using PEARLS software showed that in heparin bound intermediate state, Arg 47 and Arg 13 residues make hydrogen bonds with heparin but in the activated conformation Lys 11 and Lys 114 have more hydrogen bond interactions. In the protease bound-antithrombin-pentasaccharide complex Lys 114, Pro 12 and Lys 125 form important hydrogen bonding interactions. The results showed that A-helix and N-terminal end residues are more important in the initial interactions but D-helix is more important during the latter stage of conformational activation and during the process of protease inhibition. We carried out the residue wise Accessible Surface Area (ASA) analysis of alpha and beta ATIII native states and the results indicated major differences in burial of residues from Ser 112 to Ser 116 towards the N-terminal end. This region is involved in the P-helix formation on account of heparin binding. A cavity analysis showed a progressively larger cavity formation during activation in the region just adjacent to the heparin binding site towards the C-terminal end. We hypothesize that during the process of conformational change after heparin binding beta form of antithrombin has low energy barrier to form D-helix extension toward N and C-terminal end as compared to alpha isoform.     

Interplay between structural rigidity and electrostatic interactions in the ligand binding domain of GluR2

Using molecular dynamics (MD) simulations, computational protein modifications, and a novel theoretical methodology that determines structural rigidity/flexibility (the FIRST algorithm), we investigate how molecular structure and dynamics of the glutamate receptor ligand binding domain (GluR2 S1S2) facilitate its conformational transition. S1S2 is a two-lobe protein, which undergoes a cleft closure conformational transition upon binding an agonist in the cleft between the two lobes; hence it is expected that the mechanism of this conformational transition can be characterized as a hinge-type. However, in the rigidity analysis one lobe of the protein is identified as a single rigid cluster while the other one is structurally flexible, inconsistent with a presumed mechanical hinge mechanism. Instead, we characterize the cleft-closing transition as a load and lock mechanism. We find that when two cross-cleft hydrogen bonds are disrupted the protein undergoes a rapid cleft opening transition. At the same time, the dynamical behavior of the cleft in the presence of the glutamate ligand is only weakly affected by the S652 peptide bond in its flipped conformation observed in the crystal structure. The residue E705 plays significant role in stabilization of the closed conformation via electrostatic interactions. The presence of the E705-K730 salt bridge seems to correlate strongly withthe cleft opening transition in the MD simulations.     

Putative helix F contributes to regioselectivity of hydroxylation in mitochondrial cytochrome P450 27A1.

On the basis of alignment with structurally characterized cytochromes P450 (P450s), we have identified the putative F and G helices of mitochondrial P450s 27A1 and 11A. We introduced substitutions at Phe-207, Ile-211, and Phe-215 within putative helix F and at Trp-235 and Tyr-238 within putative helix G in P450 27A1 and compared wild type and mutants with respect to catalytic activity, product pattern, substrate binding, formation of hydrogen peroxide, and interaction with redox partner. Results indicate that the mutated residues are important for delivery of the correctly oriented substrate to the P450 active site. The I211K and F215K mutations, for example, affected the regioselectivity of P450 27A1-dependent hydroxylation reactions and conferred the P450 capacity to cleave the C-C bond of the substrate during the catalytic cycle. Studies of P450 11A1 indicate that Phe-202 has functions similar to those of its counterpart in P450 27A1 (Phe-215). We propose that putative helices F and G form the sides of the substrate-access channel, thus providing the additional mechanism to control regioselectivity of hydroxylation in mitochondrial P450s.     

Structural basis for calcium-induced inhibition of rhodopsin kinase by recoverin

Recoverin, a member of the neuronal calcium sensor branch of the EF-hand superfamily, serves as a calcium sensor that regulates rhodopsin kinase (RK) activity in retinal rod cells. We report here the NMR structure of Ca(2+)-bound recoverin bound to a functional N-terminal fragment of rhodopsin kinase (residues 1-25, called RK25). The overall main-chain structure of recoverin in the complex is similar to structures of Ca(2+)-bound recoverin in the absence of target (<1.8A root-mean-square deviation). The first eight residues of recoverin at the N terminus are solvent-exposed, enabling the N-terminal myristoyl group to interact with target membranes, and Ca(2+) is bound at the second and third EF-hands of the protein. RK25 in the complex forms an amphipathic helix (residues 4-16). The hydrophobic face of the RK25 helix (Val-9, Val-10, Ala-11, Ala-14, and Phe-15) interacts with an exposed hydrophobic groove on the surface of recoverin lined by side-chain atoms of Trp-31, Phe-35, Phe-49, Ile-52, Tyr-53, Phe-56, Phe-57, Tyr-86, and Leu-90. Residues of recoverin that contact RK25 are highly conserved, suggesting a similar target binding site structure in all neuronal calcium sensor proteins. Site-specific mutagenesis and deletion analysis confirm that the hydrophobic residues at the interface are necessary and sufficient for binding. The recoverin-RK25 complex exhibits Ca(2+)-induced binding to rhodopsin immobilized on concanavalin-A resin. We propose that Ca(2+)-bound recoverin is bound between rhodopsin and RK in a ternary complex on rod outer segment disk membranes, thereby blocking RK interaction with rhodopsin at high Ca(2+).     

Critical amino acids in phosphodiesterase-5 catalytic site that provide for high-affinity interaction with cyclic guanosine monophosphate and inhibitors

The molecular bases for phosphodiesterase 5 (PDE5) catalytic-site affinity for cyclic guanosine monophosphate (cGMP) and potency of inhibitors are poorly understood. Cocrystal structures of PDE5 catalytic (C) domain with inhibitors reveal a hydrogen bond and hydrophobic interactions with Tyr-612, hydrogen bonds with Gln-817, a hydrophobic clamp formed by Phe-820 and Val-782, and contacts with His-613, Leu-765, and Phe-786 [Sung et al. (2003) Nature 425, 98-102; Huai et al. (2004) J. Biol. Chem. 279, 13095-13101]. Present results of point mutations of full-length PDE5 showed that maximum catalysis was decreased 2650-fold in H613A and 55-fold in F820A. Catalytic-site affinities for cGMP, vardenafil, sildenafil, tadalafil, or 3-isobutyl-1-methylxanthine (IBMX) were respectively weakened 14-, 123-, 30-, 51-, and 43-fold for Y612A; 63-, 511-, 43-, 95- and 61-fold for Q817A; and 59-, 448-, 71-, 137-, and 93-fold for F820A. The data indicate that these three amino acids are major determinants of affinity for cGMP and potency of selective and nonselective inhibitors, and that higher vardenafil potency over sildenafil and tadalafil results from stronger contacts with Tyr-612, Gln-817, and Phe-820. Affinity of V782A for cGMP, vardenafil, sildenafil, tadalafil, or IBMX was reduced 5.5-, 23-, 10-, 3-, and 12-fold, respectively. Change in affinity for cGMP, vardenafil, sildenafil, or IBMX in Y612F, H613A, L765A, or F786A was less, but affinity of H613A or F786A for tadalafil was weakened 37- and 17-fold, respectively. The results quantify the role of PDE5 catalytic-site residues for cGMP and inhibitors, indicate that Tyr-612, Gln-817, and Phe-820 are the most important cGMP or inhibitor contacts studied, and identify residues that contribute to selectivity among different classes of inhibitors.     

The aromatic residues of bovine pancreatic ribonuclease studied by 1H nuclear magnetic resonance.

1. The aromatic proton resonances in the 360-MHz 1H nuclear magnetic resonance (NMR) spectrum of bovine pancreatic ribonuclease were divided into histidine, tyrosine and phenylalanine resonances by means of pH titrations and double resonance experiments. 2. Photochemically induced dynamic nuclear polarization spectra showed that one histidine (His-119) and two tyrosines are accessibly to photo-excited flavin. This permitted the identification of the C-4 proton resonance of His-119. 3. The resonances of the ring protons of Tyr-25, Tyr-76 and Tyr-115 and the C-4 proton of His-12 were identified by comparison with subtilisin-modified and nitrated ribonucleases. Other resonances were assigned tentatively to Tyr-73, Tyr-92 and Phe-46. 4. On addition of active-site inhibitors, all phenylalanine resonances broadened or disappeared. The resonance that was most affected was assigned tentatively to Phe-120. 5. Four of the six tyrosines of bovine RNase, identified as Tyr-76, Tyr-115 and, tentatively, Tyr-73 and Tyr-92, are titratable above pH 9. The rings of Tyr-73 and Tyr-115 are rapidly rotating or flipping by 180 degrees about their C beta--C gamma bond and are accessible to flavin in photochemically induced dynamic nuclear polarization experiments. Tyr-25 is involved in a pH-dependent conformational transition, together with Asp-14 and His-48. A scheme for this transition is proposed. 6. Binding of active-site inhibitors to bovine RNase only influences the active site and its immediate surroundings. These conformational changes are probably not connected with the pH-dependent transition in the region of Asp-14, Tyr-25 and His-48. 7. In NMR spectra of RNase A at elevated temperatures, no local unfolding below the temperature of the thermal denaturation was observed. NMR spectra of thermally unfolded RNase A indicated that the deviations from a random coil are small and might be caused by interactions between neighbouring residues.     

Molecular modelling of protein-carbohydrate interactions. Docking of monosaccharides in the binding site of concanavalin A.

A general procedure is described for addressing the computer simulation of protein-carbohydrate interactions. First, a molecular mechanical force field capable of performing conformational analysis of oligosaccharides has been derived by the addition of new parameters to the Tripos force field; it is also compatible with the simulation of protein. Second, a docking procedure which allows for a systematic exploration of the orientations and positions of a ligand into a protein cavity has been designed. This so-called 'crankshaft' method uses rotations and variations about/of virtual bonds connecting, via dummy atoms, the ligand to the protein binding site. Third, calculation of the relative stability of protein ligand complexes is performed. This strategy has been applied to search for all favourable interactions occurring between a lectin [concanavalin A (ConA)] and methyl alpha-D-mannopyranoside or methyl alpha-D-glucopyranoside. For each monosaccharide, different stable orientations and positions within the binding site can be distinguished. Among them, one corresponds to very favourable interactions, not only in terms of hydrogen bonding, but also in terms of van der Waals interactions. It corresponds precisely to the binding mode of methyl alpha-D-mannopyranoside into ConA as revealed by the 2.9 A resolution of the crystalline complex (Derewenda et al., 1989). Some implications of the present modelling study with respect to the molecular basis of the specificity of the interaction of lectins with various monosaccharides are presented.