Structural analysis of a HAMP domain: the linker region of the phototransducer in complex with sensory rhodopsin II

Sensory rhodopsin II, the photophobic receptor from Natronomonas pharaonis (NpSRII)5, forms a 2:2 complex with its cognate transducer (N. pharaonis halobacterial transducer of rhodopsins II (NpHtrII)) in lipid membranes. Light activation of NpSRII leads to a displacement of helix F, which in turn triggers a rotation/screw-like motion of TM2 in NpHtrII. This conformational change is thought to be transmitted through the membrane adjacent conserved signal transduction domain in histidine kinases, adenylyl cyclases, methyl-accepting chemotaxis proteins, and phosphatases (HAMP domain) to the cytoplasmic signaling domain of the transducer. The architecture and function of the HAMP domain are still unknown. In order to obtain information on the structure and dynamics of this region, EPR experiments on a truncated transducer (NpHtrII(157)) and NpSRII, site-directed spin-labeled and reconstituted into purple membrane lipids, have been carried out. A nitroxide scanning involving residues in the transducer helix TM2, in the predicted AS-1 region, and at selected positions in the following connector and AS-2 regions of the HAMP domain has been performed. Accessibility and dynamics data allowed us to identify a helical region up to residue Ala(94) in the AS-1 amphipathic sequence, followed by a highly dynamic domain protruding into the water phase. Additionally, transducer-transducer and transducer-receptor proximity relations revealed the overall architecture of the AS-1 sequences in the 2:2 complex, which are suggested to form a molten globular type of a coiled-coil bundle.     

Signaling and Adaptation Modulate the Dynamics of the Photosensoric Complex of Natronomonas pharaonis

Motile bacteria and archaea respond to chemical and physical stimuli seeking optimal conditions for survival. To this end transmembrane chemo- and photoreceptors organized in large arrays initiate signaling cascades and ultimately regulate the rotation of flagellar motors. To unravel the molecular mechanism of signaling in an archaeal phototaxis complex we performed coarse-grained molecular dynamics simulations of a trimer of receptor/transducer dimers, namely NpSRII/NpHtrII from Natronomonas pharaonis. Signaling is regulated by a reversible methylation mechanism called adaptation, which also influences the level of basal receptor activation. Mimicking two extreme methylation states in our simulations we found conformational changes for the transmembrane region of NpSRII/NpHtrII which resemble experimentally observed light-induced changes. Further downstream in the cytoplasmic domain of the transducer the signal propagates via distinct changes in the dynamics of HAMP1, HAMP2, the adaptation domain and the binding region for the kinase CheA, where conformational rearrangements were found to be subtle. Overall these observations suggest a signaling mechanism based on dynamic allostery resembling models previously proposed for E. coli chemoreceptors, indicating similar properties of signal transduction for archaeal photoreceptors and bacterial chemoreceptors.     

Transmembrane signaling of chemotaxis receptor tar: insights from molecular dynamics simulation studies

Transmembrane signaling of chemotaxis receptors has long been studied, but how the conformational change induced by ligand binding is transmitted across the bilayer membrane is still elusive at the molecular level. To tackle this problem, we carried out a total of 600-ns comparative molecular dynamics simulations (including model-building simulations) of the chemotaxis aspartate receptor Tar (a part of the periplasmic domain/transmembrane domain/HAMP domain) in explicit lipid bilayers. These simulations reveal valuable insights into the mechanistic picture of Tar transmembrane signaling. The piston-like movement of a transmembrane helix induced by ligand binding on the periplasmic side is transformed into a combination of both longitudinal and transversal movements of the helix on the cytoplasmic side as a result of different protein-lipid interactions in the ligand-off and ligand-on states of the receptor. This conformational change alters the dynamics and conformation of the HAMP domain, which is presumably a mechanism to deliver the signal from the transmembrane domain to the cytoplasmic domain. The current results are consistent with the previously suggested dynamic bundle model in which the HAMP dynamics change is a key to the signaling. The simulations provide further insights into the conformational changes relevant to the HAMP dynamics changes in atomic detail.     

Salt-driven equilibrium between two conformations in the HAMP domain from Natronomonas pharaonis: the language of signal transfer?

