LSD1/CoREST reversible opening-closing dynamics: discovery of a nanoscale clamp for chromatin and protein binding

LSD1 associated with its corepressor protein CoREST is an exceptionally relevant target for epigenetic drugs. Hypotheses for the role of LSD1/CoREST as a multidocking site for chromatin and protein binding would require significant molecular flexibility, and LSD1/CoREST large-amplitude conformational dynamics is currently unknown. Here, molecular dynamics simulation reveals that the LSD1/CoREST complex in solution functions as a reversible nanoscale binding clamp. We show that the H3 histone tail binding pocket is a potential allosteric site for regulation of the rotation of SWIRM/SANT2 domains around the Tower domain. Thus, targeting this site and including receptor flexibility are crucial strategies for future drug discovery.     

Thermodynamic characterization of the binding interaction between the histone demethylase LSD1/KDM1 and CoREST

Flavin-dependent histone demethylases catalyze the posttranslational oxidative demethylation of mono- and dimethylated lysine residues, producing formaldehyde and hydrogen peroxide in addition to the corresponding demethylated protein. In vivo, histone demethylase LSD1 (KDM1; BCH110) is a component of the multiprotein complex that includes histone deacetylases (HDAC 1 and 2) and the scaffolding protein CoREST. Although little is known about the affinities of or the structural basis for the interaction between CoREST and HDACs, the structure of CoREST(286-482) bound to an α-helical coiled-coil tower domain within LSD1 has recently been reported. Given the significance of CoREST in directing demethylation to specific nucleosomal substrates, insight into the molecular basis of the interaction between CoREST and LSD1 may suggest a new means of inhibiting LSD1 activity by misdirecting the enzyme away from nucleosomal substrates. Toward this end, isothermal titration calorimetry studies were conducted to determine the affinity and thermodynamic parameters characterizing the binding interaction between LSD1 and CoREST(286-482). The proteins tightly interact in a 1:1 stoichiometry with a dissociation constant (K(d)) of 15.9 ± 2.07 nM, and their binding interaction is characterized by a favorable enthalpic contribution near room temperature with a smaller entropic penalty at pH 7.4. Additionally, one proton is transferred from the buffer to the heterodimeric complex at pH 7.4. From the temperature dependence of the enthalpy change of interaction, a constant-pressure heat capacity change (ΔC(p)) of the interaction was determined to be -0.80 ± 0.01 kcal mol(-1) K(-1). Notably, structure-driven truncation of CoREST revealed that the central binding determinant lies within the segment of residues 293-380, also known as the CoREST "linker" region, which is a central isolated helix that interacts with the LSD1 coiled-coil tower domain to create a triple-helical bundle. Thermodynamic parameters obtained from the binding between LSD1 and the linker region of CoREST are similar to those obtained from the interaction between LSD1 and CoREST(286-482). These results provide a framework for understanding the molecular basis of protein-protein interactions that govern nucleosomal demethylation.     

Regulation of LSD1 histone demethylase activity by its associated factors

LSD1 is a recently identified human lysine (K)-specific histone demethylase. LSD1 is associated with HDAC1/2; CoREST, a SANT domain-containing corepressor; and BHC80, a PHD domain-containing protein, among others. We show that CoREST endows LSD1 with the ability to demethylate nucleosomal substrates and that it protects LSD1 from proteasomal degradation in vivo. We find hyperacetylated nucleosomes less susceptible to CoREST/LSD1-mediated demethylation, suggesting that hypoacetylated nucleosomes may be the preferred physiological substrates. This raises the possibility that histone deacetylases and LSD1 may collaborate to generate a repressive chromatin environment. Consistent with this model, TSA treatment results in derepression of LSD1 target genes. While CoREST positively regulates LSD1 function, BHC80 inhibits CoREST/LSD1-mediated demethylation in vitro and may therefore confer negative regulation. Taken together, these findings suggest that LSD1-mediated histone demethylation is regulated dynamically in vivo. This is expected to have profound effects on gene expression under both physiological and pathological conditions.     

