Enzymatic Activity versus Structural Dynamics: The Case of Acetylcholinesterase Tetramer

The function of many proteins, such as enzymes, is modulated by structural fluctuations. This is especially the case in gated diffusion-controlled reactions (where the rates of the initial diffusional encounter and of structural fluctuations determine the overall rate of the reaction) and in oligomeric proteins (where function often requires a coordinated movement of individual subunits). A classic example of a diffusion-controlled biological reaction catalyzed by an oligomeric enzyme is the hydrolysis of synaptic acetylcholine (ACh) by tetrameric acetylcholinesterase (AChEt). Despite decades of efforts, the extent to which enzymatic efficiency of AChEt (or any other enzyme) is modulated by flexibility is not fully determined. This article attempts to determine the correlation between the dynamics of AChEt and the rate of reaction between AChEt and ACh. We employed equilibrium and nonequilibrium electro-diffusion models to compute rate coefficients for an ensemble of structures generated by molecular dynamics simulation. We found that, for the static initial model, the average reaction rate per active site is ∼22-30% slower in the tetramer than in the monomer. However, this effect of tetramerization is modulated by the intersubunit motions in the tetramer such that a complex interplay of steric and electrostatic effects either guides or blocks the substrate into or from each of the four active sites. As a result, the rate per active site calculated for some of the tetramer structures is only ∼15% smaller than the rate in the monomer. We conclude that structural dynamics minimizes the adverse effect of tetramerization, allowing the enzyme to maintain similar enzymatic efficiency in different oligomerization states.     

The association of tetrameric acetylcholinesterase with ColQ tail: a block normal mode analysis

Acetylcholinesterase (AChE) rapidly hydrolyzes acetylcholine in the neuromuscular junctions and other cholinergic synapses to terminate the neuronal signal. In physiological conditions, AChE exists as tetramers associated with the proline-rich attachment domain (PRAD) of either collagen-like Q subunit (ColQ) or proline-rich membrane-anchoring protein. Crystallographic studies have revealed that different tetramer forms may be present, and it is not clear whether one or both are relevant under physiological conditions. Recently, the crystal structure of the tryptophan amphiphilic tetramerization (WAT) domain of AChE associated with PRAD ([WAT]₄PRAD), which mimics the interface between ColQ and AChE tetramer, became available. In this study we built a complete tetrameric mouse [AChE_T]₄–ColQ atomic structure model, based on the crystal structure of the [WAT]₄PRAD complex. The structure was optimized using energy minimization. Block normal mode analysis was done to investigate the low-frequency motions of the complex and to correlate the structure model with the two known crystal structures of AChE tetramer. Significant low-frequency motions among the catalytic domains of the four AChE subunits were observed, while the [WAT]₄PRAD part held the complex together. Normal mode involvement analysis revealed that the two lowest frequency modes were primarily involved in the conformational changes leading to the two crystal structures. The first 30 normal modes can account for more than 75% of the conformational changes in both cases. The evidence further supports the idea of a flexible tetramer model for AChE. This model can be used to study the implications of the association of AChE with ColQ.     

A four-to-one association between peptide motifs: four C-terminal domains from cholinesterase assemble with one proline-rich attachment domain (PRAD) in the secretory pathway.

The major type of acetylcholinesterase in vertebrates (AChET) is characterized by the presence of a short C-terminal domain of 40 residues, the 'tryptophan amphiphilic tetramerization' (WAT) domain. The presence of this domain is not necessary for catalytic activity but is responsible for hydrophobic interactions and for the capacity of AChET subunits to form quaternary associations with anchoring proteins, thereby conditioning their functional localization. In the collagen tail of asymmetric forms, we characterized a small conserved region that is sufficient for binding an AChET tetramer, the proline-rich attachment domain (PRAD). We show that the WAT domain alone is sufficient for association with the PRAD, and that it can attach foreign proteins (alkaline phosphatase, GFP) to a PRAD-containing construct with a glycophosphatidylinositol anchor (GPI), and thus anchor them to the cell surface. Furthermore, we show that isolated WAT domains, or proteins containing a WAT domain, can replace individual AChET subunits in PRAD-linked tetramers. This suggests that the four WAT domains interact with the PRAD in a similar manner. These quaternary interactions can form without intercatenary disulfide bonds. The common catalytic domains of AChE are not necessary for tetrameric assembly, although they may contribute to the stability of the tetramer.     