HAMP domains (conserved in histidine kinases, adenylyl cyclases, methyl-accepting chemotaxis proteins, and phosphatases) perform their putative function as signal transducing units in diversified environments in a variety of protein families. Here the conformational changes induced by environmental agents, namely salt and temperature, on the structure and function of a HAMP domain of the phototransducer from Natronomonas pharaonis (NpHtrII) in complex with sensory rhodopsin II (NpSRII) were investigated by site-directed spin labeling electron paramagnetic resonance. A series of spin labeled mutants were engineered in NpHtrII157, a truncated analog containing only the first HAMP domain following the transmembrane helix 2. This truncated transducer is shown to be a valid model system for a signal transduction domain anchored to the transmembrane light sensor NpSRII. The HAMP domain is found to be engaged in a "two-state" equilibrium between a highly dynamic (dHAMP) and a more compact (cHAMP) conformation. The structural properties of the cHAMP as proven by mobility, accessibility, and intra-transducer-dimer distance data are in agreement with the four helical bundle NMR model of the HAMP domain from Archaeoglobus fulgidus.     

HAMP domain signal relay mechanism in a sensory rhodopsin-transducer complex

The phototaxis receptor complex composed of sensory rhodopsin II (SRII) and the transducer subunit HtrII mediates photorepellent responses in haloarchaea. Light-activated SRII transmits a signal through two HAMP switch domains (HAMP1 and HAMP2) in HtrII that bridge the photoreceptive membrane domain of the complex and the cytoplasmic output kinase-modulating domain. HAMP domains, widespread signal relay modules in prokaryotic sensors, consist of four-helix bundles composed of two helices, AS1 and AS2, from each of two dimerized transducer subunits. To examine their molecular motion during signal transmission, we incorporated SRII-HtrII dimeric complexes in nanodiscs to allow unrestricted probe access to the cytoplasmic side HAMP domains. Spin-spin dipolar coupling measurements confirmed that in the nanodiscs, SRII photoactivation induces helix movement in the HtrII membrane domain diagnostic of transducer activation. Labeling kinetics of a fluorescein probe in monocysteine-substituted HAMP1 mutants revealed a light-induced shift of AS2 against AS1 by one-half α-helix turn with minimal other changes. An opposite shift of AS2 against AS1 in HAMP2 at the corresponding positions supports the proposal from x-ray crystal structures by Airola et al. (Airola, M. V., Watts, K. J., Bilwes, A. M., and Crane, B. R. (2010) Structure 18, 436-448) that poly-HAMP chains undergo alternating opposite interconversions to relay the signal. Moreover, we found that haloarchaeal cells expressing a HAMP2-deleted SRII-HtrII exhibit attractant phototaxis, opposite from the repellent phototaxis mediated by the wild-type di-HAMP SRII-HtrII complex. The opposite conformational changes and corresponding opposite output signals of HAMP1 and HAMP2 imply a signal transmission mechanism entailing small shifts in helical register between AS1 and AS2 alternately in opposite directions in adjacent HAMPs.     

Mechanism of regulation of receptor histidine kinases

Bacterial transmembrane receptors regulate an intracellular catalytic output in response to extracellular sensory input. To investigate the conformational changes that relay the regulatory signal, we have studied the HAMP domain, a ubiquitous intracellular module connecting input to output domains. HAMP forms a parallel, dimeric, four-helical coiled coil, and rational substitutions in our model domain (Af1503 HAMP) induce a transition in its interhelical packing, characterized by axial rotation of all four helices (the gearbox signaling model). We now illustrate how these conformational changes are propagated to a downstream domain by fusing Af1503 HAMP variants to the DHp domain of EnvZ, a bacterial histidine kinase. Structures of wild-type and mutant constructs are correlated with ligand response in vivo, clearly associating them with distinct signaling states. We propose that altered recognition of the catalytic domain by DHp, rather than a shift in position of the phospho-accepting histidine, forms the basis for regulation of kinase activity.     

Light-induced switching of HAMP domain conformation and dynamics revealed by time-resolved EPR spectroscopy

HAMP domains are widely abundant signaling modules. The putative mechanism of their function comprises switching between two distinct states. To unravel these conformational transitions, we apply site-directed spin labeling and time-resolved EPR spectroscopy to the phototactic receptor/transducer complex NpSRII/NpHtrII. We characterize the kinetic coupling of NpHtrII to NpSRII along with the activation period of the transducer and follow the transient conformational signal. The observed transient shift towards a more compact state of the HAMP domain upon light-activation agrees with structure-based calculations. It thereby validates the two modeled signaling states and integrates the domain’s dynamics into the current model.     

