PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “molecular dynamic simulation assay”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

179 records · Page 10Linked to original sources

A bicyclic and hsst2 selective somatostatin analogue: design, synthesis, conformational analysis and binding.

A backbone bridged and disulfide bridged bicyclic somatostatin analogue, compound 1 (PTR-3205), was designed and synthesized by solid-phase methodology. The binding of compound 1 to the five different somatostatin receptors, expressed in CHO or COS-7 cells, indicate a high degree of selectivity towards hsstr2. The three-dimensional structure of this compound has been determined in DMSO-d(6) and in water by 1H NMR and by molecular dynamics simulations. Similar backbone conformations were observed in both solvents. We have established direct evidence that the backbone of this bicyclic somatostatin analogue assumes a 'folded' conformation in solution, where the lactam ring extends roughly in the plane of the beta-turn. The pharmacophoric region Phe-(D)-Trp-Lys-Thr of compound 1 is in accord with that of both the Veber compound L-363,301 (Merck) and sandostatin. We believe that the enhanced selectivity towards the hsst2 receptor, in comparison with other analogues, is due to its large hydrophobic region, composed of the lactam ring and the Phe side chains at positions 1 and 8.

Animals↗

A novel endogenous inhibitor of phenoloxidase from Musca domestica has a cystine motif commonly found in snail and spider toxins.

Phenoloxidase inhibitor (POI), found in the hemolymph of housefly pupae, is a novel dopa-containing and cystine-rich peptide that competitively inhibits phenoloxidase with a Ki in the nanomolar range. [Tyr32]POI is a potential precursor molecule also found in the hemolymph that may be posttranslationally oxidized to the dopa-containing peptide after creation of a rigid structure. By employing both a solid-phase peptide synthesis system based on a 9-fluorenylmethoxycarbonyl strategy and a specific air oxidation technique to ensure correct folding, we have been able to synthesize [Tyr32]POI. The synthetic [Tyr32]POI was confirmed to be identical to the native [Tyr32]POI by coelution high-performance liquid chromatography analysis and by enzymatic analysis using the phenoloxidase inhibition assay. To determine the disulfide pairings within the peptides, a series of enzyme hydrolyses and partial reduction/alkylation steps were performed. Three cystine pairs (Cys11-Cys25, Cys18-Cys29, and Cys24-Cys36) were determined by identification of the resulting peptides. The disulfide pairings of the two adjacent Cys residues (Cys11-Cys25 and Cys24-Cys36) were unambiguously assigned by comparing the derived fragments with the two possible isomers synthesized through a novel disulfide-linking technique. The arrangement of the disulfide bridges in POI was found to be topologically identical to those found for several peptides within the inhibitor cystine knot structural family. Although these peptides share a low primary sequence homology and display a diversity of biological functions, they nonetheless share similarities in their cystine motifs and tertiary structure. The tertiary structure model of POI, which was derived through molecular dynamics and energy minimization studies using restraints with determined disulfide connectivities, suggests that POI is a new class member of the inhibitor cystine-knot structural family.

Amino Acid Sequence↗

The N-terminus moiety of the cystatin SmCys from Schistosoma mansoni regulates its inhibitory activity in vitro and in vivo.

The complete sequence of SmCys, a cystatin expressed by Schistosoma mansoni, was obtained. Constructs of SmCys consisting of deletions of 10 and 20 amino acid residues from the N-terminal of the full length recombinant protein, were cloned in the pQE-30 vector, expressed in Escherichia coli and assayed for inhibitory activity against papain. Kinetic analysis showed that SmCys -10 and SmCys -20 had K(i) values of 0.7391 and 4.9154, respectively, as compared to 0.0647, displayed by the full length recombinant. Protease inhibition by SmCys was also observed in vivo. When the recombinant products were incubated during 7 days with live schistosomula in the presence of red blood cells, only the full length product could completely inhibit the formation of haemozoin, a dark pigment formed as a by-product of haemoglobin digestion. The sequence data of the recombinant SmCys proteins were used for the construction of molecular models, which were then subjected to molecular dynamics for 2ns. In comparison to the full length, the models corresponding to the truncated constructs, showed a distinctive change on the surface charge distribution. This parameter was more pronounced in SmCys -20, which also displayed a significant displacement of the inhibitory domain, a result which could explain the kinetic data in terms of the loss of attachment sites. These changes correlated well with the progressive lack of inhibition observed for the recombinant deletion constructs, in vitro and in vivo.

