Dimerization Binding Research Articles

Non-Covalent Interactions in Molecular Crystals: Exploring the Accuracy of the Exchange-Hole Dipole Moment Model with Local Orbitals.

We present the first implementation of the exchange-hole dipole moment (XDM) model in combination with a numerical finite-support local orbital method (the SIESTA method) for the modeling of non-covalent interactions in periodic solids. The XDM model is parametrized for both the B86bPBE and PBE functionals using double-ζ- and triple-ζ-quality basis sets (DZP and TZP). The use of finite-support local orbitals is shown to have minimal impact on the computed dispersion coefficients for van der Waals molecular dimers and small molecular solids. However, the quality of the basis set affects the accuracy of calculated dimer binding energies and molecular-crystal lattice energies quite significantly; the size of the counterpoise correction indicates that this is caused by basis-set incompleteness error. In the case of the DZP basis set, its performance for weakly bound gas-phase dimers is similar to that of a double-ζ Gaussian basis set without diffuse functions. The new XDM implementation was tested on graphite and phosphorene exfoliation, and on the X23 benchmark set of molecular-crystal lattice energies. Our results indicate that lattice energies similar to plane-wave calculations can be obtained only if the counterpoise correction is applied. Alternatively, the calculated equilibrium geometries are reasonably close to the plane-wave equivalents, and composite approaches in which a single-point plane-wave calculation is used at the XDM/DZP equilibrium geometry yield good accuracy at a significantly lower computational cost.

A single XLF dimer bridges DNA ends during nonhomologous end joining.

Non-homologous end joining (NHEJ) is the primary pathway of DNA double-strand break repair in vertebrate cells, yet it remains unclear how NHEJ factors assemble a synaptic complex that bridges DNA ends. To address the role of XRCC4-like factor (XLF) in synaptic complex assembly, we employed single-molecule fluorescence imaging in Xenopus laevis egg extract, a system that efficiently joins DNA ends. We find that a single XLF dimer binds to DNA substrates just prior to formation of a ligation-competent synaptic complex between DNA ends. The interaction of both globular head domains of the XLF dimer with XRCC4 is required for efficient formation of this synaptic complex. In contrast to a model in which filaments of XLF and XRCC4 bridge DNA ends, our results indicate that binding of a single XLF dimer facilitates the assembly of a stoichiometrically well-defined synaptic complex.

Highlights from Recent Cancer Literature

Highlights from Recent Cancer Literature

A massively parallel reporter assay reveals context-dependent activity of homeodomain binding sites in vivo.

Cone-rod homeobox (CRX) is a paired-like homeodomain transcription factor (TF) and a master regulator of photoreceptor development in vertebrates. The in vitro DNA binding preferences of CRX have been described in detail, but the degree to which in vitro binding affinity is correlated with in vivo enhancer activity is not known. In addition, paired-class homeodomain TFs can bind DNA cooperatively as both homodimers and heterodimers at inverted TAAT half-sites separated by 2 or 3 nucleotides. This dimeric configuration is thought to mediate target specificity, but whether monomeric and dimeric sites encode distinct levels of activity is not known. Here, we used a massively parallel reporter assay to determine how local sequence context shapes the regulatory activity of CRX binding sites in mouse photoreceptors. We assayed inactivating mutations in more than 1700 TF binding sites and found that dimeric CRX binding sites act as stronger enhancers than monomeric CRX binding sites. Furthermore, the activity of dimeric half-sites is cooperative, dependent on a strict 3-bp spacing, and tuned by the identity of the spacer nucleotides. Saturating single-nucleotide mutagenesis of 195 CRX binding sites showed that, on average, changes in TF binding site affinity are correlated with changes in regulatory activity, but this relationship is obscured when considering mutations across multiple cis-regulatory elements (CREs). Taken together, these results demonstrate that the activity of CRX binding sites is highly dependent on sequence context, providing insight into photoreceptor gene regulation and illustrating functional principles of homeodomain binding sites that may be conserved in other cell types.

Dimerization of Coronavirus nsp9 with Diverse Modes Enhances Its Nucleic Acid Binding Affinity.