Mutational analysis of the control cable that mediates transmembrane signaling in the Escherichia coli serine chemoreceptor

During transmembrane signaling by Escherichia coli Tsr, changes in ligand occupancy in the periplasmic serine-binding domain promote asymmetric motions in a four-helix transmembrane bundle. Piston displacements of the signaling TM2 helix in turn modulate the HAMP bundle on the cytoplasmic side of the membrane to control receptor output signals to the flagellar motors. A five-residue control cable joins TM2 to the HAMP AS1 helix and mediates conformational interactions between them. To explore control cable structural features important for signal transmission, we constructed and characterized all possible single amino acid replacements at the Tsr control cable residues. Only a few lesions abolished Tsr function, indicating that the chemical nature and size of the control cable side chains are not individually critical for signal control. Charged replacements at I214 mimicked the signaling consequences of attractant or repellent stimuli, most likely through aberrant structural interactions of the mutant side chains with the membrane interfacial environment. Prolines at residues 214 to 217 also caused signaling defects, suggesting that the control cable has helical character. However, proline did not disrupt function at G213, the first control cable residue, which might serve as a structural transition between the TM2 and AS1 helix registers. Hydrophobic amino acids at S217, the last control cable residue, produced attractant-mimic effects, most likely by contributing to packing interactions within the HAMP bundle. These results suggest a helix extension mechanism of Tsr transmembrane signaling in which TM2 piston motions influence HAMP stability by modulating the helicity of the control cable segment.     

The mechanisms of HAMP-mediated signaling in transmembrane receptors

HAMP domains mediate signal transduction in over 7500 enzyme-coupled receptors represented in all kingdoms of life. The HAMP domain of the putative archaeal receptor Af1503 has a parallel, dimeric, four-helical coiled coil structure, but with unusual core packing, related to canonical packing by concerted axial rotation of the helices. This has led to the gearbox model for signal transduction, whereby the alternate packing modes correspond to signaling states. Here we present structures of a series of Af1503 HAMP variants. We show that substitution of a conserved small side chain within the domain core (A291) for larger residues induces a gradual transition in packing mode, involving both changes in helix rotation and bundle shape, which are most prominent at the C-terminal, output end of the domain. These are correlated with activity and ligand response in vitro and in vivo by incorporating Af1503 HAMP into mycobacterial adenylyl cyclase assay systems.     

Not too loose, not too tight--just right. Biphasic control of the Tsr HAMP domain

HAMP domains communicate between input and output signalling modules in a wide variety of bacterial sensor proteins. In the Tsr chemoreceptor, they convert a signal initiated by binding of serine to the periplasmic domain of the protein into regulation of receptor control of the CheA kinase, and ultimately of the direction of flagellar rotation. In this issue, Zhou et al. report an extensive mutational analysis of the Tsr HAMP domain that shows that it can assume a number of different signalling states, which presumably correspond to a variety of different conformations. The two conformational extremes of a tightly packed and a loosely packed HAMP four‐helix bundle support only low levels of CheA activity. Thus, Tsr HAMP does not function as a simple on‐off, two‐state device but rather as a dynamic structure with biphasic control. The normal physiological operating range of Tsr is proposed to be at intermediate degrees of packing of the HAMP four‐helix bundle, but HAMP domains in other proteins could occupy different portions of the conformational spectrum.     

HAMP domain structural determinants for signalling and sensory adaptation in Tsr,the Escherichia coli serine chemoreceptor

