Bound conformations for ligands for G-protein coupled receptors

The conformation of the C-terminus of the -subunit of transducin, the G-protein of vision, has been determined by transfer NOE when bound to activated (MII) rhodopsin. One hundred three new NOE constraints are apparent when light is shown on a mixture of rhodopsin bilayers and the undecapeptide. Analogs of the -peptide with covalent constraints were designed restricting the bound conformation; they stabilize MII thus supporting the deduced structure. The NMR structure of a complex of the intracellular loops of rhodopsin facilitates docking of the -peptide and also shows proximity of residues known by mutational analysis to interact to generate the activated rhodopsin-transducin interface. This constrains the location of transmembrane helices in the structure of activated rhodopsin. Methods for the prediction of affinity have been used to estimate the relative binding constants of peptide analogs with the loop complex and show strong correlation with experimental data. Various models of the rhodopsin-transmembrane helical segments have been computationally fused with distance geometry to determine the overall model which best fits the experimental data on the rhodopsin-transducin interface.     

Modeling of the complex between transducin and photoactivated rhodopsin, a prototypical G-protein-coupled receptor

Obtaining a reliable 3D model for the complex formed by photoactivated rhodopsin (R*) and its G-protein, transducin (Gtalphabetagamma), would significantly benefit the entire field of structural biology of G-protein-coupled receptors (GPCRs). In this study, we have performed extensive configurational sampling for the isolated C-terminal fragment of the alpha-subunit of transducin, Gtalpha 340-350, within cavities of photoactivated rhodopsin formed by different energetically feasible conformations of the intracellular loops. Our results suggested a new 3D model of the rhodopsin-transducin complex that fully satisfied all available experimental data on site-directed mutagenesis of rhodopsin and Gtalphabetagamma as well as data from disulfide-linking experiments. Importantly, the experimental data were not used as a priori constraints in model building. We performed a thorough comparison of existing computational models of the rhodopsin-transducin complex with each other and with current experimental data. It was found that different models suggest interactions with different molecules in the rhodopsin oligomer, that providing valuable guidance in design of specific novel experimental studies of the R*-Gtalphabetagamma complex. Finally, we demonstrated that the isolated Gtalpha 340-350 fragment does not necessarily bind rhodopsin in the same binding mode as the same segment in intact Gtalpha.     

Interaction between the main components of the vertebrate retinal rod phototransduction system in solutions of detergent n-nonyl-β-D-glucoside

cGMP-Specific phosphodiesterase (PDE6) is the key enzyme of the phototransduction system of vertebrate retinal rod outer segments (ROS). The properties of PDE in extracts prepared by solubilization of bovine ROS in a high concentration (0.5% w/v) of detergent n-nonyl-β-D-glucoside (NG) and following centrifugation (ROS-NG) have been studied. Basal PDE activity of the preparations was low, but it greatly (>50-fold) increased (up to ∼20 μmol cGMP hydrolyzed/min per mg rhodopsin (R)) in the presence of trypsin. In bleached GTPγS-containing preparations the specific PDE activity was dependent on ROS-NG concentration and was half-maximal at about 0.8 μM of ROS G protein transducin (Gt). In dark-adapted GTPγS-containing ROS-NG preparations bleaching of 0.2% of the rhodopsin resulted in half-maximal PDE activation. The same result was obtained when PDE in dark-adapted ROS-NG preparations was activated by addition of a highly purified bleached rhodopsin solubilized by 0.5% solution of NG. The results demonstrate that the presence of NG has no significant influence either on the properties of the main ROS phototrans-duction system elements (R, Gt and PDE) or on the interaction between photoactivated R and Gt and suggest that the detergent NG can be used for crystallization of the rhodopsin-transducin complex.     

A transferred NOE study of a tricyclic analog of acyclovir bound to thymidine kinase

Summary A purine derivative with an acyclic sugar analog, 3,9-dihydro-3-[(2-hydroxyethoxy)methyl]-6-ethyl-9-oxo-5H-imidazo[1,2-a]purine, was studied in the free state and in complex with herpes simplex virus thymidine kinase (HSV1 TK). Transferred NOE experiments, combined with a full relaxation matrix analysis of the substrate's spin system, resulted in a set of distance constraints for all proton pairs. These constraints were used in structure determination procedures based on simulated annealing and molecular dynamics simulations to obtain a family of structures compatible with the experimental NMR data. The results indicate that, although in both states the chains have the syn orientation with respect to the aromatic rings, in the free state the substrate's acyclic moiety is relatively disordered, while in the bound state only one specific conformation is preferred. Fluctuations can only be seen in the case of the terminal hydroxyl group, for which no NOE was recorded and hence no constraints were available.     