Coronaviruses pose serious health threats to humans and other animals. Understanding the mechanisms of their replication has important implications for global health and economic stability. Nonstructural protein 9 (nsp9) is an essential RNA binding protein for coronavirus replication. However, the mechanisms of the dimerization and nucleic acid binding of nsp9 remain elusive. Here, we report four crystal structures, including wild-type porcine delta coronavirus (PDCoV) nsp9, PDCoV nsp9-ΔN7 (N-terminal 7 amino acids deleted), wild-type porcine epidemic diarrhea virus (PEDV) nsp9, and PEDV nsp9-C59A mutant. These structures reveal the diverse dimerization forms of coronavirus nsp9. We first found that the N-finger of nsp9 from PDCoV plays a critical role in dimerization. Meanwhile, PEDV nsp9 is distinguished by the presence of a disulfide bond in the dimer interface. Interestingly, size exclusion chromatography and analytical ultracentrifugation analyses indicate that the PDCoV nsp9-ΔN7 and PEDV nsp9-C59A mutants are monomeric in solution. In addition, electrophoretic mobility shift assays and microscale thermophoresis analysis indicate that the monomeric forms of PDCoV nsp9 and PEDV nsp9 still have nucleic acid binding affinity, although it is lower than that of the wild type. Our results show that the diverse dimerization forms of coronavirus nsp9 proteins enhance their nucleic acid binding affinity.IMPORTANCE Coronaviruses cause widespread respiratory, gastrointestinal, and central nervous system diseases in humans and other animals, threatening human health and causing economic loss. Coronavirus nsp9, a member of the replication complex, is an important RNA binding subunit in the RNA-synthesizing machinery of all coronaviruses. However, the mechanisms of the dimerization and nucleic acid binding of nsp9 remain elusive. In this study we determined the nsp9 crystal structures of PDCoV and PEDV. We first found that the N-finger of nsp9 from PDCoV plays a critical role in dimerization. Meanwhile, PEDV nsp9 is distinguished by the presence of a disulfide bond in the dimer interface. This study provides a structural and functional basis for understanding the mechanism of dimerization and shows that the diverse dimerization modes of coronavirus nsp9 proteins enhance their nucleic acid binding affinity. Importantly, these findings may provide a new insight for antiviral drug development.

Torsional Potentials of Glyoxal, Oxalyl Halides, and Their Thiocarbonyl Derivatives: Challenges for Popular Density Functional Approximations.

The reliability of popular density functionals was studied for the description of torsional profiles of 36 molecules: glyoxal, oxalyl halides, and their thiocarbonyl derivatives. HF and 18 functionals of varying complexity, from local density to range-separated hybrid approximations and double-hybrid, have been considered and benchmarked against CCSD(T)-level rotational profiles. For molecules containing heavy halogens, most functionals fail to reproduce barrier heights accurately and a number of functionals introduce spurious minima. Dispersion corrections show no improvement. Calibrated torsion-corrected atom-centered potentials rectify the shortcomings of PBE and also improve on σ-hole based intermolecular binding in dimers and crystals.

Molecular model of J-aggregated pseudoisocyanine fibers.

Aggregated cyanines form ordered supramolecular structures with the potential to transport energy efficiently over long distances, a hallmark of photosynthetic light-harvesting complexes. In concentrated aqueous solution, pseudoisocyanine (PIC) spontaneously forms fibers with a chiral J-band red-shifted 1600 cm-1 from the monomeric 0-0 transition. A cryogenic transmission electron microscopy analysis of these fibers show an average fiber width of 2.89 nm, although the molecular-level structure of the aggregate is currently unknown. To determine a molecular model for these PIC fibers, the calculated spectra and dynamics using a Frenkel exciton model are compared to experiment. A chiral aggregate model in which the PIC monomers are neither parallel nor orthogonal to the long axis of the fiber is shown to replicate the experimental spectra most closely. This model can be physically realized by the sequential binding of PIC dimers and monomers to the ends of the fiber. These insights into the molecular aggregation model for aqueous PIC can also be applied to other similar cyanine-based supramolecular complexes with the potential for long-range energy transport, a key building block for the rational design of novel excitonic systems.

The mitochondrial single-stranded DNA binding protein from S. cerevisiae, Rim1, does not form stable homo-tetramers and binds DNA as a dimer of dimers.

