nCoVDock2: a docking server to predict the binding modes between COVID-19 targets and its potential ligands.

The rapid emergence of SARS-CoV-2 variants with multi-sites mutations is considered as a major obstacle for the development of drugs and vaccines. Although most of the functional proteins essential for SARS-CoV-2 have been determined, the understanding of the COVID-19 target-ligand interactions remains a key challenge. The old version of this COVID-19 docking server was built in 2020, and free and open to all users. Here, we present nCoVDock2, a new docking server to predict the binding modes for targets from SARS-CoV-2. First, the new server supports more targets. We replaced the modeled structures with newly resolved structures and added more potential targets of COVID-19, especially for the variants. Second, for small molecule docking, Autodock Vina was upgraded to the latest version 1.2.0, and a new scoring function was added for peptide or antibody docking. Third, the input interface and molecular visualization were updated for a better user experience. The web server, together with an extensive help and tutorial, are freely available at: https://ncovdock2.schanglab.org.cn.

Tumor Homing Chimeric Peptide Rhomboids to Improve Photodynamic Performance by Inhibiting Therapy-Upregulated Cyclooxygenase-2.

Negative therapeutic feedback of inflammation would extensively attenuate the antitumor effect of photodynamic therapy (PDT). In this work, tumor homing chimeric peptide rhomboids (designated as NP-Mel) are fabricated to improve photodynamic performance by inhibiting PDT-upregulated cyclooxygenase-2 (COX-2). The hydrophobic photosensitizer of protoporphyrin IX (PpIX) and palmitic acid are conjugated onto the neuropilin receptors (NRPs) targeting peptide motif (CGNKRTR) to obtain tumor homing chimeric peptide (Palmitic-K(PpIX)CGNKRTR), which can encapsulate the COX-2 inhibitor of meloxicam. The well dispersed NP-Mel not only improves the drug stability and reactive oxygen species (ROS) production ability, but also increase the breast cancer targeted drug delivery to intensify the PDT effect. In vitro and in vivo studies verify that NP-Mel will decrease the secretion of prostaglandin E2 (PGE2) after PDT treatment, inducing the downregulation of IL-6 and TNF-α expressions to suppress PDT induced inflammation. Ultimately, an improved PDT performance of NP-Mel is achieved without inducing obvious systemic toxicity, which might inspire the development of sophisticated nanomedicine in consideration of the feedback induced therapeutic resistance.

Field and Laboratory Evidence That Chlorpyrifos Exposure Reduced the Population Density of a Freshwater Snail by Increasing Juvenile Mortality.

Pesticides have been frequently detected in global freshwater ecosystems, but attempts to document changes in population dynamics of organisms upon exposure to pesticides, establish a causal relationship between exposure and population effects, and identify the key toxic events within individuals under natural field conditions remain rare. Here, we used a field survey, a reciprocal cross-transplant experiment, and a laboratory toxicity experiment to build a compelling case that exposure to the insecticide chlorpyrifos was responsible for differences in snail (Bellamya aeruginosa) densities in eastern (ELL) and western basins of Liangzi Lake in China. Our field survey and reciprocal cross-transplant experiment revealed significant differences in snail densities, juvenile percentage, survival, and relative telomere length (RTL) in the two basins. The insecticide chlorpyrifos detected in snail tissues was negatively correlated with snail densities, the percentage of juvenile snails, and RTL and had an extremely high risk quotient in ELL. In the laboratory experiment, tissue concentrations of chlorpyrifos detected in ELL were associated with reduced RTL and increased juvenile mortality in B. aeruginosa. These results support the hypothesis that chlorpyrifos exposure in ELL reduced the density of snails by reducing juvenile survival and, consequently, recruitment to the adult population.

Carrier-Free Nano-PROTACs to Amplify Photodynamic Therapy Induced DNA Damage through BRD4 Degradation.