A simple spectral-driven procedure for the refinement of DNA structures by NMR spectroscopy.

We have developed a simple and quantitative procedure (SPEDREF) for the refinement of DNA structures using experimental two-dimensional nuclear Overhauser effect (2D NOE) data. The procedure calculates the simulated 2D NOE spectrum using the full matrix relaxation method on the basis of a molecular model. The volume of all NOE peaks is measured and compared between the experimental and the calculated spectra. The difference of the experimental and simulated volumes is minimized by a conjugated gradient procedure to adjust the interproton distances in the model. An agreement factor (analogous to the crystallographic R-factor) is used to monitor the progress of the refinement. The procedure is an The agreement is considered to be complete when several parameters, including the R-factor, the energy associated with the molecule, the local conformation (as judged by the sugar pseudorotation), and the global conformation (as judged by the helical x-displacement), are refined to their respective convergence. With the B-DNA structure of d(CGATCG) as an example, we show that DNA structure may be refined to produce calculated NOE spectra that are in excellent agreement with the experimental 2D NOE spectra. This is judged to be effective by the low R-factor of approximately 15%. Moreover, we demonstrate that not only are NOE data very powerful in providing details of the local structure but, with appropriate weighting of the NOE constraints, the global structure of the DNA double helix can also be determined, even when starting with a grossly different model. The reliability and limitations of a DNA structure as determined by NMR spectroscopy are discussed.     

Three-dimensional structure of α-conotoxin EI determined by1H NMR spectroscopy

Summary α-conotoxin EI is an 18-residue peptide (RDOCCYHPTCNMSNPQIC; 4–10, 5–18) isolated from the venom ofConus ermineus, the only fish-hunting cone snail of the Atlantic Ocean. This peptide targets specifically the nicotinic acetylcholine receptor (nAChR) found in mammalian skeletal muscle and the electric organTorpedo, showing a novel selectivity profile when compared to other α-conotoxins. The 3D structure of EI has been determined by 2D-NMR methods in combination with dynamical simulated annealing protocols. A total of 133 NOE-derived distances were used to produce 13 structures with minimum energy that complied with the NOE restraints. The structure of EI is characterized by a helical loop between THr⁹ and Met¹² that is stabilized by the Cys⁴-Cys¹⁰ disulfide bond and turns involving Cys⁴-Cys⁵ and Asn¹⁴-Pro¹⁵. Other regions of the peptide appear to be flexible. The overall fold of EI is similar to that of other α4/7-conotoxins (PnIA/B, MII, EpI). However, unlike these other α4/7-conotoxins, EI targets the muscular type nAChR. The differences in selectivity can be attributed to differences in the surface charge distribution among these α4/7-conotoxins. The implications for binding of EI to the muscular nAChR are discussed with respect to the current NMR structure of EI. Supplementary material available:¹H resonance assignments of α-conotoxin EI.     

Sequential nearest-neighbor effects on computed <Superscript>13</Superscript>C<Superscript>α</Superscript> chemical shifts

To evaluate sequential nearest-neighbor effects on quantum-chemical calculations of ¹³C^α chemical shifts, we selected the structure of the nucleic acid binding (NAB) protein from the SARS coronavirus determined by NMR in solution (PDB id 2K87). NAB is a 116-residue α/β protein, which contains 9 prolines and has 50% of its residues located in loops and turns. Overall, the results presented here show that sizeable nearest-neighbor effects are seen only for residues preceding proline, where Pro introduces an overestimation, on average, of 1.73 ppm in the computed ¹³C^α chemical shifts. A new ensemble of 20 conformers representing the NMR structure of the NAB, which was calculated with an input containing backbone torsion angle constraints derived from the theoretical ¹³C^α chemical shifts as supplementary data to the NOE distance constraints, exhibits very similar topology and comparable agreement with the NOE constraints as the published NMR structure. However, the two structures differ in the patterns of differences between observed and computed ¹³C^α chemical shifts, Δ _ca,i , for the individual residues along the sequence. This indicates that the Δ _ca,i -values for the NAB protein are primarily a consequence of the limited sampling by the bundles of 20 conformers used, as in common practice, to represent the two NMR structures, rather than of local flaws in the structures.     

