Ubiquitin: Molecular modeling and simulations

Influence of Aqueous Arginine Solution on Regulating Conformational Stability and Hydration Properties of the Secondary Structural Segments of a Protein at Elevated Temperatures: A Molecular Dynamics Study

The effects of aqueous arginine solution on the conformational stability of the secondary structural segments of a globular protein, ubiquitin, and the structure and dynamics of the surrounding water and arginine were examined by performing atomistic molecular dynamics (MD) simulations. Attempts have been made to identify the osmolytic efficacy of arginine solution, and its influence in guiding the hydration properties of the protein at an elevated temperature of 450 K. The similar properties of the protein in pure water at elevated temperatures were computed and compared. Replica exchange MD simulation was performed to explore the arginine solution's sensitivity in stabilizing the protein conformations for a wide range of temperatures (300-450 K). It was observed that although all the helices and strands of the protein undergo unfolding at elevated temperature in pure water, they exhibited native-like conformational dynamics in the presence of arginine at both ambient and elevated temperatures. We find that the higher free energy barrier between the folded native and unfolded states of the protein primarily arises from the structural transformation of α-helix, relative to the strands. Our study revealed that the water structure around the secondary segments depends on the nature of amino acid compositions of the helices and strands. The reorientation of water dipoles around the helices and strands was found hindered due to the presence of arginine in the solution; such hindrance reduces the possibility of exchange of hydrogen bonds that formed between the secondary segments of protein and water (PW), and as a result, PW hydrogen bonds take longer time to relax than in pure water. On the other hand, the origin of slow relaxation of protein-arginine (PA) hydrogen bonds was identified to be due to the presence of different types of protein-bound arginine molecules, where arginine interacts with the secondary structural segments of the protein through multiple/bifurcated hydrogen bonds. These protein-bound arginine formed different kinds of bridged PA hydrogen bonds between amino acid residues of the same secondary segments or among multiple bonds and helped protein to conserve its native folded form firmly.

A case study of poly (aryl ether sulfone) hemodialysis membrane interactions with human blood: Molecular dynamics simulation and experimental analyses

Patients with end-stage renal diseases (ESRD) require specific health cares as the accumulation of toxins due to the lack of kidney functionality would affect their lives. However, the mortality rate is still high due to cardiovascular diseases, socks, etc. A majority of patients with chronic kidney disease (CKD) require hemodialysis services. Blood purifying membranes, as the main component of hemodialysis setups, however, still suffer from lack of optimum biocompatibility, which results in morbidity and mortality of hemodialysis service receiving patients. The goal of the present case study is to have an in-depth understanding of the current blood-hemodialysis membrane interactions occurring during hemodialysis sessions using poly (aryl ether sulfone)-poly (vinyl pyrrolidone) (PAES-PVP) membrane. Attenuated total reflectance-Fourier transmission infrared (ATR-FTIR) spectroscopy, Raman spectroscopy, and solid-state nuclear magnetic resonance (SSNMR) spectroscopy were used to assess the initial chemical structure of the PAES-PVP membrane along with the variations after with the infections with human blood. Furthermore, scanning electron microscopy (SEM) and Transition electron microscopy (TEM) were used to visualize the structural variation of the membrane, blood aggregations, and blood clots on the membrane surface. Besides, Molecular dynamics (MD) simulation was used to assess the interaction of PAES-PVP with major human blood proteins, in terms of interaction energy, which is a novel contribution to the area. The macromolecules (human serum albumin (HSA), human serum transferrin (TRF), and human fibrinogen (HFG)) were chosen from the plasma protein component. These protein structures were chosen based on their different molecular size. Three advanced spectroscopy techniques and two advanced visualization techniques were used for the assessment of the membranes. Spectroscopy studies revealed amine related peak displacement and intensity shifts as indices for attachment of biological species to the polymeric membrane surfaces. Raman peaks around 370, 798, and 1299 cm^-1, which experienced significant shifts that were related to carbon-nitrogen and sulfur-oxygen bonds due to protein adhesion. Visualization techniques illustrated blood protein fouling patterns and extracellular vesicles' presence in the pore structures into membranes. The findings highlight the importance of whole structure biocompatibility improvement, rather than only focusing on surface modifications of hemodialysis membranes. Molecular dynamics simulation assessment showed various interaction behaviors for different proteins suggesting molecular weight and active residues of the protein macromolecules play an important role in interacting with polymeric structure. FB had the highest interaction (4,274,749.07 kcal/mol) and binding (10,370.90 kcal/mol) energy with the PAES-PVP structure. TRF owned the lowest interaction energy with respect to its lower molecular weight and fewer active residue count.