Rim1 is the mitochondrial single-stranded DNA binding protein in Saccharomyces cerevisiae and functions to coordinate replication and maintenance of mtDNA. Rim1 can form homo-tetramers in solution and this species has been assumed to be solely responsible for ssDNA binding. We solved structures of tetrameric Rim1 in two crystals forms which differ in the relative orientation of the dimers within the tetramer. In testing whether the different arrangement of the dimers was due to formation of unstable tetramers, we discovered that while Rim1 forms tetramers at high protein concentration, it dissociates into a smaller oligomeric species at low protein concentrations. A single point mutation at the dimer–dimer interface generates stable dimers and provides support for a dimer–tetramer oligomerization model. The presence of Rim1 dimers in solution becomes evident in DNA binding studies using short ssDNA substrates. However, binding of the first Rim1 dimer is followed by binding of a second dimer, whose affinity depends on the length of the ssDNA. We propose a model where binding of DNA to a dimer of Rim1 induces tetramerization, modulated by the ability of the second dimer to interact with ssDNA.

DENR-MCTS1 heterodimerization and tRNA recruitment are required for translation reinitiation.

The succession of molecular events leading to eukaryotic translation reinitiation—whereby ribosomes terminate translation of a short open reading frame (ORF), resume scanning, and then translate a second ORF on the same mRNA—is not well understood. Density-regulated reinitiation and release factor (DENR) and multiple copies in T-cell lymphoma-1 (MCTS1) are implicated in promoting translation reinitiation both in vitro in translation extracts and in vivo. We present here the crystal structure of MCTS1 bound to a fragment of DENR. Based on this structure, we identify and experimentally validate that DENR residues Glu42, Tyr43, and Tyr46 are important for MCTS1 binding and that MCTS1 residue Phe104 is important for tRNA binding. Mutation of these residues reveals that DENR-MCTS1 dimerization and tRNA binding are both necessary for DENR and MCTS1 to promote translation reinitiation in human cells. These findings thereby link individual residues of DENR and MCTS1 to specific molecular functions of the complex. Since DENR–MCTS1 can bind tRNA in the absence of the ribosome, this suggests the DENR–MCTS1 complex could recruit tRNA to the ribosome during reinitiation analogously to the eukaryotic initiation factor 2 (eIF2) complex in cap-dependent translation.

Truncating Variants in NAA15 Are Associated with Variable Levels of Intellectual Disability, Autism Spectrum Disorder, and Congenital Anomalies

N-alpha-acetylation is a common co-translational protein modification that is essential for normal cell function in humans. We previously identified the genetic basis of an X-linked infantile lethal Mendelian disorder involving a c.109T>C (p.Ser37Pro) missense variant in NAA10, which encodes the catalytic subunit of the N-terminal acetyltransferase A (NatA) complex. The auxiliary subunit of the NatA complex, NAA15, is the dimeric binding partner for NAA10. Through a genotype-first approach with whole-exome or genome sequencing (WES/WGS) and targeted sequencing analysis, we identified and phenotypically characterized 38 individuals from 33 unrelated families with 25 different de novo or inherited, dominantly acting likely gene disrupting (LGD) variants in NAA15. Clinical features of affected individuals with LGD variants in NAA15 include variable levels of intellectual disability, delayed speech and motor milestones, and autism spectrum disorder. Additionally, mild craniofacial dysmorphology, congenital cardiac anomalies, and seizures are present in some subjects. RNA analysis in cell lines from two individuals showed degradation of the transcripts with LGD variants, probably as a result of nonsense-mediated decay. Functional assays in yeast confirmed a deleterious effect for two of the LGD variants in NAA15. Further supporting a mechanism of haploinsufficiency, individuals with copy-number variant (CNV) deletions involving NAA15 and surrounding genes can present with mild intellectual disability, mild dysmorphic features, motor delays, and decreased growth. We propose that defects in NatA-mediated N-terminal acetylation (NTA) lead to variable levels of neurodevelopmental disorders in humans, supporting the importance of the NatA complex in normal human development.

The Hcp-like protein HilE inhibits homodimerization and DNA binding of the virulence-associated transcriptional regulator HilD in Salmonella

