Assessment of thermal effects in a model of the human head implanted with a wireless active microvalve for the treatment of glaucoma creating a filtering bleb.

A 3D anatomical computational model is developed to assess thermal effects due to exposure to the electromagnetic field required to power a new investigational active implantable microvalve for the treatment of glaucoma. Such a device, located in the temporal superior eye quadrant, produces a filtering bleb, which is included in the geometry of the model, together with the relevant ocular structures. The electromagnetic field source-a planar coil-as well as the microvalve antenna and casing are also included. Exposure to the electromagnetic field source of an implanted and a non-implanted subject are simulated by solving a magnetic potential formulation, using the finite element method. The maximum SAR10 is reached in the eyebrow and remains within the limits suggested by the IEEE and ICNIRP standards. The anterior chamber, filtering bleb, iris and ciliary body are the ocular structures where more absorption occurs. The temperature rise distribution is also obtained by solving the bioheat equation with the finite element method. The numerical results are compared with the in vivo measurements obtained from four rabbits implanted with the microvalve and exposed to the electromagnetic field source.

A Model Chain Application to Estimate Mixing Layer Height Related to PM10 Dispersion Processes.

The mixing layer height (MLH) is a crucial parameter in order to investigate the near surface concentrations of air pollutants. The MLH can be estimated by measurements of some atmospheric variables, by indirect estimates based on trace gases concentration or aerosol, or by numerical models. Here, a modelling approach is proposed. The developed modelling system is based on the models WRF-ARW and CALMET. This system is applied on Firenze-Prato-Pistoia area (Central Italy), during 2010, and it is compared with in situ measurements. The aim of this work is to evaluate the use of MLH model estimates to characterize the critical episodes for PM10 in a limited area. In order to find out the meteorological conditions predisposing accumulation of PM10 in the atmosphere's lower level, some indicators are used: daily mean wind speed, cumulated rainfall, and mean MLH estimates from CALMET model. This indicator is linked to orography, which has important consequences on local weather dynamics. However, during critical events the local emission sources are crucial to the determination of threshold exceeding of PM10. Results show that the modelled MLH, together with cumulative rainfall and wind speed, can identify the meteorological conditions predisposing accumulation of air pollutant at ground level.

Characterisation and simulation of an active microvalve for glaucoma.

Glaucoma drainage device (GDD) has the potential to eliminate hypotony but still suffers from poor flow control and fibrosis. The ideal shunt should change its hydraulic resistance to achieve the desired intraocular pressure (IOP). In this study, the characterisation of a preliminary design of a new GDD is presented. This is activated by means of a diaphragm, which is actuated by conducting polymers. The valve can be manufactured employing microelectromechanical system technology by soft lithography. The characterisation process is performed by numerical simulation using the finite element method, considering the coupling between the fluid and the structure (diaphragm) obtaining the hydraulic resistance for several positions of the diaphragm. To analyse the hydraulic system of the microvalve implanted in a human eye, an equivalent circuit model was used. The parameters of the equivalent circuit model were obtained from numerical simulation. The hydraulic resistance of the designed GDD varies in the range of 13.08-0.36 mmHg min/µl compared with 3.38-0.43 mmHg min/µl for the Ahmed valve. The maximum displacement of the diaphragm in the vertical direction is 18.9 µm, and the strain in the plane is 2%. The proposed preliminary design allows to control the IOP by varying the hydraulic resistance in a greater range than the existing passive valves, and the numerical simulation facilitates the characterisation and the improvement of the design before its construction, reducing time and costs.

[Giant and bilateral adrenal myelolipoma. Case report]. / Mielolipoma surrenalico bilaterale gigante. Caso clinico.