Therapy-induced DNA damage is the most common strategy to inhibit tumor cell proliferation, but the therapeutic efficacy is limited by DNA repair machinery. Carrier-free nanoproteolysis targeting chimeras (PROTACs), designed as SDNpros, have been developed to enhance photodynamic therapy (PDT) by blocking the DNA damage repair pathway through BRD4 degradation. Specifically, SDNpros are constructed through noncovalent interactions between the photosensitizer of chlorine e6 (Ce6) and PROTACs of BRD4 degrader (dBET57) via self-assembly. SDNpro has favorable dispersibility and a uniform nanosize distribution without drug excipients. Upon light irradiation, SDNpro produces abundant reactive oxygen species (ROS) to induce DNA oxidative damage. Meanwhile, the DNA repair pathway would be interrupted by the concurrent degradation of BRD4, which could intensify the oxidative DNA damage and elevate PDT efficiency. Beneficially, SDNpro suppresses tumor growth and avoids systemic side effects, providing a promising strategy to promote the clinical translation of PROTACs for tumor treatment.

CoDock-Ligand: combined template-based docking and CNN-based scoring in ligand binding prediction.

For ligand binding prediction, it is crucial for molecular docking programs to integrate template-based modeling with a precise scoring function. Here, we proposed the CoDock-Ligand docking method that combines template-based modeling and the GNINA scoring function, a Convolutional Neural Network-based scoring function, for the ligand binding prediction in CASP15. Among the 21 targets, we obtained successful predictions in top 5 submissions for 14 targets and partially successful predictions for 4 targets. In particular, for the most complicated target, H1114, which contains 56 metal cofactors and small molecules, our docking method successfully predicted the binding of most ligands. Analysis of the failed systems showed that the predicted receptor protein presented conformational changes in the backbone and side chains of the binding site residues, which may cause large structural deviations in the ligand binding prediction. In summary, our hybrid docking scheme was efficiently adapted to the ligand binding prediction challenges in CASP15.

Impact of AlphaFold on structure prediction of protein complexes: The CASP15-CAPRI experiment.

We present the results for CAPRI Round 54, the 5th joint CASP-CAPRI protein assembly prediction challenge. The Round offered 37 targets, including 14 homodimers, 3 homo-trimers, 13 heterodimers including 3 antibody-antigen complexes, and 7 large assemblies. On average ~70 CASP and CAPRI predictor groups, including more than 20 automatics servers, submitted models for each target. A total of 21 941 models submitted by these groups and by 15 CAPRI scorer groups were evaluated using the CAPRI model quality measures and the DockQ score consolidating these measures. The prediction performance was quantified by a weighted score based on the number of models of acceptable quality or higher submitted by each group among their five best models. Results show substantial progress achieved across a significant fraction of the 60+ participating groups. High-quality models were produced for about 40% of the targets compared to 8% two years earlier. This remarkable improvement is due to the wide use of the AlphaFold2 and AlphaFold2-Multimer software and the confidence metrics they provide. Notably, expanded sampling of candidate solutions by manipulating these deep learning inference engines, enriching multiple sequence alignments, or integration of advanced modeling tools, enabled top performing groups to exceed the performance of a standard AlphaFold2-Multimer version used as a yard stick. This notwithstanding, performance remained poor for complexes with antibodies and nanobodies, where evolutionary relationships between the binding partners are lacking, and for complexes featuring conformational flexibility, clearly indicating that the prediction of protein complexes remains a challenging problem.

Cascade Immune Activation of Self-Delivery Biomedicine for Photodynamic Immunotherapy Against Metastatic Tumor.

