Curvature Footprints of Transmembrane Proteins in Simulations with the Martini Force Field.

Membranes play essential roles in biological systems and are tremendously diverse in the topologies and chemical and elastic properties that define their functions. In many cases, a given membrane may display considerable heterogeneity, with localized clusters of lipids and proteins exhibiting distinct characteristics compared to adjoining regions. These lipid-protein assemblies can span nanometers to micrometers and are associated with cellular processes such as transport and signaling. While lipid-protein assemblages are dynamic, they can be stabilized by coupling between local membrane composition and shape. Due to the inherent difficulty in resolving atomistic details of membrane proteins in their native lipid environments, these complexes are notoriously challenging to study experimentally; however, molecular dynamics (MD) simulations might be a viable alternative. Here, we aim to assess the utility of coarse-grained (CG) MD simulations with the Martini force field for studying membrane curvature induced by transmembrane (TM) proteins that are reported to generate local curvature. The direction and magnitude of curvature induced by five different TM proteins, as well as certain lipid-protein and protein-protein interactions, were found to be in good agreement with available reference data.

The C-terminal mutation beyond the catalytic site of brown spider phospholipase D significantly impacts its biological activities.

Loxosceles spider envenomation results in dermonecrosis, principally due to phospholipases D (PLDs) present in the venom. These enzymes have a strongly conserved sequence, 273ATXXDNPW280, in the C-terminal region (SMD-tail) that make contact with ß-sheets of the TIM barrel, in which the amino acids Asp277 and Trp280 establish the energetically strongest contacts. The SMD-tail is conserved in PLDs from different species but absent in the non-toxic PLD ancestral glycerophosphodiester phosphodiesterases (GDPDs). This work aims to understand the role of the C-terminal region in the structural stability and/or function of phospholipases D. Through site-directed mutagenesis of the rLiD1 protein (recombinant Loxosceles intermedia dermonecrotic protein 1), we produced two mutants: rLiD1D277A and rLiD1W280A (both with sphingomyelinase activity), in which Asp277 and Trp280 were replaced by alanine. rLiD1D277A showed similar sphingomyelinase activity but at least 2 times more dermonecrotic activity than rLiD1 (wild-type protein). Conversely, while the rLiD1W280A displayed a slight increase in sphingomyelinase activity, its biological activity was similar or lower compared to rLiD1, potentially due to its decreased thermostability and formation of amyloid aggregates. In conclusion, these new findings provide evidence that SMD-tail mutants impact the structure and function of these proteins and point out that residues outside the active site can even increase the function of these enzymes.

Water Leakage Pathway Leads to Internal Hydration of the p53 Core Domain.

The gene encoding the p53 tumor suppressor protein is the most frequently mutated oncogene in cancer patients; yet, generalized strategies for rescuing the function of different p53 mutants remain elusive. This work investigates factors that may contribute to the low inherent stability of the wild-type p53 core domain (p53C) and structurally compromised Y220C mutant. Pressure-induced unfolding of p53C was compared to p63C, the p53 family member with the highest stability, the engineered superstable p53C hexamutant (p53C HM), and lower stability p53C Y220C cancer-associated mutant. The following pressure unfolding values (P50% bar) were obtained: p53C 3346, p53C Y220C 2217, p53C HM 3943, and p63C 4326. Molecular dynamics (MD) simulations revealed that p53C Y220C was most prone to water infiltration, followed by p53C, whereas the interiors of p53C HM and p63C remained comparably dry. A strong correlation (r2 = 0.92) between P50% and extent of interior hydration was observed. The pathways of individual water molecule entry and exit were mapped and analyzed, revealing a common route preserved across the p53 family involving a previously reported pocket, along with a novel surface cleft, both of which appear to be targetable by small molecules. Potential determinants of propensity to water incursion were assessed, including backbone hydrogen bond protection and combined sequence and structure similarity. Collectively, our results indicate that p53C has an intrinsic susceptibility to water leakage, which is exacerbated in a structural class mutant, suggesting that there may be a common avenue for rescuing p53 function.