The synaptic acetylcholinesterase tetramer assembles around a polyproline II helix

Functional localization of acetylcholinesterase (AChE) in vertebrate muscle and brain depends on interaction of the tryptophan amphiphilic tetramerization (WAT) sequence, at the C-terminus of its major splice variant (T), with a proline-rich attachment domain (PRAD), of the anchoring proteins, collagenous (ColQ) and proline-rich membrane anchor. The crystal structure of the WAT/PRAD complex reveals a novel supercoil structure in which four parallel WAT chains form a left-handed superhelix around an antiparallel left-handed PRAD helix resembling polyproline II. The WAT coiled coils possess a WWW motif making repetitive hydrophobic stacking and hydrogen-bond interactions with the PRAD. The WAT chains are related by an approximately 4-fold screw axis around the PRAD. Each WAT makes similar but unique interactions, consistent with an asymmetric pattern of disulfide linkages between the AChE tetramer subunits and ColQ. The P59Q mutation in ColQ, which causes congenital endplate AChE deficiency, and is located within the PRAD, disrupts crucial WAT-WAT and WAT-PRAD interactions. A model is proposed for the synaptic AChE(T) tetramer.     

Isoflurane alters the structure and dynamics of GLIC

Pentameric ligand-gated ion channels are targets of general anesthetics. Although the search for discrete anesthetic binding sites has achieved some degree of success, little is known regarding how anesthetics work after the events of binding. Using the crystal structures of the bacterial Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC), which is sensitive to a variety of general anesthetics, we performed multiple molecular dynamics simulations in the presence and absence of the general anesthetic isoflurane. Isoflurane bound to several locations within GLIC, including the transmembrane pocket identified crystallographically, the extracellular (EC) domain, and the interface of the EC and transmembrane domains. Isoflurane also entered the channel after the pore was dehydrated in one of the simulations. Isoflurane disrupted the quaternary structure of GLIC, as evidenced in a striking association between the binding and breakage of intersubunit salt bridges in the EC domain. The pore-lining helix experienced lateral and inward radial tilting motion that contributed to the channel closure. Isoflurane binding introduced strong anticorrelated motions between different subunits of GLIC. The demonstrated structural and dynamical modulations by isoflurane aid in the understanding of the underlying mechanism of anesthetic inhibition of GLIC and possibly other homologous pentameric ligand-gated ion channels.     

Allosteric regulation of the biotin-dependent enzyme pyruvate carboxylase by acetyl-CoA

The activity of the biotin-dependent enzyme pyruvate carboxylase from many organisms is highly regulated by the allosteric activator acetyl-CoA. A number of X-ray crystallographic structures of the native pyruvate carboxylase tetramer are now available for the enzyme from Rhizobium etli and Staphylococcus aureus. Although all of these structures show that intersubunit catalysis occurs, in the case of the R. etli enzyme, only two of the four subunits have the allosteric activator bound to them and are optimally configured for catalysis of the overall reaction. However, it is apparent that acetyl-CoA binding does not induce the observed asymmetrical tetramer conformation and it is likely that, under normal reaction conditions, all of the subunits have acetyl-CoA bound to them. Thus the activation of the enzyme by acetyl-CoA involves more subtle structural effects, one of which may be to facilitate the correct positioning of Arg353 and biotin in the biotin carboxylase domain active site, thereby promoting biotin carboxylation and, at the same time, preventing abortive decarboxylation of carboxybiotin. It is also apparent from the crystal structures that there are allosteric interactions induced by acetyl-CoA binding in the pair of subunits not optimally configured for catalysis of the overall reaction.     