An examination of the binding behavior of histidine-containing peptides with immobilized metal complexes derived from the macrocyclic ligand, 1,4,7-triazacyclononane

In this study, two different experimental approaches have been employed to examine the binding behavior of histidine-containing peptides with metal ion complexes derived from the macrocyclic ligand 1,4,7-triazacyclononane (tacn). Firstly, a molecular modeling approach has been employed to derive the strain energies for test peptide sequences that have a predicted propensity to readily adopt an α-helical conformation. To this end, binuclear metal complexes were examined with peptides containing two histidine residues in different locations in a pair of peptides of the same composition but different sequence. These modeling results indicate that there are no energetic constraints for two-point binding to occur with dicopper(II) binuclear complexes when two histidine residues are appropriately placed in an α-helical conformation. Secondly, binding experiments were carried out to establish the effect of one or more histidine residues within a peptide sequence on the affinity of a peptide for these Cu(II)–tacn derived binuclear complexes when immobilized onto a chromatographic support material. The results confirm that for all chelating systems, higher affinity is achieved as the histidine number in the peptide structure increases, although the relative location of the histidine residues in these small peptides did not introduce a significant constraint to the conformation on interacting with the immobilized Cu(II) binuclear complexes.     

Solution conformation of an immunogenic peptide from HRV2: comparison with the conformation found in a complex with a Fab fragment of an anti-HRV2 neutralizing antibody

The conformation of a [15]-peptide (H-VKAETRLNPDLQPTE-NH₂) from VP2 of rhinovirus HRV2 complexed with a Fab fragment was previously shown by X-ray crystallographic studies to be similar to the one found in the corresponding region of HRV1A. Antibodies raised against this peptide bind to and neutralize HRV2. In order to identify structural features preserved in solution that may explain the ability of this short peptide to mimic the structure of the protein surface, the peptide has been studied by NMR in aqueous solution as well as under denaturing conditions. The peptide is shown to be a random coil in solution. However, the sequence forming a 3₁₀ helix in the complex is biased into a helical conformation according to NOE intensity data as well as from urea and pH titrations. This sequence adopts the same conformation in an unrelated protein. NOE data suggest that a β-turn found in the complex may be sampled in solution. Also, Glu4, interacting with Arg6 in the crystal, has a reduced pK_a value in solution. It is concluded that the local structure present in the random coil state of VP2(156–170) contains enough information to direct the production of antibodies that bind to and neutralize HRV2. © 1998 European Peptide Society and John Wiley & Sons, Ltd.     

Refinement of the NMR structures for acyl carrier protein with scalar coupling data

Structure determination of small proteins using NMR data is most commonly pursued by combining NOE derived distance constraints with inherent constraints based on chemical bonding. Ideally, one would make use of a variety of experimental observations, not just distance constraints. Here, coupling constant constraints have been added to molecular mechanics and molecular dynamics protocols for structure determination in the form of a psuedoenergy function that is minimized in a search for an optimum molecular conformation. Application is made to refinement of a structure for a 77 amino acid protein involved in fatty acid synthesis, Escherichia coli acyl carrier protein (ACP). 54 3JHN alpha coupling constants, 12 coupling constants for stereospecifically assigned side chain protons, and 450 NOE distance constraints were used to calculate the 3-D structure of ACP. A three-step protocol for a molecular dynamics calculation is described, in analogy to the protocol previously used in molecular mechanics calculations. The structures calculated with the molecular mechanics approach and the molecular dynamics approach using a rigid model for the protein show similar molecular energies and similar agreement with experimental distance and coupling constant constraints. The molecular dynamics approach shows some advantage in overcoming local minimum problems, but only when a two-state averaging model for the protein was used, did molecular energies drop significantly.     

Interactions of metarhodopsin II. Arrestin peptides compete with arrestin and transducin