HAMP domains mediate input–output transactions in many bacterial signalling proteins. To clarify the mechanistic logic of HAMP signalling, we constructed Tsr‐HAMP deletion derivatives and characterized their steady‐state signal outputs and sensory adaptation properties with flagellar rotation and receptor methylation assays. Tsr molecules lacking the entire HAMP domain or just the HAMP‐AS2 helix generated clockwise output signals, confirming that kinase activation is the default output state of the chemoreceptor signalling domain and that attractant stimuli shift HAMP to an overriding kinase‐off signalling state to elicit counter‐clockwise flagellar responses. Receptors with deletions of the AS1 helices, which free the AS2 helices from bundle‐packing constraints, exhibited kinase‐off signalling behaviour that depended on three C‐terminal hydrophobic residues of AS2. We conclude that AS2/AS2′ packing interactions alone can play an important role in controlling output kinase activity. Neither kinase‐on nor kinase‐off HAMP deletion outputs responded to sensory adaptation control, implying that out‐of‐range conformations or bundle‐packing stabilities of their methylation helices prevent substrate recognition by the adaptation enzymes. These observations support the previously proposed biphasic, dynamic‐bundle mechanism of HAMP signalling and additionally show that the structural interplay of helix‐packing interactions between HAMP and the adjoining methylation helices is critical for sensory adaptation control of receptor output.     

The HAMP domain structure implies helix rotation in transmembrane signaling

HAMP domains connect extracellular sensory with intracellular signaling domains in over 7500 proteins, including histidine kinases, adenylyl cyclases, chemotaxis receptors, and phosphatases. The solution structure of an archaeal HAMP domain shows a homodimeric, four-helical, parallel coiled coil with unusual interhelical packing, related to the canonical packing by rotation of the helices. This suggests a model for the mechanism of signal transduction, in which HAMP alternates between the observed conformation and a canonical coiled coil. We explored this mechanism in vitro and in vivo using HAMP domain fusions with a mycobacterial adenylyl cyclase and an E. coli chemotaxis receptor. Structural and functional studies show that the equilibrium between the two forms is dependent on the side-chain size of residue 291, which is alanine in the wild-type protein.     

Delineating PAS‐HAMP interaction surfaces and signalling‐associated changes in the aerotaxis receptor Aer

The Escherichia coli aerotaxis receptor, Aer, monitors cellular oxygen and redox potential via FAD bound to a cytosolic PAS domain. Here, we show that Aer‐PAS controls aerotaxis through direct, lateral interactions with a HAMP domain. This contrasts with most chemoreceptors where signals propagate along the protein backbone from an N‐terminal sensor to HAMP. We mapped the interaction surfaces of the Aer PAS, HAMP and proximal signalling domains in the kinase‐off state by probing the solvent accessibility of 129 cysteine substitutions. Inaccessible PAS‐HAMP surfaces overlapped with a cluster of PAS kinase‐on lesions and with cysteine substitutions that crosslinked the PAS β ‐scaffold to the HAMP AS‐2 helix. A refined Aer PAS‐HAMP interaction model is presented. Compared to the kinase‐off state, the kinase‐on state increased the accessibility of HAMP residues (apparently relaxing PAS‐HAMP interactions), but decreased the accessibility of proximal signalling domain residues. These data are consistent with an alternating static‐dynamic model in which oxidized Aer‐PAS interacts directly with HAMP AS‐2, enforcing a static HAMP domain that in turn promotes a dynamic proximal signalling domain, resulting in a kinase‐off output. When PAS‐FAD is reduced, PAS interaction with HAMP is relaxed and a dynamic HAMP and static proximal signalling domain convey a kinase‐on output.     

Mutational analysis of the connector segment in the HAMP domain of Tsr, the Escherichia coli serine chemoreceptor

HAMP domains are approximately 50-residue motifs, found in many bacterial signaling proteins, that consist of two amphiphilic helices joined by a nonhelical connector segment. The HAMP domain of Tsr, the serine chemoreceptor of Escherichia coli, receives transmembrane input signals from the periplasmic serine binding domain and in turn modulates output signals from the Tsr kinase control domain to elicit chemotactic responses. We created random amino acid replacements at each of the 14 connector residues of Tsr-HAMP to identify those that are critical for Tsr function. In all, we surveyed 179 connector missense mutants and identified three critical residues (G235, L237, and I241) at which most replacements destroyed Tsr function and another important residue (G245) at which most replacements impaired Tsr function. The region surrounding G245 tolerated 1-residue deletions and insertions of up to 10 glycines, suggesting a role as a relatively nonspecific, flexible linker. The critical connector residues are consistent with a structural model of the Tsr-HAMP domain based on the solution structure of an isolated thermophile HAMP domain (M. Hulko, F. Berndt, M. Gruber, J. U. Linder, V. Truffault, A. Schultz, J. Martin, J. E. Schultz, A. N. Lupas, and M. Coles, Cell 126:929-940, 2006) in which G235 defines a critical turn at the C terminus of the first helix and L237 and I241 pack against the helices, perhaps to stabilize alternative HAMP signaling conformations. Most I241 lesions locked Tsr signal output in the kinase-on mode, implying that this residue is responsible mainly for stabilizing the kinase-off signaling state. In contrast, lesions at L237 resulted in a variety of aberrant output patterns, suggesting a role in toggling output between signaling states.     

