A combinatorial scoring function for protein-RNA docking.

Protein-RNA docking is still an open question. One of the main challenges is to develop an effective scoring function that can discriminate near-native structures from the incorrect ones. To solve the problem, we have constructed a knowledge-based residue-nucleotide pairwise potential with secondary structure information considered for nonribosomal protein-RNA docking. Here we developed a weighted combined scoring function RpveScore that consists of the pairwise potential and six physics-based energy terms. The weights were optimized using the multiple linear regression method by fitting the scoring function to L_rmsd for the bound docking decoys from Benchmark II. The scoring functions were tested on 35 unbound docking cases. The results show that the scoring function RpveScore including all terms performs best. Also RpveScore was compared with the statistical mechanics-based method derived potential ITScore-PR, and the united atom-based statistical potentials QUASI-RNP and DARS-RNP. The success rate of RpveScore is 71.6% for the top 1000 structures and the number of cases where a near-native structure is ranked in top 30 is 25 out of 35 cases. For 32 systems (91.4%), RpveScore can find the binding mode in top 5 that has no lower than 50% native interface residues on protein and nucleotides on RNA. Additionally, it was found that the long-range electrostatic attractive energy plays an important role in distinguishing near-native structures from the incorrect ones. This work can be helpful for the development of protein-RNA docking methods and for the understanding of protein-RNA interactions. RpveScore program is available to the public at http://life.bjut.edu.cn/kxyj/kycg/2017116/14845362285362368_1.html Proteins 2017; 85:741-752. © 2016 Wiley Periodicals, Inc.

Allosteric transitions of ATP-binding cassette transporter MsbA studied by the adaptive anisotropic network model.

The transporter MsbA is a kind of multidrug resistance ATP-binding cassette transporter that can transport lipid A, lipopolysaccharides, and some amphipathic drugs from the cytoplasmic to the periplasmic side of the inner membrane. In this work, we explored the allosteric pathway of MsbA from the inward- to outward-facing states during the substrate transport process with the adaptive anisotropic network model. The results suggest that the allosteric transitions proceed in a coupled way. The large-scale closing motions of the nucleotide-binding domains occur first, accompanied with a twisting motion at the same time, which becomes more obvious in middle and later stages, especially for the later. This twisting motion plays an important role for the rearrangement of transmembrane helices and the opening of transmembrane domains on the periplasmic side that mainly take place in middle and later stages respectively. The topological structure plays an important role in the motion correlations above. The conformational changes of nucleotide-binding domains are propagated to the transmembrane domains via the intracellular helices IH1 and IH2. Additionally, the movement of the transmembrane domains proceeds in a nonrigid body, and the two monomers move in a symmetrical way, which is consistent with the symmetrical structure of MsbA. These results are helpful for understanding the transport mechanism of the ATP-binding cassette exporters.

Allosteric transitions of the maltose transporter studied by an elastic network model.

The maltose transporter from Escherichia coli is one of the ATP-binding cassette (ABC) transporters that utilize the energy from ATP hydrolysis to translocate substrates across cellular membranes. Until 2011, three crystal structures have been determined for maltose transporter at different states in the process of transportation. Here, based on these crystal structures, the allosteric pathway from the resting state (inward-facing) to the catalytic intermediate state (outward-facing) is studied by applying an adaptive anisotropic network model. The results suggest that the allosteric transitions proceed in a coupled way. The closing of the nucleotide-binding domains occurs first, and subsequently this conformational change is propagated to the transmembrane domains (TMD) via the EAA and EAS loops, and then to the maltose-binding protein, which facilitates the translocation of the maltose. It is also found that there exist nonrigid-body and asymmetric movements in the TMD. The cytoplasmic gate may only play the role of allosteric propagation during the transition from the pretranslocation to outward-facing states. In addition, the results show that the movment of the helical subdomain towards the RecA-like subdomain mainly occurs in the earlier stages of the transition. These results can provide some insights into the understanding of the mechanism of ABC transporters.

Detection of persistent organic pollutants binding modes with androgen receptor ligand binding domain by docking and molecular dynamics.