Do Go Chasing Waterfalls: Enoyl Reductase (FabI) in Complex with Inhibitors Stabilizes the Tetrameric Structure and Opens Water Channels.

The enzyme enoyl-ACP reductase (FabI) is the limiting step of the membrane's fatty acid biosynthesis in bacteria and a druggable target for novel antibacterial agents. The FabI active form is a homotetramer, which displays the highest affinity to inhibitors. Herein, molecular dynamics studies were carried out using the structure of FabI in complex with known inhibitors to investigate their effects on tetramerization. Our results suggest that multimerization is essential for the integrity of the catalytic site and that inhibitor binding enables the multimerization by stabilizing the substrate binding loop (SBL, L:195-200) coupled with changes in the H4/5 (QR interface). We also observed that AFN-1252 (naphtpyridinone derivative) promotes unique conformational changes affecting monomer-monomer interfaces. These changes are induced by AFN-1252 interaction with key residues in the binding sites (Ala95, Tyr146, and Tyr156). In addition, the analysis of water trajectories indicated that AFN-1252 complexes allow more water molecules to enter the binding site than triclosan and MUT056399 complexes. FabI-AFN-1252 simulations show accumulation of water molecules near the Tyr146/147 pocket, which can become a hotspot to the design of novel FabI inhibitors.

Curvature-based sorting of eight lipid types in asymmetric buckled plasma membrane models.

Curvature is a fundamental property of biological membranes and has essential roles in cellular function. Bending of membranes can be induced by their lipid and protein compositions, as well as peripheral proteins, such as those that make up the cytoskeleton. An important aspect of membrane function is the grouping of lipid species into microdomains, or rafts, which serve as platforms for specific biochemical processes. The fluid mosaic model of membranes has evolved to recognize the importance of curvature and leaflet asymmetry, and there are efforts toward evaluating their functional roles. This work investigates the effect of curvature on the sorting of lipids in buckled asymmetric bilayers containing eight lipid types, approximating an average mammalian plasma membrane, through coarse-grained (CG) molecular dynamics (MD) simulations with the Martini force field. The simulations reveal that 1) leaflet compositional asymmetry can induce curvature asymmetry, 2) lipids are sorted by curvature to different extents, and 3) curvature-based partitioning trends show moderate to strong correlations with lipid molecular volumes and head to tail bead ratios, respectively. The findings provide unique insights into the role of curvature in membrane organization, and the curvature-based sorting trends should be useful references for later investigations and potentially interpreting the functional roles of specific lipids.

Lipid distributions and transleaflet cholesterol migration near heterogeneous surfaces in asymmetric bilayers.

Specific and nonspecific protein-lipid interactions in cell membranes have important roles in an abundance of biological functions. We have used coarse-grained (CG) molecular dynamics (MD) simulations to assess lipid distributions and cholesterol flipping dynamics around surfaces in a model asymmetric plasma membrane containing one of six structurally distinct entities: aquaporin-1 (AQP1), the bacterial ß-barrel outer membrane proteins OmpF and OmpX, the KcsA potassium channel, the WALP23 peptide and a carbon nanotube (CNT). Our findings revealed varied lipid partitioning and cholesterol flipping times around the different solutes and putative cholesterol binding sites in AQP1 and KcsA. The results suggest that protein-lipid interactions can be highly variable, and that surface-dependent lipid profiles are effectively manifested in CG simulations with the Martini force field.

Anomalous structural dynamics of minimally frustrated residues in cardiac troponin C triggers hypertrophic cardiomyopathy.

