Weizu Chen

Thermal Stability and Unfolding Pathways of Sso7d and its Mutant F31A: Insight from Molecular Dynamics Simulation

Abstract The thermo-stability and unfolding behaviors of a small hyperthermophilic protein Sso7d as well as its single-point mutation F31A are studied by molecular dynamics simulation at temperatures of 300 K, 371 K and 500 K. Simulations at 300 K show that the F31A mutant displays a much larger flexibility than the wild type, which implies that the mutation obviously decreases the proteins stability. In the simulations at 371 K, although larger fluctuations were observed, both of these two maintain their stable conformations. High temperature simulations at 500 K suggest that the unfolding of these two proteins evolves along different pathways. For the wild-type protein, the C-terminal alpha-helix is melted at the early unfolding stage, whereas it is destroyed much later in the unfolding process of the F31A mutant. The results also show that the mutant unfolds much faster than its parent protein. The deeply buried aromatic cluster in the F31A mutant dissociates quickly relative to the wild-type protein at high temperature. Besides, it is found that the triple-stranded antiparallel β-sheet in the wild-type protein plays an important role in maintaining the stability of the entire structure.

Protein–protein docking with binding site patch prediction and network-based terms enhanced combinatorial scoring†

Protein–protein docking has made much progress in recent years, but challenges still exist. Here we present the application of our docking approach HoDock in CAPRI. In this approach, a binding site prediction is implemented to reduce docking sampling space and filter out unreasonable docked structures, and a network‐based enhanced combinatorial scoring function HPNCscore is used to evaluate the decoys. The experimental information was combined with the predicted binding site to pick out the most likely key binding site residues. We applied the HoDock method in the recent rounds of the CAPRI experiments, and got good results as predictors on targets 39, 40, and 41. We also got good results as scorers on targets 35, 37, 40, and 41. This indicates that our docking approach can contribute to the progress of protein–protein docking methods and to the understanding of the mechanism of protein–protein interactions. Proteins 2010.

A novel high-throughput format assay for HIV-1 integrase strand transfer reaction using magnetic beads

AbstractAim:To develop a novel high-throughput format assay to monitor the integrase (IN) strand transfer (ST) reaction in vitro and apply it to a reaction character study and the identification of antiviral drugs.Methods:The donor DNA duplex, with a sequence identical to the U5 end of HIV-1 long terminal repeats, is labeled at its 5′ end with biotin (BIO). The target DNA duplex is labeled at its 3′ end with digoxin (DIG). IN mediates the integration of donor DNA into target DNA and results in a 5′ BIO and 3′ DIG-labeled duplex DNA product. Streptavidin-coated magnetic beads were used to capture the product, and the amount of DIG was measured as the ST reaction product. The assay was optimized in 96-well microplate format for high-throughput screening purpose. Moreover, the assay was applied in a ST reaction character study, and the efficiency of the assay in the identification of antiviral compounds was tested.Results:The end-point values, measured as absorbance at 405 nm was approximately 1.5 for the IN-mediated ST reaction as compared with no more than 0.05 of background readings. The ST reaction character and the half maximal inhibitory concentration (IC50) values of 2 known IN inhibitors obtained in our assay were similar to previously reported results using other assays. The evaluation parameter Z′ factor for this assay ranged from 0.6 to 0.9.Conclusion:The assay presented here has been proven to be rapid, sensitive, and specific for the detection of IN ST activity, the reaction character study, as well as for the identification of antiviral drugs targeting IN.

Amino acid network and its scoring application in protein-protein docking

Protein-protein complex, composed of hydrophobic and hydrophilic residues, can be divided into hydrophobic and hydrophilic amino acid network structures respectively. In this paper, we are interested in analyzing these two different types of networks and find that these networks are of small-world properties. Due to the characteristic complementarity of the complex interfaces, protein-protein docking can be viewed as a particular network rewiring. These networks of correct docked complex conformations have much more increase of the degree values and decay of the clustering coefficients than those of the incorrect ones. Therefore, two scoring terms based on the network parameters are proposed, in which the geometric complementarity, hydrophobic-hydrophobic and polar-polar interactions are taken into account. Compared with a two-term energy function, a simple scoring function HPNet which includes the two network-based scoring terms shows advantages in two aspects, not relying on energy considerations and better discrimination. Furthermore, combing the network-based scoring terms with some other energy terms, a new multi-term scoring function HPNet-combine can also make some improvements to the scoring function of RosettaDock.

High-throughput real-time assay based on molecular beacons for HIV-1 integrase 3'-processing reaction.

