Generalized-ensemble method study: A helix-mimetic compound inhibits protein-protein interaction by long-range and short-range intermolecular interactions

Dynamics and interactions of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are widespread in eukaryotes and participate in a variety of important cellular processes. Numerous studies using state-of-the-art experimental and theoretical methods have advanced our understanding of IDPs and revealed that disordered regions engage in a large repertoire of intra- and intermolecular interactions through their conformational dynamics, thereby regulating many intracellular functions in concert with folded domains. The mechanisms by which IDPs interact with their partners are diverse, depending on their conformational propensities, and include induced fit, conformational selection, and their mixtures. In addition, IDPs are implicated in many diseases, and progress has been made in designing inhibitors of IDP-mediated interactions. Here we review these recent advances with a focus on the dynamics and interactions of IDPs.

REST Is Not Resting: REST/NRSF in Health and Disease

Chromatin modifications play a crucial role in the regulation of gene expression. The repressor element-1 (RE1) silencing transcription factor (REST), also known as neuron-restrictive silencer factor (NRSF) and X2 box repressor (XBR), was found to regulate gene transcription by binding to chromatin and recruiting chromatin-modifying enzymes. Earlier studies revealed that REST plays an important role in the development and disease of the nervous system, mainly by repressing the transcription of neuron-specific genes. Subsequently, REST was found to be critical in other tissues, such as the heart, pancreas, skin, eye, and vascular. Dysregulation of REST was also found in nervous and non-nervous system cancers. In parallel, multiple strategies to target REST have been developed. In this paper, we provide a comprehensive summary of the research progress made over the past 28 years since the discovery of REST, encompassing both physiological and pathological aspects. These insights into the effects and mechanisms of REST contribute to an in-depth understanding of the transcriptional regulatory mechanisms of genes and their roles in the development and progression of disease, with a view to discovering potential therapeutic targets and intervention strategies for various related diseases.

Molecular Mechanisms of Functional Modulation of Transcriptional Coactivator PC4 via Phosphorylation on Its Intrinsically Disordered Region

To investigate the effects of phosphorylation on the function of the human positive cofactor 4 (PC4), an enhanced molecular dynamics (MD) simulation was performed. The simulation system consists of the N-terminal intrinsic disordered region (IDR) of PC4 and a complex that comprises the C-terminal acidic activation domain of a herpes simplex virion protein 16 (VP16ad) and a homodimer of the C-terminal structured core domain of PC4 (PC4ctd). An earlier report of an experimental study reported that the PC4-VP16ad interaction is modulated by incremental phosphorylation of the IDR. That report also proposed a dynamic model where phosphorylated serine residues of a segment (SEAC) in the IDR contact positively charged residues (lysin and arginine) of another segment (K-rich) in the IDR. This contact formation induced by the phosphorylation results in variation of PC4-VP16ad interaction. However, this contact formation has not yet been measured directly because it is transiently formed and because the SEAC and K-rich segments are unstructured with high flexibility. We performed two simulations to mimic the incremental phosphorylation: the IDR was not phosphorylated in one simulation and only partially phosphorylated in the other. Our simulation results indicate that the phosphorylation weakens the IDR-VP16ad contact considerably with the induction of a compact structure in the IDR. This structure was stabilized by electrostatic interactions between the phosphorylated serine residues of a segment and lysine or arginine residues of another segment in the IDR, but the conformational fluctuation of this compact structure was considerably large. Consequently, the present study supports the experimentally proposed dynamic model. Results of this study can be important for computational elucidation of the functional modulation of PC4.

Computer simulation of molecular recognition in biomolecular system: from in silico screening to generalized ensembles

Prediction of ligand-receptor complex structure is important in both the basic science and the industry such as drug discovery. We report various computation molecular docking methods: fundamental in silico (virtual) screening, ensemble docking, enhanced sampling (generalized ensemble) methods, and other methods to improve the accuracy of the complex structure. We explain not only the merits of these methods but also their limits of application and discuss some interaction terms which are not considered in the in silico methods. In silico screening and ensemble docking are useful when one focuses on obtaining the native complex structure (the most thermodynamically stable complex). Generalized ensemble method provides a free-energy landscape, which shows the distribution of the most stable complex structure and semi-stable ones in a conformational space. Also, barriers separating those stable structures are identified. A researcher should select one of the methods according to the research aim and depending on complexity of the molecular system to be studied.

Fly casting with ligand sliding and orientational selection supporting complex formation of a GPCR and a middle sized flexible molecule

A GA-guided multidimensional virtual-system coupled molecular dynamics (GA-mD-VcMD) simulation was conducted to elucidate binding mechanisms of a middle-sized flexible molecule, bosentan, to a GPCR protein, human endothelin receptor type B (hETB). GA-mD-VcMD is a generalized ensemble method that produces a free-energy landscape of the ligand-receptor binding by searching large-scale motions accompanied with stable maintenance of the fragile cell-membrane structure. All molecular components (bosentan, hETB, membrane, and solvent) were represented with an all-atom model. Then sampling was conducted from conformations where bosentan was distant from the binding site in the hETB binding pocket. The deepest basin in the resultant free-energy landscape was assigned to native-like complex conformation. The following binding mechanism was inferred. First, bosentan fluctuating randomly in solution is captured using a tip region of the flexible N-terminal tail of hETB via nonspecific attractive interactions (fly casting). Bosentan then slides occasionally from the tip to the root of the N-terminal tail (ligand–sliding). During this sliding, bosentan passes the gate of the binding pocket from outside to inside of the pocket with an accompanying rapid reduction of the molecular orientational variety of bosentan (orientational selection). Last, in the pocket, ligand–receptor attractive native contacts are formed. Eventually, the native-like complex is completed. The bosentan-captured conformations by the tip-region and root-region of the N-terminal tail correspond to two basins in the free-energy landscape. The ligand-sliding corresponds to overcoming of a free-energy barrier between the basins.

Visualized and Quantitative Conformational Analysis of Peptidomimetics

Protein-protein interactions (PPIs) are fundamentally important and challenging drug targets. Peptidomimetic molecules of various types have been developed to modulate PPIs. A particularly promising drug discovery strategy, structural peptidomimetics, was designed based on special mimicking of side-chain C_α-C_β bonds. It is simple and versatile. Nevertheless, no quantitative method has been established to evaluate its similarity to a target peptide motif. We developed two methods that enable visual, comprehensive, and quantitative analysis of peptidomimetics: peptide conformation distribution (PCD) plot and peptidomimetic analysis (PMA) map. These methods specifically examine multiple side-chain C_α-C_β bonds of a peptide fragment motif and their corresponding bonds (pseudo-C_α-C_β bonds) in a mimetic molecule instead of φ and ψ angles of a single amino acid in the traditional Ramachandran plot. The PCD plot is an alignment-free method, whereas the PMA map is an alignment-based method providing distinctive and complementary analysis. Results obtained from analysis using these two methods indicate our multifacial α-helix mimetic scaffold 12 as an excellent peptidomimetic that can precisely mimic the spatial positioning of side-chain functional groups of α-helix. These methods are useful for visualized and quantified evaluation of peptidomimetics and for the rational design of new mimetic scaffolds.