Cross-linking reconsidered: binding and cross-linking fields and the cellular response.

We analyze a model for the reversible cross-linking of cell surface receptors by a collection of bivalent ligands with different affinities for the receptor as would be found in a polyclonal anti-receptor serum. We assume that the amount of cross-linking determines, via a monotonic function, the rate at which cells become activated and divide. In addition to the density of receptors on the cell surface, two quantities, the binding field and the cross-linking field, are needed to characterize the cross-linking curve, i.e., the equilibrium concentration of cross-linked receptors plotted as a function of the total ligand site concentration. The binding field is the sum of all ligand site concentrations weighted by their respective binding affinities, and the cross-linking field is the sum of all ligand site concentrations weighted by the product of their respective binding and cross-linking affinity and the total receptor density. Assuming that the cross-linking affinity decreases if the binding affinity decreases, we find that the height of the cross-linking curve decreases, its width narrows, and its center shifts to higher ligand site concentrations as the affinities decrease. Moreover, when we consider cross-linking-induced proliferation, we find that there is a minimum cross-linking affinity that must be surpassed before a clone can expand. We also show that under many circumstances a polyclonal antiserum would be more likely than a monoclonal antibody to lead to cross-linking-induced proliferation.     

Notochord-derived Shh concentrates in close association with the apically positioned basal body in neural target cells and forms a dynamic gradient during neural patterning

Sonic hedgehog (Shh) ligand secreted by the notochord induces distinct ventral cell identities in the adjacent neural tube by a concentration-dependent mechanism. To study this process, we genetically engineered mice that produce bioactive, fluorescently labeled Shh from the endogenous locus. We show that Shh ligand concentrates in close association with the apically positioned basal body of neural target cells, forming a dynamic, punctate gradient in the ventral neural tube. Both ligand lipidation and target field response influence the gradient profile, but not the ability of Shh to concentrate around the basal body. Further, subcellular analysis suggests that Shh from the notochord might traffic into the neural target field by means of an apical-to-basal-oriented microtubule scaffold. This study, in which we directly observe, measure, localize and modify notochord-derived Shh ligand in the context of neural patterning, provides several new insights into mechanisms of Shh morphogen action.     

Ion coordination in the amphotericin B channel.

The antifungal polyene antibiotic amphotericin B forms channels in lipid membranes that are permeable to ions, water, and nonelectrolytes. Anion, cation, and ion pair coordination in the water-filled pore of the "barrel" unit of the channels was studied by molecular dynamics simulations. Unlike the case of the gramicidin A channel, the water molecules do not create a single-file configuration in the pore, and some cross sections of the channel contain three or four water molecules. Both the anion and cation are strongly bound to ligand groups and water molecules located in the channel. The coordination number of the ions is about six. The chloride has two binding sites in the pore. The binding with water is dominant; more than four water molecules are localized in the anion coordination sphere. Three motifs of the ion coordination were monitored. The dominant motif occurs when the anion is bound to one ligand group. The ion is bound to two or three ligand groups in the less favorable configurations. The strong affinity of cations to the channel is determined by the negatively charged ligand oxygens, whose electrostatic field dominates over the field of the hydrogens. The ligand contribution to the coordination number of the sodium ion is noticeably higher than in the case of the anion. As in the case of the anion, there are three motifs of the cation coordination. The favorable one occurs when the cation is bound to two ligand oxygens. In the less favorable cases, the cation is bound to three or four oxygens. In the contact ion pair, the cation and anion are bound to two ligand oxygens and one ligand hydrogen, respectively. There exist intermediate solvent-shared states of the ion pair. The average distances between ions in these states are twice as large as that of the contact ion pair. The stability of the solvent-shared state is defined by the water molecule oriented along the electrostatic field of both ions.     

Novel heme ligand displacement by CO in the soluble hemophore HasA and its proximal ligand mutants: implications for heme uptake and release

HasASM, a hemophore secreted by the Gram-negative bacteria Serratia marcescens, extracts heme from host hemoproteins and shuttles it to HasRSM, a specific hemophore outer membrane receptor. Heme iron in HasASM is in a six-coordinate ferric state. It is linked to the protein by the heretofore uncommon axial ligand set, His32 and Tyr75. A third residue of the heme pocket, His83, plays a crucial role in heme ligation through hydrogen bonding to Tyr75. The vibrational frequencies of coordinated carbon monoxide constitute a sensitive probe of trans ligand field, FeCO structure, and electrostatic landscape of the distal heme pockets of heme proteins. In this study, carbonyl complexes of wild-type (WT) HasASM and its heme pocket mutants His32Ala, Tyr75Ala, and His83Ala were characterized by resonance Raman spectroscopy. The CO complexes of WT HasASM, HasASM(His32Ala), and HasASM(His83Ala) exhibit similar spectral features and fall above the line that correlates nuFe-CO and nuC-O for proteins having a proximal imidazole ligand. This suggests that the proximal ligand field in these CO adducts is weaker than that for heme-CO proteins bearing a histidine axial ligand. In contrast, the CO complex of HasASM(Tyr75Ala) has resonance Raman signatures consistent with ImH-Fe-CO ligation. These results reveal that in WT HasASM, the axial ImH side chain of His32 is displaced by CO. This is in contrast to other heme proteins known to have the His/Tyr axial ligand set, wherein the phenolic side chain of the Tyr ligand dissociates upon CO addition. The displacement of His32 and its stabilization in an unbound state is postulated to be relevant to heme uptake and/or release.     