The adrenal myelolipoma is a relatively rare benign tumour of adipose cell and bone marrow elements, non functioning and asymptomatic. Giant and bilateral adrenal myelolipoma is quite rare. The Authors report a case of bilateral adrenal myelolipoma, a giant one (> 15 cm) on the left side and a small one (> 4 cm) on the right with constant pain in a 57-year-old man, shown by computerized tomography. The hormonal blood tests were normal. Surgical excision was performed for large left symptomatic mass, by open laparotomy, and biopsy for right minor adrenal lesion. Histology confirmed diagnosis of myelolipoma for both masses. Follow-up to 6-12 months did'nt show any change of the right myelolipoma. The authors agree with the need to remove the giant adrenal myelolipoma, because the lesion > 10 cm have a high risk of cancer and hemorrhagic complication, while for small myelolipoma (< 6 cm) 6-12 months follow-up is the appropriate choice.

[Acetate free biofiltration (AFB) and ionic dialysance]. / AFB e dialysance ionica.

Among various methods and models to evaluate the dialysis dose, the on-line measure of ionic dialysance gives good results. Compared to measured KT/V, ionic dialysance does not require blood sampling and, consequently, it reduces the risks and the costs. Functional evaluations during acetate free biofiltration (AFB treatments) confirm the usefulness of the on-line continuous measure by biosensors, especially when they are integrated in the dialysis machinery.

Activation of the cannabinoid CB1 receptor may involve a W6 48/F3 36 rotamer toggle switch.

The cannabinoid CB1 receptor, a member of the Rhodopsin (Rho) family of G protein coupled receptors (GPCRs), exhibits high levels of constitutive activity. In contrast, Rho exhibits an exquisite lack of constitutive activity. In Rho, W6.48(265) on transmembrane helix 6 (TMH6) is flanked by aromatic residues at positions i-4 (F6.44) and i + 3 (Y6.51), while in CB1 the residues i-4 and i + 3 to W6.48 are leucines (L6.44 and L6.51). Based upon spectroscopic evidence, W6.48 has been proposed to undergo a rotamer switch (chi1 g+ -->trans) upon activation of Rho. In the work reported here, the biased Monte Carlo method, Conformational Memories (CM) was used to test the hypothesis that the high constitutive activity exhibited by CB1 may be due, in part, to the lack of aromatic residues i-4 and i + 3 from W6.48. In this work, the W6.48 rotamer shift (chi1 g+ -->trans) was used as the criterion for activation. Conformational Memories (CM) calculations on WT CB1 TMH6 and L6.44F and L6.51Y mutant TMH6s revealed that an aromatic residue at 6.44 tends to disfavor the W6.48 chi1 g+ -->trans transition and an aromatic residue at 6.51 would require a concomitant movement of the Y6.51 chi1 from trans-->g+ when the W6.48 chi1 undergoes a g+ -->trans shift. In contrast, CM calculations on WT CB1 TMH6 revealed that the presence of leucines at 6.44 and 6.51 provide W6.48 with greater conformational mobility, with a W6.48 transchi1 preferred. Conformational Memories calculations also revealed that the W6.48 chi1 g+ -->trans transition in WT CB1 TMH6 is correlated with the degree of kinking in TMH6. The average proline kink angles for TMH6 were higher for helices with a W6.48 g+ chi1 than for those with a W6.48 transchi1. These results are consistent with experimental evidence that TMH6 straightens during activation. Transmembrane helix (TMH) bundle models of the inactive (R) and active (R*) states of CB1 were then probed for interactions that may constrain W6.48 in the inactive state of CB1. These studies revealed that F3.36 (transchi1) helps to constrain W6.48 in a g+ chi1 in the inactive (R) state of CB1. In the R* state, these studies suggest that F3.36 must assume a g+ chi1 in order to allow W6.48 to shift to a transchi1. These results suggest that the W6.48/F3.36 interaction may act as the 'toggle switch' for CB1 activation, with W6.48 chi1 g+/F3.36 chi1 trans representing the inactive (R) and W6.48 chi1 trans/F3.36 chi1 g+ representing the active (R*) state of CB1.

Agonist-induced conformational changes at the cytoplasmic side of transmembrane segment 6 in the beta 2 adrenergic receptor mapped by site-selective fluorescent labeling.