Cardiac TnC (cTnC) is highly conserved among mammals, and genetic variants can result in disease by perturbing Ca2+-regulation of myocardial contraction. Here, we report the molecular basis of a human mutation in cTnC's αD-helix (TNNC1-p.C84Y) that impacts conformational dynamics of the D/E central-linker and sampling of discrete states in the N-domain, favoring the "primed" state associated with Ca2+ binding. We demonstrate cTnC's αD-helix normally functions as a central hub that controls minimally frustrated interactions, maintaining evolutionarily conserved rigidity of the N-domain. αD-helix perturbation remotely alters conformational dynamics of the N-domain, compromising its structural rigidity. Transgenic mice carrying this cTnC mutation exhibit altered dynamics of sarcomere function and hypertrophic cardiomyopathy. Together, our data suggest that disruption of evolutionary conserved molecular frustration networks by a myofilament protein mutation may ultimately compromise contractile performance and trigger hypertrophic cardiomyopathy.

PyAutoFEP: An Automated Free Energy Perturbation Workflow for GROMACS Integrating Enhanced Sampling Methods.

Free energy perturbation (FEP) calculations are now routinely used in drug discovery to estimate the relative FEB (RFEB) of small molecules to a biomolecular target of interest. Using enhanced sampling can improve the correlation between predictions and experimental data, especially in systems with conformational changes. Due to the large number of perturbations required in drug discovery campaigns, the manual setup of FEP calculations is no longer viable. Here, we introduce PyAutoFEP, a flexible and open-source tool to aid the setup of RFEB FEP. PyAutoFEP is written in Python3, and automates the generation of perturbation maps, dual topologies, system building and molecular dynamics (MD), and analysis. PyAutoFEP supports multiple force fields, incorporates replica exchange with solute tempering (REST) and replica exchange with solute scaling (REST2) enhanced sampling methods, and allows flexible λ values along perturbation windows. To validate PyAutoFEP, it was applied to a set of 14 Farnesoid X receptor ligands, a system included in the drug design data resource grand challenge 2. An 88% mean correct sign prediction was achieved, and 75% of the predictions had an error below 1.5 kcal/mol. Results using Amber03/GAFF, CHARMM36m/CGenFF, and OPLS-AA/M/LigParGen had Pearson's r values of 0.71 ± 0.13, 0.30 ± 0.27, and 0.66 ± 0.20, respectively. The Amber03/GAFF and OPLS-AA/M/LigParGen results were on par with the top grand challenge 2 submissions. Applying REST2 improved the results using CHARMM36m/CGenFF (Pearson's r = 0.43 ± 0.21) but had little impact on the other force fields. CHARMM36-YF and CHARMM36-WYF modifications did not yield improved predictions compared to CHARMM36m. Finally, we estimated the probability of finding a molecule 1 pKi better than a lead when using PyAutoFEP to screen 10 or 100 analogues. The probabilities, when compared to random sampling, increased up to sevenfold when 100 molecules were to be screened, suggesting that PyAutoFEP would likely be useful for lead optimization. PyAutoFEP is available on GitHub at https://github.com/lmmpf/PyAutoFEP.

Conformational dynamics of Tetracenomycin aromatase/cyclase regulate polyketide binding and enzyme aggregation propensity.

BACKGROUND: The N-terminal domain of Tetracenomycin aromatase/cyclase (TcmN), an enzyme derived from Streptomyces glaucescens, is involved in polyketide cyclization, aromatization, and folding. Polyketides are a diverse class of secondary metabolites produced by certain groups of bacteria, fungi, and plants with various pharmaceutical applications. Examples include antibiotics, such as tetracycline, and anticancer drugs, such as doxorubicin. Because TcmN is a promising enzyme for in vitro production of polyketides, it is important to identify conditions that enhance its thermal resistance and optimize its function. METHODS: TcmN unfolding, stability, and dynamics were evaluated by fluorescence spectroscopy, circular dichroism, nuclear magnetic resonance 15N relaxation experiments, and microsecond molecular dynamics (MD) simulations. RESULTS: TcmN thermal resistance was enhanced at low protein and high salt concentrations, was pH-dependent, and denaturation was irreversible. Conformational dynamics on the µs-ms timescale were detected for residues in the substrate-binding cavity, and two predominant conformers representing opened and closed cavity states were observed in the MD simulations. CONCLUSION: Based on the results, a mechanism was proposed in which the thermodynamics and kinetics of the TcmN conformational equilibrium modulate enzyme function by favoring ligand binding and avoiding aggregation. GENERAL SIGNIFICANCE: Understanding the principles underlying TcmN stability and dynamics may help in designing mutants with optimal properties for biotechnological applications.