AbstractAim:To develop a high-throughput real-time assay based on molecular beacons to monitor the integrase 3′-processing reaction in vitro and apply it to inhibitor screening.Methods:The recombinant human immunodeficiency virus (HIV)-1 integrase (IN) is incubated with a 38 mer oligonucleotide substrate, a sequence identical to the U5 end of HIV-1 long terminal repeats (LTR). Based on the fluorescence properties of molecular beacons, the substrate is designed to form a stem-loop structure labeled with a fluorophore at the 5′ end and a quencher at the 3′ end. IN cleaves the terminal 3′-dinucleotide containing the quencher, resulting in an increase in fluorescence which can be monitored on a spectrofluorometer. To optimize this assay, tests were performed to investigate the effects of substrates, enzyme and the metal ion concentrations on the IN activity and optimal parameters were obtained. Moreover, 2 IN inhibitors were employed to test the performance of this assay in antiviral compound screening.Results:The fluorescent intensity of the reaction mixture varies linearly with time and is proportional to the velocity of the 3′-processing reaction. Tests were performed and the results showed that the optimal rate was obtained for a reaction mixture containing 50 mg/L recombinant HIV-1 IN, 400 nmol/L substrate, and 10 mmol/L Mn2+. The IN 3′-processing reaction under the optimal conditions showed a more than 18-fold increase in the fluorescence intensity compared to the enzyme-free control. The IC50 values of the IN inhibitors obtained in our assay were similar to the values obtained from a radiolabeled substrate assay.Conclusion:Our results demonstrated that this is a fast, reliable, and sensitive method to monitor HIV IN 3′-processing reaction and that it can be used for inhibitor screening.

Uncovering the Properties of Energy-Weighted Conformation Space Networks with a Hydrophobic-Hydrophilic Model

The conformation spaces generated by short hydrophobic-hydrophilic (HP) lattice chains are mapped to conformation space networks (CSNs). The vertices (nodes) of the network are the conformations and the links are the transitions between them. It has been found that these networks have “small-world” properties without considering the interaction energy of the monomers in the chain, i. e. the hydrophobic or hydrophilic amino acids inside the chain. When the weight based on the interaction energy of the monomers in the chain is added to the CSNs, it is found that the weighted networks show the “scale-free” characteristic. In addition, it reveals that there is a connection between the scale-free property of the weighted CSN and the folding dynamics of the chain by investigating the relationship between the scale-free structure of the weighted CSN and the noted parameter Z score. Moreover, the modular (community) structure of weighted CSNs is also studied. These results are helpful to understand the topological properties of the CSN and the underlying free-energy landscapes.

Molecular dynamics simulation of the transmembrane subunit of BtuCD in the lipid bilayer

Based on the crystal structure of the vitamin B12 transporter protein of Escherichia coli (BtuCD) a system consisting of the BtuCD transmembrane domain (BtuC) and the palmitoyloleoyl phosphatidylcholine (POPC) lipid bilayer was constructed in silica, and a more-than-57-nanosecond molecular dynamics (MD) simulation was performed on it to reveal the intrinsic functional motions of BtuC. The results showed that a stable protein-lipid bilayer was obtained and the POPC lipid bilayer was able to adjust its thickness to match the embedded BtuC which underwent relatively complicated motions. These results may help to understand the mechanism of transmembrane substrate transport at the atomic level.

Studies on binding free energies and the binding mode by docking and MM-PBSA in gp41-ligand complex

Despite the synthetic peptides inhibit HIV-1 entry, its application may be limited due to the high cost of the peptide production and lack of its oral availability. Thus, it is necessary to identify the small molecule inhibitors reacting with the same or overlapping target sites on gp41 recognizing the antiviral peptides. In this paper, a small inhibitor (TP1) is docked into the hydrophobic grooves of gp41 by using Autodock software, resulting in five alternative energetically favorable models. The molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) method is applied to calculate the binding free energies. The data from other studies were used to define our preferred models. We found that only one binding mode is supported by the experimental evidence. The model could be used to design more effective HIV-1 inhibitors targeted to the HIV-1 gp41 core structure.

A holistic molecular docking approach for predicting protein-protein complex structure

A holistic protein-protein molecular docking approach, HoDock, was established, composed of such steps as binding site prediction, initial complex structure sampling, refined complex structure sampling, structure clustering, scoring and final structure selection. This article explains the detailed steps and applications for CAPRI Target 39. The CAPRI result showed that three predicted binding site residues, A191HIS, B512ARG and B531ARG, were correct, and there were five submitted structures with a high fraction of correct receptor-ligand interface residues, indicating that this docking approach may improve prediction accuracy for protein-protein complex structures.

Peptides Design Based on the Interfacial Helix of Integrase Dimer

HIV-I integrase (IN) plays a crucial role in the retroviral life cycle. The peptides derived from the helix of IN were reported to have the potency of inhibition. We designed a series of peptides based on interface helices alpha1 and alpha5 with the aim of increasing their inhibitory activity. The helix-forming tendency and the affinity with IN were essential for interfacial peptide inhibitors. The MD simulation and AGADIR prediction both showed favorable results for the designed peptides. The binding mode and binding free energy of peptide and IN were investigated subsequently to test our design. The improvement in binding free energy compared with that of alpha1 and alpha5 indicates that some of the designed peptides may have a higher potency for inhibiting the dimerization of IN. This study provides some useful information for rational design of IN peptide inhibitor