Active and inactive state structures of unliganded Lactobacillus casei allosteric L‐lactate dehydrogenase

Lactobacillus casei L ‐lactate dehydrogenase (LCLDH) is activated through the homotropic and heterotropic activation effects of pyruvate and fructose 1,6‐bisphosphate (FBP), respectively, and exhibits unusually high pH‐dependence in the allosteric effects of these ligands. The active (R) and inactive (T) state structures of unliganded LCLDH were determined at 2.5 and 2.6 Å resolution, respectively. In the catalytic site, the structural rearrangements are concerned mostly in switching of the orientation of Arg171 through the flexible intersubunit contact at the Q‐axis subunit interface. The distorted orientation of Arg171 in the T state is stabilized by a unique intra‐helix salt bridge between Arg171 and Glu178, which is in striking contrast to the multiple intersubunit salt bridges in Lactobacillus pentosus nonallosteric L ‐lactate dehydrogenase. In the backbone structure, major structural rearrangements of LCLDH are focused in two mobile regions of the catalytic domain. The two regions form an intersubunit linkage through contact at the P‐axis subunit interface involving Arg185, replacement of which with Gln severely decreases the homotropic and hetertropic activation effects on the enzyme. These two regions form another intersubunit linkage in the Q‐axis related dimer through the rigid NAD‐binding domain, and thus constitute a pivotal frame of the intersubunit linkage for the allosteric motion, which is coupled with the concerted structural change of the four subunits in a tetramer, and of the binding sites for pyruvate and FBP. The unique intersubunit salt bridges, which are observed only in the R state structure, are likely involved in the pH‐dependent allosteric equilibrium. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Normal‐mode‐based modeling of allosteric couplings that underlie cyclic conformational transition in F1 ATPase

F₁ ATPase, a rotary motor comprised of a central stalk ( γ subunit) enclosed by three α and β subunits alternately arranged in a hexamer, features highly cooperative binding and hydrolysis of ATP. Despite steady progress in biophysical, biochemical, and computational studies of this fascinating motor, the structural basis for cooperative ATPase involving its three catalytic sites remains not fully understood. To illuminate this key mechanistic puzzle, we have employed a coarse‐grained elastic network model to probe the allosteric couplings underlying the cyclic conformational transition in F₁ ATPase at a residue level of detail. We will elucidate how ATP binding and product (ADP and phosphate) release at two catalytic sites are coupled with the rotation of γ subunit via various domain motions in α ₃ β ₃ hexamer (including intrasubunit hinge‐bending motions in β subunits and intersubunit rigid‐body rotations between adjacent α and β subunits). To this end, we have used a normal‐mode‐based correlation analysis to quantify the allosteric couplings of these domain motions to local motions at catalytic sites and the rotation of γ subunit. We have then identified key amino acid residues involved in the above couplings, some of which have been validated against past studies of mutated and γ ‐truncated F₁ ATPase. Our finding strongly supports a binding change mechanism where ATP binding to the empty catalytic site triggers a series of intra‐ and intersubunit domain motions leading to ATP hydrolysis and product release at the other two closed catalytic sites. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Modeling Neuronal Nicotinic and GABA Receptors: Important Interface Salt-Links and Protein Dynamics