Amino Acid Sequence↗

Structural analysis of a biologically active echistatin analogue des(46-49)-[Ala8,37]-echistatin gamma with three disulfide bonds by 2D-NMR and computer graphics.

An echistatin analogue, designated as des(46-49)-[Ala8,37]-echistatin gamma, was synthesized chemically by solid-phase peptide synthesis. The analogue was made by replacing Cys8 and Cys37 residues with two alanines and the deletion of C-terminal peptide 46-49 of echistatin gamma, resulting in an artificial polypeptide of 45 amino acids with three disulfide bonds. In the platelet aggregation assay, the analogue exhibits almost the same activity as echistatin gamma, indicating that the linear sequence of des(46-49)-[Ala8,37]-echistatin gamma contains all of the primary-structure information that is required for proper folding of this synthetic polypeptide. The tertiary structure of the analogue, as determined from high-resolution nuclear magnetic resonance (NMR) coupled with dynamic simulated annealing, is very similar to that of echistatin alpha1 which differs from echistatin gamma by 8 residues. In particular the two important sites of the Arg-Gly-Asp (RGD) loop and the C-terminal Lys45, both of which show some degree of disorder, are maintained in similar spatial orientation and proximity as those in echistatin alpha 1 even without the constraint provided by the disulfide bond of the (Cys8-Cys37) pair. These results provide new insights in further defining distinct structural features of echistatin gamma, which are involved in supporting the active polypeptide conformation to achieve biological activity in the absence of one pair of disulfide bonds.

Amino Acid Sequence↗

Genetic network inference: from co-expression clustering to reverse engineering.

MOTIVATION: Advances in molecular biological, analytical and computational technologies are enabling us to systematically investigate the complex molecular processes underlying biological systems. In particular, using high-throughput gene expression assays, we are able to measure the output of the gene regulatory network. We aim here to review datamining and modeling approaches for conceptualizing and unraveling the functional relationships implicit in these datasets. Clustering of co-expression profiles allows us to infer shared regulatory inputs and functional pathways. We discuss various aspects of clustering, ranging from distance measures to clustering algorithms and multiple-cluster memberships. More advanced analysis aims to infer causal connections between genes directly, i.e. who is regulating whom and how. We discuss several approaches to the problem of reverse engineering of genetic networks, from discrete Boolean networks, to continuous linear and non-linear models. We conclude that the combination of predictive modeling with systematic experimental verification will be required to gain a deeper insight into living organisms, therapeutic targeting and bioengineering.

Animals↗

Glucocorticoid receptor interactions with glucocorticoids: evaluation by molecular modeling and functional analysis of glucocorticoid receptor mutants.

In the treatment of inflammatory skin diseases, there are some glucocorticoid (GC) double esters combining pronounced antiinflammatory activity and minor atrophogenic side effects. The reason, however, is only poorly understood. To investigate interactions of GCs with the ligand-binding domain of the glucocorticoid receptor (GR), we measured receptor-binding potency of a series of GC esters including their metabolites and performed a molecular modeling study using progesterone receptor crystal structure data. Ligand docking to the GR-binding pocket showed good fitting of GC 17-esters corresponding to their high receptor-binding affinity, and unfavorable sterical interactions for GC 21-esters with substituents larger than propionate. Molecular dynamics simulations served to visualize induced fit procedures. Ligand docked GC conformations after dynamics simulations were used for generation of a 3D quantitative structure-activity relationship model. Using a set of 11 steroids, this model showed a correlation coefficient (r(2)) of 0.98, a leave-one-out cross validation (q(2)) of 0.79 and was able to predict binding affinity of further six ligands with a standard error of prediction of 0.33. Moreover, interactions of Asn-564 and Met-639 with the steroids were investigated by studying GR mutants of these amino acids. Met-639 participates in hydrophobic interactions mainly with GC side chains, while Asn-564 forms a hydrogen bond to the C11-OH group of the steroid. Asn-564 is shown to be very important for ligand binding and even more for target gene activation and transcription factor repression.