The environmentally sensitive, sulfhydryl-reactive, fluorescent probe N,N'-dimethyl-N-(iodoacetyl)-N'-(7-nitrobenz-2-oxa-1,3-diazol-4-yl) ethylene-diamine (IANBD) was used as a molecular reporter of agonist-induced conformational changes in the beta(2) adrenergic receptor, a prototype hormone-activated G protein-coupled receptor. In the background of a mutant beta(2) adrenergic receptor, with a minimal number of endogenous cysteine residues, new cysteines were introduced in positions 269(6.31), 270(6.32), 271(6.33), and 272(6.34) at the cytoplasmic side of transmembrane segment (TM) 6. The resulting mutant receptors were fully functional and bound both agonists and antagonist with high affinities also upon IANBD labeling. Fluorescence spectroscopy analysis of the purified and site-selectively IANBD-labeled mutants suggested that the covalently attached fluorophore was exposed to a less polar environment at all four positions upon agonist binding. Whereas evidence for only a minor change in the molecular environment was obtained for positions 269(6.31) and 270(6.32), the full agonist isoproterenol caused clear dose-dependent and reversible increases in fluorescence emission at positions 271(6.33) and 272(6.34). The data suggest that activation of G protein-coupled receptors, which are activated by "diffusible" ligands, involves a structural rearrangement corresponding to the cytoplasmic part of TM 6. The preferred conformations of the IANBD moiety attached to the inserted cysteines were predicted by employing a computational method that incorporated the complex hydrophobic/hydrophilic environment in which the cysteines reside. Based on these preferred conformations, it is suggested that the spectral changes reflect an agonist-promoted movement of the cytoplasmic part of TM 6 away from the receptor core and upwards toward the membrane bilayer.

Targeting antigen-specific T cells by genetically engineered antigen presenting cells. A strategy for specific immunotherapy of autoimmune disease.

We describe a strategy for specific immunotherapy of autoimmune disease based on targeting the antigen-specific T cells in an experimental model of myasthenia gravis. To address the problem of heterogeneity of the T cell repertoire, we have genetically engineered antigen presenting cells (APCs) to process and present epitopes of the autoantigen, acetylcholine receptor (AChR), to the entire spectrum of AChR-specific syngeneic T cells. APCs derived from BALB/c mice were stably transfected with cDNA for the key immunogenic domain of the AChR alpha-subunit, flanked by sequences of the lysosome-associated membrane protein (LAMP) that direct APCs to process and present the antigen via the MHC Class II pathway. Transfected APCs strongly stimulated AChR-specific T cells from BALB/c mice. Fas ligand, or antibody to Fas, abrogated the T cell response, by inducing apoptosis of the APC-stimulated T cells. The new results of this investigation are (1) that autoreactive T cells can be effectively targeted by autologous APCs that are engineered to present the relevant autoantigen, and (2) that these specifically targeted and activated T cells can be profoundly inhibited by agents that trigger the Fas-mediated apoptosis pathway. The present findings suggest that engineering APCs for simultaneous presentation of the autoantigen and delivery of FasL will provide a powerful strategy for the elimination of autoreactive T cells.

A self-consistent, microenvironment modulated screened coulomb potential approximation to calculate pH-dependent electrostatic effects in proteins.

An improved approach is presented for calculating pH-dependent electrostatic effects in proteins using sigmoidally screened Coulomb potentials (SCP). It is hypothesized that a key determinant of seemingly aberrant behavior in pKa shifts is due to the properties of the unique microenvironment around each residue. To help demonstrate this proposal, an approach is developed to characterize the microenvironments using the local hydrophobicity/hydrophilicity around each residue of the protein. The quantitative characterization of the microenvironments shows that the protein is a complex mosaic of differing dielectric regions that provides a physical basis for modifying the dielectric screening functions: in more hydrophobic microenvironments the screening decreases whereas the converse applies to more hydrophilic regions. The approach was applied to seven proteins providing more than 100 measured pKa values and yielded a root mean square deviation of 0.5 between calculated and experimental values. The incorporation of the local hydrophobicity characteristics into the algorithm allowed the resolution of some of the more intractable problems in the calculation of pKa. Thus, the divergent shifts of the pKa of Glu-35 and Asp-66 in hen egg white lysozyme, which are both about 90% buried, was correctly predicted. Mechanistically, the divergence occurs because Glu-35 is in a hydrophobic microenvironment, while Asp-66 is in a hydrophilic microenvironment. Furthermore, because the calculation of the microenvironmental effects takes very little CPU time, the computational speed of the SCP formulation is conserved. Finally, results from different crystal structures of a given protein were compared, and it is shown that the reliability of the calculated pKa values is sufficient to allow identification of conformations that may be more relevant for the solution structure.

