Folding Simulations of a Nuclear Receptor Box-Containing Peptide Demonstrate the Structural Persistence of the LxxLL Motif Even in the Absence of Its Cognate Receptor

A deep intronic splice-altering AIRE variant causes APECED syndrome through antisense oligonucleotide-targetable pseudoexon inclusion

Autoimmune polyendocrinopathy-candidiasis-ectodermal dystrophy (APECED) is a life-threatening monogenic autoimmune disorder primarily caused by biallelic deleterious variants in the autoimmune regulator (AIRE) gene. We prospectively evaluated 104 patients with clinically diagnosed APECED syndrome and identified 17 patients (16%) from 14 kindreds lacking biallelic AIRE variants in exons or flanking intronic regions; 15 had Puerto Rican ancestry. Through whole-genome sequencing, we identified a deep intronic AIRE variant (c.1504-818 G>A) cosegregating with the disease in all 17 patients. We developed a culture system of AIRE-expressing primary patient monocyte-derived dendritic cells and demonstrated that c.1504-818 G>A creates a cryptic splice site and activates inclusion of a 109-base pair frame-shifting pseudoexon. We also found low-level AIRE expression in patient-derived lymphoblastoid cell lines (LCLs) and confirmed pseudoexon inclusion in independent extrathymic AIRE-expressing cell lines. Through protein modeling and transcriptomic analyses of AIRE-transfected human embryonic kidney 293 and thymic epithelial cell 4D6 cells, we showed that this variant alters the carboxyl terminus of the AIRE protein, abrogating its function. Last, we developed an antisense oligonucleotide (ASO) that reversed pseudoexon inclusion and restored the normal AIRE transcript sequence in LCLs. Thus, our findings revealed c.1504-818 G>A as a founder APECED-causing AIRE variant in the Puerto Rican population and uncovered pseudoexon inclusion as an ASO-reversible genetic mechanism underlying APECED.

The Curious Case of A31P, a Topology-Switching Mutant of the Repressor of Primer Protein: A Molecular Dynamics Study of Its Folding and Misfolding

The effect of mutations on protein structures is usually rather localized and minor. Finding a mutation that can single-handedly change the fold and/or topology of a protein structure is a rare exception. The A31P mutant of the homodimeric Repressor of primer (Rop) protein is one such exception: This single mutation ─and as demonstrated by two independent crystal structure determinations─ can convert the canonical (left-handed/all-antiparallel) 4-α-helical bundle of Rop to a new form (right-handed/mixed parallel and antiparallel bundle) displaying a previously unobserved "bisecting U" topology. The main problem with understanding the dramatic effect of this mutation on the folding of Rop is to understand its very existence: Most computational methods appear to agree that the mutation should have had no appreciable effect, with the majority of energy minimization methods and protein structure prediction protocols indicating that this mutation is fully consistent with the native Rop structure, requiring only a local and minor change at the mutation site. Here we use two long (10 μs each) molecular dynamics simulations to compare the stability and dynamics of the native Rop versus a hypothetical structure that is identical with the native Rop but is carrying this single Alanine31 to Proline mutation. Comparative analysis of the two trajectories convincingly shows that, in contrast to the indications from energy minimization ─but in agreement with the experimental data─, this hypothetical native-like A31P structure is unstable, with its turn regions almost completely unfolding, even under the relatively mild 320 K NpT simulations that we have used for this study. We discuss the implication of these findings for the folding of the A31P mutant, especially with respect to the proposed model of a double-funneled energy landscape.

A proof-of-concept study of the secondary structure of influenza A, B M2 and MERS- and SARS-CoV E transmembrane peptides using folding molecular dynamics simulations in a membrane mimetic solvent