Intermolecular insights into allosteric inhibition of histone lysine-specific demethylase 1

BackgroundHistone lysine-specific demethylase 1 (LSD1) has become a potential anticancer target for the novel drug discovery. Recent reports have shown that SP2509 and its derivatives strongly inhibit LSD1 as allosteric inhibitors. However, the binding mechanism of these allosteric inhibitors in the allosteric site of LSD1 is not known yet.MethodsThe stability and binding mechanism of allosteric inhibitors in the binding site of LSD1 were evaluated by molecular docking, ligand-based pharmacophore, molecular dynamics (MD) simulations, molecular mechanics generalized born surface area (MM/GBSA) analysis, quantum mechanics/molecular mechanics (QM/MM) calculation and Hirshfeld surface analysis.ResultsThe conformational geometry and the intermolecular interactions of allosteric inhibitors showed high binding affinity towards allosteric site of LSD1 with the neighboring amino acids (Gly358, Cys360, Leu362, Asp375 and Glu379). Meanwhile, MD simulations and MM/GBSA analysis were performed on selected allosteric inhibitors in complex with LSD1 protein, which confirmed the high stability and binding affinity of these inhibitors in the allosteric site of LSD1.ConclusionThe simulation results revealed the crucial factors accounting for allosteric inhibitors of LSD1, including different protein–ligand interactions, the positions and conformations of key residues, and the ligands flexibilities. Meanwhile, a halogen bond interaction between chlorine atom of ligand and key residues Trp531 and His532 was recurrent in our analysis confirming its importance.General significanceOverall, our research analyzed in depth the binding modes of allosteric inhibitors with LSD1 and could provide useful information for the design of novel allosteric inhibitors.     

Modeling of DNA binding to the condensin hinge domain using molecular dynamics simulations guided by atomic force microscopy

The condensin protein complex compacts chromatin during mitosis using its DNA-loop extrusion activity. Previous studies proposed scrunching and loop-capture models as molecular mechanisms for the loop extrusion process, both of which assume the binding of double-strand (ds) DNA to the hinge domain formed at the interface of the condensin subunits Smc2 and Smc4. However, how the hinge domain contacts dsDNA has remained unknown. Here, we conducted atomic force microscopy imaging of the budding yeast condensin holo-complex and used this data as basis for coarse-grained molecular dynamics simulations to model the hinge structure in a transient open conformation. We then simulated the dsDNA binding to open and closed hinge conformations, predicting that dsDNA binds to the outside surface when closed and to the outside and inside surfaces when open. Our simulations also suggested that the hinge can close around dsDNA bound to the inside surface. Based on these simulation results, we speculate that the conformational change of the hinge domain might be essential for the dsDNA binding regulation and play roles in condensin-mediated DNA-loop extrusion.     

Catalytic mechanism investigation of lysine-specific demethylase 1 (LSD1): a computational study

Lysine-specific demethylase 1 (LSD1), the first identified histone demethylase, is a flavin-dependent amine oxidase which specifically demethylates mono- or dimethylated H3K4 and H3K9 via a redox process. It participates in a broad spectrum of biological processes and is of high importance in cell proliferation, adipogenesis, spermatogenesis, chromosome segregation and embryonic development. To date, as a potential drug target for discovering anti-tumor drugs, the medical significance of LSD1 has been greatly appreciated. However, the catalytic mechanism for the rate-limiting reductive half-reaction in demethylation remains controversial. By employing a combined computational approach including molecular modeling, molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations, the catalytic mechanism of dimethylated H3K4 demethylation by LSD1 was characterized in details. The three-dimensional (3D) model of the complex was composed of LSD1, CoREST, and histone substrate. A 30-ns MD simulation of the model highlights the pivotal role of the conserved Tyr761 and lysine-water-flavin motif in properly orienting flavin adenine dinucleotide (FAD) with respect to substrate. The synergy of the two factors effectively stabilizes the catalytic environment and facilitated the demethylation reaction. On the basis of the reasonable consistence between simulation results and available mutagenesis data, QM/MM strategy was further employed to probe the catalytic mechanism of the reductive half-reaction in demethylation. The characteristics of the demethylation pathway determined by the potential energy surface and charge distribution analysis indicates that this reaction belongs to the direct hydride transfer mechanism. Our study provides insights into the LSD1 mechanism of reductive half-reaction in demethylation and has important implications for the discovery of regulators against LSD1 enzymes.     