Exploration of biologically relevant conformations of anandamide, 2-arachidonylglycerol, and their analogues using conformational memories.

The endogenous cannabinoid anandamide (N-arachidonoylethanolamide) has been shown to possess higher affinity for the cannabinoid CB1 receptor than for the CB2 receptor. Carrier-mediated transport has been proposed to be essential for the termination of the biological effects of anandamide, while hydrolysis of anandamide is performed by a membrane-bound amidohydrolase, fatty acid amidohydrolase (FAAH). As interaction of anandamide with each of these targets occurs in different environments, the conformations of anandamide for interaction with each target may differ. To ascertain what conformations of anandamide, a highly flexible molecule, are favored in polar and nonpolar environments, the new method of Conformational Memories (CM) was used. CM has been shown to achieve complete conformational sampling of highly flexible ligands, to converge in a very practical number of steps, and to be capable of overcoming energy barriers very efficiently (Guarnieri et al. J. Am. Chem. Soc. 1996, 118, 5580). The generalized Born/surface area (GB/SA) continuum solvation models for chloroform and for water were used in the CM calculations. As a means of validation, CM was first applied to arachidonic acid because both experimental and theoretical conformational studies of arachidonic acid have been reported. CM was also applied to sn-2-arachidonylglycerol (2-AG), another endogenous CB ligand; to a 1,1-dimethylheptyl derivative of anandamide, an analogue with higher CB1 affinity than anandamide; and to N-(2-hydroxyethyl)prostaglandin-B2-ethanolamide (PGB2-EA), a prostanoid ligand which does not bind to CB1. Consistent with the literature, arachidonic acid was found to exist in an extended, angle-iron shape and in back-folded conformations which were U, J, or helical in shape. The angle-iron and U-shapes were both highly populated conformations with the angle-iron preferred in CHCl3 and the U-shape preferred in H2O. Results for anandamide and 2-AG paralleled those for arachidonic acid with the exception that anandamide in water does not adopt a pure extended conformation but, rather, favors a hybrid-extended/U-shape. For the dimethyl-heptyl derivative of anandamide, the U-shape was found to be predominant in both environments (42% in CHCl3, 38% in H2O), but the population of the angle-iron shape was still significant (25% in CHCl3, 29% in H2O). For all of these ligands, J-shaped conformers constituted from 7% to 17% of the conformer population, while the helical shape was less than 5%. In both CHCl3 and H2O, the presence of the five-membered ring attenuates the ability of PGB2-EA to adopt an extended conformation. PGB2-EA was found instead to exist predominantly in an L-shape (i.e., distorted U-shape). The low probability of PGB2-EA adopting an extended conformation may be why PGB2-EA shows such low affinity for the CB1 receptor. The conformational information obtained here for anandamide and 2-AG may be useful in the design of rigid analogues which mimic the preferred molecular conformations (shapes) of these ligands. Such rigid analogues may be useful in deducing the bioactive conformation of these endogenous cannabinoids, not only at the CB receptors but also at the FAAH enzyme active site and possibly at the binding site(s) on the newly proposed anandamide transporter.

A proposed structure for transmembrane segment 7 of G protein-coupled receptors incorporating an asn-Pro/Asp-Pro motif.