Sequence similarity in 3D for comparison of protein families.

Homologous proteins are often compared by pairwise sequence alignment, and structure superposition if the atomic coordinates are available. Unification of sequence and structure data is an important task in structural biology. Here, we present the Sequence Similarity 3D (SS3D) method of integrating sequence and structure information. SS3D is a distance and substitution matrix-based method for straightforward visualization of regions of similarity and difference between homologous proteins. This work details the SS3D approach, and demonstrates its utility through case studies comparing members of several protein families. The examples show that SS3D can effectively highlight biologically important regions of similarity and dissimilarity. We anticipate that the method will be useful for numerous structural biology applications, including, but not limited to, studies of binding specificity, structure-function relationships, and evolutionary pathways. SS3D is available with a manual and tutorial at https://github.com/0x462e41/SS3D/.

The Status of p53 Oligomeric and Aggregation States in Cancer.

Despite being referred to as the guardian of the genome, when impacted by mutations, p53 can lose its protective functions and become a renegade. The malignant transformation of p53 occurs on multiple levels, such as altered DNA binding properties, acquisition of novel cellular partners, or associating into different oligomeric states. The consequences of these transformations can be catastrophic. Ongoing studies have implicated different oligomeric p53 species as having a central role in cancer biology; however, the correlation between p53 oligomerization status and oncogenic activities in cancer progression remains an open conundrum. In this review, we summarize the roles of different p53 oligomeric states in cancer and discuss potential research directions for overcoming aberrant p53 function associated with them. We address how misfolding and prion-like amyloid aggregation of p53 seem to play a crucial role in cancer development. The misfolded and aggregated states of mutant p53 are prospective targets for the development of novel therapeutic strategies against tumoral diseases.

The synthetic oleanane triterpenoid CDDO-Me binds and inhibits pyruvate kinase M2.

BACKGROUND: The M2 isoform of the glycolytic enzyme pyruvate kinase (PKM2) is one of the key components in the Warburg effect, and an important regulator of cancer cell metabolism. Elevated PKM2 expression is a hallmark of numerous tumor types, making it a promising target for cancer therapy. METHODS: Migration of H1299 lung tumor cells treated with synthetic oleanane triterpenoid derivatives CDDO-Me and CDDO-Im was monitored using scratch and transwell assays. Direct binding and inhibition of PKM2 activity by CDDO-Me was demonstrated by pull-down and activity assays. PKM2 localization in the absence and presence of CDDO-Me or CDDO-Im was determined by subcellular fractionation and immunofluorescence microscopy. Involvement of PKM2 in tumor cell migration was assessed using a stable PKM2 knockdown cell line. RESULTS: We demonstrate that migration of H1299 lung tumor cells is inhibited by CDDO-Me and CDDO-Im in scratch and transwell assays. CDDO-Me binds directly and specifically to recombinant PKM2, leading to a reduction of its catalytic activity. PKM2 knockdown cells exhibit significantly lower migration compared to control cells when subjected to glucose and oxygen deprivation, but not under regular conditions. CONCLUSIONS: The results suggest that PKM2 expression in a tumor-like environment contributes to cell migration, and that PKM2 activity can be down regulated by synthetic triterpenoid derivatives.

Biophysical characterization of p53 core domain aggregates.

