A highly ordered, nonprotective MALAT1 ENE structure is adopted prior to triplex formation.

The 3' end of the ∼7 kb lncRNA MALAT1 contains an evolutionarily and structurally conserved element for nuclear expression (ENE) which confers protection from cellular degradation pathways. Formation of an ENE triple helix is required to support transcript accumulation, leading to persistent oncogenic activity of MALAT1 in multiple cancer types. Though the specific mechanism of triplex-mediated protection remains unknown, the MALAT1 ENE triplex has been identified as a promising target for therapeutic intervention. Interestingly, a maturation step of the nascent lncRNA 3' end is required prior to triplex formation. We hypothesize that disruption of the maturation or folding process may be a viable mechanism of inhibition. To assess putative cotranscriptional ENE conformations prior to triplex formation, we perform microsecond MD simulations of a partially folded ENE conformation and the ENE triplex. We identify a highly ordered ENE structure prior to triplex formation. Extensive formation of U•U base pairs within the large U-rich internal loops produces a global rod-like architecture. We present a three-dimensional structure of the isolated ENE motif, the global features of which are consistent with small angle X-ray scattering (SAXS) experiments. Our structural model represents a nonprotective conformation of the MALAT1 ENE, providing a molecular description useful for future mechanistic and inhibition studies. We anticipate that targeting stretches of U•U pairs within the ENE motif will prove advantageous for the design of therapeutics targeting this oncogenic lncRNA.

Selective Small-Molecule Targeting of a Triple Helix Encoded by the Long Noncoding RNA, MALAT1.

Metastasis-associated lung adenocarcinoma transcript 1 ( Malat1/ MALAT1, mouse/human), a highly conserved long noncoding (lnc) RNA, has been linked with several physiological processes, including the alternative splicing, nuclear organization, and epigenetic modulation of gene expression. MALAT1 has also been implicated in metastasis and tumor proliferation in multiple cancer types. The 3' terminal stability element for nuclear expression (ENE) assumes a triple-helical configuration that promotes its nuclear accumulation and persistent function. Utilizing a novel small molecule microarray strategy, we identified multiple Malat1 ENE triplex-binding chemotypes, among which compounds 5 and 16 reduced Malat1 RNA levels and branching morphogenesis in a mammary tumor organoid model. Computational modeling and Förster resonance energy transfer experiments demonstrate distinct binding modes for each chemotype, conferring opposing structural changes to the triplex. Compound 5 modulates Malat1 downstream genes without affecting Neat1, a nuclear lncRNA encoded in the same chromosomal region as Malat1 with a structurally similar ENE triplex. Supporting this observation, the specificity of compound 5 for Malat1 over Neat1 and a virus-coded ENE was demonstrated by nuclear magnetic resonance spectroscopy. Small molecules specifically targeting the MALAT1 ENE triplex lay the foundation for new classes of anticancer therapeutics and molecular probes for the treatment and investigation of MALAT1-driven cancers.

Configurational Preference of the Glutamate Receptor Ligand Binding Domain Dimers.

Ionotropic glutamate receptors are a family of tetrameric ion channels with functional states consisting of nonconducting, conducting, and desensitized states that are starting to become well characterized by electrophysiological and biophysical studies. However, the structure and relative energetics of these states beyond the general structure of the receptor are still not well understood. It is known that the interface between monomeric subunits of the tetramer plays a major role in distinguishing these functional states. We have used umbrella sampling and multimicrosecond molecular dynamics simulations of the GluA2 AMPA subtype glutamate receptor ligand-binding domain (LBD) dimers to characterize a natural propensity of the LBD dimers for various configurational states. Our results show a proposed desensitized conformation of the LBD dimer is a highly preferable conformation of the LBD dimer without the influence of other receptor domains or crystallographic conditions. This has been demonstrated by both free protein simulations of 5 µs duration, as well as by computed free energy difference between the active and desensitized states. At the same time, the simulations performed using the same protocols revealed that for the LBD mutant L483Y, known to lack desensitization, the postulated active state of the LBD dimer is indeed the preferred configurational state, which remained stable in the simulations. Our findings pave the path for developing more detailed hypotheses of the full receptor activation mechanism. Combined with the energetics of glutamate binding to the LBD and the energy required to open the transmembrane pore helices, our results strongly support a hypothesis that the low absolute free-energy state is the desensitized state of the intact AMPA receptor.

All atom NMDA receptor transmembrane domain model development and simulations in lipid bilayers and water.

