Induced fit or conformational selection for RNA/U1A folding

Chemoproteomic capture of RNA binding activity in living cells

Proteomic methods for RNA interactome capture (RIC) rely principally on crosslinking native or labeled cellular RNA to enrich and investigate RNA-binding protein (RBP) composition and function in cells. The ability to measure RBP activity at individual binding sites by RIC, however, has been more challenging due to the heterogenous nature of peptide adducts derived from the RNA-protein crosslinked site. Here, we present an orthogonal strategy that utilizes clickable electrophilic purines to directly quantify protein-RNA interactions on proteins through photoaffinity competition with 4-thiouridine (4SU)-labeled RNA in cells. Our photo-activatable-competition and chemoproteomic enrichment (PACCE) method facilitated detection of >5500 cysteine sites across ~3000 proteins displaying RNA-sensitive alterations in probe binding. Importantly, PACCE enabled functional profiling of canonical RNA-binding domains as well as discovery of moonlighting RNA binding activity in the human proteome. Collectively, we present a chemoproteomic platform for global quantification of protein-RNA binding activity in living cells.

Cytotoxic Activity of Phytoconstituents Isolated from Monotheca buxifolia against Hepatocellular Carcinoma Cell Line HepG2: In Vitro and Molecular Docking Studies

Natural products and conventional chemotherapeutic drugs are believed to enhance anticancer treatment efficacy while lowering toxicity. The current study investigates the cytotoxic and apoptogenic effects of Monotheca buxifolia bioactive compounds on HepG2 cell lines. MTT assay was used to assess the effect on the viability of HepG2 cells. Morphological changes were investigated. Annexin-V-FITC/PI was used to demonstrate apoptotic activity. A molecular dynamics simulation study was carried out to investigate the compound binding pattern in the active site of the PPRAδ protein. MTT and annexin V-FITC/PI assays revealed that the isolated compounds lauric acid, oleanolic acid, and bis(2-ethylhexyl) phthalate inhibited the growth of hepatocellular cancer cells. The IC50 value for lauric acid was 56.46 ± 1.20 μg/mL, 31.94 ± 1.03 μg/mL for oleanolic acid, and 83.80 ± 2.18 μg/mL for bis(2-ethylhexyl) phthalate. Apoptosis was observed in 29.5, 52.1 and 22.4% of HepG2 cells treated with lauric acid, oleanolic acid, and bis(2-ethylhexyl) phthalate, respectively, after 24 h of treatment. Morphological assays and Hoechst staining microscopy revealed that the treatment caused morphological changes in the cell membrane and nuclear condensation. The high fluctuation indicates that various interactions were highly potent and widely adopted, and vice versa. Oleanolic acid displayed high residue fluctuation, remaining stable in the active site of the PPRAδ protein and involved in various interactions while remaining locally fluctuating in the binding sites of the other two compounds. These findings concluded that lauric acid, oleanolic acid, and bis(2-ethylhexyl) phthalate have a significant apoptogenic effect against HepG2 cells in inducing apoptosis. Our findings suggest that these bioactive compounds could be used as adjuvant therapies.

Binding to the conserved and stably folded guide RNA pseudoknot induces Cas12a conformational changes during ribonucleoprotein assembly

Ribonucleoproteins (RNPs) comprise one or more RNA and protein molecules that interact to form a stable complex, which commonly involves conformational changes in the more flexible RNA components. Here, we propose that Cas12a RNP assembly with its cognate CRISPR RNA (crRNA) guide instead proceeds primarily through Cas12a conformational changes during binding to more stable, prefolded crRNA 5' pseudoknot handles. Phylogenetic reconstructions and sequence and structure alignments revealed that the Cas12a proteins are divergent in sequence and structure while the crRNA 5' repeat region, which folds into a pseudoknot and anchors binding to Cas12a, is highly conserved. Molecular dynamics simulations of three Cas12a proteins and their cognate guides revealed substantial flexibility for unbound apo-Cas12a. In contrast, crRNA 5' pseudoknots were predicted to be stable and independently folded. Limited trypsin hydrolysis, differential scanning fluorimetry, thermal denaturation, and CD analyses supported conformational changes of Cas12a during RNP assembly and an independently folded crRNA 5' pseudoknot. This RNP assembly mechanism may be rationalized by evolutionary pressure to conserve CRISPR loci repeat sequence, and therefore guide RNA structure, to maintain function across all phases of the CRISPR defense mechanism.

Personal Precise Force Field for Intrinsically Disordered and Ordered Proteins Based on Deep Learning

Intrinsically disordered proteins (IDPs) are proteins without a fixed three-dimensional (3D) structure under physiological conditions and are associated with Parkinson's disease, Alzheimer's disease, cancer, cardiovascular disease, amyloidosis, diabetes, and other diseases. Experimental methods can hardly capture the ensemble of diverse conformations for IDPs. Molecular dynamics (MD) simulations can sample continuous conformations that might provide a valuable complement to experimental data. However, the accuracy of MD simulations depends on the quality of force field. In particular, the evolutionary conservation and coevolution of IDPs introduce that current force fields could not precisely reproduce the conformation of IDPs. In order to improve the performance of force field, deep learning and reweighting methods were used to automatically generate personal force field parameters for intrinsically disordered and ordered proteins. At first, the deep learning method predicted more accuracy φ/ψ dihedral of residue than the previous method. Then, reweighting optimized the personal force field parameters for each residue. Finally, typical representative systems such as IDPs, structure protein, and fast-folding protein were used to evaluate this force field. The results indicate that two personal force field parameters (named PPFF1 and PPFF1_af2) could better reproduce the experimental observables than ff03CMAP force field. In summary, this strategy will provide feasibility for the development of precise personal force fields.