Protein motions in the Cys-loop ligand-gated ion receptors that govern the gating mechanism are still not well understood. The details as to how motions in the ligand-binding domain are translated to the transmembrane domain and how subunit rotations are linked to bring about the cooperative movements involved in gating are under investigation. Homology models of the α4β2 nicotinic acetylcholine (nACh) and β2α1γ2 GABA receptors were constructed based on the torpedo neuromuscular-like nicotinic receptor structure. The template constructed for the full electron microscopy structure must be considered more reliable for structure-function studies due to the preservation of the E45-R209 salt-link. Many other salt-links are seen to transiently form, including switching off of the E45-R209 link, within a network of potential salt-links at the binding domain to the transmembrane domain interface region. Several potentially important intersubunit salt-links form in both the nAChR and GABAR structures during the simulation and appear conserved across many subunit combinations, such as the salt-link between α4.E262 and β2.K255 in nAChR (β2.E262 and α1.K263 in GABAR), at the top of the pore-lining M2 helices, and the intersubunit link of R210 on the M1-linker to E168 on the β8-sheet of the adjacent subunit in the GABA receptor (E175-K46 being the structurally equivalent link in the nAChR, with reversed polarity). A network of other salt-links may be vital for transmitting the cooperative gating motions between subunits that become biased upon ligand binding. The changes seen in the simulations suggest that this network of salt-links helps to set limits and specific states for the conformational changes involved in gating of the receptor. We hope that these hypotheses will be tested experimentally in the near future.     

The C-terminal t peptide of acetylcholinesterase forms an alpha helix that supports homomeric and heteromeric interactions.

Acetylcholinesterase subunits of type T (AChET) possess an alternatively spliced C-terminal peptide (t peptide) which endows them with amphiphilic properties, the capacity to form various homo-oligomers and to associate, as a tetramer, with anchoring proteins containing a proline rich attachment domain (PRAD). The t peptide contains seven conserved aromatic residues. By spectroscopic analyses of the synthetic peptides covering part or all of the t peptide of Torpedo AChET, we show that the region containing the aromatic residues adopts an alpha helical structure, which is favored in the presence of lipids and detergent micelles: these residues therefore form a hydrophobic cluster in a sector of the helix. We also analyzed the formation of disulfide bonds between two different AChET subunits, and between AChET subunits and a PRAD-containing protein [the N-terminal fragment of the ColQ protein (QN)] possessing two cysteines upstream or downstream of the PRAD. This shows that, in the complex formed by four T subunits with QN (T4-QN), the t peptides are not folded on themselves as hairpins but instead are all oriented in the same direction, antiparallel to that of the PRAD. The formation of disulfide bonds between various pairs of cysteines, introduced by mutagenesis at various positions in the t peptides, indicates that this complex possesses a surprising flexibility.     

Open interface and large quaternary structure movements in 3D domain swapped proteins: insights from molecular dynamics simulations of the C-terminal swapped dimer of ribonuclease A

Bovine pancreatic ribonuclease (RNase A) forms two three-dimensional (3D) domain swapped dimers. Crystallographic investigations have revealed that these dimers display completely different quaternary structures: one dimer (N-dimer), which presents the swapping of the N-terminal helix, is characterized by a compact structure, whereas the other (C-dimer), which is stabilized by the exchange of the C-terminal end, shows a rather loose assembly of the two subunits. The dynamic properties of monomeric RNase A and of the N-dimer have been extensively characterized. Here, we report a molecular dynamics investigation carried out on the C-dimer. This computational experiment indicates that the quaternary structure of the C-dimer undergoes large fluctuations. These motions do not perturb the proper folding of the two subunits, which retain the dynamic properties of RNase A and the N-dimer. Indeed, the individual subunits of the C-dimer display the breathing motion of the beta-sheet structure, which is important for the enzymatic activity of pancreatic-like ribonucleases. In contrast to what has been observed for the N-dimer, the breathing motion of the two subunits of the C-dimer is not coupled. This finding suggests that the intersubunit communications in a 3D domain swapped dimer strongly rely on the extent of the interchain interface. Furthermore, the observation that the C-dimer is endowed with a high intrinsic flexibility holds interesting implications for the specific properties of 3D domain swapped dimers. Indeed, a survey of the quaternary structures of the other 3D domain swapped dimers shows that large variations are often observed when the structural determinations are conducted in different experimental conditions. The 3D domain swapping phenomenon coupled with the high flexibility of the quaternary structure may be relevant for protein-protein recognition, and in particular for the pathological aggregations.     

The essential dynamics of Cu, Zn superoxide dismutase: suggestion of intersubunit communication.