A 3-D model of ligand transport in a deforming extracellular space

Cells communicate through shed or secreted ligands that traffic through the interstitium. Force-induced changes in interstitial geometry can initiate mechanotransduction responses through changes in local ligand concentrations. To gain insight into the temporal and spatial evolution of such mechanotransduction responses, we developed a 3-D computational model that couples geometric changes observed in the lateral intercellular space (LIS) of mechanically loaded airway epithelial cells to the diffusion-convection equations that govern ligand transport. By solving the 3-D fluid field under changing boundary geometries, and then coupling the fluid velocities to the ligand transport equations, we calculated the temporal changes in the 3-D ligand concentration field. Our results illustrate the steady-state heterogeneities in ligand distribution that arise from local variations in interstitial geometry, and demonstrate that highly localized changes in ligand concentration can be induced by mechanical loading, depending on both local deformations and ligand convection effects. The occurrence of inhomogeneities at steady state and in response to mechanical loading suggest that local variations in ligand concentration may have important effects on cell-to-cell variations in basal signaling state and localized mechanotransduction responses.     

Ligand-induced rearrangement of the dimeric metabotropic glutamate receptor 1alpha

The extracellular domain of the metabotropic glutamate receptor 1alpha (mGluR1alpha) forms a dimer and the ligand, glutamate, induces a structural rearrangement in this domain. However, the conformational change in the cytoplasmic domain, which is critical for mGluR1alpha's interaction with G proteins, remains unclear. Here we investigated the ligand-induced conformational changes in the cytoplasmic domain by fluorescence resonance energy transfer (FRET) analysis of mGluR1alpha labeled with fluorescent protein(s) under total internal reflection field microscopy. Upon ligand binding, the intersubunit FRET efficiency between the second loops increased, whereas that between first loops decreased. In contrast, the intrasubunit FRET did not change clearly. These results show that ligand binding does not change the structure of each subunit, but does change the dimeric allocation of the cytoplasmic regions, which may underlie downstream signaling.     

Studies on the substrate-induced spectral change of cytochrome P-450 in liver microsomes.

The spectral changes of cytochrome P-450 caused by the addition of small molecules to liver microsomes were investigated precisely and the following conclusions were reached. 1. The Type I spectral change was entirely due to the interaction of the cytochrome with a hydrocarbon residue in a ligand. To induce the modified Type II spectral change, the presence of a hydroxyl group in a ligand was required. Compounds which contain a basic amino group induced the Type II spectral change. 2. The Type I spectral change was caused by the interaction of a ligand with the 419-nm form of cytochrome P-450, with its concomitant conversion to the 394-nm form. Whereas, compounds inducing modified Type II spectral change interacted with the 394nm form of the cytochrome. In this case, however, the 394-nm form was not converted back to the 419-nm form but was converted to a new state showing an absorption peak at 416 nm. The Type II spectral change-inducing interaction of a ligand with the cytochrome could occur with all forms of the cytochrome. 3. Both Type II and modified Type II compounds bound to the cytochrome at heme iron, and converted the cytochrome into modified ferrihemochromes. On the other hand, the Type I interaction occurred ina protein moiety of the cytochrome, and probably caused a conformational change of the cytochrome accompanied either by weakening of the internal ligand interaction or by displacement of the ligand with another one having a weaker field at the heme iron. 4. Type I and each of other two types of binding of compounds with cytochrome P-450 could occur simultaneously.     

Determination of the binding mode and interacting amino-acids for dibasic H3 receptor antagonists

Due to its involvement in major CNS functions, the histamine H3 receptor (H3R) is the subject of intensive medicinal chemistry investigation, supported by the range of modern drug discovery tools, such as receptor modeling and ligand docking. Although the receptor models described to date share a majority of common traits, they display discrete alternatives in amino-acid conformation, rendering ligand binding modes quite different. Such variations impede structure-based drug design in the H3R field. In the present study, we used a combination of medicinal chemistry, receptor-guided and ligand-based methods to elucidate the binding mode of antagonists. The approaches converged towards a ligand orientation perpendicular to the membrane plane, bridging Glu206 of the transmembrane helix 5 to acidic amino acids of the extracellular loops. This consensus will help future structure-based drug design for H3R ligands.     