PTMint database of experimentally verified PTM regulation on protein-protein interaction

MOTIVATION

Post-translational modification (PTM) is an important biochemical process. which includes six most well-studied types: phosphorylation, acetylation, methylation, sumoylation, ubiquitylation and glycosylation. PTM is involved in various cell signaling pathways and biological processes. Abnormal PTM status is closely associated with severe diseases (such as cancer and neurologic diseases) by regulating protein functions, such as protein-protein interactions (PPIs). A set of databases was constructed separately for PTM sites and PPI; however, the resource of regulation for PTM on PPI is still unsolved.

RESULTS

Here, we firstly constructed a public accessible database of PTMint (PTMs that are associated with PPIs) (https://ptmint.sjtu.edu.cn/) that contains manually curated complete experimental evidence of the PTM regulation on PPIs in multiple organisms, including Homo sapiens, Arabidopsis thaliana, Caenorhabditis elegans, Drosophila melanogaster, Saccharomyces cerevisiae and Schizosaccharomyces pombe. Currently, the first version of PTMint encompassed 2477 non-redundant PTM sites in 1169 proteins affecting 2371 protein-protein pairs involving 357 diseases. Various annotations were systematically integrated, such as protein sequence, structure properties and protein complex analysis. PTMint database can help to insight into disease mechanism, disease diagnosis and drug discovery associated with PTM and PPI.

AVAILABILITY AND IMPLEMENTATION

PTMint is freely available at: https://ptmint.sjtu.edu.cn/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Molecular Dynamics Simulations of Protein RNA Complexes by Using an Advanced Electrostatic Model

Protein-RNA interactions are integral to the biological functions of RNA. It is well recognized that molecular dynamics (MD) simulations of protein-RNA complexes are more challenging than those of each component. The difficulty arises from the strong electrostatic interactions and the delicate balance between various types of physical forces at the interface. Previously, MD simulations of protein-RNA complexes have predominantly employed fixed-charge force fields. Although force field modifications have been developed to address problems identified in the simulations, some protein-RNA structures are still hard to reproduce by simulations. Here, we present MD simulations of two representative protein-RNA complexes using the AMOEBA polarizable force field. The van der Waals parameters were refined to reproduce accurate quantum-mechanical data of base-base and base-amino acid interactions. It was found that the refined parameters produced a more stable hydrogen-bond network in the interface. One of the complexes remained stable during the short simulations, whereas it could quickly break down in previous simulations using fixed-charge force fields. There was reversible breaking and formation of hydrogen bonds that are observed in the crystal structure, which may indicate the difference in solution and crystal structures. While further improvement and validation of the force fields are still needed, this work demonstrates that polarizable force fields are promising for the study of protein-RNA complexes.

Understanding the binding specificities of mRNA targets by the mammalian Quaking protein

Mammalian Quaking (QKI) protein, a member of STAR family of proteins is a mRNA-binding protein, which post-transcriptionally modulates the target RNA. QKI protein possesses a maxi-KH domain composed of single heterogeneous nuclear ribonucleoprotein K homology (KH) domain and C-terminal QUA2 domain, that binds a sequence-specific QKI RNA recognition element (QRE), CUAAC. To understand the binding specificities for different mRNA sequences of the KH-QUA2 domain of QKI protein, we introduced point mutations at different positions in the QRE resulting in twelve different mRNA sequences with single nucleotide change. We carried out long unbiased molecular dynamics simulations using two different sets of recently updated forcefield parameters: AMBERff14SB+RNAχOL3 and CHARMM36 (with CMAP correction). We analyzed the changes in intermolecular dynamics as a result of mutation. Our results show that AMBER forcefields performed better to model the interactions between mRNA and protein. We also calculated the binding affinities of different mRNA sequences and found that the relative order correlates to the reported experimental studies. Our study shows that the favorable binding with the formation of stable complex will occur when there is an increase of the total intermolecular contacts between mRNA and protein, but without the loss of native contacts within the KH-QUA domain.

Interpreting the Dynamics of Binding Interactions of snRNA and U1A Using a Coarse-Grained Model

The binding interactions of small nuclear RNAs (snRNA) and the associated protein factors are critical to the function of spliceosomes in alternatively splicing primary RNA transcripts. Although molecular dynamics simulations are a powerful tool to interpret the mechanism of biological processes, the atomic-level simulations are, however, too expensive and with limited accuracy for the large-size systems, such as snRNA-protein complexes. We extend the coarse-grained Gaussian network model, which models the RNA-protein complexes as a harmonic chain of Cα, P, and O4' atoms, to investigating the impact of the snRNA-binding interaction on the dynamic stability of the human U1A protein, which is a major component of the spliceosomal U1 small nuclear ribonucleoprotein particle. The results reveal that the first and third loops and the C-terminal helix regions of the U1A domain undergo a significant loss of flexibility upon the RNA binding due to the forming of mostly electrostatic and hydrogen bond interactions with RNA 5' stem and loop. By examining the residues whose mutations significantly change the binding free energy between U1A and snRNA, the Gaussian network model-based calculations show that not only the residues at the binding sites that are traditionally considered to play a major role in U1A-RNA association but also those residues that are far away from the RNA-binding interface can participate in the long-range allosteric signal transmission; these calculations are quantitatively consistent with the data observed in the recent snRNA binding experiments. The study demonstrates a useful avenue to utilize the simplified elastic network model to investigate the dynamics characteristics of the biologically important macromolecular interactions.

Order-disorder transition of intrinsically disordered kinase inducible transactivation domain of CREB

Transcription factor cyclic Adenosine monophosphate response-element binding protein plays a critical role in the cyclic AMP response pathway via its intrinsically disordered kinase inducible transactivation domain (KID). KID is one of the most studied intrinsically disordered proteins (IDPs), although most previous studies focus on characterizing its disordered state structures. An interesting question that remains to be answered is how the order-disorder transition occurs at experimental conditions. Thanks to the newly developed IDP-specific force field ff14IDPSFF, the quality of conformer sampling for IDPs has been dramatically improved. In this study, molecular dynamics (MD) simulations were used to study the order-to-disorder transition kinetics of KID based on the good agreement with the experiment on its disordered-state properties. Specifically, we tested four force fields, ff99SBildn, ff99IDPs, ff14IDPSFF, and ff14IDPs in the simulations of KID and found that ff14IDPSFF can generate more diversified disordered conformers and also reproduce more accurate experimental secondary chemical shifts. Kinetics analysis of MD simulations demonstrates that the order-disorder transition of KID obeys the first-order kinetics, and the transition nucleus is I127/L128/L141. The possible transition pathways from the nucleus to the last folded residues were identified as I127-R125-L138-L141-S143-A145 and L128-R125-L138-L141-S143-A145 based on a residue-level dynamical network analysis. These computational studies not only provide testable prediction/hypothesis on the order-disorder transition of KID but also confirm that the ff14IDPSFF force field can be used to explore the correlation between the structure and function of IDPs.

Chimeric Guides Probe and Enhance Cas9 Biochemical Activity

DNA substitutions in RNA can probe the importance of A-form structure, 2'-hydroxyl contacts, and conformational constraints within RNA-guided enzymes. Using this approach, we found that Cas9 biochemical activity tolerated significant substitution with DNA nucleotides in the clustered regularly interspaced short palindromic repeat RNA (crRNA). Only minimal RNA content was needed in or near the seed region. Simultaneous substitution at all positions with predicted crRNA-Cas9 2'-hydroxyl contacts had no effect on enzyme activity. The trans-activating crRNA (tracrRNA) also tolerated >50% substitution with DNA. DNA substitutions in the tracrRNA-pairing region of crRNA consistently enhanced cleavage activity while maintaining or improving target specificity. Together, results point to a prominent role for guide:target A-form-like helical structure and a possible regulatory role for the crRNA-tracrRNA pairing motif. A model chimeric crRNA with high activity did not significantly alter RNP assembly or target binding but did reduce Cas9 ribonucleoprotein stability, suggesting effects through conformation or dynamics. Cas9 directed by chimeric RNA-DNA guides may represent a cost-effective synthetic or molecular biology tool for robust and specific DNA cleavage.

Synergistic regulation mechanism of iperoxo and LY2119620 for muscarinic acetylcholine M2 receptor

Muscarinic acetylcholine receptors are GPCRs that regulate the activity of a diverse array of central and peripheral functions in the human body, including the parasympathetic actions of acetylcholine. The M2 muscarinic receptor subtype plays a key role in modulating cardiac function and many important central processes. The orthosteric agonist and allosteric modulator can bind the pocket of M2. However, the detailed relationship between orthosteric agonist and allosteric modulator of M2 is still unclear. In this study, we intend to elucidate the residue-level regulation mechanism and pathway via a combined approach of dynamical correlation network and molecular dynamics simulation. Specifically computational residue-level fluctuation correlation data was analyzed to reveal detailed dynamics signatures in the regulation process. A hypothesis of "synergistic regulation" is proposed to reveal the cooperation affection between the orthosteric agonist and allosteric modulator, which is subsequently validated by perturbation and mutation analyses. Two possible synergistic regulation pathways of 2CU-I178-Y403-W400-F396-L114-Y440-Nb9 and IXO-V111-F396-L114-Y440-Nb9 were identified by the shortest path algorithm and were confirmed by the mutation of junction node. Furthermore, the efficiency of information transfer of bound M2 is significant higher than any single binding system. Our study shows that targeting the synergistic regulation pathways may better regulate the calcium channel of M2. The knowledge gained in this study may help develop drugs for diseases of the central nervous system and metabolic disorders.

Structural Changes of RNA in Complex with Proteins in the SRP

The structural flexibility of RNA allows it to exist in several shapes and sizes. Thus, RNA is functionally diverse and is known to be involved in processes such as catalysis, ligand binding, and most importantly, protein recognition. RNA can adopt different structures, which can often dictate its functionality. When RNA binds onto protein to form a ribonucleoprotein complex (RNP), multiple interactions and conformational changes occur with the RNA and protein. However, there is the question of whether there is a specific pattern for these changes to occur upon recognition. In particular when RNP complexity increases with the addition of multiple proteins/RNA, it becomes difficult to structurally characterize the overall changes using the current structural determination techniques. Hence, there is a need to use a combination of biochemical, structural and computational modeling to achieve a better understanding of the processes that RNPs are involved. Nevertheless, there are well-characterized systems that are evolutionarily conserved [such as the signal recognition particle (SRP)] that give us important information on the structural changes of RNA and protein upon complex formation.

Mechanism of mRNA-STAR domain interaction: Molecular dynamics simulations of Mammalian Quaking STAR protein

STAR proteins are evolutionary conserved mRNA-binding proteins that post-transcriptionally regulate gene expression at all stages of RNA metabolism. These proteins possess conserved STAR domain that recognizes identical RNA regulatory elements as YUAAY. Recently reported crystal structures show that STAR domain is composed of N-terminal QUA1, K-homology domain (KH) and C-terminal QUA2, and mRNA binding is mediated by KH-QUA2 domain. Here, we present simulation studies done to investigate binding of mRNA to STAR protein, mammalian Quaking protein (QKI). We carried out conventional MD simulations of STAR domain in presence and absence of mRNA, and studied the impact of mRNA on the stability, dynamics and underlying allosteric mechanism of STAR domain. Our unbiased simulations results show that presence of mRNA stabilizes the overall STAR domain by reducing the structural deviations, correlating the ‘within-domain’ motions, and maintaining the native contacts information. Absence of mRNA not only influenced the essential modes of motion of STAR domain, but also affected the connectivity of networks within STAR domain. We further explored the dissociation of mRNA from STAR domain using umbrella sampling simulations, and the results suggest that mRNA binding to STAR domain occurs in multi-step: first conformational selection of mRNA backbone conformations, followed by induced fit mechanism as nucleobases interact with STAR domain.

Allosteric mechanism of cyclopropylindolobenzazepine inhibitors for HCV NS5B RdRp via dynamic correlation network analysis

HCV RNA dependent RNA polymerase (RdRp) nonstructural protein 5B (NS5B) is a major target against hepatitis C virus (HCV) for antiviral therapy. Recently discovered cyclopropylindolobenzazepine derivatives have been considered as the most potent for their ability to bind the thumb site 1 domain and allosterically inhibit HCV NS5B RdRp activity. However, the allosteric mechanism for these derivatives has not been clarified at the molecular level. In this study, fluctuation correlation networks were constructed based on all-atom molecular dynamics simulations to elucidate the allosteric mechanism. The fluctuation correlation networks between free and M2 bound NS5B are significantly different. Information can better transfer from the allosteric site to the catalytic site for bound NS5B than for free NS5B. Thus, the hypothesis of "binding induced allosteric regulation" is proposed to link the enzyme activation and inhibitor binding and then confirmed by the mutant network. Finally, one possible allosteric pathway was identified with the shortest path and evaluated by the perturbation of the network. These methods will be helpful to identify the allosteric pathway of other proteins and to design new drugs targeting the pathway.

Conformation Dynamics of the Intrinsically Disordered Protein c-Myb with the ff99IDPs Force Field

The intrinsically disordered protein c-Myb plays a critical role in cellular proliferation and differentiation. Loss of c-myb function results in embryonic lethality due to failure of fetal hepatic hematopoiesis. The conformation dynamics of the intrinsically disordered c-Myb are still unknown. Here, molecular dynamics (MD) simulations with the intrinsically disordered protein force field ff99IDPs were used to study the conformation dynamics. In comparison with ff99SBildn, ff99IDPs can reproduce more diverse disordered conformers of c-Myb. The predicted secondary chemical shift under ff99IDPs is more close to that of experiment data than that under ff99SBildn. Therefore, ff99IDPs can sample native molten globule, native pre-molten globule and native coil conformers for c-Myb. These results are consistent with those of other intrinsically disordered proteins. Kinetic analysis of MD simulations shows that c-Myb folds via a two-state process and indicates that c-Myb folds in the order of tertiary folding and helical folding. The folding nucleus of KEL plays an essential role in stabilizing the folding state with dynamic correlation networks. The influences of solvent models for TIP3P, TIP4P-EW and TIP5P were also investigated and it was found that TIP3P and ff99IDPs are the best combination to research the conformer sampling of c-Myb. These results reveal the conformation dynamics of c-Myb and confirm that the ff99IDPs force field can be used to research the relationship between structure and function of other intrinsically disordered proteins.

Crystal Structure of StnA for the Biosynthesis of Antitumor Drug Streptonigrin Reveals a Unique Substrate Binding Mode

Streptonigrin methylesterase A (StnA) is one of the tailoring enzymes that modify the aminoquinone skeleton in the biosynthesis pathway of Streptomyces species. Although StnA has no significant sequence homology with the reported α/β-fold hydrolases, it shows typical hydrolytic activity in vivo and in vitro. In order to reveal its functional characteristics, the crystal structures of the selenomethionine substituted StnA (SeMet-StnA) and the complex (S185A mutant) with its substrate were resolved to the resolution of 2.71 Å and 2.90 Å, respectively. The overall structure of StnA can be described as an α-helix cap domain on top of a common α/β hydrolase domain. The substrate methyl ester of 10'-demethoxystreptonigrin binds in a hydrophobic pocket that mainly consists of cap domain residues and is close to the catalytic triad Ser185-His349-Asp308. The transition state is stabilized by an oxyanion hole formed by the backbone amides of Ala102 and Leu186. The substrate binding appears to be dominated by interactions with several specific hydrophobic contacts and hydrogen bonds in the cap domain. The molecular dynamics simulation and site-directed mutagenesis confirmed the important roles of the key interacting residues in the cap domain. Structural alignment and phylogenetic tree analysis indicate that StnA represents a new subfamily of lipolytic enzymes with the specific binding pocket located at the cap domain instead of the interface between the two domains.

Protein-RNA interactions: structural biology and computational modeling techniques

RNA-binding proteins are functionally diverse within cells, being involved in RNA-metabolism, translation, DNA damage repair, and gene regulation at both the transcriptional and post-transcriptional levels. Much has been learnt about their interactions with RNAs through structure determination techniques and computational modeling. This review gives an overview of the structural data currently available for protein-RNA complexes, and discusses the technical issues facing structural biologists working to solve their structures. The review focuses on three techniques used to solve the 3-dimensional structure of protein-RNA complexes at atomic resolution, namely X-ray crystallography, solution nuclear magnetic resonance (NMR) and cryo-electron microscopy (cryo-EM). The review then focuses on the main computational modeling techniques that use these atomic resolution data: discussing the prediction of RNA-binding sites on unbound proteins, docking proteins, and RNAs, and modeling the molecular dynamics of the systems. In conclusion, the review looks at the future directions this field of research might take.

Dynamics Correlation Network for Allosteric Switching of PreQ1 Riboswitch

Riboswitches are a class of metabolism control elements mostly found in bacteria. Due to their fundamental importance in bacteria gene regulation, riboswitches have been proposed as antibacterial drug targets. Prequeuosine (preQ₁) is the last free precursor in the biosynthetic pathway of queuosine that is crucial for translation efficiency and fidelity. However, the regulation mechanism for the preQ₁ riboswitch remains unclear. Here we constructed fluctuation correlation network based on all-atom molecular dynamics simulations to reveal the regulation mechanism. The results suggest that the correlation network in the bound riboswitch is distinctly different from that in the apo riboswitch. The community network indicates that the information freely transfers from the binding site of preQ₁ to the expression platform of the P3 helix in the bound riboswitch and the P3 helix is a bottleneck in the apo riboswitch. Thus, a hypothesis of “preQ₁-binding induced allosteric switching” is proposed to link riboswitch and translation regulation. The community networks of mutants support this hypothesis. Finally, a possible allosteric pathway of A50-A51-A52-U10-A11-G12-G56 was also identified based on the shortest path algorithm and confirmed by mutations and network perturbation. The novel fluctuation network analysis method can be used as a general strategy in studies of riboswitch structure-function relationship.

Synergistic Allosteric Mechanism of Fructose-1,6-bisphosphate and Serine for Pyruvate Kinase M2 via Dynamics Fluctuation Network Analysis

Pyruvate kinase M2 (PKM2) plays a key role in tumor metabolism and regulates the rate-limiting final step of glycolysis. In tumor cells, there are two allosteric effectors for PKM2: fructose-1,6-bisphosphate (FBP) and serine. However, the relationship between FBP and serine for allosteric regulation of PKM2 is unknown. Here we constructed residue/residue fluctuation correlation network based on all-atom molecular dynamics simulations to reveal the regulation mechanism. The results suggest that the correlation network in bound PKM2 is distinctly different from that in the free state, FBP/PKM2, or Ser/PKM2. The community network analysis indicates that the information can freely transfer from the allosteric sites of FBP and serine to the substrate site in bound PKM2, while there exists a bottleneck for information transfer in the network of the free state. Furthermore, the binding free energy between the substrate and PKM2 for bound PKM2 is significantly lower than either of FBP/PKM2 or Ser/PKM2. Thus, a hypothesis of "synergistic allosteric mechanism" is proposed for the allosteric regulation of FBP and serine. This hypothesis was further confirmed by the perturbational and mutational analyses of community networks and binding free energies. Finally, two possible synergistic allosteric pathways of FBP-K433-T459-R461-A109-V71-R73-MG2-OXL and Ser-I47-C49-R73-MG2-OXL were identified based on the shortest path algorithm and were confirmed by the network perturbation analysis. Interestingly, no similar pathways could be found in the free state. The process targeting on the allosteric pathways can better regulate the glycolysis of PKM2 and significantly inhibit the progression of tumor.

Selectivity Mechanism of ATP-Competitive Inhibitors for PKB and PKA

Protein kinase B (PKB) acts as a central node on the PI3K kinase pathway. Constitutive activation and overexpression of PKB have been identified to involve in various cancers. However, protein kinase A (PKA) sharing high homology with PKB is essential for metabolic regulation. Therefore, specific targeting on PKB is crucial strategy in drug design and development for antitumor. Here, we had revealed the selectivity mechanism for PKB inhibitors with molecular dynamics simulation and 3D-QSAR methods. Selective inhibitors of PKB could form more hydrogen bonds and hydrophobic contacts with PKB than those with PKA. This could explain that selective inhibitor M128 is more potent to PKB than to PKA. Then, 3D-QSAR models were constructed for these selective inhibitors and evaluated by test set compounds. 3D-QSAR model comparison of PKB inhibitors and PKA inhibitors reveals possible methods to improve the selectivity of inhibitors. These models can be used to design new chemical entities and make quantitative prediction of the specific selective inhibitors before resorting to in vitro and in vivo experiment.

Combined mechanism of conformational selection and induced fit in U1A-RNA molecular recognition

In this study, we demonstrate that U1A-RNA molecular recognition is mediated by a combined mechanism of conformational selection and induced fit. The binding of U1A to RNA has been discussed in the context of induced fit that involves the reorientation of the α-helix in the C-terminal region (Helix-C) of U1A to permit RNA access only when U1A correctly recognizes RNA. However, according to our molecular dynamics simulations, even in the absence of RNA, Helix-C spontaneously reoriented to permit RNA access. Nonetheless, such a conformational change was still incomplete. Helix-C was often partially or even fully unfolded and in an infrequent RNA-accessible conformation, which can be detected using state-of-the-art nuclear magnetic resonance methodology. These results suggest that the formation of an energetically stabilized complex is promoted by specific interactions between U1A and RNA. In conclusion, in the recognition of RNA by U1A protein, we propose a combined mechanism that requires the reorientation of Helix-C and the subsequent contact with RNA through conformational selection, although the stabilization of the U1A-RNA complex is caused by induced fit. We further propose a modification to the conventional assumption regarding the mechanism of U1A-RNA molecular recognition.

Kink turn sRNA folding upon L7Ae binding using molecular dynamics simulations

The kink-turn sRNA motif in archaea, whose combination with protein L7Ae initializes the assembly of small ribonucleoprotein particles (sRNPs), plays a key role in ribosome maturation and the translation process. Although many studies have been reported on this motif, the mechanism of sRNA folding coupled with protein binding is still poorly understood. Here, room and high temperature molecular dynamics (MD) simulations were performed on the complex of 25-nt kink-turn sRNA and L7Ae. The average RMSD values between the bound and corresponding apo structures and Kolmogorov-Smirnov P test analysis indicate that sRNA may follow an induced fit mechanism upon binding with L7Ae, both locally and globally. These conclusions are further supported by high-temperature unfolding kinetic analysis. Principal component analysis (PCA) found both closing and opening motions of the kink-turn sRNA. This might play a key role in the sRNP assembly and methylation catalysis. These combined computational methods can be used to study the specific recognition of other sRNAs and proteins.

Insights into the Binding of Intrinsically Disordered Proteins from Molecular Dynamics Simulation

Intrinsically disordered proteins (IDPs) are a class of protein that, in the native state, possess no well-defined secondary or tertiary structure, existing instead as dynamic ensembles of conformations. They are biologically important, with approximately 20% of all eukaryotic proteins disordered, and found at the heart of many biochemical networks. To fulfil their biological roles, many IDPs need to bind to proteins and/or nucleic acids. And while unstructured in solution, IDPs typically fold into a well-defined three-dimensional structure upon interaction with a binding partner. The flexibility and structural diversity inherent to IDPs makes this coupled folding and binding difficult to study at atomic resolution by experiment alone, and computer simulation currently offers perhaps the best opportunity to understand this process. But simulation of coupled folding and binding is itself extremely challenging; these molecules are large and highly flexible, and their binding partners, such as DNA or cyclins, are also often large. Therefore, their study requires either or both simplified representations and advanced enhanced sampling schemes. It is not always clear that existing simulation techniques, optimized for studying folded proteins, are well-suited to IDPs. In this article, we examine the progress that has been made in the study of coupled folding and binding using molecular dynamics simulation. We summarise what has been learnt, and examine the state of the art in terms of both methodologies and models. We also consider the lessons to be learnt from advances in other areas of simulation and highlight the issues that remain of be addressed.

Global conformational selection and local induced fit for the recognition between intrinsic disordered p53 and CBP

The transactivation domain (TAD) of tumor suppressor p53 can bind with the nuclear coactivator binding domain (NCBD) of cyclic-AMP response element binding protein (CBP) and activate transcription. NMR experiments demonstrate that both apo-NCBD and TAD are intrinsic disordered and bound NCBD/TAD undergoes a transition to well folded. The recognition mechanism between intrinsic disordered proteins is still hotly debated. Molecular dynamics (MD) simulations in explicit solvent are used to study the recognition mechanism between intrinsic disordered TAD and NCBD. The average RMSD values between bound and corresponding apo states and Kolmogorov-Smirnov P test analysis indicate that TAD and NCBD may follow an induced fit mechanism. Quantitative analysis indicates there is also a global conformational selection. In summary, the recognition of TAD and NCBD might obey a local induced fit and global conformational selection. These conclusions are further supported by high-temperature unbinding kinetics and room temperature landscape analysis. These methods can be used to study the recognition mechanism of other intrinsic disordered proteins.

Single-molecule observation of the induction of k-turn RNA structure on binding L7Ae protein

The k-turn is a commonly occurring structural motif that introduces a tight kink into duplex RNA. In free solution, it can exist in an extended form, or by folding into the kinked structure. Binding of proteins including the L7Ae family can induce the formation of the kinked geometry, raising the question of whether this occurs by passive selection of the kinked structure, or a more active process in which the protein manipulates the RNA structure. We have devised a single-molecule experiment whereby immobilized L7Ae protein binds Cy3-Cy5-labeled RNA from free solution. We find that all bound RNA is in the kinked geometry, with no evidence for transitions to an extended form at the millisecond timescale of the camera. Furthermore, real-time binding experiments provide no evidence for a more extended intermediate even at the earliest times, at a time resolution of 16 ms. The data support a passive conformational selection model by which the protein selects a fraction of RNA that is already in the kinked conformation, thereby drawing the equilibrium into this form.

Intrinsic flexibility of snRNA hairpin loops facilitates protein binding

Stem-loop II of U1 snRNA and Stem-loop IV of U2 snRNA typically have 10 or 11 nucleotides in their loops. The fluorescent nucleobase 2-aminopurine was used as a substitute for the adenines in each loop to probe the local and global structures and dynamics of these unusually long loops. Using steady-state and time-resolved fluorescence, we find that, while the bases in the loops are stacked, they are able to undergo significant local motion on the picosecond/nanosecond timescale. In addition, the loops have a global conformational change at low temperatures that occurs on the microsecond timescale, as determined using laser T-jump experiments. Nucleobase and loop motions are present at temperatures far below the melting temperature of the hairpin stem, which may facilitate the conformational change required for specific protein binding to these RNA loops.

Atomistic mechanism of microRNA translation upregulation via molecular dynamics simulations

MicroRNAs are endogenous 23–25 nt RNAs that play important gene-regulatory roles in animals and plants. Recently, miR369-3 was found to upregulate translation of TNFα mRNA in quiescent (G0) mammalian cell lines. Knock down and immunofluorescence experiments suggest that microRNA-protein complexes (with FXR1 and AGO2) are necessary for the translation upregulation. However the molecular mechanism of microRNA translation activation is poorly understood. In this study we constructed the microRNA-mRNA-AGO2-FXR1 quadruple complex by bioinformatics and molecular modeling, followed with all atom molecular dynamics simulations in explicit solvent to investigate the interaction mechanisms for the complex. A combined analysis of experimental and computational data suggests that AGO2-FXR1 complex relocalize microRNA:mRNA duplex to polysomes in G0. The two strands of dsRNA are then separated upon binding of AGO2 and FXR1. Finally, polysomes may improve the translation efficiency of mRNA. The mutation research confirms the stability of microRNA-mRNA-FXR1 and illustrates importance of key residue of Ile304. This possible mechanism can shed more light on the microRNA-dependent upregulation of translation.

Enhancing the promiscuous phosphotriesterase activity of a thermostable lactonase (GkaP) for the efficient degradation of organophosphate pesticides

The phosphotriesterase-like lactonase (PLL) enzymes in the amidohydrolase superfamily hydrolyze various lactones and exhibit latent phosphotriesterase activities. These enzymes serve as attractive templates for in vitro evolution of neurotoxic organophosphates (OPs) with hydrolytic capabilities that can be used as bioremediation tools. Here, a thermostable PLL from Geobacillus kaustophilus HTA426 (GkaP) was targeted for joint laboratory evolution with the aim of enhancing its catalytic efficiency against OP pesticides. By a combination of site saturation mutagenesis and whole-gene error-prone PCR approaches, several improved variants were isolated. The most active variant, 26A8C, accumulated eight amino acid substitutions and demonstrated a 232-fold improvement over the wild-type enzyme in reactivity (k(cat)/K(m)) for the OP pesticide ethyl-paraoxon. Concomitantly, this variant showed a 767-fold decrease in lactonase activity with δ-decanolactone, imparting a specificity switch of 1.8 × 10(5)-fold. 26A8C also exhibited high hydrolytic activities (19- to 497-fold) for several OP pesticides, including parathion, diazinon, and chlorpyrifos. Analysis of the mutagenesis sites on the GkaP structure revealed that most mutations are located in loop 8, which determines substrate specificity in the amidohydrolase superfamily. Molecular dynamics simulation shed light on why 26A8C lost its native lactonase activity and improved the promiscuous phosphotriesterase activity. These results permit us to obtain further insights into the divergent evolution of promiscuous enzymes and suggest that laboratory evolution of GkaP may lead to potential biological solutions for the efficient decontamination of neurotoxic OP compounds.

Insight into the stability of cross-β amyloid fibril from VEALYL short peptide with molecular dynamics simulation

Amyloid fibrils are found in many fatal neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, type II diabetes, and prion disease. The VEALYL short peptide from insulin has been confirmed to aggregate amyloid-like fibrils. However, the aggregation mechanism of amyloid fibril is poorly understood. Here, we utilized molecular dynamics simulation to analyse the stability of VEALYL hexamer. The statistical results indicate that hydrophobic residues play key roles in stabilizing VEALYL hexamer. Single point and two linkage mutants confirmed that Val1, Leu4, and Tyr5 of VEALYL are key residues. The consistency of the results for the VEALYL oligomer suggests that the intermediate states might be trimer (3-0) and pentamer(3-2). These results can help us to obtain an insight into the aggregation mechanism of amyloid fibril. These methods can be used to study the stability of amyloid fibril from other short peptides.

Revealing interaction mode between HIV-1 protease and mannitol analog inhibitor

HIV protease is a key enzyme to play a key role in the HIV-1 replication cycle and control the maturation from HIV viruses to an infectious virion. HIV-1 protease has become an important target for anti-HIV-1 drug development. Here, we used molecular dynamics simulation to study the binding mode between mannitol derivatives and HIV-1 protease. The results suggest that the most active compound (M35) has more stable hydrogen bonds and stable native contacts than the less active one (M17). These mannitol derivatives might have similar interaction mode with HIV-1 protease. Then, 3D-QSAR was used to construct quantitative structure-activity models. The cross-validated q(2) values are found as 0.728 and 0.611 for CoMFA and CoMSIA, respectively. And the non-cross-validated r(2) values are 0.973 and 0.950. Nine test set compounds validate the model. The results show that this model possesses better prediction ability than the previous work. This model can be used to design new chemical entities and make quantitative prediction of the bioactivities for HIV-1 protease inhibitors before resorting to in vitro and in vivo experiment.

RNA folding pathways and the self-assembly of ribosomes

Many RNAs do not directly code proteins but are nonetheless indispensable to cellular function. These strands fold into intricate three-dimensional shapes that are essential structures in protein synthesis, splicing, and many other processes of gene regulation and expression. A variety of biophysical and biochemical methods are now showing, in real time, how ribosomal subunits and other ribonucleoprotein complexes assemble from their molecular components. Footprinting methods are particularly useful for studying the folding of long RNAs: they provide quantitative information about the conformational state of each residue and require little material. Data from footprinting complement the global information available from small-angle X-ray scattering or cryo-electron microscopy, as well as the dynamic information derived from single-molecule Förster resonance energy transfer (FRET) and NMR methods. In this Account, I discuss how we have used hydroxyl radical footprinting and other experimental methods to study pathways of RNA folding and 30S ribosome assembly. Hydroxyl radical footprinting probes the solvent accessibility of the RNA backbone at each residue in as little as 10 ms, providing detailed views of RNA folding pathways in real time. In conjunction with other methods such as solution scattering and single-molecule FRET, time-resolved footprinting of ribozymes showed that stable domains of RNA tertiary structure fold in less than 1 s. However, the free energy landscapes for RNA folding are rugged, and individual molecules kinetically partition into folding pathways that lead through metastable intermediates, stalling the folding or assembly process. Time-resolved footprinting was used to follow the formation of tertiary structure and protein interactions in the 16S ribosomal RNA (rRNA) during the assembly of 30S ribosomes. As previously observed in much simpler ribozymes, assembly occurs in stages, with individual molecules taking different routes to the final complex. Interactions occur concurrently in all domains of the 16S rRNA, and multistage protection of binding sites of individual proteins suggests that initial encounter complexes between the rRNA and ribosomal proteins are remodeled during assembly. Equilibrium footprinting experiments showed that one primary binding protein was sufficient to stabilize the tertiary structure of the entire 16S 5'-domain. The rich detail available from the footprinting data showed that the secondary assembly protein S16 suppresses non-native structures in the 16S 5'-domain. In doing so, S16 enables a conformational switch distant from its own binding site, which may play a role in establishing interactions with other domains of the 30S subunit. Together, the footprinting results show how protein-induced changes in RNA structure are communicated over long distances, ensuring cooperative assembly of even very large RNA-protein complexes such as the ribosome.

Male germ cell apoptosis and epigenetic histone modification induced by Tripterygium wilfordii Hook F

Multiglycosides of Tripterygium wilfordii Hook f (GTW), a Chinese herb-derived medicine used as a remedy for rheumatoid arthritis, are considered to be a reversible anti-fertility drug affecting the mammalian spermatids. However, the mechanism behind this effect is still unknown. To study the possible mechanism behind the impact of GTW on spermatogenesis, we administered 4 groups of 4-week-old male mice with different doses of GTW. We found a dose-dependent decrease in the number of germ cells after 40 days of GTW treatment, and an increase in apoptotic cells from the low-dose to the high-dose group. During this same period the dimethylated level of histone H3 lysine 9 (H3K9me2) in GTW-treated testes germ cells declined. Additionally, spermatogonial stem cells (SSCs) from 6-day-old mice were isolated to evaluate the possible effect of GTW or triptolide on development of SSCs. We found a significantly higher incidence of apoptosis and lower dimethylation level of H3K9me2 in the SSCs of GTW or triptolide treatment than in controls. Thus, these data suggest that the GTW-induced apoptosis might be responsible for the fertility impairment in mice. This damage could be traced back to the early stages of spermatogenesis. GTW also affected the epigenetic modification of H3K9 in spermatogenesis. Molecular dynamics simulation suggested that triptolide and dimethylated or trimethylated H3K9 might have similar interaction mechanisms with EED (embryonic ectoderm development). These candidate activation mechanisms provide the first glimpse into the pathway of GTW-induced gonad toxicity, which is crucial for further research and clinical application.

In silico elucidation of the recognition dynamics of ubiquitin

Elucidation of the mechanism of biomacromolecular recognition events has been a topic of intense interest over the past century. The inherent dynamic nature of both protein and ligand molecules along with the continuous reshaping of the energy landscape during the binding process renders it difficult to characterize this process at atomic detail. Here, we investigate the recognition dynamics of ubiquitin via microsecond all-atom molecular dynamics simulation providing both thermodynamic and kinetic information. The high-level of consistency found with respect to experimental NMR data lends support to the accuracy of the in silico representation of the conformational substates and their interconversions of free ubiquitin. Using an energy-based reweighting approach, the statistical distribution of conformational states of ubiquitin is monitored as a function of the distance between ubiquitin and its binding partner Hrs-UIM. It is found that extensive and dense sampling of conformational space afforded by the µs MD trajectory is essential for the elucidation of the binding mechanism as is Boltzmann sampling, overcoming inherent limitations of sparsely sampled empirical ensembles. The results reveal a population redistribution mechanism that takes effect when the ligand is at intermediate range of 1–2 nm from ubiquitin. This mechanism, which may be depicted as a superposition of the conformational selection and induced fit mechanisms, also applies to other binding partners of ubiquitin, such as the GGA3 GAT domain.

Molecular recognition plays a central role in many biological processes, ensuring specific and efficient interaction between binding partners. Various models for describing the mechanisms of molecular recognition have been proposed, but the validation of these models has been traditionally difficult due to the transient and complex nature of the dynamic recognition process. In the present study, we aim at visually characterizing the mutual interplay between human ubiquitin and its ligands via microsecond time scale molecular dynamics simulation, which is validated rigorously against experimental NMR data. Taking advantage of Boltzmann sampling of molecular dynamics snapshots, we statistically reweight the populations of ubiquitin in the presence of its ligand molecule at intermediate distance range (1–2 nm) to examine the population redistribution mechanisms. These results offer new atomistic insights into this vital protein-protein recognition event.

Multistep kinetics of the U1A-SL2 RNA complex dissociation

The U1A-SL2 RNA complex is a model system for studying interactions between RNA and the RNA recognition motif (RRM), which is one of the most common RNA binding domains. We report here kinetic studies of dissociation of the U1A-SL2 RNA complex, using laser temperature jump and stopped-flow fluorescence methods with U1A proteins labeled with the intrinsic chromophore tryptophan. An analysis of the kinetic data suggests three phases of dissociation with time scales of ∼100 μs, ∼50 ms, and ∼2 s. We propose that the first step of dissociation is a fast rearrangement of the complex to form a loosely bound complex. The intermediate step is assigned to be the dissociation of the U1A-SL2 RNA complex, and the final step is assigned to a reorganization of the U1A protein structure into the conformation of the free protein. These assignments are consistent with previous proposals based on thermodynamic, NMR, and surface plasmon resonance experiments and molecular dynamics simulations. Together, these results begin to build a comprehensive model of the complex dynamic processes involved in the formation and dissociation of an RRM-RNA complex.

Drug resistant mechanism of diaryltriazine analog inhibitors of HIV-1 reverse transcriptase using molecular dynamics simulation and 3D-QSAR

Diaryltriazine inhibitors have highly potent and effective bioactivities for the wild type of HIV-1 reverse transcription. To design new drug of antimutant HIV-1 reverse transcriptase, the mechanism of drug resistance for four types of mutants was revealed. Molecular dynamics simulations suggest that Lys101, Leu100, Lys103, Tyr181, and Tyr188 are key residues. Different mutants of key residues may have different interaction modes and lead to different drug resistances. Then, CoMFA and CoMSIA methods were employed to construct 3D quantitative structure-activity relationship models. These models were evaluated by test set compounds. These models can be used to make quantitative prediction of their bioactivities for lead compounds before resorting to in vitro and in vivo experimentation.