Different conformations of the kinase-on and kinase-off signaling states in the Aer HAMP domain

HAMP domains are sensory transduction modules that connect input and output domains in diverse signaling proteins from archaea, bacteria, and lower eukaryotes. Here, we employed in vivo disulfide cross-linking to explore the structure of the HAMP domain in the Escherichia coli aerotaxis receptor Aer. Using an Aer HAMP model based on the structure of Archaeoglobus fulgidus Af1503-HAMP, the closest residue pairs at the interface of the HAMP AS-1 and AS-2' helices were determined and then replaced with cysteines and cross-linked in vivo. Except for a unique discontinuity in AS-2, the data suggest that the Aer HAMP domain forms a parallel four-helix bundle that is similar to the structure of Af1503. The HAMP discontinuity was associated with a segment of AS-2 that was recently shown to interact with the Aer-PAS sensing domain. The four-helix HAMP bundle and its discontinuity were maintained in both the kinase-on and kinase-off states of Aer, although differences in the rates of disulfide formation also indicated the existence of different HAMP conformations in the kinase-on and kinase-off states. In particular, the kinase-on state was accompanied by significantly increased disulfide formation rates at the distal end of the HAMP four-helix bundle. This indicates that HAMP signaling may be associated with a tilting of the AS-1 and AS-2' helices, which may be the signal that is transmitted to the kinase control region of Aer.     

Segmental Helical Motions and Dynamical Asymmetry Modulate Histidine Kinase Autophosphorylation

Histidine kinases (HKs) are dimeric receptors that participate in most adaptive responses to environmental changes in prokaryotes. Although it is well established that stimulus perception triggers autophosphorylation in many HKs, little is known on how the input signal propagates through the HAMP domain to control the transient interaction between the histidine-containing and ATP-binding domains during the catalytic reaction. Here we report crystal structures of the full cytoplasmic region of CpxA, a prototypical HK involved in Escherichia coli response to envelope stress. The structural ensemble, which includes the Michaelis complex, unveils HK activation as a highly dynamic process, in which HAMP modulates the segmental mobility of the central HK α-helices to promote a strong conformational and dynamical asymmetry that characterizes the kinase-active state. A mechanical model based on our structural and biochemical data provides insights into HAMP-mediated signal transduction, the autophosphorylation reaction mechanism, and the symmetry-dependent control of HK kinase/phosphatase functional states.     

Role of the HAMP domain region of sensory rhodopsin transducers in signal transduction

Archaea are able to sense light via the complexes of sensory rhodopsins I and II and their corresponding chemoreceptor-like transducers HtrI and HtrII. Though generation of the signal has been studied in detail, the mechanism of its propagation to the cytoplasm remains obscured. The cytoplasmic part of the transducer consists of adaptation and kinase activity modulating regions, connected to transmembrane helices via two HAMP (histidine kinases, adenylyl cyclases, methyl-accepting chemotaxis proteins, phosphatases) domains. The inter-HAMP region of Natronomonas pharaonis HtrII (NpHtrII) was found to be α-helical [Hayashi, K., et al. (2007) Biochemistry 46, 14380-14390]. We studied the inter-HAMP regions of NpHtrII and other phototactic signal transducers by means of molecular dynamics. Their structure is found to be a bistable asymmetric coiled coil, in which the protomers are longitudinally shifted by ~1.3 ?. The free energy penalty for the symmetric structure is estimated to be 1.2-1.5 kcal/mol depending on the molarity of the solvent. Both flanking HAMP domains are mechanistically coupled to the inter-HAMP region and are asymmetric. The longitudinal shift in the inter-HAMP region is coupled with the in-plane displacement of the cytoplasmic part by 8.6 ? relative to the transmembrane part. The established properties suggest that (1) the signal may be transduced through the inter-HAMP domain switching and (2) the inter-HAMP region may allow cytoplasmic parts of the transducers to come sufficiently close to each other to form oligomers.     