Transmembrane segment (TMS) 7 has been shown to play an important role in the signal transduction function of G-protein-coupled receptors (GPCRs). Although transmembrane segments are most likely to adopt a helical structure, results from a variety of experimental studies involving TMS 7 are inconsistent with it being an ideal alpha-helix. Using results from a search of the structure database and extensive simulated annealing Monte Carlo runs with the new Conformational Memories method, we have identified the conserved (N/D)PxxY region of TMS 7 as the major determinant for deviation of TMS 7 from ideal helicity. The perturbation consists of an Asx turn and a flexible "hinge" region. The Conformational Memories procedure yielded a model structure of TMS 7 which, unlike an ideal alpha-helix, is capable of accommodating all of the experimentally derived geometrical criteria for the interactions of TMS 7 in the transmembrane bundle of GPCRs. In the context of the entire structure of a transmembrane bundle model for the 5HT2a receptor, the specific perturbation of TMS 7 by the NP sequence suggests a structural hypothesis for the pattern of amino acid conservation observed in TMS 1, 2, and 7 of GPCRs. The structure resulting from the incorporation of the (N/D)P motif satisfies fully the H-bonding capabilities of the 100% conserved polar residues in these TMSs, in agreement with results from mutagenesis experiments. The flexibility introduced by the specific structural perturbation produced by the (NP/DP) motif in TMS 7 is proposed to have a significant role in receptor activation.

Functional microdomains in G-protein-coupled receptors. The conserved arginine-cage motif in the gonadotropin-releasing hormone receptor.

An Arg present in the third transmembrane domain of all rhodopsin-like G-protein-coupled receptors is required for efficient signal transduction. Mutation of this Arg in the gonadotropin-releasing hormone receptor to Gln, His, or Lys abolished or severely impaired agonist-stimulated inositol phosphate generation, consistent with Arg having a role in receptor activation. To investigate the contribution of the surrounding structural domain in the actions of the conserved Arg, an integrated microdomain modeling and mutagenesis approach has been utilized. Two conserved residues that constrain the Arg side chain to a limited number of conformations have been identified. In the inactive wild-type receptor, the Arg side chain is proposed to form an ionic interaction with Asp3.49(138). Experimental results for the Asp3. 49(138) --> Asn mutant receptor show a modestly enhanced receptor efficiency, consistent with the hypothesis that weakening the Asp3. 49(138)-Arg3.50(139) interaction by protonation of the Asp or by the mutation to Asn favors activation. With activation, the Asp3. 49(138)-Arg3.50(139) ionic bond would break, and the unrestrained Arg would be prevented from orienting itself toward the water phase by a steric clash with Ile3.54(143). The mutation Ile3.54(143) --> Ala, which eliminates this clash in simulations, causes a marked reduction in measured receptor signaling efficiency, implying that solvation of Arg3.50(139) prevents it from functioning in the activation of the receptor. These data are consistent with residues Asp3.49(138) and Ile3.54(143) forming a structural motif, which helps position Arg in its appropriate inactive and active receptor conformations.

Computer simulation studies of the fully solvated wild-type and mutated GnRH in extended and beta-turn conformations.

The conformational preference of the gonadotropin-releasing hormone (GnRH) and its Lys-8 mutant, studied earlier with a continuum model, was revisited using an explicit solvent model and thermodynamic integration to calculate the solvents contribution to the conformation-dependence of its free energy. In addition, the Proximity Criterion was used to further analyze the effects of conformational changes.

The enhanced immune response to the HIV gp160/LAMP chimeric gene product targeted to the lysosome membrane protein trafficking pathway.

The lysosome-associated membrane proteins (LAMP), found in the outer membrane of lysosomes and also in a multilaminar compartment that contains major histocompatibility complex class II (MHC II) proteins, are directed to their localization by a cytoplasmic carboxyl-terminal sequence. Our studies of the immune response to LAMP-targeted proteins has led to the application of a HIV-1 gp160/LAMP chimeric gene as a novel means to enhance the MHC II presentation of gp160. Immunofluorescence microscopy confirmed that the gp160/LAMP protein had a cellular localization corresponding to that of lysosomes. Pulse-chase analysis confirmed that the rates of synthesis of gp160/LAMP and wild type gp160 were comparable and that both proteins were processed to gp120 at similar rates. However, the gp160/LAMP was degraded more rapidly than the wild type gp160. MHC II-mediated T cell proliferation assays performed with cloned human cell lines showed that gp160/LAMP stimulated greater responses than did the wild type gp160. Moreover, mice vaccinated with recombinant vaccinia expressing gp160/LAMP had greater gp160-specific lymphoproliferation responses and higher titers of anti-V3 loop antibodies than mice vaccinated with recombinant vaccinia expressing wild type gp160.