Carbon monoxide binding to the heme group at the dimeric interface modulates structure and copper accessibility in the Cu,Zn superoxide dismutase from Haemophilus ducreyi: in silico and in vitro evidences

X-ray absorption near-edge structure (XANES) spectroscopy and molecular dynamics (MD) simulations have been jointly applied to the study of the Cu,Zn superoxide dismutase from Haemophilus ducreyi (HdSOD) in interaction with the carbon monoxide molecule. The configurational flexibility of the Fe(II)-heme group, intercalated between the two subunits, has been sampled by MD simulations and included in the XANES data analysis without optimization in the structural parameter space. Our results provide an interpretation of the observed discrepancy in the Fe-heme distances as detected by extended X-ray absorption fine structure (EXAFS) spectroscopy and the classical XANES analysis, in which the structural parameters are optimized in a unique structure. Moreover, binding of the CO molecule to the heme induces a long range effect on the Cu,Zn active site, as evidenced by both MD simulations and in vitro experiments. MD simulation of the CO bound system, in fact, highlighted a structural rearrangement of the protein-protein hydrogen bond network in the region of the Cu,Zn active site, correlated with an increase in water accessibility at short distance from the copper atom. In line, in vitro experiments evidenced an increase of copper accessibility to a chelating agent when the CO molecule binds to the heme group, as compared to a heme deprived HdSOD. Altogether, our results support the hypothesis that the HdSOD is a heme-sensor protein, in which binding to small gaseous molecules modulates the enzyme superoxide activity as an adaptive response to the bacterial environment.     

Interplay between structural rigidity and electrostatic interactions in the ligand binding domain of GluR2

Using molecular dynamics (MD) simulations, computational protein modifications, and a novel theoretical methodology that determines structural rigidity/flexibility (the FIRST algorithm), we investigate how molecular structure and dynamics of the glutamate receptor ligand binding domain (GluR2 S1S2) facilitate its conformational transition. S1S2 is a two-lobe protein, which undergoes a cleft closure conformational transition upon binding an agonist in the cleft between the two lobes; hence it is expected that the mechanism of this conformational transition can be characterized as a hinge-type. However, in the rigidity analysis one lobe of the protein is identified as a single rigid cluster while the other one is structurally flexible, inconsistent with a presumed mechanical hinge mechanism. Instead, we characterize the cleft-closing transition as a load and lock mechanism. We find that when two cross-cleft hydrogen bonds are disrupted the protein undergoes a rapid cleft opening transition. At the same time, the dynamical behavior of the cleft in the presence of the glutamate ligand is only weakly affected by the S652 peptide bond in its flipped conformation observed in the crystal structure. The residue E705 plays significant role in stabilization of the closed conformation via electrostatic interactions. The presence of the E705-K730 salt bridge seems to correlate strongly withthe cleft opening transition in the MD simulations.     

Computational investigation of the HIV-1 Rev multimerization using molecular dynamics simulations and binding free energy calculations