Photodynamic therapy (PDT) can generate reactive oxygen species (ROS) to cause cell apoptosis and induce immunogenic cell death (ICD) to activate immune response, becoming a promising antitumor modality. However, the overexpressions of indoleamine 2,3-dioxygenase (IDO) and programmed cell death ligand 1 (PD-L1) on tumor cells would reduce cytotoxic T cells infiltration and inhibit the immune activation. In this paper, a simple but effective nanosystem is developed to solve these issues for enhanced photodynamic immunotherapy. Specifically, it has been constructed a self-delivery biomedicine (CeNB) based on photosensitizer chlorine e6 (Ce6), IDO inhibitor (NLG919), and PD1/PDL1 blocker (BMS-1) without the need for extra excipients. Of note, CeNB possesses fairly high drug content (nearly 100%), favorable stability, and uniform morphology. More importantly, CeNB-mediated IDO inhibition and PD1/PDL1 blockade greatly improve the immunosuppressive tumor microenvironments to promote immune activation. The PDT of CeNB not only inhibits tumor proliferation but also induces ICD response to activate immunological cascade. Ultimately, self-delivery CeNB tremendously suppresses the tumor growth and metastasis while leads to a minimized side effect. Such simple and effective antitumor strategy overcomes the therapeutic resistance against PDT-initiated immunotherapy, suggesting a potential for metastatic tumor treatment in clinic.

In vitro testing and computational analysis of specific phytochemicals with antiviral activities considering their possible applications against COVID-19.

The purpose of this study was to investigate the reservoir of natural products against the SARS-CoV-2 virus and to identify suitable candidates in order to recommend appropriate phytotherapy. Adequately prepared 65 molecules from traditional Chinese medicine with proven antiviral properties were subjected to docking analysis using AutoDock Vina 4 software with the aim to investigate binding affinity and interactions of compounds with Mpro from the SARS-CoV-2 virus. Biflavonoids and tannins show best docking scores with -9,80 kcal/mol for biflavonoids and -9,00 kcal/mol for tannins. Biflavonoids: amentoflavone, agathistaflavone, robustaflavone, hinokiflavone and rhusflavanone were tested for their radical scavenging activity. Partition coefficients were examined by RP-HPLC. Evaluation of drug-likeness properties of investigated biflavonoids suggested rhusflavanone as a molecule with the best ADMET characteristics. Anti-inflammatory activity of rhusflavanone was investigated in LPS stimulated RAW264.7 macrophages. Tested biflavonoids exibit beneficial effects against inflammation by scavenging free radicals and by suppressing the production of proinflammatory mediators by macrophages. Both predictions of affinity spectra for substances (PASS) and in vitro testing showed promising biological activity of investigated biflavonoids. A Quantum chemical study was performed in order to calculate the thermodynamic, molecular orbital, and electrostatic potential of selected molecules and to compare their biological and chemical features. Our results highlighted antioxidant, anti-inflammatory and antiviral properties of investigated compounds, emphasizing the significance of biflavonoid moiety to selected characteristics, which encourage further investigational strategies against COVID-19.

Prediction of protein assemblies, the next frontier: The CASP14-CAPRI experiment.

We present the results for CAPRI Round 50, the fourth joint CASP-CAPRI protein assembly prediction challenge. The Round comprised a total of twelve targets, including six dimers, three trimers, and three higher-order oligomers. Four of these were easy targets, for which good structural templates were available either for the full assembly, or for the main interfaces (of the higher-order oligomers). Eight were difficult targets for which only distantly related templates were found for the individual subunits. Twenty-five CAPRI groups including eight automatic servers submitted ~1250 models per target. Twenty groups including six servers participated in the CAPRI scoring challenge submitted ~190 models per target. The accuracy of the predicted models was evaluated using the classical CAPRI criteria. The prediction performance was measured by a weighted scoring scheme that takes into account the number of models of acceptable quality or higher submitted by each group as part of their five top-ranking models. Compared to the previous CASP-CAPRI challenge, top performing groups submitted such models for a larger fraction (70-75%) of the targets in this Round, but fewer of these models were of high accuracy. Scorer groups achieved stronger performance with more groups submitting correct models for 70-80% of the targets or achieving high accuracy predictions. Servers performed less well in general, except for the MDOCKPP and LZERD servers, who performed on par with human groups. In addition to these results, major advances in methodology are discussed, providing an informative overview of where the prediction of protein assemblies currently stands.