HilD is an AraC-like transcriptional regulator that plays a central role in Salmonella virulence. HilD controls the expression of the genes within the Salmonella pathogenicity island 1 (SPI-1) and of several genes located outside SPI-1, which are mainly required for Salmonella invasion of host cells. The expression, amount, and activity of HilD are tightly controlled by the activities of several factors. The HilE protein represses the expression of the SPI-1 genes through its interaction with HilD; however, the mechanism by which HilE affects HilD is unknown. In this study, we used genetic and biochemical assays revealing how HilE controls the transcriptional activity of HilD. We found that HilD needs to assemble in homodimers to induce expression of its target genes. Our results further indicated that HilE individually interacts with each the central and the C-terminal HilD regions, mediating dimerization and DNA binding, respectively. We also observed that these interactions consistently inhibit HilD dimerization and DNA binding. Interestingly, a computational analysis revealed that HilE shares sequence and structural similarities with Hcp proteins, which act as structural components of type 6 secretion systems in Gram-negative bacteria. In conclusion, our results uncover the molecular mechanism by which the Hcp-like protein HilE controls dimerization and DNA binding of the virulence-promoting transcriptional regulator HilD. Our findings may indicate that HilE's activity represents a functional adaptation during the evolution of Salmonella pathogenicity.

Patient‐Specific Mechanisms of Lafora Disease Mutations in the Human Glycogen Phosphatase

Lafora disease (LD) is a fatal neurodegenerative epilepsy caused by defects in glycogen metabolism. Classically, LD manifests in adolescence with progressive myoclonus epilepsy, rapid neurodegeneration and death 10 years after onset. Insoluble inclusions of abnormal glycogen known as Lafora bodies (LBs) are found in the brain and other tissues of Lafora patients. Genetic ablation of glycogen synthesis eliminates LBs and rescues LD in mice, indicating LBs are the pathological cause of LD.50% of LD cases are caused by recessive mutations in the Epilepsy Progressive Myoclonus 2A (EPM2A) gene encoding laforin, the glycogen phosphatase. Laforin is a bimodular protein containing a carbohydrate binding module (CBM) and a dual specificity phosphatase domain. In the absence of laforin, glycogen becomes hyperphosphorylated and abberantly branched, leading to LB formation. >90 distinct missense, nonsense, and frameshift mutations have been described in LD patients. We previously determined the crystal structure of laforin and defined the structural basis of its function.The objective of this study is to define the molecular defect(s) of >30 pathogenic mutations from LD patients and determine why some are linked to milder phenotypes.Using the laforin crystal structure and a variety of cellular, biochemical, and biophysical tools, we discovered a hierarchy of mutations that differentially affect aspects of laforin function. Specifically, we discovered that a large group of pathogenic mutations cluster in the CBM, leading to destabilization, domain decoupling, and reduced protein‐protein interactions. Further, a mutation associated with an extremely mild case of LD impairs dimerization and preferential binding to abnormal glucan substrates.Our study reveals the importance of the CBM for both stability and glucan binding, defines the atomic‐level effects of patient mutations, and provides a biochemical explanation for late‐onset phenotypes. Our work not only provides a framework for defining novel LD mutations and understanding patient‐specific disease progression, but also deepens our understanding of the enzymatic regulation of glycogen.Support or Funding InformationThis work was supported by NIH Grants R01 NS070899 (M.S.G.), P01 NS097197 (M.S.G.), and F31 NS093892 (M.K.B.); NSF Grant IIA‐1355438 (M.S.G.); Mizutani Foundation for Glycoscience Award (M.S.G.).This abstract is from the Experimental Biology 2018 Meeting. There is no full text article associated with this abstract published in The FASEB Journal.

Mechanism of agonism and antagonism of the Pseudomonas aeruginosa quorum sensing regulator QscR with non‐native ligands

Pseudomonas aeruginosa is an opportunistic pathogen that uses the process of quorum sensing (QS) to coordinate the expression of many virulence genes. During quorum sensing, N-acyl-homoserine lactone (AHL) signaling molecules regulate the activity of three LuxR-type transcription factors, LasR, RhlR and QscR. To better understand P. aeruginosa QS signal reception, we examined the mechanism underlying the response of QscR to synthetic agonists and antagonists using biophysical and structural approaches. The structure of QscR bound to a synthetic agonist reveals a novel mode of ligand binding supporting a general mechanism for agonist activity. In turn, antagonists of QscR with partial agonist activity were found to destabilize and greatly impair QscR dimerization and DNA binding. These results highlight the diversity of LuxR-type receptor responses to small molecule agonists and antagonists and demonstrate the potential for chemical strategies for the selective targeting of individual QS systems.

The binding of DCC-P3 motif and FAK-FAT domain mediates the initial step of netrin-1/DCC signaling for axon attraction