Amino Acid Substitution↗

Analysis of structure-activity relationships for MT2 selective antagonists by melatonin MT1 and MT2 receptor models.

Three-dimensional homology models of human MT(1) and MT(2) melatonin receptors were built with the aim to investigate the structure-activity relationships (SARs) of MT(2) selective antagonists. A common interaction pattern was proposed for a series of structurally different MT(2) selective antagonists, which were positioned within the binding site by docking and simulated annealing. The proposed antagonist binding mode to the MT(2) receptor is characterized by the accommodation of the out-of-plane substituents in a hydrophobic pocket, which resulted as being fundamental for the explanation of the antagonist behavior and the MT(2) receptor selectivity. Moreover, to assess the ability of the MT(2) receptor model to reproduce the SARs of MT(2) antagonists, three new derivatives of the MT(2) selective antagonist N-[1-(4-chloro-benzyl)-4-methoxy-1H-indol-2-ylmethyl]-propionamide (7) were synthesized and tested for their receptor affinity and intrinsic activity. These compounds were docked into the MT(2) receptor model and were submitted to molecular dynamics studies, providing results in qualitative agreement with the experimental data. These results confirm the importance of the out-of-plane group in receptor binding and selectivity and provide a partial validation of the proposed G protein-coupled receptor model.

Amino Acid Sequence↗

Statistical kinetics of macromolecular dynamics.

Fluctuations in biochemical processes can provide insights into the underlying kinetics beyond what can be gleaned from studies of average rates alone. Historically, analysis of fluctuating transmembrane currents supplied information about ion channel conductance states and lifetimes before single-channel recording techniques emerged. More recently, fluctuation analysis has helped to define mechanochemical pathways and coupling ratios for the motor protein kinesin as well as to probe the contributions of static and dynamic disorder to the kinetics of single enzymes. As growing numbers of assays are developed for enzymatic or folding behaviors of single macromolecules, the range of applications for fluctuation analysis increases. To evaluate specific biochemical models against experimental data, one needs to predict analytically the distribution of times required for completion of each reaction pathway. Unfortunately, using traditional methods, such calculations can be challenging for pathways of even modest complexity. Here, we derive an exact expression for the distribution of completion times for an arbitrary pathway with a finite number of states, using a recursive method to solve algebraically for the appropriate moment-generating function. To facilitate comparisons with experiments on processive motor proteins, we develop a theoretical formalism for the randomness parameter, a dimensionless measure of the variance in motor output. We derive the randomness for motors that take steps of variable sizes or that move on heterogeneous substrates, and then discuss possible applications to enzymes such as RNA polymerase, which transcribes varying DNA sequences, and to myosin V and cytoplasmic dynein, which may advance by variable increments.

Algorithms↗

A fractional programming approach to efficient DNA melting temperature calculation.