Here, we have carried out a proof-of-concept molecular dynamics (MD) simulation with adaptive tempering in a membrane mimetic environment to study the folding of single-pass membrane peptides. We tested the influenza A M2 viroporin, influenza B M2 viroporin, and protein E from coronaviruses MERS-Cov-2 and SARS-CoV-2 peptides with known experimental secondary structures in membrane bilayers. The two influenza-derived peptides are significantly different in the peptide sequence and secondary structure and more polar than the two coronavirus-derived peptides. Through a total of more than 50 μs of simulation time that could be accomplished in trifluoroethanol (TFE), as a membrane model, we characterized comparatively the folding behavior, helical stability, and helical propensity of these transmembrane peptides that match perfectly their experimental secondary structures, and we identified common motifs that reflect their quaternary organization and known (or not) biochemical function. We showed that BM2 is organized into two structurally distinct parts: a significantly more stable N-terminal half, and a fast-converting C-terminal half that continuously folds and unfolds between α-helical structures and non-canonical structures, which are mostly turns. In AM2, both the N-terminal half and C-terminal half are very flexible. In contrast, the two coronavirus-derived transmembrane peptides are much more stable and fast helix-formers when compared with the influenza ones. In particular, the SARS-derived peptide E appears to be the fastest and most stable helix-former of all the four viral peptides studied, with a helical structure that persists almost without disruption for the whole of its 10 μs simulation. By comparing the results with experimental observations, we benchmarked TFE in studying the conformation of membrane and hydrophobic peptides. This work provided accurate results suggesting a methodology to run long MD simulations and predict structural properties of biologically important membrane peptides.

Folding molecular dynamics simulation of T-peptide, a HIV viral entry inhibitor: Structure, dynamics, and comparison with the experimental data

Peptide T is a synthetic octapeptide fragment, which corresponds to the region 185-192 of the gp120 HIV coat protein and functions as a viral entry inhibitor. In this work, a folding molecular dynamics simulation of peptide T in a membrane-mimicking (DMSO) solution was performed with the aim of characterizing the peptide's structural and dynamical properties. We show that peptide T is highly flexible and dynamic. The main structural characteristics observed were rapidly interconverting short helical stretches and turns, with a notable preference for the formation of β-turns. The simulation also indicated that the C-terminal part appears to be more stable than the rest of the peptide, with the most preferred conformation for residues 5-8 being a β-turn. In order to validate the accuracy of the simulations, we compared our results with the experimental NMR data obtained for the T-peptide in the same solvent. In agreement with the simulation, the NMR data indicated the presence of a preferred structure in solution that was consistent with a β-turn comprising the four C-terminal residues. An additional comparison between the experimental and simulation-derived chemical shifts also showed a reasonable agreement between experiment and simulation, further validating the simulation-derived structural characterization of the T-peptide. We conclude that peptide folding simulations produce physically relevant results even when performed with organic solvents that were not part of the force field parameterization procedure.

The balance between side-chain and backbone-driven association in folding of the α-helical influenza A transmembrane peptide

The correct balance between attractive, repulsive and peptide hydrogen bonding interactions must be attained for proteins to fold correctly. To investigate these important contributors, we sought a comparison of the folding between two 25-residues peptides, the influenza A M2 protein transmembrane domain (M2TM) and the 25-Ala (Ala₂₅ ). M2TM forms a stable α-helix as is shown by circular dichroism (CD) experiments. Molecular dynamics (MD) simulations with adaptive tempering show that M2TM monomer is more dynamic in nature and quickly interconverts between an ensemble of various α-helical structures, and less frequently turns and coils, compared to one α-helix for Ala₂₅ . DFT calculations suggest that folding from the extended structure to the α-helical structure is favored for M2TM compared with Ala₂₅ . This is due to CH⋯O attractive interactions which favor folding to the M2TM α-helix, and cannot be described accurately with a force field. Using natural bond orbital (NBO) analysis and quantum theory atoms in molecules (QTAIM) calculations, 26 CH⋯O interactions and 22 NH⋯O hydrogen bonds are calculated for M2TM. The calculations show that CH⋯O hydrogen bonds, although individually weaker, have a cumulative effect that cannot be ignored and may contribute as much as half of the total hydrogen bonding energy, when compared to NH⋯O, to the stabilization of the α-helix in M2TM. Further, a strengthening of NH⋯O hydrogen bonding interactions is calculated for M2TM compared to Ala₂₅ . Additionally, these weak CH⋯O interactions can dissociate and associate easily leading to the ensemble of folded structures for M2TM observed in folding MD simulations.

CLPTM1L induces estrogen receptor β signaling-mediated radioresistance in non-small cell lung cancer cells

INTRODUCTION