A 300-ps molecular dynamics simulation of the whole Cu, Zn superoxide dismutase dimer has been carried out in water, and the trajectory has been analyzed by the essential dynamics method. The results indicate that the motion is defined by few preferred directions identified by the first four to six eigenvectors and that the motion of the two monomers at each instant is not symmetrical. The vectors symmetrical to the eigenvectors are significantly sampled, suggesting that, on average, the motions of the two subunits will exchange. Large intra- and intersubunit motions involving different subdomains of the protein are observed. A mechanical coupling between the two subunits is also suggested, because displacements of the loops surrounding the active site in one monomer are correlated with the motion of parts of the second toward the intersubunit interface.     

Acetylcholinesterase: C-terminal domains,molecular forms and functional localization

Acetylcholinesterase (AChE) possesses short C-terminal peptides that are not necessary for catalytic activity. These peptides belong to different classes (R, H, T, S) and define the post-translational processing and targeting of the enzyme. In vertebrates, subunits of type H (AChE_H) and of type T (AChE_T) are the most important: AChE_H subunits produced glycolipid (GPI)-anchored dimers and AChE_T subunits produce hetero-oligometric forms such as membrane-bound tetramer in the mammalian brain (containing a 20 kDa hydrophobic protein) and asymmetric collagen-tailed forms in neuromuscular junctions (containing a specific collagen, ColQ). The T peptide allows the formation of tetrameric assemblies with a proline-rich attachment domain (PRAD) of collagen ColQ. These complex molecular structures condition the functional localization of the enzyme in the supramolecular architecture of cholinegic synapses.     

Propofol Binding to the Resting State of the Gloeobacter violaceus Ligand-gated Ion Channel (GLIC) Induces Structural Changes in the Inter- and Intrasubunit Transmembrane Domain (TMD) Cavities

General anesthetics exert many of their CNS actions by binding to and modulating membrane-embedded pentameric ligand-gated ion channels (pLGICs). The structural mechanisms underlying how anesthetics modulate pLGIC function remain largely unknown. GLIC, a prokaryotic pLGIC homologue, is inhibited by general anesthetics, suggesting anesthetics stabilize a closed channel state, but in anesthetic-bound GLIC crystal structures the channel appears open. Here, using functional GLIC channels expressed in oocytes, we examined whether propofol induces structural rearrangements in the GLIC transmembrane domain (TMD). Residues in the GLIC TMD that frame intrasubunit and intersubunit water-accessible cavities were individually mutated to cysteine. We measured and compared the rates of modification of the introduced cysteines by sulfhydryl-reactive reagents in the absence and presence of propofol. Propofol slowed the rate of modification of L240C (intersubunit) and increased the rate of modification of T254C (intrasubunit), indicating that propofol binding induces structural rearrangements in these cavities that alter the local environment near these residues. Propofol acceleration of T254C modification suggests that in the resting state propofol does not bind in the TMD intrasubunit cavity as observed in the crystal structure of GLIC with bound propofol (Nury, H., Van Renterghem, C., Weng, Y., Tran, A., Baaden, M., Dufresne, V., Changeux, J. P., Sonner, J. M., Delarue, M., and Corringer, P. J. (2011) Nature 469, 428–431). In silico docking using a GLIC closed channel homology model suggests propofol binds to intersubunit sites in the TMD in the resting state. Propofol-induced motions in the intersubunit cavity were distinct from motions associated with channel activation, indicating propofol stabilizes a novel closed state.     

hERG potassium channel gating is mediated by N- and C-terminal region interactions