Carrier Free Photodynamic Synergists for Oxidative Damage Amplified Tumor Therapy.

Tumor cells adapt to excessive oxidative stress by actuating reactive oxygen species (ROS)-defensing system, leading to a resistance to oxidation therapy. In this work, self-delivery photodynamic synergists (designated as PhotoSyn) are developed for oxidative damage amplified tumor therapy. Specifically, PhotoSyn are fabricated by the self-assembly of chlorine e6 (Ce6) and TH588 through π-π stacking and hydrophobic interactions. Without additional carriers, nanoscale PhotoSyn possess an extremely high drug loading rate (up to 100%) and they are found to be fairly stable in aqueous phase with a uniform size distribution. Intravenously injected PhotoSyn prefer to accumulate at tumor sites for effective cellular uptake. More importantly, TH588-mediated MTH1 inhibition could destroy the ROS-defensing system of tumor cells by preventing the elimination of 8-oxo-2'-deoxyguanosine triphosphate (8-oxo-dG), thereby exacerbating the oxidative DNA damage induced by the photodynamic therapy (PDT) of Ce6 under light irradiation. As a consequence, PhotoSyn exhibit enhanced photo toxicity and a significant antitumor effect. This amplified oxidative damage strategy improves the PDT efficiency with a reduced side effect by increasing the lethality of ROS without generating superabundant ROS, which would provide a new insight for developing self-delivery nanoplatforms in photodynamic tumor therapy in clinic.

COVID-19 Docking Server: a meta server for docking small molecules, peptides and antibodies against potential targets of COVID-19.

MOTIVATION: The coronavirus disease 2019 (COVID-19) caused by a new type of coronavirus has been emerging from China and led to thousands of death globally since December 2019. Despite many groups have engaged in studying the newly emerged virus and searching for the treatment of COVID-19, the understanding of the COVID-19 target-ligand interactions represents a key challenge. Herein, we introduce COVID-19 Docking Server, a web server that predicts the binding modes between COVID-19 targets and the ligands including small molecules, peptides and antibodies. RESULTS: Structures of proteins involved in the virus life cycle were collected or constructed based on the homologs of coronavirus, and prepared ready for docking. The meta-platform provides a free and interactive tool for the prediction of COVID-19 target-ligand interactions and following drug discovery for COVID-19. AVAILABILITY AND IMPLEMENTATION: http://ncov.schanglab.org.cn. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Novel STAT3 small-molecule inhibitors identified by structure-based virtual ligand screening incorporating SH2 domain flexibility.

Efforts to develop STAT3 inhibitors have focused on its SH2 domain starting with short phosphotyrosylated peptides based on STAT3 binding motifs, e.g. pY905LPQTV within gp130. Despite binding to STAT3 with high affinity, issues regarding stability, bioavailability, and membrane permeability of these peptides, as well as peptidomimetics such as CJ-887, have limited their further clinical development and led to interest in small-molecule inhibitors. Some small molecule STAT3 inhibitors, identified using structure-based virtual ligand screening (SB-VLS); while having favorable drug-like properties, suffer from weak binding affinities, possibly due to the high flexibility of the target domain. We conducted molecular dynamic (MD) simulations of the SH2 domain in complex with CJ-887, and used an averaged structure from this MD trajectory as an "induced-active site" receptor model for SB-VLS of 110,000 compounds within the SPEC database. Screening was followed by re-docking and re-scoring of the top 30% of hits, selection for hit compounds that directly interact with pY + 0 binding pocket residues R609 and S613, and testing for STAT3 targeting in vitro, which identified two lead hits with good activity and favorable drug-like properties. Unlike most small-molecule STAT3 inhibitors previously identified, which contain negatively-charged moieties that mediate binding to the pY + 0 binding pocket, these compounds are uncharged and likely will serve as better candidates for anti-STAT3 drug development. IMPLICATIONS: SB-VLS, using an averaged structure from molecular dynamics (MD) simulations of STAT3 SH2 domain in a complex with CJ-887, a known peptidomimetic binder, identify two highly potent, neutral, low-molecular weight STAT3-inhibitors with favorable drug-like properties.