N-methyl-d-aspartate receptors (NMDARs) are members of the ionotropic glutamate receptor family that mediate excitatory synaptic transmission in the central nervous system. The channels of NMDARs are permeable to Ca2+ but blocked by Mg2+, distinctive properties that underlie essential brain processes such as induction of synaptic plasticity. However, due to limited structural information about the NMDAR transmembrane ion channel forming domain, the mechanism of divalent cation permeation and block is understood poorly. In this paper we developed an atomistic model of the transmembrane domain (TMD) of NMDARs composed of GluN1 and GluN2A subunits (GluN1/2A receptors). The model was generated using (a) a homology model based on the structure of the NaK channel and a partially resolved structure of an AMPA receptor (AMPAR), and (b) a partially resolved X-ray structure of GluN1/2B NMDARs. Refinement and extensive Molecular Dynamics (MD) simulations of the NMDAR TMD model were performed in explicit lipid bilayer membrane and water. Targeted MD with simulated annealing was introduced to promote structure refinement. Putative positions of the Mg2+ and Ca2+ ions in the ion channel divalent cation binding site are proposed. Differences in the structural and dynamic behavior of the channel protein in the presence of Mg2+ or Ca2+ are analyzed. NMDAR protein conformational flexibility was similar with no ion bound to the divalent cation binding site and with Ca2+ bound, whereas Mg2+ binding reduced protein fluctuations. While bound at the binding site both ions retained their preferred ligand coordination numbers: 6 for Mg2+, and 7-8 for Ca2+. Four asparagine side chain oxygens, a back-bone oxygen, and a water molecule participated in binding a Mg2+ ion. The Ca2+ ion first coordination shell ligands typically included four to five side-chain oxygen atoms of the binding site asparagine residues, two water molecules and zero to two backbone oxygens of the GluN2B subunits. These results demonstrate the importance of high-resolution channel structures for elucidation of mechanisms of NMDAR permeation and block.

The Hydrophobic Effect Contributes to the Closed State of a Simplified Ion Channel through a Conserved Hydrophobic Patch at the Pore-Helix Crossing.

Ion selectivity-filter structures are strikingly similar throughout the large family of K(++) channels and other p-loop-like receptors (i.e., glutamate receptors). At the same time, the triggers for opening these channels, or gating, are diverse. Two questions that remain unanswered regarding these channels are: (1) what force(s) stabilize the closed non-conducting channel-pore conformation? And (2) what is the free energy associated with transitioning from a closed (non-conducting) to an open (conducting) channel-pore conformation? The effects of charge and hydrophobicity on the conformational states of a model tetrameric biological ion channel are shown utilizing the amino acid sequence from the K(+) channel KcsA as the model "channel". Its widely conserved hydrophobic bundle crossing located adjacent to the lipid head-groups at the intracellular side of the membrane was calculated to have a 5 kcal/mol free energy difference between modeled open and closed conformations. Simulated mutants of amino acids within the hydrophobic region significantly contribute to the size of this difference. Specifically for KcsA, these residues are part of the pH sensor important for channel gating and our results are in agreement with published electrophysiology data. Our simulations support the idea that the hydrophobic effect contributes significantly to the stability of the closed conformation in tetrameric ion channels with a hydrophobic bundle crossing positioned in proximity to the lipid head groups of the biological membrane.

Crosslinking constraints and computational models as complementary tools in modeling the extracellular domain of the glycine receptor.

The glycine receptor (GlyR), a member of the pentameric ligand-gated ion channel superfamily, is the major inhibitory neurotransmitter-gated receptor in the spinal cord and brainstem. In these receptors, the extracellular domain binds agonists, antagonists and various other modulatory ligands that act allosterically to modulate receptor function. The structures of homologous receptors and binding proteins provide templates for modeling of the ligand-binding domain of GlyR, but limitations in sequence homology and structure resolution impact on modeling studies. The determination of distance constraints via chemical crosslinking studies coupled with mass spectrometry can provide additional structural information to aid in model refinement, however it is critical to be able to distinguish between intra- and inter-subunit constraints. In this report we model the structure of GlyBP, a structural and functional homolog of the extracellular domain of human homomeric α1 GlyR. We then show that intra- and intersubunit Lys-Lys crosslinks in trypsinized samples of purified monomeric and oligomeric protein bands from SDS-polyacrylamide gels may be identified and differentiated by MALDI-TOF MS studies of limited resolution. Thus, broadly available MS platforms are capable of providing distance constraints that may be utilized in characterizing large complexes that may be less amenable to NMR and crystallographic studies. Systematic studies of state-dependent chemical crosslinking and mass spectrometric identification of crosslinked sites has the potential to complement computational modeling efforts by providing constraints that can validate and refine allosteric models.

Characterizing the energetic states of the GluR2 ligand binding domain core-dimer.