MOTIVATION: In a wide range of experimental techniques in biology, there is a need for an efficient method to calculate the melting temperature of pairings of two single DNA strands. Avoiding cross-hybridization when choosing primers for the polymerase chain reaction or selecting probes for large-scale DNA assays are examples where the exact determination of melting temperatures is important. Beyond being exact, the method has to be efficient, as these techniques often require the simultaneous calculation of melting temperatures of up to millions of possible pairings. The problem is to simultaneously determine the most stable alignment of two sequences, including potential loops and bulges, and calculate the corresponding melting temperature. RESULTS: As the melting temperature can be expressed as a fraction in terms of enthalpy and entropy differences of the corresponding annealing reaction, we propose to use a fractional programming algorithm, the Dinkelbach algorithm, to solve the problem. To calculate the required differences of enthalpy and entropy, the Nearest Neighbor model is applied. Using this model, the substeps of the Dinkelbach algorithm in our problem setting turn out to be calculations of alignments which optimize an additive score function. Thus, the usual dynamic programming techniques can be applied. The result is an efficient algorithm to determine melting temperatures of two DNA strands, suitable for large-scale applications such as primer or probe design. AVAILABILITY: The software is available for academic purposes from the authors. A web interface is provided at http://www.zaik.uni-koeln.de/bioinformatik/fptm.html

Algorithms↗

Procedures for reliable estimation of viral fitness from time-series data.