Tris(1,3-dichloro-2-propyl)phosphate Reduces Growth Hormone Expression via Binding to Growth Hormone Releasing Hormone Receptors and Inhibits the Growth of Crucian Carp.

Tris(1,3-dichloro-2-propyl)phosphate (TDCIPP) has commonly been used as an additive flame retardant and frequently detected in the aquatic environment and in biological samples worldwide. Recently, it was found that exposure to TDCIPP inhibited the growth of zebrafish, but the relevant molecular mechanisms remained unclear. In this study, 5 day-old crucian carp (Carassius auratus) larvae were treated with 0.5, 5, or 50 µg/L TDCIPP for 90 days; the effect on growth was evaluated; and related molecular mechanisms were explored. Results demonstrated that 5 or 50 µg/L TDCIPP treatment significantly inhibited the growth of crucian carp and downregulated the expression of growth hormones (ghs), growth hormone receptor (ghr), and insulin-like growth factor 1 (igf1). Molecular docking, dual-luciferase reporter gene assay, and in vitro experiments demonstrated that TDCIPP could bind to the growth hormone releasing hormone receptor protein of crucian carp and disturb the stimulation of growth hormone releasing hormone to the expression of ghs, resulting in the decrease of the mRNA level of gh1 and gh2 in pituitary cells. Our findings provide new perceptions into the molecular mechanisms of developmental toxicity of TDCIPP in fish.

Template-based modeling and ab-initio docking using CoDock in CAPRI.

Integration of template-based modeling, global sampling and precise scoring is crucial for the development of molecular docking programs with improved accuracy. We combined template-based modeling and ab-initio docking protocol as hybrid docking strategy called CoDock for the docking and scoring experiments of the seventh CAPRI edition. For CAPRI rounds 38-45, we obtained acceptable or better models in the top 10 submissions for eight out of the 16 evaluated targets as predictors, nine out of the 16 targets as scorers. Especially, we submitted acceptable models for all of the evaluated protein-oligosaccharide targets. For the CASP13-CAPRI experiment (round 46), we obtained acceptable or better models in the top 5 submissions for 10 out of the 20 evaluated targets as predictors, 11 out of the 20 targets as scorers. The failed cases for our group were mainly the difficult targets and the protein-peptide systems in CAPRI and CASP13-CAPRI experiments. In summary, this CAPRI edition showed that our hybrid docking strategy can be efficiently adapted to the increasing variety of challenges in the field of molecular interactions.

Effects of triphenyl phosphate on ciliate protozoa Tetrahymena thermophila following acute exposure and sub-chronic exposure.

Triphenyl phosphate (TPHP) is one of the most widely used organophosphate flame retardants (OPFRs) and is frequently detected in a variety of environmental media. Previous studies reported that TPHP had toxic effects on vertebrates, but little toxic information was available in lower trophic aquatic organisms which were more sensitive to the exposure of many toxic substances. In this study, protozoa Tetrahymena thermophila (T. thermophila) were exposed to 0, 0.01, 0.17 or 2.35 mg/L TPHP for 5 days to study the effects of sub-chronic exposure on theoretical population, cell viability, cell size and number of cilia. Additionally, the effects of TPHP on gene transcription were assessed by transcriptome sequencing technique (RNA-Seq). Cell viability and number of cilia were significantly reduced in all TPHP exposure groups compared with the control. In addition, exposure to 0.17 or 2.35 mg/L TPHP significantly reduced the theoretical population, circumference and body width, and there was a significant decrease in body length in the 2.35 mg/L exposure group. Comparative transcriptome sequencing identified a total of 4105 up- and 4487 down-regulated genes after exposure to 2.35 mg/L TPHP for 5 days compared with the control. KEGG analysis showed that dysfunction of pathways associated with ribosome, spliceosome, phagosome, proteasome and protein processing in endoplasmic reticulum in this study might be responsible for the toxicity of T. thermophila caused by TPHP. In general, the results indicated that TPHP had an adverse effect on the protozoa T. thermophila.