The HIV Rev protein mediates the nuclear export of viral mRNA, and is thereby essential for the production of late viral proteins in the replication cycle. Rev forms a large organized multimeric protein-protein complex for proper functioning. Recently, the three-dimensional structures of a Rev dimer and tetramer have been resolved and provide the basis for a thorough structural analysis of the binding interaction. Here, molecular dynamics (MD) and binding free energy calculations were performed to elucidate the forces thriving dimerization and higher order multimerization of the Rev protein. It is found that despite the structural differences between each crystal structure, both display a similar behavior according to our calculations. Our analysis based on a molecular mechanics-generalized Born surface area (MM/GBSA) and a configurational entropy approach demonstrates that the higher order multimerization site is much weaker than the dimerization site. In addition, a quantitative hot spot analysis combined with a mutational analysis reveals the most contributing amino acid residues for protein interactions in agreement with experimental results. Additional residues were found in each interface, which are important for the protein interaction. The investigation of the thermodynamics of the Rev multimerization interactions performed here could be a further step in the development of novel antiretrovirals using structure based drug design. Moreover, the variability of the angle between each Rev monomer as measured during the MD simulations suggests a role of the Rev protein in allowing flexibility of the arginine rich domain (ARM) to accommodate RNA binding.     

Molecular dynamics simulations of the Apo-, Holo-, and acyl-forms of Escherichia coli acyl carrier protein

Acyl carrier protein (ACP) is an essential co-factor protein in fatty acid biosynthesis that shuttles covalently bound fatty acyl intermediates in its hydrophobic pocket to various enzyme partners. To characterize acyl chain-ACP interactions and their influence on enzyme interactions, we performed 19 molecular dynamics (MD) simulations of Escherichia coli apo-, holo-, and acyl-ACPs. The simulations were started with the acyl chain in either a solvent-exposed or a buried conformation. All four short-chain (< or = C10) and one long-chain (C16) unbiased acyl-ACP MD simulation show the transition of the solvent-exposed acyl chain into the hydrophobic pocket of ACP, revealing its pathway of acyl chain binding. Although the acyl chain resides inside the pocket, Thr-39 and Glu-60 at the entrance stabilize the phosphopantetheine linker through hydrogen bonding. Comparisons of the different ACP forms indicate that the loop region between helices II and III and the prosthetic linker may aid in substrate recognition by enzymes of fatty acid synthase systems. The MD simulations consistently show that the hydrophobic binding pocket of ACP is best suited to accommodate an octanoyl group and is capable of adjusting in size to accommodate chain lengths as long as decanoic acid. The simulations also reveal a second, novel binding mode of the acyl chains inside the hydrophobic binding pocket directed toward helix I. This study provides a detailed dynamic picture of acyl-ACPs that is in excellent agreement with available experimental data and, thereby, provides a new understanding of enzyme-ACP interactions.     

Achieving Peptide Binding Specificity and Promiscuity by Loops: Case of the Forkhead-Associated Domain

The regulation of a series of cellular events requires specific protein–protein interactions, which are usually mediated by modular domains to precisely select a particular sequence from diverse partners. However, most signaling domains can bind to more than one peptide sequence. How do proteins create promiscuity from precision? Moreover, these complex interactions typically occur at the interface of a well-defined secondary structure, α helix and β sheet. However, the molecular recognition primarily controlled by loop architecture is not fully understood. To gain a deep understanding of binding selectivity and promiscuity by the conformation of loops, we chose the forkhead-associated (FHA) domain as our model system. The domain can bind to diverse peptides via various loops but only interact with sequences containing phosphothreonine (pThr). We applied molecular dynamics (MD) simulations for multiple free and bound FHA domains to study the changes in conformations and dynamics. Generally, FHA domains share a similar folding structure whereby the backbone holds the overall geometry and the variety of sidechain atoms of multiple loops creates a binding surface to target a specific partner. FHA domains determine the specificity of pThr by well-organized binding loops, which are rigid to define a phospho recognition site. The broad range of peptide recognition can be attributed to different arrangements of the loop interaction network. The moderate flexibility of the loop conformation can help access or exclude binding partners. Our work provides insights into molecular recognition in terms of binding specificity and promiscuity and helpful clues for further peptide design.