Machine Learning Driven Analysis of Large Scale Simulations Reveals Conformational Characteristics of Ubiquitin Chains

Understanding the conformational characteristics of protein complexes in solution is crucial for a deeper insight in their biological function. Molecular dynamics simulations performed on high performance computing plants and with modern simulation techniques can be used to obtain large data sets that contain conformational and thermodynamic information about biomolecular systems. While this can in principle give a detailed picture of protein-protein interactions in solution and therefore complement experimental data, it also raises the challenge of processing exceedingly large high-dimensional data sets with several million samples. Here we present a novel method for the characterization of protein-protein interactions, which combines a neural network based dimensionality reduction technique to obtain a two-dimensional representation of the conformational space with a density based clustering algorithm for state detection and a metric which assesses the (dis)similarity between different conformational spaces. This method is highly scalable and therefore makes the analysis of massive data sets computationally tractable. We demonstrate the power of this approach to large scale data analysis by characterizing highly dynamic conformational phase spaces of differently linked ubiquitin (Ub) oligomers from coarse-grained simulations. We are able to extract a protein-protein interaction model for two unlinked Ub proteins which is then used to determine how the Ub-Ub interaction pattern is altered in Ub oligomers by the introduction of a covalent linkage. We find that the Ub chain conformational ensemble depends highly on the linkage type and for some cases also on the Ub chain length. By this, we obtain insight into the conformational characteristics of different Ub chains and how this may contribute to linkage type and chain length specific recognition.

Increased expression of Triad1 is associated with neuronal apoptosis after intracerebral hemorrhage in adult rats

Objective: It has been demonstrated that Triad1 (2 RING fingers and double RING finger linked 1) negatively regulates myeloid cell growth and induces cell apoptosis. However, its functions in intracerebral hemorrhage (ICH) disease have not been conducted. In this study, the role of Triad1 in rat model of ICH was explored.Methods: We observe an increasing expression of Triad1 in areas adjacent to hematoma after ICH. Immunofluorescence shows that Triad1 is colocalized with neurons, while not microglia or astrocyte, indicates its correlation with neuronal activities following ICH.Results: As neuronal apoptosis is the most crucial event in ICH disease, the expression of active caspase-3 and p53 is also enhanced around the hematoma, which is consistent with Triad1 in expression tendency. In turn, Triad1 depletion in primary cortical neurons decreased the apoptosis of neurons after using Triad1 shRNA.Conclusion: We conclude that inhibition of Triad1 expression might protect the brain from secondary damage following ICH.

Regulation of Cancer Immune Checkpoint: Mono- and Poly-Ubiquitination: Tags for Fate