Blind prediction of homo- and hetero-protein complexes: The CASP13-CAPRI experiment.

We present the results for CAPRI Round 46, the third joint CASP-CAPRI protein assembly prediction challenge. The Round comprised a total of 20 targets including 14 homo-oligomers and 6 heterocomplexes. Eight of the homo-oligomer targets and one heterodimer comprised proteins that could be readily modeled using templates from the Protein Data Bank, often available for the full assembly. The remaining 11 targets comprised 5 homodimers, 3 heterodimers, and two higher-order assemblies. These were more difficult to model, as their prediction mainly involved "ab-initio" docking of subunit models derived from distantly related templates. A total of ~30 CAPRI groups, including 9 automatic servers, submitted on average ~2000 models per target. About 17 groups participated in the CAPRI scoring rounds, offered for most targets, submitting ~170 models per target. The prediction performance, measured by the fraction of models of acceptable quality or higher submitted across all predictors groups, was very good to excellent for the nine easy targets. Poorer performance was achieved by predictors for the 11 difficult targets, with medium and high quality models submitted for only 3 of these targets. A similar performance "gap" was displayed by scorer groups, highlighting yet again the unmet challenge of modeling the conformational changes of the protein components that occur upon binding or that must be accounted for in template-based modeling. Our analysis also indicates that residues in binding interfaces were less well predicted in this set of targets than in previous Rounds, providing useful insights for directions of future improvements.

Design, synthesis of novel purin-6-one derivatives as phosphodiesterase 2 (PDE2) inhibitors: The neuroprotective and anxiolytic-like effects.

Phosphodiesterase 2 (PDE2) has received much attention for the potential treatment of the central nervous system (CNS) disorders. Herein, based on the existing PDE2 inhibitors and their binding modes, a series of purin-6-one derivatives were designed, synthesized and evaluated for PDE2 inhibitory activities, which led to the discovery of the best compounds 6p and 6s with significant inhibitory potency (IC50: 72 and 81 nM, respectively). Docking simulation was performed to insert compound 6s into the crystal structure of PDE2 at the active site to determine the binding mode. Furthermore, compound 6s significantly protected HT-22 cells against corticosterone-induced cytotoxicity and rescued corticosterone-induced decreases in cAMP and cGMP levels. It also produced anxiolytic-like effect in the elevated plus-maze test and exhibited favorable pharmacokinetic properties in vivo. These results might bring significant instruction for further development of potent PDE2 inhibitors.

CoDockPP: A Multistage Approach for Global and Site-Specific Protein-Protein Docking.

Protein-protein docking technology is an effective approach to study the molecular mechanism of essential biological processes mediated by complex protein-protein interactions. The fast Fourier transform (FFT) correlation approach makes a good balance between the exhaustive global sampling and the computational efficiency for protein-protein docking. However, it is difficult to integrate the precise knowledge-based scoring function and site constraint information into the FFT-based approach. New docking strategies with the capability of combining both global sampling and precise scoring are strongly needed. We propose a multistage protein-protein docking strategy called CoDockPP. This program takes full advantage of the sampling efficiency of the FFT-based method to choose the valid ligand protein poses with good surface complementarity. The retained poses are transformed to the real Cartesian space for the implementation of site constraints and atomic scoring. Site constraints and a rapid table lookup scoring are applied to gradually reduce the candidate poses to a tractable number. To enhance the accuracy of docking prediction, the best fast-scoring states are expanded the local sampling points and then these neighbor poses are further evaluated by the precise knowledge-based scoring function. By testing on protein-protein docking benchmark 5.0, CoDockPP remarkably improves the success rate and hit count in both ab initio docking and site-specific docking, especially in difficult cases. The server is free and open to all users with no login requirement at http://codockpp.schanglab.org.cn .

