Molecular dynamics simulations of the adsorption of an intrinsically disordered protein: Force field and water model evaluation in comparison with experiments

Streamlining NMR Chemical Shift Predictions for Intrinsically Disordered Proteins: Design of Ensembles with Dimensionality Reduction and Clustering

By merging advanced dimensionality reduction (DR) and clustering algorithm (CA) techniques, our study advances the sampling procedure for predicting NMR chemical shifts (CS) in intrinsically disordered proteins (IDPs), making a significant leap forward in the field of protein analysis/modeling. We enhance NMR CS sampling by generating clustered ensembles that accurately reflect the different properties and phenomena encapsulated by the IDP trajectories. This investigation critically assessed different rapid CS predictors, both neural network (e.g., Sparta+ and ShiftX2) and database-driven (ProCS-15), and highlighted the need for more advanced quantum calculations and the subsequent need for more tractable-sized conformational ensembles. Although neural network CS predictors outperformed ProCS-15 for all atoms, all tools showed poor agreement with HN CSs, and the neural network CS predictors were unable to capture the influence of phosphorylated residues, highly relevant for IDPs. This study also addressed the limitations of using direct clustering with collective variables, such as the widespread implementation of the GROMOS algorithm. Clustered ensembles (CEs) produced by this algorithm showed poor performance with chemical shifts compared to sequential ensembles (SEs) of similar size. Instead, we implement a multiscale DR and CA approach and explore the challenges and limitations of applying these algorithms to obtain more robust and tractable CEs. The novel feature of this investigation is the use of solvent-accessible surface area (SASA) as one of the fingerprints for DR alongside previously investigated α carbon distance/angles or ϕ/ψ dihedral angles. The ensembles produced with SASA tSNE DR produced CEs better aligned with the experimental CS of between 0.17 and 0.36 r2 (0.18-0.26 ppm) depending on the system and replicate. Furthermore, this technique produced CEs with better agreement than traditional SEs in 85.7% of all ensemble sizes. This study investigates the quality of ensembles produced based on different input features, comparing latent spaces produced by linear vs nonlinear DR techniques and a novel integrated silhouette score scanning protocol for tSNE DR.

Exploring the Functional Landscape of the p53 Regulatory Domain: The Stabilizing Role of Post-Translational Modifications

This study focuses on the intrinsically disordered regulatory domain of p53 and the impact of post-translational modifications. Through fully atomistic explicit water molecular dynamics simulations, we show the wealth of information and detailed understanding that can be obtained by varying the number of phosphorylated amino acids and implementing a restriction in the conformational entropy of the N-termini of that intrinsically disordered region. The take-home message for the reader is to achieve a detailed understanding of the impact of phosphorylation with respect to (1) the conformational dynamics and flexibility, (2) structural effects, (3) protein interactivity, and (4) energy landscapes and conformational ensembles. Although our model system is the regulatory domain p53 of the tumor suppressor protein p53, this study contributes to understanding the general effects of intrinsically disordered phosphorylated proteins and the impact of phosphorylated groups, more specifically, how minor changes in the primary sequence can affect the properties mentioned above.

Developing tardigrade-inspired material: Track membranes functionalized with Dsup protein for cell-free DNA isolation

When developing functionalized biomaterials, the proteins from extremophilic organisms, in particular unique tardigrade disordered proteins, are of great value. The damage suppressor protein (Dsup), initially discovered in the tardigrade Ramazzottius varieornatus and found to be an efficient DNA protector under oxidative and irradiation stress, has been hypothesized to possess a good potential for the development of the material, which can isolate cell-free DNA. With this in mind, DNA-nonadsorbing polyethylene terephthalate track membranes have been functionalized using the Dsup protein via covalent bonding with glutaraldehyde. The filtration experiments have verified the ability of track membranes with the immobilized Dsup protein to adsorb cell-free DNA, with an accumulation capacity of 70 ± 19 mg m-2. The resulting track membrane-based biomaterial might be used in various devices for filtration and separation of cell-free DNA molecules from biological solutions and environmental samples, and also for their accumulation, storage, and further manipulation.