BACKGROUND: Persistent organic pollutants (POPs) are persistent in the environment after release from industrial compounds, combustion productions or pesticides. The exposure of POPs has been related to various reproductive disturbances, such as reduced semen quality, testicular cancer, and imbalanced sex ratio. Among POPs, dichlorodiphenyldichloroethylene (4,4'-DDE) and polychlorinated biphenyls (PCBs) are the most widespread and well-studied compounds. Recent studies have revealed that 4,4'-DDE is an antagonist of androgen receptor (AR). However, the mechanism of the inhibition remains elusive. CB-153 is the most common congener of PCBs, while the action of CB-153 on AR is still under debate. RESULTS: Molecular docking and molecular dynamics (MD) approaches have been employed to study binding modes and inhibition mechanism of 4,4'-DDE and CB-153 against AR ligand binding domain (LBD). Several potential binding sites have been detected and analyzed. One possible binding site is the same binding site of AR natural ligand androgen 5α-dihydrotestosterone (DHT). Another one is on the ligand-dependent transcriptional activation function (AF2) region, which is crucial for the co-activators recruitment. Besides, a novel possible binding site was observed for POPs with low binding free energy with the receptor. Detailed interactions between ligands and the receptor have been represented. The disrupting mechanism of POPs against AR has also been discussed. CONCLUSIONS: POPs disrupt the function of AR through binding to three possible biding sites on AR/LBD. One of them shares the same binding site of natural ligand of AR. Another one is on AF2 region. The third one is in a cleft near N-terminal of the receptor. Significantly, values of binding free energy of POPs with AR/LBD are comparable to that of natural ligand androgen DHT.

Identification of functionally key residues in AMPA receptor with a thermodynamic method.

AMPA receptor mediates the fast excitatory synaptic transmission in the central nervous system, and it is activated by the binding of glutamate that results in the opening of the transmembrane ion channel. In the present work, the thermodynamic method developed by our group was improved and then applied to identify the functionally key residues that regulate the glutamate-binding affinity of AMPA receptor. In our method, the key residues are identified as those whose perturbation largely changes the ligand binding free energy of the protein. It is found that besides the ligand binding sites, other residues distant from the binding cleft can also influence the glutamate binding affinity through a long-range allosteric regulation. These allosteric sites include the hinge region of the ligand binding cleft, the dimer interface of the ligand binding domain, the linkers between the ligand binding domain and the transmembrane domain, and the interface between the N-terminal domain and the ligand binding domain. Our calculation results are consistent with the available experimental data. The results are helpful for our understanding of the mechanism of long-range allosteric communication in the AMPA receptor and the mechanism of channel opening triggered by glutamate binding.

The current status and challenges in the development of fusion inhibitors as therapeutics for HIV-1 infection.

HIV-1 membrane fusion as a part of the process of viral entry in the target cells is facilitated by gp41 and gp120, which are encoded by Env gene of HIV-1. Based on the structure and the mechanism researches, new treatment options targeting HIV-1 entry process have been proposed. Enfuvirtide, which mimics amino acid sequences of viral envelope glycoprotein gp41, is the first HIV-1 fusion inhibitor approved by FDA. Although it fulfills vital functions by binding to gp41 and abolishing the membrane fusion reaction when used in combination, it could induce drug resistant virus variants. Currently, a number of design and modification schemes have been presented, a large number of prospective fusion peptides have emerged. For these fusion inhibitors, multiple mutations in gp41 have been associated with the loss of susceptibility to agents. This review reported the current developments and innovative designs of HIV-1 membrane fusion inhibitors.

The interactions and recognition of cyclic peptide mimetics of Tat with HIV-1 TAR RNA: a molecular dynamics simulation study.

The interaction of HIV-1 trans-activator protein Tat with its cognate trans-activation response element (TAR) RNA is critical for viral transcription and replication. Therefore, it has long been considered as an attractive target for the development of antiviral compounds. Recently, the conformationally constrained cyclic peptide mimetics of Tat have been tested to be a promising family of lead peptides. Here, we focused on two representative cyclic peptides termed as L-22 and KP-Z-41, both of which exhibit excellent inhibitory potency against Tat and TAR interaction. By means of molecular dynamics simulations, we obtained a detailed picture of the interactions between them and HIV-1 TAR RNA. In results, it is found that the binding modes of the two cyclic peptides to TAR RNA are almost identical at or near the bulge regions, whereas the binding interfaces at the apical loop exhibit large conformational heterogeneity. In addition, it is revealed that electrostatic interaction energy contributes much more to KP-Z-41 complex formation than to L-22 complex, which is the main source of energy that results in a higher binding affinity of KP-Z-41 over-22 for TAR RNA. Furthermore, we identified a conserved motif RRK (Arg-Arg-Lys) that is shown to be essential for specific binding of this class of cyclic peptides to TAR RNA. This work can provide a useful insight into the design and modification of cyclic peptide inhibitors targeting the association of HIV-1 Tat and TAR RNA.

A new residue-nucleotide propensity potential with structural information considered for discriminating protein-RNA docking decoys.

Understanding the key factors that influence the preferences of residue-nucleotide interactions in specific protein-RNA interactions has remained a research focus. We propose an effective approach to derive residue-nucleotide propensity potentials through considering both the types of residues and nucleotides, and secondary structure information of proteins and RNAs from the currently largest nonredundant and nonribosomal protein-RNA interaction database. To test the validity of the potentials, we used them to select near-native structures from protein-RNA docking poses. The results show that considering secondary structure information, especially for RNAs, greatly improves the predictive power of pair potentials. The success rate is raised from 50.7 to 65.5% for the top 2000 structures, and the number of cases in which a near-native structure is ranked in top 50 is increased from 7 to 13 out of 17 cases. Furthermore, the exclusion of ribosomes from the database contributes 8.3% to the success rate. In addition, some very interesting findings follow: (i) the protein secondary structure element π-helix is strongly associated with RNA-binding sites; (ii) the nucleotide uracil occurs frequently in the most preferred pairs in which the unpaired and non-Watson-Crick paired uracils are predominant, which is probably significant in evolution. The new residue-nucleotide potentials can be helpful for the progress of protein-RNA docking methods, and for understanding the mechanisms of protein-RNA interactions.

Identification of key residues for protein conformational transition using elastic network model.

Proteins usually undergo conformational transitions between structurally disparate states to fulfill their functions. The large-scale allosteric conformational transitions are believed to involve some key residues that mediate the conformational movements between different regions of the protein. In the present work, a thermodynamic method based on the elastic network model is proposed to predict the key residues involved in protein conformational transitions. In our method, the key functional sites are identified as the residues whose perturbations largely influence the free energy difference between the protein states before and after transition. Two proteins, nucleotide binding domain of the heat shock protein 70 and human/rat DNA polymerase ß, are used as case studies to identify the critical residues responsible for their open-closed conformational transitions. The results show that the functionally important residues mainly locate at the following regions for these two proteins: (1) the bridging point at the interface between the subdomains that control the opening and closure of the binding cleft; (2) the hinge region between different subdomains, which mediates the cooperative motions between the corresponding subdomains; and (3) the substrate binding sites. The similarity in the positions of the key residues for these two proteins may indicate a common mechanism in their conformational transitions.

Insight into the inhibitory mechanism and binding mode between D77 and HIV-1 integrase by molecular modeling methods.

Integrase is an essential enzyme in the life cycle of Human immunoficiency virus type 1 (HIV-1) and also an important target for designing integrase inhibitors. In this paper, the binding modes between the wild type integrase core domain (ICD) and the W131A mutant ICD with the benzoic acid derivative--D77 were investigated using the molecular docking combined with molecular dynamics (MD) simulations. The result of MD simulations showed that the W131A substitution affected the flexibility of the region 150-167 in both the monomer A and B of the mutant type ICD. In principle, D77 interacted with the residues around the Lens Epithelium-Derived Growth Factor (LEDGF/p75) binding site which is nearby the HIV-1 integrase dimer interface. However, the specific binding modes for D77-wild type integrase and D77-mutant integrase systems are various. According to the binding mode of D77 with the wild type ICD, D77 can effectively intervene with the binding of LEDGF/p75 to integrase due to a steric hindrance effect around the LEDGF/p75 binding site. In addition, we found that D77 might also affect its inhibitory action by reducing the flexibility of the region 150-167 of integrase. Through energy decomposition calculated with the Molecular Mechanics Generalized Born Surface Area approach to estimate the binding affinity, it seems likely that W131 and E170 are indispensable for the ligand binding, as characterized by the largest binding affinity. All the above results are consistent with the experimental data, providing us with some helpful information not only for the understanding of the mechanism of this kind of inhibitor but also for the rational drug design.

An analysis of the influence of protein intrinsic dynamical properties on its thermal unfolding behavior.

The influence of the protein topology-encoded dynamical properties on its thermal unfolding motions was studied in the present work. The intrinsic dynamics of protein topology was obtained by the anisotropic network model (ANM). The ANM has been largely used to investigate protein collective functional motions, but it is not well elucidated if this model can also reveal the preferred large-scale motions during protein unfolding. A small protein barnase is used as a typical case study to explore the relationship between protein topology-encoded dynamics and its unfolding motions. Three thermal unfolding simulations at 500 K were performed for barnase and the entire unfolding trajectories were sampled and partitioned into several windows. For each window, the preferred unfolding motions were investigated by essential dynamics analysis, and then associated with the intrinsic dynamical properties of the starting conformation in this window, which is detected by ANM. The results show that only a few slow normal modes imposed by protein structure are sufficient to give a significant overlap with the preferred unfolding motions. Especially, the large amplitude unfolding movements, which imply that the protein jumps out of a local energy basin, can be well described by a single or several ANM slow modes. Besides the global motions, it is also found that the local residual fluctuations encoded in protein structure are highly correlated with those in the protein unfolding process. Furthermore, we also investigated the relationship between protein intrinsic flexibility and its unfolding events. The results show that the intrinsic flexible regions tend to unfold early. Several early unfolding events can be predicted by analysis of protein structural flexibility. These results imply that protein structure-encoded dynamical properties have significant influences on protein unfolding motions.

Computer-aided design, synthesis, and biological activity evaluation of potent fusion inhibitors targeting HIV-1 gp41.

This discovered and optimized several novel HIV-1 fusion inhibitors and further evaluated the inhibitory activities of these compounds in vitro. Here, we have reported the computer-aided design, synthesis, and biological evaluation of a series of small molecule fusion inhibitors targeting HIV-1 gp41. Based on the structure of inhibitor (NB2), we carried out de novo design and screened out a series of novel structure molecules by using Leapfrog and Autodock programs. Our structure-based modification obtained a potent fusion inhibitor (IC50 = 41.1 µg/mL). Several novel compounds were discovered as fusion inhibitors, which suggested that our design methodology is reliable, paving the way for de novo design of novel small-molecule HIV inhibitors targeting gp41.

[Effect of pre-electroacupuncture at "Zusanli" (ST 36) on DA and 5-HT contents and their ratio in hypothalamus and striatum in exercise rats].

OBJECTIVE: To observe anti-fatigue effect and mechanisms of pre-electroacupuncture (EA) at "Zusanli" (ST 36) in rats undergoing acute treadmill running. METHODS: Fifty male SD rats were randomly divided into three groups: a quiet group (group Q, n = 10), a model group (group M, n = 20) and an EA preconditioning group (group EAP, n = 20). After adaptation for undergoing treadmill running, all the rats in group M and group EAP were trained on acute treadmill running. Besides, EA with continuous waves, 2 Hz in frequency and 2 mA in intensity was applied at bilateral "Zusanli" (ST 36) for 30 min, which was applied once daily for continuous 6 days before treadmill running for the rats in Group EAP. Plasma lactate contents were measured immediately and 3 hours after treadmill running, respectively. Changes of dopamine (DA) and serotonin (5-HT) contents obtained immediately and 3 hours after treadmill running, respectively, in hypothalamus and striatum, were detected and compared, and DA/5-HT ratios were calculated. RESULTS: Compared with group Q, the levels of blood lactate and hypothalamic 5-HT tented to increase in rats of group M, and the contents of hypothalamic DA increased significantly (P < 0.01), while the contents of striatal DA and 5-HT in group M decreased significantly (both P < 0.01) at 3 h after treadmill running. Immediately after treadmill running, the contents of DA and 5-HT increased significantly in hypothalamus (both P < 0.01), but decreased significantly in striatum (both P < 0.01) in group EAP, compared with those in group M. Moreover, EA pretreatment markedly decreased the levels of blood lactate (P < 0.05) and hypothalamic 5-HT (P < 0.01), and obviously elevated the ratio of DA/5-HT in the hypothalamus (P < 0.01) at 3 h after treadmill running. CONCLUSION: Preventive EA at "Zusanli" (ST 36) can accelerate recovery from fatigue, which may be related to its reducing accumulation of blood lactate, elevating DA/ 5-HT ratio in the hypothalamus of the rats undergoing treadmill running.

Role of electrostatic interactions for the stability and folding behavior of cold shock protein.

The impacts of three charged-residue-involved mutations, E46A, R3E, and R3E/L66E, on the thermostability and folding behavior of the cold shock protein from the themophile Bacillus caldolyticus (Bc-Csp) were investigated by using a modified Go-like model, in which the nonspecific electrostatic interactions of charged residues were taken into account. Our simulation results show that the wild-type Bc-Csp and its three mutants are all two-sate folders, which is consistent with the experimental observations. It is found that these three mutations all lead to a decrease of protein thermodynamical stability, and the effect of R3E mutation is the strongest. The lower stability of these three mutants is due to the increase of the enthalpy of the folded state and the entropy of the unfolded state. Using this model, we also studied the folding kinetics and the folding/unfolding pathway of the wild-type Bc-Csp as well as its three mutants and then discussed the effects of electrostatic interactions on the folding kinetics. The results indicate that the substitutions at positions 3 and 46 largely decrease the folding kinetics, whereas the mutation of residue 66 only slightly decreases the folding rate. This result agrees well with the experimental observations. It is also found that these mutations have little effects on the folding transition state and the folding pathway, in which the N-terminal beta sheet folds earlier than the C-terminal region. We also investigated the detailed unfolding pathway and found that it is really the reverse of the folding pathway, providing the validity of our simulation results.

Design, synthesis and anti-HIV integrase evaluation of N-(5-chloro-8-hydroxy-2-styrylquinolin-7-yl)benzenesulfonamide derivatives.

Styrylquinoline derivatives are demonstrated to be HIV-1 integrase inhibitors. On the basis of our previous CoMFA analysis of a series of styrylquinoline derivatives, N-[(2-substituted-styryl)-5-chloro-8-hydroxyquinolin-7-yl]-benzenesulfonamide derivatives were designed and synthesized,and their possible HIV IN inhibitory activity was evaluated.

Study on the binding mode of the integrase with DNA via steered molecular dynamics simulation.

Human immunodeficiency virus type 1 (HIV-1) integrase (IN) is an essential enzyme in the lifecycle of this virus and also an important target for the study of anti-HIV drugs. In the current work, a model for the active site of IN and viral DNA was built by combining experimental data with the results of steered molecular dynamics simulation. The model was then taken into a series automatic molecular docking calculations with two groups of inhibitors. According to the results of molecular docking, the inhibitors of the second group share a similar binding model with those of the first group, though they have no common scaffold. The newly built model of the IN-DNA complex is helpful for our subsequent research on the design of IN inhibitors.

Study on the inhibitory mechanism and binding mode of the hydroxycoumarin compound NSC158393 to HIV-1 integrase by molecular modeling.

Human immunodeficiency virus type 1 integrase (IN) is an essential enzyme in the life cycle of this virus and also an important target for the study of anti-HIV drugs. In this work, the binding modes of the wild type IN core domain and the two mutants, that is, W132G and C130S, with the 4-hydroxycoumarin compound NSC158393 were evaluated by using the "relaxed complex" molecular docking approach combined with molecular dynamics (MD) simulations. Based on the monomer MD simulations, both of the two substitutions affect not only the stability of the 128-136 peptides, but also the flexibility of the functional 140s loop. In principle, NSC158393 binds the 128-136 peptides of IN; however, the specific binding modes for the three systems are various. According to the binding mode of NSC158393 with WT, NSC158393 can effectively interfere with the stability of the IN dimer by causing a steric hindrance around the monomer interface. Additionally, through the comparative analysis of the MD trajectories of the wild type IN and the IN-NSC158393 complex, we found that NSC15893 may also exert its inhibitory function by diminishing the mobility of the function loop of IN. Three key binding residues, that is, W131, K136, and G134, were discovered by energy decomposition calculated with the Molecular Mechanics Generalized Born Surface Area method. Characterized by the largest binding affinity, W131 is likely to be indispensable for the ligand binding. All the above results are consistent with experiment data, providing us some helpful information for understanding the mechanism of the coumarin-based inhibitors.

Molecular dynamics simulations of the bacterial periplasmic heme binding proteins ShuT and PhuT.

ShuT and PhuT are two periplasmic heme binding proteins that shuttle heme between the outer and inner membranes of the Gram-negative bacteria. Periplasmic binding proteins (PBPs) generally exhibit considerable conformational changes during the ligand binding process, whereas ShuT and PhuT belong to a class of PBPs that do not show such behavior based on their apo and holo crystal structures. By employing a series of molecular dynamic simulations on the ShuT and the PhuT, the dynamics and functions of the two PBPs were investigated. Through monitoring the distance changes between the two conserved glutamates of ShuT and PhuT, it was found the two PBPs were more flexible than previously assumed, exhibiting obvious opening-closing motions which were more remarkable in the apo runs of ShuT. Based on the results of the domain motion analysis, large scale conformational transitions were found in all apo runs of ShuT and PhuT, hinting that the domain motions of the two PBPs may be intrinsic. On the basis of the results of the principle component analysis, distinct opening-closing and twisting motion tendencies were observed not only in the apo, but also in the holo simulations of the two PBPs. The Gaussian network model was applied in order to analyze the hinge bending regions. The most important bending regions of ShuT and PhuT are located around the midpoints of their respective connecting helixes. Finally, the flexibilities and the details of the simulations of ShuT and PhuT were discussed. Characterized by the remarkably large flexibilities, the loop constituted by Ala 169, Gly170 and Gly171 of ShuT and the beta-turn constituted by Ala176, Gly177 and Gly178 of PhuT may be important for the functions of the two PBPs. Furthermore, the Asn254 of ShuT and the Arg228 of PhuT may be indispensable for the binding or unbinding of heme, since it is involved in the important hydrogen bonding to the propionate side-chains of heme.

[Effects of transcutaneous electric acupoint stimulation on plasma SOD and MDA in rats with sports fatigue].

OBJECTIVE: To observe the effect of transcutaneous electric acupoint stimulation (TEAS) on plasma superoxide dismutase (SOD) and malondialdehyde (MDA) in rats with sports fatigue so as to explore its mechanisms in resisting exercise-induced fatigue. METHODS: Twenty-seven male adult SD rats were randomly divided into control group (n 9), model group (n=9) and TEAS group (n=9). Sports fatigue model was established by using treadmill method, i.e. forcing the rat to run 10 min at a speed of 10 m/min, 15 m/min, 20 m/min, 24 m/min and 28 m/min, once daily for 6 days. TEAS (continuous waves, 2 Hz, 5 mA) was applied to unilateral "Zusanli" (ST 36) for 30 min, once per day for 7 days. The exhausted exercise time from starting running to exhaustion was recorded on the 7th day and the lactate levels, SOD activity and MDA contents in plasma were measured by lactate oxidase method, xanthine oxidase method and thiobarbituric acid method respectively. RESULTS: The duration of exhausted exercise in model group and TEAS group were (56.00 +/- 12.27) min and (70.88 +/- 13.74) min respectively, displaying significant increase in exercise tolerance after TEAS (P<0.05). The lactate level and MDA content in model group were markedly higher than those in control group (P<0.05); while compared with model group, lactate and MDA contents in TEAS group were significantly lower (P<0.05), and plasma SOD activity in TEAS group was significantly higher (P<0.01). No significant differences were found in plasma lactate content between control and TEAS groups and in SOD activity between model group and control group (P>0.05). SOD/MDA of model group was significantly lower than those of control group and TEAS group (P<0.05, 0.01). CONCLUSION: TEAS at "Zusanli"(ST 36) can effectively postpone exercise-induced fatigue by reducing accumulation of blood lactate, improving anti-oxidative ability and relieving lipid peroxidation in the rat.