Gain‐of‐function mutations cluster in distinct regions associated with the signalling pathway in the PAS domain of the aerotaxis receptor,Aer

The Aer receptor monitors internal energy (redox) levels in Escherichia coli with an FAD‐containing PAS domain. Here, we randomly mutagenized the region encoding residues 14–119 of the PAS domain and found 72 aerotaxis‐defective mutants, 24 of which were gain‐of‐function, signal‐on mutants. The mutations were mapped onto an Aer homology model based on the structure of the PAS–FAD domain in NifL from Azotobacter vinlandii. Signal‐on lesions clustered in the FAD binding pocket, the β‐scaffolding and in the N‐cap loop. We suggest that the signal‐on lesions mimic the ‘signal‐on’ state of the PAS domain, and therefore may be markers for the signal‐in and signal‐out regions of this domain. We propose that the reduction of FAD rearranges the FAD binding pocket in a way that repositions the β‐scaffolding and the N‐cap loop. The resulting conformational changes are likely to be conveyed directly to the HAMP domain, and on to the kinase control module. In support of this hypothesis, we demonstrated disulphide band formation between cysteines substituted at residues N98C or I114C in the PAS β‐scaffold and residue Q248C in the HAMP AS‐2 helix.     

Genetic analysis of the HAMP domain of the Aer aerotaxis sensor localizes flavin adenine dinucleotide-binding determinants to the AS-2 helix

HAMP domains are signal transduction domains typically located between the membrane anchor and cytoplasmic signaling domain of the proteins in which they occur. The prototypical structure consists of two helical amphipathic sequences (AS-1 and AS-2) connected by a region of undetermined structure. The Escherichia coli aerotaxis receptor, Aer, has a HAMP domain and a PAS domain with a flavin adenine dinucleotide (FAD) cofactor that senses the intracellular energy level. Previous studies reported mutations in the HAMP domain that abolished FAD binding to the PAS domain. In this study, using random and site-directed mutagenesis, we identified the distal helix, AS-2, as the component of the HAMP domain that stabilizes FAD binding. AS-2 in Aer is not amphipathic and is predicted to be buried. Mutations in the sequence coding for the contiguous proximal signaling domain altered signaling by Aer but did not affect FAD binding. The V264M residue replacement in this region resulted in an inverted response in which E. coli cells expressing the mutant Aer protein were repelled by oxygen. Bioinformatics analysis of aligned HAMP domains indicated that the proximal signaling domain is conserved in other HAMP domains that are not involved in chemotaxis or aerotaxis. Only one null mutation was found in the coding sequence for the HAMP AS-1 and connector regions, suggesting that these are not active signal transduction sites. We consider a model in which the signal from FAD is transmitted across a PAS-HAMP interface to AS-2 or the proximal signaling domain.     

Interactions between the PAS and HAMP domains of the Escherichia coli aerotaxis receptor Aer

The Escherichia coli energy-sensing Aer protein initiates aerotaxis towards environments supporting optimal cellular energy. The Aer sensor is an N-terminal, FAD-binding, PAS domain. The PAS domain is linked by an F1 region to a membrane anchor, and in the C-terminal half of Aer, a HAMP domain links the membrane anchor to the signaling domain. The F1 region, membrane anchor, and HAMP domain are required for FAD binding. Presumably, alterations in the redox potential of FAD induce conformational changes in the PAS domain that are transmitted to the HAMP and C-terminal signaling domains. In this study we used random mutagenesis and intragenic pseudoreversion analysis to examine functional interactions between the HAMP domain and the N-terminal half of Aer. Missense mutations in the HAMP domain clustered in the AS-2 alpha-helix and abolished FAD binding to Aer, as previously reported. Three amino acid replacements in the Aer-PAS domain, S28G, A65V, and A99V, restored FAD binding and aerotaxis to the HAMP mutants. These suppressors are predicted to surround a cleft in the PAS domain that may bind FAD. On the other hand, suppression of an Aer-C253R HAMP mutant was specific to an N34D substitution with a predicted location on the PAS surface, suggesting that residues C253 and N34 interact or are in close proximity. No suppressor mutations were identified in the F1 region or membrane anchor. We propose that functional interactions between the PAS domain and the HAMP AS-2 helix are required for FAD binding and aerotactic signaling by Aer.