Human ether-á-go-go-related gene (hERG) potassium channels have voltage-dependent closing (deactivation) kinetics that are unusually slow. A Per-Arnt-Sim (PAS) domain in the cytoplasmic N-terminal region of hERG regulates slow deactivation by making a direct interaction with another part of the hERG channel. The mechanism for slow deactivation is unclear, however, because the other regions of the channel that participate in regulation of deactivation are not known. To identify other functional determinants of slow deactivation, we generated hERG channels with deletions of the cytoplasmic C-terminal regions. We report that hERG channels with deletions of the cyclic nucleotide-binding domain (CNBD) had accelerated deactivation kinetics that were similar to those seen in hERG channels lacking the PAS domain. Channels with dual deletions of the PAS domain and the CNBD did not show further acceleration in deactivation, indicating that the PAS domain and the CNBD regulate deactivation by a convergent mechanism. A recombinant PAS domain that we previously showed could directly regulate PAS domain-deleted channels did not regulate channels with dual deletions of the PAS domain and CNBD, suggesting that the PAS domain did not interact with CNBD-deleted channels. Biochemical protein interaction assays showed that glutathione S-transferase (GST)-PAS (but not GST) bound to a CNBD-containing fusion protein. Coexpression of PAS domain-deleted subunits (with intact C-terminal regions) and CNBD-deleted subunits (with intact N-terminal regions) resulted in channels with partially restored slow deactivation kinetics, suggesting regulatory intersubunit interactions between PAS domains and CNBDs. Together, these data suggest that the mechanism for regulation of slow deactivation in hERG channels is an interaction between the N-terminal PAS domain and the C-terminal CNBD.     

Letter: functional form of RNA synthesis termination factor rho

The ability of rho factor to effect in vitro termination of RNA synthesis requires the association of rho monomers into an oligomer consisting of at least four subunits. It was found that (1) rho factor activity has a sigmoidal dependence upon concentration, and (2) rho factor's sedimentation coefficient decreases with decreasing concentration, proceeding through 8.5 S (tetramer), 6.2 S (dimer) and finally 4.3 S (monomer) forms as the concentration approaches the apparent equilibrium binding constant. Rho factor may function at specific sites in the DNA template through the co-operative binding of subunits into an oligomer which surrounds the DNA helix.     

Anesthetic Binding in a Pentameric Ligand-Gated Ion Channel: GLIC

Cys-loop receptors are molecular targets of general anesthetics, but the knowledge of anesthetic binding to these proteins remains limited. Here we investigate anesthetic binding to the bacterial Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC), a structural homolog of cys-loop receptors, using an experimental and computational hybrid approach. Tryptophan fluorescence quenching experiments showed halothane and thiopental binding at three tryptophan-associated sites in the extracellular (EC) domain, transmembrane (TM) domain, and EC-TM interface of GLIC. An additional binding site at the EC-TM interface was predicted by docking analysis and validated by quenching experiments on the N200W GLIC mutant. The binding affinities (K_D) of 2.3 ± 0.1 mM and 0.10 ± 0.01 mM were derived from the fluorescence quenching data of halothane and thiopental, respectively. Docking these anesthetics to the original GLIC crystal structure and the structures relaxed by molecular dynamics simulations revealed intrasubunit sites for most halothane binding and intersubunit sites for thiopental binding. Tryptophans were within reach of both intra- and intersubunit binding sites. Multiple molecular dynamics simulations on GLIC in the presence of halothane at different sites suggested that anesthetic binding at the EC-TM interface disrupted the critical interactions for channel gating, altered motion of the TM23 linker, and destabilized the open-channel conformation that can lead to inhibition of GLIC channel current. The study has not only provided insights into anesthetic binding in GLIC, but also demonstrated a successful fusion of experiments and computations for understanding anesthetic actions in complex proteins.     

Structures of D-xylose isomerase from Arthrobacter strain B3728 containing the inhibitors xylitol and D-sorbitol at 2.5 A and 2.3 A resolution, respectively