Arrestin blocks the interaction of rhodopsin with the G protein transducin (G(t)). To characterize the sites of arrestin that interact with rhodopsin, we have utilized a spectrophotometric peptide competition assay. It is based on the stabilization of the active intermediates metarhodopsin II (MII) and phosphorylated MII by G(t) and arrestin, respectively (extra MII monitor). The protocol involves native disc membranes and three sets of peptides 10-30 amino acids in length spanning the arrestin sequence. In the absence of arrestin, not one of the peptides by itself had an effect on the amount of MII formed. However, inhibition of arrestin-dependent extra MII was found for the peptides at residues 11-30 and 51-70 (IC(50) < 100 microm) and residues 231-260 (IC(50) < 200 microm). A similar pattern of inhibition by arrestin peptides was seen when arrestin was replaced by G(t) or the farnesylated G(t)gamma C-terminal peptide. Only arrestin-(11-30) inhibited MII.G(t) less (IC(50) = 300 microm) than phosphorylated MII.arrestin. We interpreted the data by competition of the arrestin peptides for interaction sites at rhodopsin, exposed in the MII conformation and specific for both arrestin and G(t). The arrestin sites are located in both the C- and N-terminal domains of the arrestin structure.     

Phosphorylation modulates the affinity of light-activated rhodopsin for G protein and arrestin

Reduced effector activity and binding of arrestin are widely accepted consequences of GPCR phosphorylation. However, the effect of receptor multiphosphorylation on G protein activation and arrestin binding parameters has not previously been quantitatively examined. We have found receptor phosphorylation to alter both G protein and arrestin binding constants for light-activated rhodopsin in proportion to phosphorylation stoichiometry. Rod disk membranes containing different average receptor phosphorylation stoichiometries were combined with G protein or arrestin, and titrated with a series of brief light flashes. Binding of G(t) or arrestin to activated rhodopsin augmented the 390 nm MII optical absorption signal by stabilizing MII as MII.G or MII.Arr. The concentration of active arrestin or G(t) and the binding constant of each to MII were determined using a nonlinear least-squares (Simplex) reaction model analysis of the titration data. The binding affinity of phosphorylated MII for G(t) decreased while that for arrestin increased with each added phosphate. G(t) binds more tightly to MII at phosphorylation levels less than or equal to two phosphates per rhodopsin; at higher phosphorylation levels, arrestin binding is favored. However, arrestin was found to bind much more slowly than G(t) at all phosphorylation levels, perhaps allowing time for phosphorylation to gradually reduce receptor-G protein interaction before arrestin capping of rhodopsin. Sensitivity of the binding constants to ionic strength suggests that a strong membrane electrostatic component is involved in both the reduction of G(t) binding and the increase of arrestin binding with increasing rhodopsin phosphorylation.     

Light-regulated localization of the beta-subunit of Gq-type G-protein in the crayfish photoreceptors

In crayfish photoreceptor cells, Gq-type G-protein plays a central role in the phototransduction pathway, and the translocation of Gqα has been proposed as one of the molecular mechanisms to control photoreceptor sensitivity. We here investigated β subunit of Gq and its localization profiles under various light conditions in the crayfish photoreceptor cells to understand the functional characteristic of visual Gq in the phototransduction pathway. An immunoprecipitation experiment was performed using an anti-Gqα antibody and a thiol-cleavable crosslinker. A 39 kDa protein was co-immunoprecipitated with Gqα, but not by irradiation, in the presence of GTPγS. The partial amino acid sequence of the 39 kDa protein was similar to Gβe in Drosophila photoreceptors, indicating that the crayfish Gβ which combines with Gqα is a Gβe homologue. Immunohistochemical and immunoblot analyses revealed that the amount of the Gβ decreased in the rhabdomeric membranes and increased in the cytoplasm in the light, compared with that in the dark. The profile of the translocation was similar to that reported for Gqα. Since both α and βγ subunits are necessary for G-proteins to be activated by rhodopsin in the rhabdom, the light-modulated translocation of a Gβe homologue possibly controls the amount of Gq which can be activated by light-stimulated rhodopsin. Accepted: 27 June 1998     