Contributions to the binding free energy of ligands to avidin and streptavidin

The free energy of binding of a ligand to a macromolecule is here formally decomposed into the (effective) energy of interaction, reorganization energy of the ligand and the macromolecule, conformational entropy change of the ligand and the macromolecule, and translational and rotational entropy loss of the ligand. Molecular dynamics simulations with implicit solvation are used to evaluate these contributions in the binding of biotin, biotin analogs, and two peptides to avidin and streptavidin. We find that the largest contribution opposing binding is the protein reorganization energy, which is calculated to be from 10 to 30 kcal/mol for the ligands considered here. The ligand reorganization energy is also significant for flexible ligands. The translational/rotational entropy is 4.5-6 kcal/mol at 1 M standard state and room temperature. The calculated binding free energies are in the correct range, but the large statistical uncertainty in the protein reorganization energy precludes precise predictions. For some complexes, the simulations show multiple binding modes, different from the one observed in the crystal structure. This finding is probably due to deficiencies in the force field but may also reflect considerable ligand flexibility.     

Model systems for the coordination chemistry of cytochrome P-450.

Studies on model complexes have supported the presence of a mercaptide as the fifth ligand of cytochrome P-450 monooxygenases. When alcohol or thiol ligands are added to the sixth coordination position of a five-coordinated 4-nitrobenzene thiolate complex of FeIII protoporphyrin IX dimethyl ester chloride low spin complexes with optical and EPR-spectra very similar to cytochrome P-450 are obtained. From a comparison with all ligands of cytochrome P-450 and the model complexes it is concluded that a hard ligands must occupy the sixth coordination position of cytochrome P-450. An imidazole group is less likely, also in view of the ligand field parameters. The significance of the fifth and sixth ligand of cytochrome P-450 is discussed with respect to the monooxygenase mechanism.     

Recognition and discrimination of gases by the oxygen-sensing signal transducer protein HemAT as revealed by FTIR spectroscopy

The determination of ligand binding properties is a key step in our understanding of gas sensing and discrimination by gas sensory proteins. HemAT is a newly discovered signal transducer heme protein that recognizes O(2) and discriminates against other gases such as CO and NO. We have used FTIR spectroscopy on CO- and NO-bound sensor domain HemAT and sensor domain distal mutants Y70F, T95A, R91A, and L92A to gain insight into the structure of the iron-bound ligand at ambient temperature. These mutations were designed to perturb the electrostatic field near the iron-bound gaseous ligand and also allow us to investigate the communication pathway between the distal residues of the protein and the heme. We show the formation of both H-bonded and non-H-bonded conformations in the CO-bound forms. In addition, we report the presence of multiple conformations in the NO-bound forms. Such distal H-bonding is crucial for ligand binding and activation by the heme. The comparison of the O(2), NO, and CO data demonstrates that Thr95 and Tyr70 are crucial for ligand recognition and discrimination and, thus, for specific sensing of gases, and L92 is crucial for controlling the conformational changes of the Thr95 and Tyr70 residues upon NO binding.     

Electron paramagnetic resonance and optical spectroscopic study of the Cu2+-tRNA system

The interaction of copper ions with tRNA has been studied by optical and EPR spectroscopies. The interaction results in two different paramagnetic complexes characterized by a tetragonal symmetry of the ligand electric field sensed by the ions. The complete set of the spin Hamiltonian parameters has been extracted by computer simulation with the Monte Carlo method. Hypotheses concerning the putative ligands are put forward.     

Evaluation of hydrogen bonds of ecdysteroids in the ligand–receptor interactions using a protein modeling system

The insect molting hormone, 20-hydroxyecdysone (20E) and its analogs (ecdysteroids) specifically bind to the ecdysone receptor. Previously, we synthesized various ecdysteroids containing the side chain moiety of ponasterone A (PonA), and measured the binding activity against Drosophila Kc cells to study the structure–activity relationship. Here we quantitatively analyzed the structure–activity relationship for the ligand binding of ecdysteroids including 20E and PonA. Since the hydrogen bonding (HB) is one of the important physicochemical properties for ligand binding to the ecdysteroid receptor, the number of possible HBs between the ligand molecule and the receptor was manually counted in the modeled ligand–receptor complex for all compounds. The construction of the ligand–receptor model was executed by the full-automatic modeling system (FAMS) in which calculation was done by simulated annealing. The binding potency of 15 ecdysteroids to Kc-cells were linearly correlated (r² = 0.63) with the number of HBs which are observed between ligand and receptor molecule. Contribution of steric and electrostatic effects on the ligand–receptor binding was also examined using a three-dimensional quantitative structure–activity relationship (3-D QSAR), comparative molecular field analysis (CoMFA).     

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.     