Elucidating the Membrane Binding Process of a Disordered Protein: Dynamic Interplay of Anionic Lipids and the Polybasic Region

Intrinsically disordered regions of proteins are responsible for many biological processes such as in the case of liver kinase B1 (LKB1)-a serine/threonine kinase relevant for cell proliferation and cell polarity. LKB1 becomes fully activated upon recruitment to the plasma membrane by binding of its disordered C-terminal polybasic motif consisting of eight lysines/arginines to phospholipids. Here, we present extensive molecular dynamics (MD) simulations of the polybasic motif interacting with a model membrane composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1-palmitoyl-2-oleyl phosphatidic acid (PA) and cell culture experiments. Protein-membrane binding effects are due to the electrostatic interactions between the polybasic amino acids and PAs. For significant binding, the first three lysines turn out to be dispensable, which was also recapitulated in cell culture using transfected GFP-LKB1 variants. LKB1-membrane binding results in nonmonotonous changes in the structure of the protein as well as the membrane, in particular, accumulation of PAs and reduced thickness at the protein-membrane contact area. The protein-lipid binding turns out to be highly dynamic due to an interplay of PA-PA repulsion and protein-PA attraction. The thermodynamics of this interplay is captured by a statistical fluctuation model, which allows the estimation of both energies. Quantification of the significance of each polar amino acid in the polybasic provides detailed insights into the molecular mechanism of protein-membrane binding of LKB1. These results can likely be transferred to other proteins, which interact by intrinsically disordered polybasic regions with anionic membranes.

Exploring the Role of Globular Domain Locations on an Intrinsically Disordered Region of p53: A Molecular Dynamics Investigation

The pre-tetramerization loop (PTL) of the human tumor suppressor protein p53 is an intrinsically disordered region (IDR) necessary for the tetramerization process, and its flexibility contributes to the essential conformational changes needed. Although the IDR can be accurately simulated in the traditional manner of molecular dynamics (MD) with the end-to-end distance (EEdist) unhindered, we sought to explore the effects of restraining the EEdist to the values predicted by electron microscopy (EM) and other distances. Simulating the PTL trajectory with a restrained EEdist , we found an increased agreement of nuclear magnetic resonance (NMR) chemical shifts with experiments. Additionally, we observed a plethora of secondary structures and contacts that only appear when the trajectory is restrained. Our findings expand the understanding of the tetramerization of p53 and provide insight into how mutations could make the protein impotent. In particular, our findings demonstrate the importance of restraining the EEdist in studying IDRs and how their conformations change under different conditions. Our results provide a better understanding of the PTL and the conformational dynamics of IDRs in general, which are useful for further studies regarding mutations and their effects on the activity of p53.

Prediction of interface between regions of varying degrees of order or disorderness in intrinsically disordered proteins from dihedral angles

Intrinsically disordered proteins (IDPs) are proteins that do not form uniquely defined three-dimensional (3-D) structures. Experimental research on IDPs is difficult since they go against the traditional protein structure-function paradigm. Although there are several predictors of disorder based on amino acid sequences, but very limited based on the 3-D structures of proteins. Dihedral angles have a significant role in predicting protein structure because they establish a protein's backbone, which, coupled with its side chain, establishes its overall shape. Here, we have carried out atomistic Molecular Dynamics (MD) simulations on four different proteins: one ordered protein (Monellin), two partially disordered proteins (p53-TAD and Amyloid beta (Aβ1-42) peptide), and one completely disordered protein (Histatin 5). The MD simulation trajectories for the corresponding four proteins were used to conduct dihedral angle (ϕ and ѱ) analysis. Then, the average dihedral angles for each of the residues were calculated and plotted against the residue index. We noticed steep rises or falls in the average ϕ value at certain locations in the plot. These sudden shifts in the average ϕ value reflect the interface between regions of varying degrees of order or disorderness in intrinsically disordered proteins. Using this method, the probable conformer of a protein with a higher degree of disorder can be found among the ensembles of structures sampled during the MD simulations. The results of our study offer new understandings on precisely identifying regions of various degrees of disorder in intrinsically disordered proteins.Communicated by Ramaswamy H. Sarma.