The structures of D-xylose isomerase from Arthrobacter strain B3728 containing the polyol inhibitors xylitol and D-sorbitol have been solved at 2.5 A and 2.3 A, respectively. The structures have been refined using restrained least-squares refinement methods. The final crystallographic R-factors for the D-sorbitol (xylitol) bound molecules, for 43,615 (32,989) reflections are 15.6 (14.7). The molecule is a tetramer and the asymmetric unit of the crystal contains a dimer, the final model of which, incorporates a total of 6086 unique protein, inhibitor and magnesium atoms together with 535 bound solvent molecules. Each subunit of the enzyme contains two domains: the main domain is a parallel-stranded alpha-beta barrel, which has been reported in 14 other enzymes. The C-terminal domain is a loop structure consisting of five helical segments and is involved in intermolecular contacts between subunits that make up the tetramer. The structures have been analysed with respect to molecular symmetry, intersubunit contacts, inhibitor binding and active site geometry. The refined model shows the two independent subunits to be similar apart from local deviations due to solvent contacts in the solvent-exposed helices. The enzyme is dependent on a divalent cation for catalytic activity. Two metal ions are required per monomer, and the high-affinity magnesium(II) site has been identified from the structural results presented here. The metal ion is complexed, at the high-affinity site, by four carboxylate side-chains of the conserved residues, Glu180, Glu216, Asp244 and Asp292. The inhibitor polyols are bound in the active site in an extended open chain conformation and complete an octahedral co-ordination shell for the magnesium cation via their oxygen atoms O-2 and O-4. The active site lies in a deep pocket near the C-terminal ends of the beta-strands of the barrel domain and includes residues from a second subunit. The tetrameric molecule can be considered to be a dimer of "active" dimers, the active sites being composed of residues from both subunits. The analysis has revealed the presence of several internal salt-bridges stabilizing the tertiary and quaternary structure. One of these, between Asp23 and Arg139, appears to play a key role in stabilizing the active dimer and is conserved in the known sequences of this enzyme.(ABSTRACT TRUNCATED AT 400 WORDS)     

MCM ring hexamerization is a prerequisite for DNA-binding

The hexameric Minichromosome Maintenance (MCM) protein complex forms a ring that unwinds DNA at the replication fork in eukaryotes and archaea. Our recent crystal structure of an archaeal MCM N-terminal domain bound to single-stranded DNA (ssDNA) revealed ssDNA associating across tight subunit interfaces but not at the loose interfaces, indicating that DNA-binding is governed not only by the DNA-binding residues of the subunits (MCM ssDNA-binding motif, MSSB) but also by the relative orientation of the subunits. We now extend these findings by showing that DNA-binding by the MCM N-terminal domain of the archaeal organism Pyrococcus furiosus occurs specifically in the hexameric oligomeric form. We show that mutants defective for hexamerization are defective in binding ssDNA despite retaining all the residues observed to interact with ssDNA in the crystal structure. One mutation that exhibits severely defective hexamerization and ssDNA-binding is at a conserved phenylalanine that aligns with the mouse Mcm4(Chaos3) mutation associated with chromosomal instability, cancer, and decreased intersubunit association.     

Access Path to the Ligand Binding Pocket May Play a Role in Xenobiotics Selection by AhR

Understanding of multidrug binding at the atomic level would facilitate drug design and strategies to modulate drug metabolism, including drug transport, oxidation, and conjugation. Therefore we explored the mechanism of promiscuous binding of small molecules by studying the ligand binding domain, the PAS-B domain of the aryl hydrocarbon receptor (AhR). Because of the low sequence identities of PAS domains to be used for homology modeling, structural features of the widely employed HIF-2α and a more recent suitable template, CLOCK were compared. These structures were used to build AhR PAS-B homology models. We performed molecular dynamics simulations to characterize dynamic properties of the PAS-B domain and the generated conformational ensembles were employed in in silico docking. In order to understand structural and ligand binding features we compared the stability and dynamics of the promiscuous AhR PAS-B to other PAS domains exhibiting specific interactions or no ligand binding function. Our exhaustive in silico binding studies, in which we dock a wide spectrum of ligand molecules to the conformational ensembles, suggest that ligand specificity and selection may be determined not only by the PAS-B domain itself, but also by other parts of AhR and its protein interacting partners. We propose that ligand binding pocket and access channels leading to the pocket play equally important roles in discrimination of endogenous molecules and xenobiotics.