The use of pseudocontact shifts to refine solution structures of paramagnetic metalloproteins: Met80Ala cyano-cytochrome c as an example

 The availability of NOE constraints and of the relative solution structure of a paramagnetic protein permits the use of pseudocontact shifts as further structural constraints. We have developed a strategy based on: (1) determination of the χ tensor anisotropy parameters from the starting structure; (2) recalculation of a new structure by using NOE and pseudocontact shift constraints simultaneously; (3) redetermination of the χ tensor anisotropy parameters from the new structure, and so on until self-consistency. The system investigated is the cyanide derivative of a variant of the oxidized Saccharomyces cerevisiae iso-1-cytochrome c containing the Met80Ala mutation. The structure has been substantially refined. It is shown that the analysis of the deviation of the experimental pseudocontact shifts from those calculated using the starting structure may be unsound, as may the simple structure refinement based on the pseudocontact shift constraints only. Received: 11 July 1995 / Accepted: 30 October 1995     

Validation of the use of intermolecular NOE constraints for obtaining docked structures of protein-ligand complexes

Summary The use of intermolecular NOEs for docking a small ligand molecule into its target protein has been investigated with the aim of determining the effectiveness and methodology of this type of NOE docking calculation. A high-resolution X-ray structure of a protein-ligand complex has been used to simulate loose distance constraints of varying degrees of quality, typical of those estimated from experimental NOE intensities. These simulated data were used to examine the effect of the number, distribution and representation of the experimental constraints on the precision and accuracy of the calculated structures. A standard simulated annealing protocol was used, as well as a more novel method based on rigid-body dynamics. The results showed some analogies with those from similar studies on complete protein NMR structure determinations, but it was found that more constraints per torsion angle are required to define docked structures of similar quality. The effectiveness of different NOE-constraint averaging methods was explored and the benefits of using R^–6 averaging rather than centre averaging with small sets of NOE constraints were shown. The starting protein structure used in docking calculations was obtained from previous X-ray or NMR structure studies on a related complex. The effects on the calculated conformations of introducing structural differences into the binding site of the initial protein structure were also considered.To whom correspondence should be addressed.     

Association of α-subunits with nucleotide-modified β-subunits induces asymmetry in the catalytic sites of the F1-ATPase α3β3

The photoaffinity spin-labeled ATP analog, 2-N₃-SL-adenosine triphosphate (ATP), was used to covalently modify isolated β-subunits from F₁-ATPase of the thermophilic bacterium PS3. Approximately 1.2 mol of the nucleotide analog bound to the isolated subunit in the dark. Irradiation leads to covalent incorporation of the nucleotide into the binding site. ESR spectra of the complex show a signal that is typical for protein-immobilized radicals. Addition of isolated α-subunits to the modified β-subunits results in ESR spectra with two new signals indicative of two distinctly different environments of the spin-label, e.g., two distinctly different conformations of the catalytic sites. The relative ratio of the signals is approx 2∶1 in favor of the more closed conformation. The data show for the first time that when nucleotides are bound to isolated β-subunits, binding of α-subunits induces asymmetry in the catalytic sites even in the absence of the γ-subunit. This work was supported by a grant from the Deutsche Forschungsgemeinschaft to PDV.     

Binding of transducin and transducin-derived peptides to rhodopsin studies by attenuated total reflection-Fourier transform infrared difference spectroscopy.

Fourier transform infrared difference spectroscopy combined with the attenuated total reflection technique allows the monitoring of the association of transducin with bovine photoreceptor membranes in the dark. Illumination causes infrared absorption changes linked to formation of the light-activated rhodopsin-transducin complex. In addition to the spectral changes normally associated with meta II formation, prominent absorption increases occur at 1735 cm-1, 1640 cm-1, 1550 cm-1, and 1517 cm-1. The D2O sensitivity of the broad carbonyl stretching band around 1735 cm-1 indicates that a carboxylic acid group becomes protonated upon formation of the activated complex. Reconstitution of rhodopsin into phosphatidylcholine vesicles has little influence on the spectral properties of the rhodopsin-transducin complex, whereas pH affects the intensity of the carbonyl stretching band. AC-terminal peptide comprising amino acids 340-350 of the transducin alpha-subunit reproduces the frequencies and isotope sensitivities of several of the transducin-induced bands between 1500 and 1800 cm-1, whereas an N-terminal peptide (aa 8-23) does not. Therefore, the transducin-induced absorption changes can be ascribed mainly to an interaction between the transducin-alpha C-terminus and rhodopsin. The 1735 cm-1 vibration is also seen in the complex with C-terminal peptides devoid of free carboxylic acid groups, indicating that the corresponding carbonyl group is located on rhodopsin.     