A surface on the G protein beta-subunit involved in interactions with adenylyl cyclases.

Receptor activation of heterotrimeric G proteins dissociates G alpha from the G betagamma complex, allowing both to regulate effectors. Little is known about the effector-interaction regions of G betagamma. We had used molecular modeling to dock a peptide encoding the region of residues 956-982 of adenylyl cyclase (AC) 2 onto Gbeta to identify residues on Gbeta that may interact with effectors. Based on predictions from the model, we synthesized peptides encoding sequences of residues 86-105 (Gbeta 86-105) and 115-135 (Gbeta 115-135) from Gbeta. The Gbeta 86-105 peptide inhibited G betagamma stimulation of AC2 and blocked G betagamma inhibition of AC1 and by itself inhibited calmodulin-stimulated AC1, thus displaying partial agonist activity. Substitution of Met-101 with Asn in this peptide resulted in the loss of both the inhibitory and partial agonist activities. Most activities of the Gbeta 115-135 peptide were similar to those of Gbeta 86-105 but Gbeta 115-135 was less efficacious in blocking G betagamma inhibition of AC1. Substitution of Tyr-124 with Val in the Gbeta 115-135 peptide diminished all of its activities. These results identify the region encoded by amino acids 84-143 of Gbeta as a surface that is involved in transmitting signals to effectors.

Making DNA add.

Recent studies have demonstrated the feasibility of using DNA-based experiments to compute solutions to combinatorial problems. However, a prerequisite for designing a computer useful in a wide range of applications is the ability to perform mathematical calculations. The development of a DNA-based algorithm for addition is presented. The DNA representation of two nonnegative binary numbers is presented in a form permitting a chain of primer extension reactions to carry out the addition operation. To demonstrate the feasibility of this algorithm, a simple example was executed biochemically.

Restricted analogues provide evidence of a biologically active conformation of thyrotropin-releasing hormone.

Thyrotropin-releasing hormone (TRH) is a tripeptide (< Glu-His-Pro-NH2) that signals through a G protein-coupled receptor. TRH is a highly flexible molecule that can assume many conformations in solution. To attempt to delineate the biologically active conformation of TRH, we synthesized a pair of conformationally restricted cyclohexyl/Ala2-TRH analogues. The diastereomeric analogues use a lactam ring to restrict two of the six free torsional angles of TRH and constrain the X-Pro-NH2 peptide bond to trans. Unrestricted cyclohexyl/Ala2-TRH exhibited a 650-fold lower affinity than TRH for TRH receptor and was 430-fold less potent than TRH in stimulating inositol phosphate second messenger formation. One diastereomer exhibited higher affinity and potency than the unrestricted analogue despite the presence of the methylene bridge and fused ring, whereas the other showed lower affinity and potency. Computer simulations predicted that the positions of the cyclohexyl/Ala2 and Pro-NH2 moieties relative to < glutamate were different in the two analogues and that the conformation of the higher affinity analogue is different from that of trans-TRH in solution but is superimposable on that of trans-TRH found in a model of the TRH/TRH receptor complex. These experimental findings identify a favored relative position of < glutamate and Pro-NH2 in the more active conformation of two diastereomeric analogues of TRH and provide independent support for the model of the TRH/TRH receptor complex.

A refined model of the thyrotropin-releasing hormone (TRH) receptor binding pocket. Experimental analysis and energy minimization of the complex between TRH and TRH receptor.