Low-temperature electron paramagnetic resonance study of the ferricytochrome c-cardiolipin complex

The electron paramagnetic resonance spectrum of the ferricytochrome c complex with cardiolipin was observed at temperatures below 20 K. For the low-spin iron(III) heme system complexed with the negatively charged lipid, the tetragonal and rhombic ligand field parameters (delta/lambda = 3.58, V/lambda = 1.82) differ significantly from those (delta/lambda = 2.53, V/lambda = 1.49) of the free ferricytochrome c sample. The g values of the complex (gx = 1.54 +/- 0.02, gy = 2.26 +/- 0.01, gz = 3.02 +/- 0.01) are compared to the values for free ferricytochrome c (gx = 1.25 +/- 0.02, gy = 2.25 +/- 0.01, gz = 3.04 +/- 0.01). Spectral alterations are interpreted in terms of the ligand field changes induced within the heme group by association with the negatively charged phosphoglyceride.     

Probing the local electronic and geometric properties of the heme iron center in a H-NOX domain

Heme-Nitric oxide and/or OXygen binding (H-NOX) proteins are a family of diatomic gas binding hemoproteins that have attracted intense research interest. Here we employ X-ray absorption near-edge structure (XANES) spectroscopy to study the nitric oxide (NO) binding site of H-NOX. This is the first time this technique has been utilized to examine the NO/H-NOX signaling pathway. XANES spectra of wildtype and a point mutant (proline 115 to alanine, P115A) of the H-NOX domain from Thermoanaerobacter tengcongensis (Tt H-NOX) were obtained and analyzed for ferrous and ferric complexes of the protein. This work provides specific structural characterization of the solution state of several Tt H-NOX ferrous complexes (− unligated, − NO, and − CO) that were previously unavailable. Our iron K-edges indicate effective charge on the iron center in the various complexes and report on the electronic environment of heme iron. We analyzed the ligand field indicator ratio (LFIR), which is extracted from XANES spectra, for each complex, providing an understanding of ligand field strength, spin state of the central iron, movement of the iron atom upon ligation, and ligand binding properties. In particular, our LFIRs indicate that the heme iron is dramatically displaced towards the distal pocket during ligand binding. Based on these results, we propose that iron displacement towards the distal heme pocket is an essential step in signal initiation in H-NOX proteins. This provides a mechanistic link between ligand binding and the changes in heme and protein conformation that have been observed for H-NOX family members during signaling.     

Cytochrome c-554 from Methylosinus trichosporium OB3b; a protein that belongs to the cytochrome c2 family and exhibits a HALS-Type EPR signal

A small soluble cytochrome c-554 purified from Methylosinus trichosporium OB3b has been purified and analyzed by amino acid sequencing, mass spectrometry, visible, CD and EPR spectroscopies. It is found to be a mono heme protein with a characteristic cytochrome c fold, thus fitting into the class of cytochrome c(2), which is the bacterial homologue of mitochondrial cytochrome c. The heme iron has a Histidine/Methionine axial ligation and exhibits a highly anisotropic/axial low spin (HALS) EPR signal, with a g(max) at 3.40, and ligand field parameters V/ξ = 0.99, Δ/ξ = 4.57. This gives the rhombicity V/Δ = 0.22. The structural basis for this HALS EPR signal in Histidine/Methionine ligated hemes is not resolved. The ligand field parameters observed for cytochrome c-554 fits the observed pattern for other cytochromes with similar ligation and EPR behaviour.     

Independent-systems mechanisms for f–f transition probabilities in lanthanide coordination compounds

The Laporte-forbidden f–f electronic transitions of lanthanide coordination compounds acquire an electric-dipole probability by two mechanisms in the general independent-systems model, where overlap between the charge distributions of the metal ion and the ligands is neglected. A first-order electric-dipole transition moment arises, either from the mixing of the f–f with f–d and f–g electron promotions under the electrostatic field of the ligands, or from transient dipoles induced in the ligand groups, by an allowed even-multipole electric moment of the f–f excitation. The electrostatic field and the ligand polarization mechanisms make complementary intensity contributions to the f–f transitions of a given Ln(III) complex, dependent upon the rank of the leading electric multipole moment. The polarization mechanism contributes principally to the intensities of the ligand-hypersensitive 2²-pole f–f transitions, whereas the electrostatic mechanism is predominant for the 2⁶-pole transition intensities, and makes the more important contribution in the 2⁴-pole cases. Applied initially to Ln(III) complexes containing monoatomic ligands, which have an effective isotropic polarizability, the ligand polarization mechanism is found to depend, on extension to the corresponding polyatomic ligand cases, upon the anisotropy of the ligand polarizability tensor in complexes belonging to the higher non-centric symmetries.