Aggregation is the cause of numerous protein conformation diseases. A common facet of these maladies is the transition of a protein from its functional native state into higher order forms, such as oligomers and amyloid fibrils. p53 is an essential tumor suppressor that is prone to such conformational transitions, resulting in its compromised ability to avert cancer. This work explores the biophysical properties of early-, mid-, and late-stage p53 core domain (p53C) aggregates. Atomistic and coarse-grained molecular dynamics (MD) simulations suggest that early- and mid-stage p53C aggregates have a polymorphic topology of antiparallel and parallel ß-sheets that localize to the core amyloidogenic sequence. Both topologies involve similar extents of interstrand mainchain hydrogen bonding, while sidechain interactions could play a role in regulating strand orientation. The free energy difference between the antiparallel and parallel states was within statistical uncertainty. Negative stain electron microscopy of mature fibrils shows a wide distribution of fiber widths, indicating that polymorphism may extend to the quaternary structure level. Circular dichroism of the fibrils was indicative of ß-sheet rich structures in atypical conformations. The Raman spectrum of aggregated p53C was consistent with a mixture of arranged ß-sheets and heterogeneous structural elements, which is compatible with the MD findings of an ordered ß-sheet nucleus flanked by disordered structure. Structural polymorphism is a common property of amyloids; however, because certain polymorphs of the same protein can be more harmful than others, going forward it will be pertinent to establish correlations between p53C aggregate structure and pathology.

Aggregation-primed molten globule conformers of the p53 core domain provide potential tools for studying p53C aggregation in cancer.

The functionality of the tumor suppressor p53 is altered in more than 50% of human cancers, and many individuals with cancer exhibit amyloid-like buildups of aggregated p53. An understanding of what triggers the pathogenic amyloid conversion of p53 is required for the further development of cancer therapies. Here, perturbation of the p53 core domain (p53C) with subdenaturing concentrations of guanidine hydrochloride and high hydrostatic pressure revealed native-like molten globule (MG) states, a subset of which were highly prone to amyloidogenic aggregation. We found that MG conformers of p53C, probably representing population-weighted averages of multiple states, have different volumetric properties, as determined by pressure perturbation and size-exclusion chromatography. We also found that they bind the fluorescent dye 4,4'-dianilino-1,1'-binaphthyl-5,5'-disulfonic acid (bis-ANS) and have a native-like tertiary structure that occludes the single Trp residue in p53. Fluorescence experiments revealed conformational changes of the single Trp and Tyr residues before p53 unfolding and the presence of MG conformers, some of which were highly prone to aggregation. p53C exhibited marginal unfolding cooperativity, which could be modulated from unfolding to aggregation pathways with chemical or physical forces. We conclude that trapping amyloid precursor states in solution is a promising approach for understanding p53 aggregation in cancer. Our findings support the use of single-Trp fluorescence as a probe for evaluating p53 stability, effects of mutations, and the efficacy of therapeutics designed to stabilize p53.

Prediction of Binding Energy of Keap1 Interaction Motifs in the Nrf2 Antioxidant Pathway and Design of Potential High-Affinity Peptides.

Nuclear factor erythroid 2-related factor 2 (Nrf2) is a transcription factor and principal regulator of the antioxidant pathway. The Kelch domain of Kelch-like ECH-associated protein 1 (Keap1) binds to motifs in the N-terminal region of Nrf2, promoting its degradation. There is interest in developing ligands that can compete with Nrf2 for binding to Kelch, thereby activating its transcriptional activities and increasing antioxidant levels. Using experimental Δ Gbind values of Kelch-binding motifs determined previously, a revised hydrophobicity-based model was developed for estimating Δ Gbind from amino acid sequence and applied to rank potential uncharacterized Kelch-binding motifs identified from interaction databases and BLAST searches. Model predictions and molecular dynamics (MD) simulations suggested that full-length MAD2A binds Kelch more favorably than a high-affinity 20-mer Nrf2 E78P peptide, but that the motif in isolation is not a particularly strong binder. Endeavoring to develop shorter peptides for activating Nrf2, new designs were created based on the E78P peptide, some of which showed considerable propensity to form binding-competent structures in MD, and were predicted to interact with Kelch more favorably than the E78P peptide. The peptides could be promising new ligands for enhancing the oxidative stress response.