Netrin-1 plays a key role in axon guidance through binding to its receptor, Deleted in Colorectal Cancer (DCC). The initial step of signaling inside the cell after netrin-1/DCC ligation is the binding of DCC cytoplasmic P3 motif to focal adhesion targeting (FAT) domain of focal adhesion kinase (FAK). Here we report the crystal structure of P3/FAT complex. The helical P3 peptide interacts with a helix-swapped FAT dimer in a 2:2 ratio. Dimeric FAT binding is P3-specific and stabilized by a calcium ion. Biochemical studies showed that DCC-P3 motif and calcium ion could facilitate FAT dimerization in solution. Axon guidance assays confirm that the DCC/FAK complex is essential for netrin-1-induced chemoattraction. We propose that netrin-1/DCC engagement creates a small cluster of P3/FAT for FAK recruitment close to the cell membrane, which exerts a concerted effect with PIP2 for FAK signaling. We also compare P3/FAT binding with paxillin/FAT binding and discuss their distinct recognition specificity on a common FAT domain for axon attraction versus integrin signaling, respectively.

Neuropathic MORC2 mutations perturb GHKL ATPase dimerization dynamics and epigenetic silencing by multiple structural mechanisms

Missense mutations in MORC2 cause neuropathies including spinal muscular atrophy and Charcot–Marie–Tooth disease. We recently identified MORC2 as an effector of epigenetic silencing by the human silencing hub (HUSH). Here we report the biochemical and cellular activities of MORC2 variants, alongside crystal structures of wild-type and neuropathic forms of a human MORC2 fragment comprising the GHKL-type ATPase module and CW-type zinc finger. This fragment dimerizes upon binding ATP and contains a hinged, functionally critical coiled-coil insertion absent in other GHKL ATPases. We find that dimerization and DNA binding of the MORC2 ATPase module transduce HUSH-dependent silencing. Disease mutations change the dynamics of dimerization by distinct structural mechanisms: destabilizing the ATPase-CW module, trapping the ATP lid, or perturbing the dimer interface. These defects lead to the modulation of HUSH function, thus providing a molecular basis for understanding MORC2-associated neuropathies.

Computationally Determined Free Energy Profiles of the Mucin-1 Transmembrane Homodimer

Determining the principles governing lipid-protein interactions, specifically those in single-pass transmembrane protein dimerization, are fundamental problems in the field of membrane biophysics. Mucin-1 (MUC1), a mucus constituent found on the apical side of epithelial cells, typically acts as a barrier against pathogens but has also been identified as an oncogene. When overexpressed, the C-terminal subunit of MUC1 is known to homodimerize and undergo nuclear translocation as well as interact with growth factors and other kinases in signaling cascades to promote metastasis, inflammation, and drug resistance. MUC1 homodimerization in the membrane occurs due to a pair of disulfide bonds formed between two adjacent monomers. Strong dimerization affinity in bacterial membranes comparable to glycophorin A has been observed experimentally, and weak affinity is seen even without the cytosolic membrane-proximal cysteine residues. Several mutations in the TMD, most notably Ala-1180-Leu, also affect dimerization. All-atom molecular dynamics simulations of the transmembrane and juxtamembrane domains of MUC1 have been used to better quantify the dimer structure and binding energetics. Umbrella sampling and free energy perturbation methods were used to generate potential of mean force (PMF) curves, which were then compared to experimental assay activity.

Multiscale Simulations of Membrane Recognition by Lipid Kinases

Peripheral membrane proteins associate with lipids in the plasma membrane in many important cell signalling processes. Type I phosphatidylinositol phosphate kinase A (PIP5K1A) is one such example, which phosphorylates the head group of the lipid phosphatidylinositol 4-phosphate (PI4P) to generate phosphatidylinositol (4,5)-bisphosphate (PIP2). The dynamic interactions of peripheral proteins with their target membranes is difficult to capture on experimentally. A recent crystal structure at 3.3 Å resolution enables the further investigation of the structure and dynamics of PIP5K1A and its interactions with lipids in the plasma membrane. Here we use coarse-grained (CG) and atomistic (AT) MD simulations of the PIP5K1A protein kinase to investigate the nature of monomer and dimer binding to the plasma membrane and their effects on lipid clustering. The AT simulations suggest that whilst dimerisation of the kinase results in an increase in the area of the binding surface, the monomers may not bind concurrently. This in turn suggests that an allosteric mechanism for stabilisation of the productive binding conformation is more likely. We have also used simulations to investigate the relationship between conformational dynamics of the activation loop and interactions with lipids in the cell membrane. Furthermore, large-scale CG simulations (1 million particles) reveal how lipid composition and membrane curvature relate to the interactions of an ensemble of 25 kinases with a model membrane. PIP kinases are implicated in many cancers as a result of their influence on cell processes. Studies in mice and human tissue samples suggest that inhibition of PIP5K1A is a potential target for treatment. Simulations of the interactions of PIP5K1A with realistic models of cell membranes will aid design of drugs targeted at these and other related lipid kinases.