Enhancing the Biopharmacological Characteristics of Asperosaponin VI: Unveiling Dynamic Self-Assembly Phase Transitions in the Gastrointestinal Environment

Purpose

Asperosaponin VI (ASP VI) as an active ingredient of Dipsacus asperoides, which has a wide range of biological and pharmacological activity. However, its development and application are restricted due to the poor gastrointestinal permeability and oral bioavailability. This investigation aims to reveal the influence of the self-assembled structure by the interaction between ASP VI and endogenous components NaTC and/or DOPC in the gastrointestinal environment on its biopharmaceutical properties, and novelty elucidated the molecular mechanism for the formation of self-assembled nanomicelles.

Methods

This change in phase state in gastrointestinal fluids is characterized by dynamic light scattering (DLS) and transmission electron microscope (TEM). UPLC-Q-TOF-MS was used to analyze the composition of phase components and the exposure of nanomicelles in vivo. Molecular dynamics simulation (MDS) was applied to preliminarily elucidate the self-assembly mechanism of ASP VI in the gastrointestinal environment. Furthermore, theS8 promoting absorption mechanism of nanomicelles were investigated through in vivo pharmacokinetic experiments, parallel artificial membrane permeability assay (PAMPA), quadruple single-pass intestinal perfusion in rats, and Caco-2 cell monolayer model.

Results

We demonstrated that the ASP VI could spontaneously form dynamic self-assembled structures with sodium taurocholate (NaTC) and dipalmitoyl phosphatidylcholine (DOPC) during gastrointestinal solubilization, which promoted the gastrointestinal absorption and permeability of ASP VI and increased its exposure in vivo, thus improving the biopharmacological characteristics of ASP VI. Moreover, ASP VI-NaTC-DOPC-self-assembled nanostructures (ASP VI-NaTC-DOPC-SAN) manifested higher cellular uptake in Caco-2 cells as evidenced by flow cytometry and confocal microscopy, and this study also preliminarily revealed the mechanism of self-assembly formation of ASP VI with endogenous components NaTC and DOPC driven by electrostatic and hydrogen bonding interactions.

Conclusion

This study provides evidence that the dynamic self-assembled phase transition may play a key role in improving the biopharmacological characteristics of insoluble or low permeability active ingredients during the gastrointestinal dissolution of Chinese medicines.

Cross-Linked α-Synuclein as Inhibitor of Amyloid Formation

The aggregation and amyloid formation of α-synuclein is associated with Parkinson's disease and other synucleinopathies. In its native, monomeric form α-synuclein is an intrinsically disordered protein represented by highly dynamic conformational ensembles. Inhibition of α-synuclein aggregation using small molecules, peptides, or proteins has been at the center of interest in recent years. Our aim was to explore the effects of cross-linking on the structure and aggregation/amyloid formation properties of α-synuclein. Comparative analysis of available high-resolution amyloid structures and representative structural models and MD trajectory of monomeric α-synuclein revealed that potential cross-links in the monomeric protein are mostly incompatible with the amyloid forms and thus might inhibit fibrillation. Monomeric α-synuclein has been intramolecularly chemically cross-linked under various conditions using different cross-linkers. We determined the location of cross-links and their frequency using mass spectrometry and found that most of them cannot be realized in the amyloid structures. The inhibitory potential of cross-linked proteins has been experimentally investigated using various methods, including thioflavin-T fluorescence and transmission electron microscopy. We found that conformational constraints applied by cross-linking fully blocked α-synuclein amyloid formation. Moreover, DTSSP-cross-linked molecules exhibited an inhibitory effect on the aggregation of unmodified α-synuclein as well.