The antagonism, stalemate and compromise between the immune system and tumor cells is closely associated with tumor development and progression. In recent years, tumor immunotherapy has made continuous breakthroughs. It has become an important approach for cancer treatment, improving the survival and prognosis of more and more tumor patients. Further investigating the mechanism of tumor immune regulation, and exploring tumor immunotherapy targets with high specificity and wide applicability will provide researchers and clinicians with favorable weapons towards cancer. Ubiquitination affects protein fate through influencing the activity, stability and location of target protein. The regulation of substrate protein fate by ubiquitination is involved in cell cycle, apoptosis, transcriptional regulation, DNA repair, immune response, protein degradation and quality control. E3 ubiquitin ligase selectively recruits specific protein substrates through specific protein-protein interactions to determine the specificity of the overall ubiquitin modification reaction. Immune-checkpoint inhibitory pathway is an important mechanism for tumor cells to evade immune killing, which can inhibit T cell activity. Blocking the immune checkpoints and activating T cells through targeting the negative regulatory factors of T cell activation and removing the "brake" of T lymphocytes can enhance T cells immune response against tumors. Therefore, blocking the immune checkpoint is one of the methods to enhance the activity of T cells, and it is also a hot target for the development of anti-tumor drugs in recent years, whose inhibitors have shown good effect in specific tumor treatment. Ubiquitination, as one of the most important posttranslational modification of proteins, also modulates the expression, intracellular trafficking, subcellular and membranous location of immune checkpoints, regulating the immune surveillance of T cells to tumors.

Molecular dynamics simulation of proteins under high pressure: Structure, function and thermodynamics

BACKGROUND

Molecular dynamics (MD) simulation is well-recognized as a powerful tool to investigate protein structure, function, and thermodynamics. MD simulation is also used to investigate high pressure effects on proteins. For conducting better MD simulation under high pressure, the main issues to be addressed are: (i) protein force fields and water models were originally developed to reproduce experimental properties obtained at ambient pressure; and (ii) the timescale to observe the pressure effect is often much longer than that of conventional MD simulations.

SCOPE OF REVIEW

First, we describe recent developments in MD simulation methodologies for studying the high-pressure structure and dynamics of protein molecules. These developments include force fields for proteins and water molecules, and enhanced simulation techniques. Then, we summarize recent studies of MD simulations of proteins in water under high pressure.

MAJOR CONCLUSIONS

Recent MD simulations of proteins in solution under pressure have reproduced various phenomena identified by experiments using high pressure, such as hydration, water penetration, conformational change, helix stabilization, and molecular stiffening.

GENERAL SIGNIFICANCE

MD simulations demonstrate differences in the properties of proteins and water molecules between ambient and high-pressure conditions. Comparing the results obtained by MD calculations with those obtained experimentally could reveal the mechanism by which biological molecular machines work well in collaboration with water molecules.

Towards a molecular basis of ubiquitin signaling: A dual-scale simulation study of ubiquitin dimers

Covalent modification of proteins by ubiquitin or ubiquitin chains is one of the most prevalent post-translational modifications in eukaryotes. Different types of ubiquitin chains are assumed to selectively signal respectively modified proteins for different fates. In support of this hypothesis, structural studies have shown that the eight possible ubiquitin dimers adopt different conformations. However, at least in some cases, these structures cannot sufficiently explain the molecular basis of the selective signaling mechanisms. This indicates that the available structures represent only a few distinct conformations within the entire conformational space adopted by a ubiquitin dimer. Here, molecular simulations on different levels of resolution can complement the structural information. We have combined exhaustive coarse grained and atomistic simulations of all eight possible ubiquitin dimers with a suitable dimensionality reduction technique and a new method to characterize protein-protein interfaces and the conformational landscape of protein conjugates. We found that ubiquitin dimers exhibit characteristic linkage type-dependent properties in solution, such as interface stability and the character of contacts between the subunits, which can be directly correlated with experimentally observed linkage-specific properties.

Post-translational modification of proteins by covalent attachment of ubiquitin is a key cellular process, regulating for example the fate and recycling of proteins. We present a new method to combine multiscale simulation with advanced analysis methods to characterize the states of ubiquitin-ubiquitin conjugates. We found that the linkage position affects the conformational space of ubiquitin dimers, determining the number and stability of relevant states, the character of subunit contacts and the nature of the surface exposed to possible binding partners.

Phylogenetic Resolution in Juglans Based on Complete Chloroplast Genomes and Nuclear DNA Sequences

