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.

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.

[Effects of LncRNA UNC5B-AS1 on adhesion, invasion and migration of lung cancer cells and its mechanism].

Objective: To investigate the effects of long-chain non-coding RNA (lncRNA) UNC5B-AS1 on the adhesion, invasion and migration of lung cancer cells by regulating the expression of miR-218-5p. Methods: Twenty specimens of lung cancer patients and corresponding paracancerous tissues were surgically removed and collected from the oncology department of Chongqing Three Gorges Central Hospital from June 2017 to June 2019. Real-time quantitative PCR (qRT-PCR) was used to detect the expressions of UNC5B-AS1 in human bronchial epithelial cells HBE and different lung cancer cells of A549, H1437, H1975, H1299 and H460. UNC5B-AS1 siRNA was transfected into lung cancer A549 cells. Adhesion assay, transwell invasion assay and scratch assay were used to detect the effect of UNC5B-AS1 on adhesion, invasion and migration of A549 cells. qRT-PCR and dual luciferase reporter gene were used for the detection and identification of UNC5B-AS1 targeting miR-218-5p. The expression of epithelial-mesenchymal transition (EMT)-related protein was detected by Western blot. Results: The expression of UNC5B-AS1 in lung cancer tissues and cells was significantly higher than that in adjacent tissues and bronchial epithelial cells (P＜0.05). The expression of UNC5B-AS1 in lung cancer A549 cells was the highest (P＜0.05). Down-regulation of UNC5B-AS1 expression inhibited adhesion, invasion and migration of A549 cells (P＜0.05). qRT-PCR and dual luciferase reporter assay experiments showed that UNC5B-AS1 targeted the regulation of miR-218-5p expression. Down-regulation of UNC5B-AS1 inhibited E-cadherin protein expression and promoted Vimentin and Twist protein expression. Conclusion: lncRNA UNC5B-AS1 promotes adhesion, invasion and migration of lung cancer cells through targeted regulation of miR-218-5p expression, and its mechanism may be related to the promotion of EMT.

Prediction of the binding mode between BMS-378806 and HIV-1 gp120 by docking and molecular dynamics simulation.

BMS-378806 is a newly discovered small molecule that effectively blocks the binding of CD4 with gp120. The binding mode of this kind of inhibitor remains unknown. In this paper, AutoDock 3.0 in conjunction with molecular dynamics simulation, accommodating the receptor's flexibility, was used to explore the binding mode between BMS-378806 and gp120. Two structures, Mode I and Mode II, with the lowest docking energy were selected as different representative binding modes. The analysis of the results from the molecular dynamics simulation indicated that the binding of BMS-348806 in Mode II is more stable. The average structure of Mode II was analyzed and compared with the experimental data. The conclusion was that BMS-378806 inserts the azaindole ring deeply into the PHE43 cavity and makes contact with a number of residues in the cavity, on the cavity and near the cavity. This study benefits the understanding of the mechanism of this kind of inhibitor and may provide useful information for rational drug design.

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.

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.

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.

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.

Exploring binding mode for styrylquinoline HIV-1 integrase inhibitors using comparative molecular field analysis and docking studies.

AIM: To understand pharmacophore properties of styrylquinoline derivatives and to design inhibitors of HIV-1 integrase. METHODS: Comparative molecular field analysis (CoMFA) was performed to analyze three-dimensional quantitative structure-activity relationship (3D-QSAR) of styrylquinoline derivatives. Thirty-eight compounds were randomly divided into a training set of 28 compounds and a test set of 10 compounds. The stability of 3D-QSAR models was proved by the analysis of cross-validated and non-cross-validated methods. Moreover, the binding mode of these compounds and integrase was constructed by AutoDock program. RESULTS: The CoMFA model of the training compounds was reasonably predicted with cross-validated coefficient (q2) and conventional (r2) values (up to 0.696 and 0.754). Then the model was validated by the test set. The resulting CoMFA maps visualized structural requirements for the biological activity of these inhibitors. Docking results showed that a carboxyl group at C-7 and a hydroxyl group at C-8 in the quinoline subunit, bound closely to the divalent metal cofactor (Mg2+) around the integrase catalytic site. Moreover, there is a linear correlation between the binding energy of the inhibitors with integrase and their inhibitory effect. CONCLUSIONS: The present study indicated that the CoMFA model together with docking results could give us helpful hints for drug design as well as interpretation of the binding affinity between these inhibitors and integrase.