Seven transmembrane (TM) spanning, G protein-coupled receptors (GPCRs) appear to bind large glycoprotein hormones predominantly within their extracellular domains, small nonpeptidic ligands within the TM helical bundle, and peptide ligands within the extracellular domains and TM bundle. The tripeptide thyrotropin-releasing hormone (TRH, pyroGlu-His-ProNH2) may bind entirely within the TM bundle of the TRH receptor (TRH-R). We have previously demonstrated direct binding contacts between the pyroGlu of TRH and two residues in TM helix 3 (TM-3) of TRH-R and proposed a model of the binding pocket of TRH-R [Perlman, J. H., Laakkonen, L., Osman, R., & Gershengorn, M. C. (1994) J. Biol. Chem. 269, 23383-23386]. Here, we provide evidence for two additional direct interactions between TRH and TRH-R. One interaction is between the aromatic ring of Tyr 282 of TM-6 and His of TRH. This is based on a large increase in the half-maximally effective concentration (EC50) of TRH for stimulation of inositol phosphate formation by Y282A TRH-R and a loss of selectivity of this mutant receptor for TRH analogs substituted at His. We provide evidence for another interaction between Arg 306 of TM-7 and the terminal carboxamide of TRH. Using four direct interactions as anchors, a refined model of the TRH-R binding pocket was constructed using geometry optimization through energy minimization. A novel method for modeling GPCRs based on Monte Carlo and stochastic dynamics simulations is presented in the accompanying paper [Laakkonen, L. J., Guarnieri, F., Perlman, J. H., Gershengorn, M. C., & Osman, R. (1996) Biochemistry 35, 7651-7663].

A refined model of the thyrotropin-releasing hormone (TRH) receptor binding pocket. Novel mixed mode Monte Carlo/stochastic dynamics simulations of the complex between TRH and TRH receptor.

Previous mutational and computational studies of the thyrotropin-releasing hormone (TRH) receptor identified several residues in its binding pocket [see accompanying paper, Perlman et al. (1996) Biochemistry 35, 7643-7650]. On the basis of the initial model constructed with standard energy minimization techniques, we have conducted 15 mixed mode Monte Carlo/stochastic dynamics (MC-SD) simulations to allow for extended sampling of the conformational states of the ligand and the receptor in the complex. A simulated annealing protocol was adopted in which the complex was cooled from 600 to 310 K in segments of 30 ps of the MC-SD simulations for each change of 100 K. Analysis of the simulation results demonstrated that the mixed mode MC-SD protocol maintained the desired temperature in the constant temperature simulation segments. The elevated temperature and the repeating simulations allowed for adequate sampling of the torsional space of the complex with successful conservation of the general structure and good helicity of the receptor. For the analysis of the interaction between TRH and the binding pocket, TRH was divided into four groups consisting of pyroGlu, His, ProNH2, and the backbone. The pairwise interaction energies of the four separate portions of TRH with the corresponding residues in the receptor provide a physicochemical basis for the understanding of ligand-receptor complexes. The interaction of pyroGlu with Tyr106 shows a bimodal distribution that represents two populations: one with a H-bond and another without it. Asp195 was shown to compete with pyroGlu for the H-bond to Tyr106. Simulations in which Asp195 was interacting with Arg283, thus removing it from the vicinity of Tyr106, resulted in a stable H-bond to pyroGlu. In all simulations His showed a van der Waals attraction to Tyr282 and a weak electrostatic repulsion from Arg 306. The ProNH2 had a strong and frequent H-bonding interaction with Arg306. The backbone carbonyls show a frequent H-bonding interaction with the OH group of Tyr282 and strong, often multiple, interactions with Arg306. Three structures, which maintained these interactions simultaneously, were selected as candidates for ligand-receptor complexes. These show persistent interactions of TRH with Ile 109 and Ile 116 in HX 3 and with Tyr310 and Ser313 in HX 7, which will be tested to refine the structure of the ligand-receptor complex. The superposition of the three structures shows the extent of structural flexibility of the receptor and the ligand in the complex. The backbone of TRH inside the receptor is in an alpha-helical conformation, suggesting that the receptor, through its interaction with the ligand, provides the energy required for the conformational change in the ligand from an extended to the folded form.