Walnuts (Juglans of the Juglandaceae) are well-known economically important resource plants for the edible nuts, high-quality wood, and medicinal use, with a distribution from tropical to temperate zones and from Asia to Europe and Americas. There are about 21 species in Juglans. Classification of Juglans at section level is problematic, because the phylogenetic position of Juglans cinerea is disputable. Lacking morphological and DNA markers severely inhibited the development of related researches. In this study, the complete chloroplast genomes and two nuclear DNA regions (the internal transcribed spacer and ubiquitin ligase gene) of 10 representative taxa of Juglans were used for comparative genomic analyses in order to deepen the understanding on the application value of genetic information for inferring the phylogenetic relationship of the genus. The Juglans chloroplast genomes possessed the typical quadripartite structure of angiosperms, consisting of a pair of inverted repeat regions separated by a large single-copy region and a small single-copy region. All the 10 chloroplast genomes possessed 112 unique genes arranged in the same order, including 78 protein-coding, 30 tRNA, and 4 rRNA genes. A combined sequence data set from two nuclear DNA regions revealed that Juglans plants could be classified into three branches: (1) section Juglans, (2) section Cardiocaryon including J. cinerea which is closer to J. mandshurica, and (3) section Rhysocaryon. However, three branches with a different phylogenetic topology were recognized in Juglans using the complete chloroplast genome sequences: (1) section Juglans, (2) section Cardiocaryon, and (3) section Rhysocaryon plus J. cinerea. The molecular taxonomy of Juglans is almost compatible to the morphological taxonomy except J. cinerea (section Trachycaryon). Based on the complete chloroplast genome sequence data, the divergence time between section Juglans and section Cardiocaryon was 44.77 Mya, while section Rhysocaryon diverged from other sections in the genus Juglans was 47.61 Mya. Eleven of the 12 small inversions in the chloroplast genomes provided valuable phylogenetic information for classification of walnut plants at section and species levels. Our results are valuable for future studies on Juglans genetic diversity and will enhance the understanding on the phylogenetic evolution of Juglandaceae.

Enhancing backbone sampling in Monte Carlo simulations using internal coordinates normal mode analysis

Normal mode methods are becoming a popular alternative to sample the conformational landscape of proteins. In this study, we describe the implementation of an internal coordinate normal mode analysis method and its application in exploring protein flexibility by using the Monte Carlo method PELE. This new method alternates two different stages, a perturbation of the backbone through the application of torsional normal modes, and a resampling of the side chains. We have evaluated the new approach using two test systems, ubiquitin and c-Src kinase, and the differences to the original ANM method are assessed by comparing both results to reference molecular dynamics simulations. The results suggest that the sampled phase space in the internal coordinate approach is closer to the molecular dynamics phase space than the one coming from a Cartesian coordinate anisotropic network model. In addition, the new method shows a great speedup (∼5-7×), making it a good candidate for future normal mode implementations in Monte Carlo methods.

A new nuclear DNA marker from ubiquitin ligase gene region for genetic diversity detection of walnut germplasm resources

Development of more sensitive nuclear DNA markers for identification of species, particularly closely allied taxa has been a challenging task that has attracted interest from scientists in fields of biotechnological development and genetic diversity detection. In this study, the sequence of the ubiquitin ligase gene (UBE3) region of nuclear DNA was tested for applicability and efficacy in revealing genetic diversity of walnut resources, with an emphasis on inter- and intra-specific levels. Analysis on genetic relationship among the taxa was conducted with the neighbor-joining (NJ) method. The number of variable bases in the UBE3 region was 20 sites. All nine taxa (species/variety/cultivars) were distinguished using the UBE3 sequence. In addition, each taxon was characterized molecularly with a unique nucleotide molecular formula using ten variable base sites derived from the nuclear DNA UBE3 gene sequence. This study presents a good complementary methodology for developing new DNA markers for identification of genus Juglans.