In order to develop a better understanding of the evolutionary dynamics of HIV drug resistance, it is necessary to quantify accurately the in vivo fitness costs of resistance mutations. However, the reliable estimation of such fitness costs is riddled with both theoretical and experimental difficulties. Experimental fitness assays typically suffer from the shortcoming that they are based on in vitro data. Fitness estimates based on the mathematical analysis of in vivo data, however, are often questionable because the underlying assumptions are not fulfilled. In particular, the assumption that the replication rate of the virus population is constant in time is frequently grossly violated. By extending recent work of Marée and colleagues, we present here a new approach that corrects for time-dependent viral replication in time-series data for growth competition of mutants. This approach allows a reliable estimation of the relative replicative capacity (with confidence intervals) of two competing virus variants growing within the same patient, using longitudinal data for the total plasma virus load, the relative frequency of the two variants and the death rate of infected cells. We assess the accuracy of our method using computer-generated data. An implementation of the developed method is freely accessible on the Web (http://www.eco.ethz.ch/fitness.html).

Anti-HIV Agents↗

Activity-structure correlations in divergent lectin evolution: fine specificity of chicken galectin CG-14 and computational analysis of flexible ligand docking for CG-14 and the closely related CG-16.

Gene duplication and sequence divergence are driving forces toward establishing protein families. To examine how sequence changes affect carbohydrate specificity, the two closely related proto-type chicken galectins CG-14 and CG-16 were selected as models. Binding properties were analyzed using a highly sensitive solid-phase assay. We tested 56 free saccharides and 34 well-defined glycoproteins. The two galectins share preference for the II (Galbeta1-4GlcNAc) versus I (Galbeta1-3GlcNAc) version of beta-galactosides. A pronounced difference is found owing to the reactivity of CG-14 with histo-blood group ABH active oligosaccharides and A/B active glycoproteins. These experimental results prompted to determine activity-structure correlations by modeling. Computational analysis included consideration of the flexibility of binding partners and the presence of water molecules. It provided a comparative description of complete carbohydrate recognition domains, which had so far not been characterized in animal galectins. The structural models assigned II, I selectivity to a region downstream of the central Trp moiety. Docking revealed that the tetrasaccharides can be accommodated in their free-state low-energy conformations. CG-14's preference for A versus B epitopes could be attributed to a contact between His124 and the N-acetyl group of GalNAc. Regarding intergalectin comparison, the Ala53/Cys51 exchange affects the interaction potential of His54/His52. Close inspection of simulated dynamic interplay revealed reorientation of His124 at the site of the His124/Glu123 substitution, with potential impact on ligand dissociation. In summary, this study identifies activity differences and provides information on their relation to structural divergence, epitomizing the value of this combined approach beyond galectins.

Amino Acid Sequence↗

The Ig superfamily cell adhesion molecule, apCAM, mediates growth cone steering by substrate-cytoskeletal coupling.

Dynamic cytoskeletal rearrangements are involved in neuronal growth cone motility and guidance. To investigate how cell surface receptors translate guidance cue recognition into these cytoskeletal changes, we developed a novel in vitro assay where beads, coated with antibodies to the immunoglobulin superfamily cell adhesion molecule apCAM or with purified native apCAM, replaced cellular substrates. These beads associated with retrograde F-actin flow, but in contrast to previous studies, were then physically restrained with a microneedle to simulate interactions with noncompliant cellular substrates. After a latency period of approximately 10 min, we observed an abrupt increase in bead-restraining tension accompanied by direct extension of the microtubule-rich central domain toward sites of apCAM bead binding. Most importantly, we found that retrograde F-actin flow was attenuated only after restraining tension had increased and only in the bead interaction axis where preferential microtubule extension occurred. These cytoskeletal and structural changes are very similar to those reported for growth cone interactions with physiological targets. Immunolocalization using an antibody against the cytoplasmic domain of apCAM revealed accumulation of the transmembrane isoform of apCAM around bead-binding sites. Our results provide direct evidence for a mechanical continuum from apCAM bead substrates through the peripheral domain to the central cytoplasmic domain. By modulating functional linkage to the underlying actin cytoskeleton, cell surface receptors such as apCAM appear to enable the application of tensioning forces to extracellular substrates, providing a mechanism for transducing retrograde flow into guided growth cone movement.

Actins↗

Study of disabling T-cell activation and inhibiting T-cell-mediated immunopathology reveals a possible inverse agonist activity of CD4 peptidomimetics.

We designed a new class of aromatically modified exocyclic peptides based on the structure of CD4 by engineering one of the cysteine residues in a peptidomimetic derived from the CDR3 region of the CD4 molecule. All three species mediate inhibition of T-cell proliferation at concentrations ranging from 10 to 100 microM. The mimetics CD4-Cys and CD4-Met bind to sCD4 with affinities ranging from 1 to 2 microM, while CD4-Ser shows poor binding in radioisotope assay. Though these mimetics have similar structures, they exhibit different biochemical and biological functions. Activation of T-cells as measured by thymidine incorporation or IL-2 production revealed that CD4-Cys and CD4-Ser mimetics behave as classical antagonists. On the other hand, the CD4-Met species inhibited T-cell proliferation with an IC(50) of 30 microM but unexpectedly increased IL-2 secretion modestly at a less than 3 microM concentration. In experimental autoimmune encephalitis (EAE), CD4-Ser and CD4-Cys mimetics reduced the severity of EAE symptoms while the CD4-Met mimetic exacerbated the conditions. We propose that CD4-Cys and CD4-Ser are classical antagonists, but CD4-Met may possess properties of an inverse agonist. The structure-activity relationship of mimetics reveals that a minor change in the net hydropathic value is enough to alter the dynamic nature of the receptor-ligand complex.

Animals↗

A new combined computational and NMR-spectroscopical strategy for the identification of additional conformational constraints of the bound ligand in an aprotic solvent.

This study documents the feasibility of switching to an aprotic medium in sugar receptor research. The solvent change offers additional insights into mechanistic details of receptor--carbohydrate ligand interactions. If a receptor retained binding capacity in an aprotic medium, solvent-exchangeable protons of the ligand would not undergo transfer and could act as additional sensors, thus improving the level of reliability in conformational analysis. To probe this possibility, we first focused on hevein, the smallest lectin found in nature. The NMR-spectroscopic measurements verified complexation, albeit with progressively reduced affinity by more than 1.5 orders of magnitude, in mixtures of up to 50% dimethyl sulfoxide (DMSO). Since hevein lacks the compact beta-strand arrangement of other sugar receptors, such a structural motif may confer enhanced resistance to solvent exchange. Two settings of solid-phase activity assays proved this assumption for three types of alpha- and/or beta-galactoside-binding proteins, that is, a human immunoglobulin G (IgG) subfraction, the mistletoe lectin, and a member of the galectin family of animal lectins. Computer-assisted calculations and NMR experiments also revealed no conspicuous impact of the solvent on the conformational properties of the tested ligands. To define all possible nuclear Overhauser effect (NOE) contacts in a certain conformation and to predict involvement of exchangeable protons, we established a new screening protocol applicable during a given molecular dynamics (MD) trajectory and calculated population densities of distinct contacts. Experimentally, transferred NOE (tr-NOE) experiments with IgG molecules and the disaccharide Gal'alpha1-3Galbeta1-R in DMSO as solvent disclosed that such an additional crosspeak, that is, Gal'OH2--GalOH4, was even detectable for the bound ligand under conditions in which spin diffusion effects are suppressed. Further measurements with the plant lectin and galectins confirmed line broadening of ligand signals and gave access to characteristic crosspeaks in the aprotic solvent and its mixtures with water. Our combined biochemical, computational, and NMR-spectroscopical strategy is expected to contribute notably to the precise elucidation of the geometry of ligands bound to compactly folded sugar receptors and of the role of water molecules in protein--ligand (carbohydrate) recognition, with relevance to areas beyond the glycosciences.

Antimicrobial Cationic Peptides↗

Xyloglucan sidechains modulate binding to cellulose during in vitro binding assays as predicted by conformational dynamics simulations.

Cross-links between cellulose microfibrils and xyloglucan (XG) molecules play a major role in defining the structural properties of plant cell walls and the regulation of growth and development of dicotyledonous plants. How these cross-links are established and how they are regulated has yet to be determined. In a previous study, preliminary data were presented which suggested that the different sidechains of XG may play a role in controlling cellulose microfibril-XG interactions. In this study, this question is addressed directly by analyzing to what extent the different sidechains of pea cell wall XG and nasturtium seed storage XG affect their binding to cellulose microfibrils. Of particular importance to this study are the chemical data indicating that pea XG possesses a trisaccharide sidechain, which is not found in nasturtium XG. To this end, conformational dynamic simulations have been used to predict whether oligosaccharides representative of pea and nasturtium XG can adopt a hypothesized cellulose-binding conformation and which of these XGs exhibits a preferential ability to bind cellulose. Extensive analysis of the conformational forms populated during 300 K and high-temperature Monte Carlo simulations established that a planar, sterically accessible, glucan backbone is essential for optimal cellulose-binding. For the trisaccharide sidechain-containing oligosaccharide as found in pea XG, sidechain orientation appeared to regulate the gradual acquisition of this hypothesized cellulose binding conformation. Thus, conformational forms were identified that included the twisted backbone (non-planar) putative solution form of XG, forms in which the trisaccharide sidechain orientation enables increased backbone planarity and steric accessibility, and finally a planar, sterically accessible, backbone. By applying these conformational requirements for cellulose binding, it has been determined that pea XG possesses a two- to threefold occurrence of the cellulose binding conformation than nasturtium XG. Based on this finding, it was predicted that pea XG would bind to cellulose at a higher rate than nasturtium XG. In vitro binding assays showed that pea XG-avicel binding does indeed occur at a twofold higher rate than nasturtium XG-avicel binding. The enhanced ability of pea cell wall XG over nasturtium seed storage XG to associate with cellulose is consistent with a structural role of the former during epicotyl growth where efficient association with cellulose is a requirement. In contrast, the relatively low ability of nasturtium XG to bind cellulose is consistent with the need to enhance the accessibility of this polymer to glycanases during germination. These findings suggest potential roles for XG sidechain substitution, enabling XG to function in a variety of different biological contexts.

Binding Sites↗

Identification of the binding sites and selectivity of sarpogrelate, a novel 5-HT2 antagonist, to human 5-HT2A, 5-HT2B and 5-HT2C receptor subtypes by molecular modeling.

The aim of the present study was to investigate the binding sites interactions and the selectivity of sarpogrelate to human 5-HT(2) receptor family (5-HT(2A), 5-HT(2B) and 5-HT(2C) receptor subtypes) using molecular modeling. Rhodopsin (RH) crystal structures were used as template to build structural models of the human serotonin-2A and -2C receptors (5-HT(2A)R, 5-HT(2C)R), whereas for 5-HT(2B)R, we used our previously published three-dimensional (3D) models based on bacteriorhodopsin (BR). Sarpogrelate, a novel 5-HT(2)R antagonist, was docked to the receptors. Molecular dynamics (MD) simulations produced the strongest interaction for 5-HT(2A)R/sarpogrelate complex. Upon binding, sarpogrelate constraints aromatic residues network (Trp(3.28), Phe(5.47), Trp(6.48), Phe(6.51), Phe(6.52) in 5-HT(2A)R; Phe(3.35), Phe(6.51), Trp(7.40) in 5-HT(2B)R; Trp(3.28), Phe(3.35), Phe(5.47), Trp(6.48), Phe(6.51), Phe(6.52) in 5-HT(2C)R) in a stacked configuration, preventing activation of the receptor. The models suggest that the structural origin of the selectivity of sarpogrelate to 5-HT(2A)R vs both 5-HT(2B)R and 5-HT(2C)R comes from the following results: (1) The tight interaction between the antagonist and the transmembrane domain (TMD) 3. Asp(3.32) neutralizes the cationic head and interacts simultaneously with carboxylic group hydrogen of the antagonist molecule. (2) Due to steric hindrance, Ser(5.46) (vs Ala(5.46) in 5HT(2B) and 5HT(2C)) prevents sarpogrelate to enter deeply inside the hydrophobic core of the helix bundle and to interact with Pro(5.50). (3) The side chain of Ile(4.56) (vs Ile(4.56) in 5HT(2B)R and Val(4.56) in 5HT(2C)R) constraints sarpogrelate to adjust its position by translating toward the strongly attractive Asp(3.32). These results are in good agreement with binding affinities (pKi) of sarpogrelate for 5-HT(2) receptor family expressed in transfected cell.

Amino Acid Sequence↗

A critical role for a tyrosine residue in the cannabinoid receptors for ligand recognition.

Previous mutation and modeling studies have identified an aromatic cluster in the transmembrane helix (TMH) 3-4-5 region as important for ligand binding at the CB(1) and CB(2) cannabinoid receptors. Through novel mixed mode Monte Carlo/Stochastic Dynamics (MC/SD) calculations, we tested the importance of aromaticity at position 5.39(275) in CB(1). MC/SD calculations were performed on wild-type (WT) CB(1) and two mutants, Y5.39(275)F and Y5.39(275)I. Results indicated that while the CB(1) Y5.39(275)F mutant is very similar to WT, the Y5.39(275)I mutant shows pronounced topology changes in the TMH 3-4-5 region. Site-directed mutagenesis studies of tyrosine 5.39 to phenylalanine (Y-->F) or isoleucine (Y-->I) in both CB(1) and CB(2) were performed to determine the functional role of this amino acid in each receptor subtype. HEK 293 cells transfected with mutant receptor cDNAs were evaluated in radioligand binding and cyclic AMP assays. The CB(1) mutant and WT receptors were also co-expressed with G-protein-coupled inwardly rectifying channels (GIRK1 and GIRK4) in Xenopus oocytes to assess functional coupling. The Y-->F mutation resulted in cannnabinoid receptors with subtle differences in WT binding and signal transduction. In contrast, the Y-->I mutations produced receptors that could not produce signal transduction or bind to multiple cannabinoid compounds. However, immunofluorescence data indicate that the Y-->I mutation was compartmentalized and expressed at a level similar to that of the WT cannabinoid receptor. These results underscore the importance of aromaticity at position CB(1) 5.39(275) and CB(2) 5.39(191) for ligand recognition in the cannabinoid receptors.

Amino Acid Sequence↗