Modular discovery of monomeric and dimeric transcription factor binding motifs for large data sets.

In some dimeric cases of transcription factor (TF) binding, the specificity of dimeric motifs has been observed to differ notably from what would be expected were the two factors to bind to DNA independently of each other. Current motif discovery methods are unable to learn monomeric and dimeric motifs in modular fashion such that deviations from the expected motif would become explicit and the noise from dimeric occurrences would not corrupt monomeric models. We propose a novel modeling technique and an expectation maximization algorithm, implemented as software tool MODER, for discovering monomeric TF binding motifs and their dimeric combinations. Given training data and seeds for monomeric motifs, the algorithm learns in the same probabilistic framework a mixture model which represents monomeric motifs as standard position-specific probability matrices (PPMs), and dimeric motifs as pairs of monomeric PPMs, with associated orientation and spacing preferences. For dimers the model represents deviations from pure modular model of two independent monomers, thus making co-operative binding effects explicit. MODER can analyze in reasonable time tens of Mbps of training data. We validated the tool on HT-SELEX and ChIP-seq data. Our findings include some TFs whose expected model has palindromic symmetry but the observed model is directional.

Impact of intracellular ionic strength on dimer binding in the NF-kB Inducing kinase.

Improper signaling of the nuclear factor-κB (NF-κB) pathway plays a critical role in many inflammatory disease states including cancer, stroke, and viral infections. Although the signaling pathways are known, how these molecular mechanisms respond to changes in the intracellular microenvironment such as pH, ionic strength, and temperature, remains elusive. Molecular dynamics simulations were employed to differentiate the structural dynamics of the NF-κB Inducing Kinase (NIK), a protein kinase responsible for invoking the non-canonical NF-κB pathway, in its native and mutant form, and in the absence and presence of salt concentration in efforts to probe whether changes in the ionic environment stabilize or destabilize the NIK dimer. Analyses of structure-activity and conformational-activity relationships indicate that the protein-protein interactions are sensitive to changes in the ionic strength. Ligand binding pockets as well as regions between the oligomer interface either compress or expand, affecting both local and distal intermolecular interactions that result in stabilization or destabilization in the protein assembly.

Molecular dynamics simulation unveils the conformational flexibility of the interdomain linker in the bacterial transcriptional regulator GabR from Bacillus subtilis bound to pyridoxal 5'-phosphate.

GabR from Bacillus subtilis is a transcriptional regulator belonging to the MocR subfamily of the GntR regulators. The structure of the MocR regulators is characterized by the presence of two domains: i) a N-terminal domain, about 60 residue long, possessing the winged-Helix-Turn-Helix (wHTH) architecture with DNA recognition and binding capability; ii) a C-terminal domain (about 350 residue) folded as the pyridoxal 5’-phosphate (PLP) dependent aspartate aminotransferase (AAT) with dimerization and effector binding functions. The two domains are linked to each other by a peptide bridge. Although structural and functional characterization of MocRs is proceeding at a fast pace, virtually nothing is know about the molecular changes induced by the effector binding and on how these modifications influence the properties of the regulator. An extensive molecular dynamics simulation on the crystallographic structure of the homodimeric B. subtilis GabR has been undertaken with the aim to envisage the role and the importance of conformational flexibility in the action of GabR. Molecular dynamics has been calculated for the apo (without PLP) and holo (with PLP bound) forms of the GabR. A comparison between the molecular dynamics trajectories calculated for the two GabR forms suggested that one of the wHTH domain detaches from the AAT-like domain in the GabR PLP-bound form. The most evident conformational change in the holo PLP-bound form is represented by the rotation and the subsequent detachment from the subunit surface of one of the wHTH domains. The movement is mediated by a rearrangement of the linker connecting the AAT domain possibly triggered by the presence of the negative charge of the PLP cofactor. This is the second most significant conformational modification. The C-terminal section of the linker docks into the “active site” pocket and establish stabilizing contacts consisting of hydrogen-bonds, salt-bridges and hydrophobic interactions.