Molecular Dynamics Simulations of Wild Type and Mutants of SAPAP in Complexed with Shank3.

Specific interactions between scaffold protein SH3 and multiple ankyrin repeat domains protein 3 (Shank3) and synapse-associated protein 90/postsynaptic density-95⁻associated protein (SAPAP) are essential for excitatory synapse development and plasticity. In a bunch of human neurological diseases, mutations on Shank3 or SAPAP are detected. To investigate the dynamical and thermodynamic properties of the specific binding between the N-terminal extended PDZ (Post-synaptic density-95/Discs large/Zonaoccludens-1) domain (N-PDZ) of Shank3 and the extended PDZ binding motif (E-PBM) of SAPAP, molecular dynamics simulation approaches were used to study the complex of N-PDZ with wild type and mutated E-PBM peptides. To compare with the experimental data, 974QTRL977 and 966IEIYI970 of E-PBM peptide were mutated to prolines to obtain the M4P and M5P system, respectively. Conformational analysis shows that the canonical PDZ domain is stable while the ßN extension presents high flexibility in all systems, especially for M5P. The high flexibility of ßN extension seems to set up a barrier for the non-specific binding in this area and provide the basis for specific molecular recognition between Shank3 and SAPAP. The wild type E-PBM tightly binds to N-PDZ during the simulation while loss of binding is observed in different segments of the mutated E-PBM peptides. Energy decomposition and hydrogen bonds analysis show that M4P mutations only disrupt the interactions with canonical PDZ domain, but the interactions with ßN1' remain. In M5P system, although the interactions with ßN1' are abolished, the binding between peptide and the canonical PDZ domain is not affected. The results indicate that the interactions in the two-binding site, the canonical PDZ domain and the ßN1' extension, contribute to the binding between E-PBM and N-PDZ independently. The binding free energies calculated by MM/GBSA (Molecular Mechanics/Generalized Born Surface Area) are in agreement with the experimental binding affinities. Most of the residues on E-PBM contribute considerably favorable energies to the binding except A963 and D964 in the N-terminal. The study provides information to understand the molecular basis of specific binding between Shank3 and SAPAP, as well as clues for design of peptide inhibitors.

Exploring the molecular basis of RNA recognition by the dimeric RNA-binding protein via molecular simulation methods.

RNA-binding protein with multiple splicing (RBPMS) is critical for axon guidance, smooth muscle plasticity, and regulation of cancer cell proliferation and migration. Recently, different states of the RNA-recognition motif (RRM) of RBPMS, one in its free form and another in complex with CAC-containing RNA, were determined by X-ray crystallography. In this article, the free RRM domain, its wild type complex and 2 mutant complex systems are studied by molecular dynamics (MD) simulations. Through comparison of free RRM domain and complex systems, it's found that the RNA binding facilitates stabilizing the RNA-binding interface of RRM domain, especially the C-terminal loop. Although both R38Q and T103A/K104A mutations reduce the binding affinity of RRM domain and RNA, the underlining mechanisms are different. Principal component analysis (PCA) and Molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) methods were used to explore the dynamical and recognition mechanisms of RRM domain and RNA. R38Q mutation is positioned on the homodimerization interface and mainly induces the large fluctuations of RRM domains. This mutation does not directly act on the RNA-binding interface, but some interfacial hydrogen bonds are weakened. In contrast, T103A/K104A mutations are located on the RNA-binding interface of RRM domain. These mutations obviously break most of high occupancy hydrogen bonds in the RNA-binding interface. Meanwhile, the key interfacial residues lose their favorable energy contributions upon RNA binding. The ranking of calculated binding energies in 3 complex systems is well consistent with that of experimental binding affinities. These results will be helpful in understanding the RNA recognition mechanisms of RRM domain.