Radioresistance is a major challenge in lung cancer radiotherapy, and new radiosensitizers are urgently needed. Estrogen receptor β (ERβ) is involved in the progression of non-small cell lung cancer (NSCLC), however, the role of ERβ in the response to radiotherapy in lung cancer remains elusive. In the present study, we investigated the mechanism underlying ERβ-mediated transcriptional activation and radioresistance of NSCLC cells.

METHODS

Quantitative real-time PCR, western blot and immunohistochemistry were used to detect the expression of CLPTM1L, ERβ and other target genes. The mechanism of CLPTM1L in modulation of radiosensitivity was investigated by chromatin immunoprecipitation assay, luciferase reporter gene assay, immunofluorescence staining, confocal microscopy, coimmunoprecipitation and GST pull-down assays. The functional role of CLPTM1L was detected by function assays in vitro and in vivo.

RESULTS

CLPTM1L expression was negatively correlated with the radiosensitivity of NSCLC cell lines, and irradiation upregulated CLPTM1L in radioresistant (A549) but not in radiosensitive (H460) NSCLC cells. Meanwhile, IR induced the translocation of CLPTM1L from the cytoplasm into the nucleus in NSCLC cells. Moreover, CLPTM1L induced radioresistance in NSCLC cells. iTRAQ-based analysis and cDNA microarray identified irradiation-related genes commonly targeted by CLPTM1L and ERβ, and CLPTM1L upregulated ERβ-induced genes CDC25A, c-Jun, and BCL2. Mechanistically, CLPTM1L coactivated ERβ by directly interacting with ERβ through the LXXLL NR (nuclear receptor)-binding motif. Functionally, ERβ silencing was sufficient to block CLPTM1L-enhanced radioresistance of NSCLC cells in vitro. CLPTM1L shRNA treatment in combination with irradiation significantly inhibited cancer cell growth in NSCLC xenograft tumors in vivo.

CONCLUSIONS

The present results indicate that CLPTM1L acts as a critical coactivator of ERβ to promote the transcription of its target genes and induce radioresistance of NSCLC cells, suggesting a new target for radiosensitization in NSCLC therapy. Video Abstract.

Folding Molecular Dynamics Simulation of a gp41-Derived Peptide Reconcile Divergent Structure Determinations

T-20 peptide is the first FDA-approved fusion inhibitor against AIDS/HIV-1 gp41 protein. Part of it, the gp41[659-671] peptide, that contains the complete epitope for the neutralizing 2F5 monoclonal antibody, has been found experimentally in a number of divergent structures. Herein, we attempt to reconcile them by using unbiased large-scale all-atom molecular dynamics folding simulations. We show that our approach can successfully capture the peptide's heterogeneity and reach each and every experimentally determined conformation in sub-angstrom accuracy, whilst preserving the peptide's disordered nature. Our analysis also unveils that the minor refinements within the AMBER99SB family of force fields can lead to appreciable differences in the predicted conformational stability arising from subtle differences in the helical- and β-region of the Ramachandran plot. Our work underlines the contribution of molecular dynamics simulation in structurally characterizing pharmacologically important peptides of ambiguous structure.

ETS-1 induces Sorafenib-resistance in hepatocellular carcinoma cells via regulating transcription factor activity of PXR

Transcription factor E26 transformation specific sequence 1 (ETS-1) is a primary regulator in the metastasis of human cancer cells, especially hepatocellular carcinoma (HCC) cells; and it would affect the prognosis of HCC patients who received chemotherapies. However, the regulatory role of ETS-1 in the resistance of HCC cells to molecular-targeting agent remains poorly understood. In the present work, we demonstrate that high ETS-1 expression correlates with poor prognosis of advanced HCC patients received Sorafenib treatment. Mechanistically, ETS-1 binds to nuclear Pregnane X receptor (PXR) directly and enhances PXR's transcription factor activity, which further leads to the induction of the PXR's downstream multi-drug resistance related genes. Overexpression of ETS-1 accelerates the metabolic clearance of Sorafenib in HCC cells and leads to the better survival and faster migration of those cells. The therapeutic studies show that ETS-1 promotes the Sorafenib-resistance of HCC tumor models and ETS-1 blockade enhances the anti-tumor capacity of Sorafenib by decreasing PXR activation. Thus, our study suggests that ETS-1 could enhance the activation of PXR and be a potential therapeutic target for overcoming Sorafenib resistance in HCC treatment.