HMG-D complexed to a bulge DNA: an NMR model

An NMR model is presented for the structure of HMG-D, one of the DROSOPHILA: counterparts of mammalian HMG1/2 proteins, bound to a particular distorted DNA structure, a dA(2) DNA bulge. The complex is in fast to intermediate exchange on the NMR chemical shift time scale and suffers substantial linebroadening for the majority of interfacial resonances. This essentially precludes determination of a high-resolution structure for the interface based on NMR data alone. However, by introducing a small number of additional constraints based on chemical shift and linewidth footprinting combined with analogies to known structures, an ensemble of model structures was generated using a computational strategy equivalent to that for a conventional NMR structure determination. We find that the base pair adjacent to the dA(2) bulge is not formed and that the protein recognizes this feature in forming the complex; intermolecular NOE enhancements are observed from the sidechain of Thr 33 to all four nucleotides of the DNA sequence step adjacent to the bulge. Our results form the first experimental demonstration that when binding to deformed DNA, non-sequence-specific HMG proteins recognize the junction between duplex and nonduplex DNA. Similarities and differences of the present structural model relative to other HMG-DNA complex structures are discussed.     

Methods of NMR structure refinement: molecular dynamics simulations improve the agreement with measured NMR data of a C-terminal peptide of GCN4-p1

The C-terminal trigger sequence is essential in the coiled-coil formation of GCN4-p1; its conformational properties are thus of importance for understanding this process at the atomic level. A solution NMR model structure of a peptide, GCN4p16–31, encompassing the GCN4-p1 trigger sequence was proposed a few years ago. Derived using a standard single-structure refinement protocol based on 172 nuclear Overhauser effect (NOE) distance restraints, 14 hydrogen-bond and 11 ϕ torsional-angle restraints, the resulting set of 20 NMR model structures exhibits regular α-helical structure. However, the set slightly violates some measured NOE bounds and does not reproduce all 15 measured ³J(H_N-H_Cα)-coupling constants, indicating that different conformers of GCN4p16–31 might be present in solution. With the aim to resolve structures compatible with all NOE upper distance bounds and ³J-coupling constants, we executed several structure refinement protocols employing unrestrained and restrained molecular dynamics (MD) simulations with two force fields. We find that only configurational ensembles obtained by applying simultaneously time-averaged NOE distance and ³J-coupling constant restraining with either force field reproduce all the experimental data. Additionally, analyses of the simulated ensembles show that the conformational variability of GCN4p16–31 in solution admitted by the available set of 187 measured NMR data is larger than represented by the set of the NMR model structures. The conformations of GCN4p16–31 in solution differ in the orientation not only of the side-chains but also of the backbone. The inconsistencies between the NMR model structures and the measured NMR data are due to the neglect of averaging effects and the inclusion of hydrogen-bond and torsional-angle restraints that have little basis in the primary, i.e. measured NMR data.     

On the use of one-step perturbation to investigate the dependence of NOE-derived atom–atom distance bound violations of peptides upon a variation of force-field parameters

The method of one-step perturbation can be used to predict from a single molecular dynamics simulation the values of observable quantities as functions of variations in the parameters of the Hamiltonian or biomolecular force field used in the simulation. The method is used to predict violations of nuclear overhauser effect (NOE) distance bounds measured in nuclear magnetic resonance (NMR) experiments by atom–atom distances of the NOE atom pairs when varying force-field parameters. Predictions of NOE distance bound violations between different versions of the GROMOS force field for a hexa-β-peptide in solution show that the technique works for rather large force-field parameter changes as well as for very different NOE bound violation patterns. The effect of changing individual force-field parameters on the NOE distance bound violations of the β-peptide and an α-peptide was investigated too. One-step perturbation, which in this case is equivalent to reweighting configurations, constitutes an efficient technique to predict many values of different quantities from a single conformational ensemble for a particular system, which makes it a powerful force-field development technique that easily reduces the number of required separate simulations by an order of magnitude.