Tetrameric ligand binding domains of the family of ionotropic glutamate receptors assemble as dimers-of-dimers. Crystallographic studies of several glutamate receptor subtype isolated core-dimers suggest a single stable dimeric conformation. A binding domain dimer has not been captured in other conformations without the aid of biochemical methods to disrupt a critical dimer interface. Molecular dynamics simulations and continuum electrostatics calculations reveal that the active glutamate bound form of the ligand-binding domain found in typical crystal structures is the preferred energetic state of the isolated core-dimer in the presence of agonist glutamate. A desensitized conformational state is a higher energy ligand-bound state of the core-dimer. The resting apo conformational state is comparatively the least energetically favored conformation and does not contain a single state but a set of energetically equivalent conformational core-dimer states. We hypothesize the energetic balance of an open versus closed transmembrane region must be included to characterize the absolute energetic states of the full receptor, which in the presence of the ligand is believed to be a desensitized state.

Energetics of the cleft closing transition and the role of electrostatic interactions in conformational rearrangements of the glutamate receptor ligand binding domain.

The ionotropic glutamate receptors are localized in the pre- and postsynaptic membrane of neurons in the brain. Activation by the principal excitatory neurotransmitter glutamate allows the ligand binding domain to change conformation, communicating opening of the channel for ion conduction. The free energy of the GluR2 S1S2 ligand binding domain (S1S2) closure transition was computed using a combination of thermodynamic integration and umbrella sampling modeling methods. A path that involves lowering the charge on E705 was chosen to clarify the role of this binding site residue. A continuum electrostatics approach in S1S2 is used to show E705, located in the ligand binding cleft, stabilizes the closed conformation of S1S2 via direct interactions with other protein residues, not through the ligand. In the closed conformation, in the absence of a ligand, S1S2 is somewhat more closed than what has been reported in X-ray structures. A semiopen conformation has been identified which is characterized by disruption of a single cross-cleft interaction and differs only slightly in energy from the fully closed S1S2. The fully open S1S2 conformation exhibits a wide energy well and shares structural similarity with the apo S1S2 crystal structure. Hybrid continuum electrostatics/MD calculations along the chosen closure transition pathway reveal solvation energies, and electrostatic interaction energies between two lobes of the protein increase the relative energetic difference between the open and closed conformational states. By analyzing the role of several cross-cleft contacts as well as other binding site residues, we demonstrate how S1S2 interactions facilitate formation of the closed conformation of the GluR2 ligand binding domain.

Ligand-dependent linkage of the ATP site to inhibition gate closure in the KATP channel.

Major advances have been made on the inhibition gate and ATP site of the K(ir)6.2 subunit of the K(ATP) channel, but little is known about conformational coupling between the two. ATP site mutations dramatically disrupt ATP-dependent gating without effect on ligand-independent gating, observed as interconversions between active burst and inactive interburst conformations in the absence of ATP. This suggests that linkage between site and gate is conditionally dependent on ATP occupancy. We studied all substitutions at position 334 of the ATP site in K(ir)6.2deltaC26 that express in Xenopus oocytes. All substitutions disrupted ATP-dependent gating by 10-fold or more. Only positive-charged arginine or lysine at 334, however, slowed ligand-independent gating from the burst, and this was in some but not all patches. Moreover, the polycationic peptide protamine reversed the slowed gating from the burst of 334R mutant channels, and speeded the slow gating from the burst of wild-type SUR1/K(ir)6.2 in the absence of ATP. Our results support a two-step ligand-dependent linkage mechanism for K(ir)6.2 channels in which ATP-occupied sites function to electrostatically dissociate COOH-terminal domains from the membrane, then as in all K(ir) channels, free COOH-terminal domains and inner M2 helices transit to a lower energy state for gate closure.

Anesthetic interaction with ketosteroid isomerase: insights from molecular dynamics simulations.

The nature and the sites of interactions between anesthetic halothane and homodimeric Delta5-3-ketosteroid isomerase (KSI) are characterized by flexible ligand docking and confirmed by 1H-15N NMR. The dynamics consequence of halothane interaction and the implication of the dynamic changes to KSI function are studied by multiple 5-ns molecular dynamics simulations in the presence and absence of halothane. Both docking and MD simulations show that halothane prefer the amphiphilic dimeric interface to the hydrophobic active site of KSI. Halothane occupancy at the dimer interface disrupted the intersubunit hydrogen bonding formed either directly through side chains of polar residues or indirectly through the mediation of the interfacial water molecules. Moreover, in the presence of halothane, the exchange rate of the bound waters with bulk water was increased. Halothane perturbation to the dimer interface affected the overall flexibility of the active site. This action is likely to contribute to the halothane-induced reduction of the KSI activity. The allosteric halothane modulation of the dynamics-function relationship of KSI without direct competition at the enzymatic active sites may be generalized to offer a unifying explanation of anesthetic action on a diverse range of multidomain neuronal proteins that are potentially relevant to clinical general anesthesia.