Investigation of the binding mode of a novel cruzain inhibitor by docking, molecular dynamics, ab initio and MM/PBSA calculations.

Chagas disease remains a major health problem in South America, and throughout the world. The two drugs clinically available for its treatment have limited efficacy and cause serious adverse effects. Cruzain is an established therapeutic target of Trypanosoma cruzi, the protozoan that causes Chagas disease. Our group recently identified a competitive cruzain inhibitor (compound 1) with an IC50 = 15 µM that is also more synthetically accessible than the previously reported lead, compound 2. Prior studies, however, did not propose a binding mode for compound 1, hindering understanding of the structure-activity relationship and optimization. Here, the cruzain binding mode of compound 1 was investigated using docking, molecular dynamics (MD) simulations with ab initio derived parameters, ab initio calculations, and MM/PBSA. Two ligand protonation states and four binding poses were evaluated. A careful ligand parameterization method was employed to derive more physically meaningful parameters than those obtained by automated tools. The poses of unprotonated 1 were unstable in MD, showing large conformational changes and diffusing away from the binding site, whereas the protonated form showed higher stability and interaction with negatively charged residues Asp161 and Cys25. MM/PBSA also suggested that these two residues contribute favorably to binding of compound 1. By combining results from MD, ab initio calculations, and MM/PBSA, a binding mode of 1 is proposed. The results also provide insights for further optimization of 1, an interesting lead compound for the development of new cruzain inhibitors.

Targeting the Prion-like Aggregation of Mutant p53 to Combat Cancer.

Prion-like behavior of several amyloidogenic proteins has been demonstrated in recent years. Despite having functional roles in some cases, irregular aggregation can have devastating consequences. The most commonly known amyloid diseases are Alzheimer's, Parkinson's, and Transmissible Spongiform Encephalopathies (TSEs). The pathophysiology of prion-like diseases involves the structural transformation of wild-type (wt) proteins to transmissible forms that can convert healthy proteins, generating aggregates. The mutant form of tumor suppressor protein, p53, has recently been shown to exhibit prion-like properties. Within the context of p53 aggregation and the search for ways to avert it, this review emphasizes discoveries, approaches, and research from our laboratory and others. Although its standard functions are strongly connected to tumor suppression, p53 mutants and aggregates are involved in cancer progression. p53 aggregates are heterogeneous assemblies composed of amorphous aggregates, oligomers, and amyloid-like fibrils. Evidence of these structures in tumor tissues, the in vitro capability for p53 mutants to coaggregate with wt protein, and the detection of cell-to-cell transmission indicate that cancer has the basic characteristics of prion and prion-like diseases. Various approaches aim to restore p53 functions in cancer. Methods include the use of small-molecule and peptide stabilizers of mutant p53, zinc administration, gene therapy, alkylating and DNA intercalators, and blockage of p53-MDM2 interaction. A primary challenge in developing small-molecule inhibitors of p53 aggregation is the large number of p53 mutations. Another issue is the inability to recover p53 function by dissociating mature fibrils. Consequently, efforts have emerged to target the intermediate species of the aggregation reaction. Φ-value analysis has been used to characterize the kinetics of the early phases of p53 aggregation. Our experiments using high hydrostatic pressure (HHP) and chemical denaturants have helped to clarify excited conformers of p53 that are prone to aggregation. Molecular dynamics (MD) and phasor analysis of single Trp fluorescence signals point toward the presence of preamyloidogenic conformations of p53, which are not observed for p63 or p73. Exploring the features of competent preamyloidogenic states of wt and different p53 mutants may provide a framework for designing personalized drugs for the restoration of p53 function. Protection of backbone hydrogen bonds (BHBs) has been shown to be an important factor for the stability of amyloidogenic proteins and was employed to identify and stabilize the structural defect resulting from the p53 Y220C mutation. Using MD simulations, we compared BHB protection factors between p53 family members to determine the donor-acceptor pairs in p53 that exhibit lower protection. The identification of structurally vulnerable sites in p53 should provide new insights into rational designs that can rapidly be screened using our experimental methodology. Through continued and combined efforts, the outlook is positive for the development of strategies for regulating p53 amyloid transformation.

Aggregation tendencies in the p53 family are modulated by backbone hydrogen bonds.

The p53 family of proteins is comprised of p53, p63 and p73. Because the p53 DNA binding domain (DBD) is naturally unstable and possesses an amyloidogenic sequence, it is prone to form amyloid fibrils, causing loss of functions. To develop p53 therapies, it is necessary to understand the molecular basis of p53 instability and aggregation. Light scattering, thioflavin T (ThT) and high hydrostatic pressure (HHP) assays showed that p53 DBD aggregates faster and to a greater extent than p63 and p73 DBDs, and was more susceptible to denaturation. The aggregation tendencies of p53, p63, and p73 DBDs were strongly correlated with their thermal stabilities. Molecular Dynamics (MD) simulations indicated specific regions of structural heterogeneity unique to p53, which may be promoted by elevated incidence of exposed backbone hydrogen bonds (BHBs). The results indicate regions of structural vulnerability in the p53 DBD, suggesting new targetable sites for modulating p53 stability and aggregation, a potential approach to cancer therapy.

Aggregation and Prion-Like Properties of Misfolded Tumor Suppressors: Is Cancer a Prion Disease?

Prion diseases are disorders that share several characteristics that are typical of many neurodegenerative diseases. Recently, several studies have extended the prion concept to pathological aggregation in malignant tumors involving misfolded p53, a tumor-suppressor protein. The aggregation of p53 and its coaggregation with p53 family members, p63 and p73, have been shown. Certain p53 mutants exert a dominant-negative regulatory effect on wild-type (WT) p53. The basis for this dominant-negative effect is that amyloid-like mutant p53 converts WT p53 into an aggregated species, leading to a gain-of-function (GoF) phenotype and the loss of its tumor-suppressor function. Recently, it was shown that p53 aggregates can be internalized by cells and can coaggregate with endogenous p53, corroborating the prion-like properties of p53 aggregates. The prion-like behavior of oncogenic p53 mutants provides an explanation for its dominant-negative and GoF properties, including the high metastatic potential of cancer cells carrying p53 mutations. The inhibition of p53 aggregation appears to represent a promising target for therapeutic intervention in patients with malignant tumors.

Characterization of the Free State Ensemble of the CoRNR Box Motif by Molecular Dynamics Simulations.

Intrinsically disordered proteins (IDPs) and regions are highly prevalent in eukaryotic proteomes, and like folded proteins, they perform essential biological functions. Interaction sites in folded proteins are generally formed by tertiary structures, whereas IDPs use short segments called linear motifs (LMs). Despite their short length and lack of stable structure, LMs may have considerable structural propensities, which often resemble bound-state conformations with targets. Structural data is crucial for understanding the molecular basis of protein interactions and development of targeted pharmaceuticals, but IDPs present considerable challenges to experimental techniques. As a result, IDPs are largely underrepresented in the Protein Data Bank. In the face of experimental challenges, molecular dynamics (MD) simulations have proven to be a useful tool for structural characterization of IDPs. Here, the free state ensemble of the nuclear receptor corepressor 1 (NCOR1) CoRNR box 3 motif, which is important for binding to nuclear receptors to control gene expression, is studied using MD simulations of a total of 8 µs. Transitions between disordered and α-helical conformations resembling a bound-state structure were observed throughout the trajectory, indicating that the motif may have a natural conformational bias toward bound-state structures. The data shows that the disordered and folded populations are separated by a low energy (4-6 kJ/mol) barrier, and the presence of off-pathway intermediates, leading to a C-terminally folded species that cannot efficiently transition into a completely folded conformation. Structural transitions and folding pathways within the free state ensemble were well-described by principal component analysis (PCA) of the peptide backbone dihedral angles, with the analysis providing insight for increasing structural homogeneity of the ensemble.