Tianyang Sun

Interaction of P-glycoprotein with anti-tumor drugs: the site, gate and pathway

Understanding the mechanism and pathway of anti-cancer drugs to be pumped out by P-glycoprotein (P-gp) in cancer cell is very important for the successful chemotherapy. P-gp is a member of ATP-binding cassette (ABC) transporters. In this study, random accelerated molecular dynamics (RAMD) simulation was used to explore the potential egress pathway of ligands from the binding pocket. This could be considered as a reverse process of drug binding. The most possible portal of drugs to dissociate is TM4/TM6, which is almost the same for different drugs, such as paclitaxel and doxorubicin. The interactions in the binding site are found to be remarkably stronger than that outside of the binding site. The results were suggested by the free energy calculation between P-gp and different drugs from metadynamics simulation. All the results indicate that the flexibility of inner residues, especially the residue Phe339, is very important for the drugs to access the binding site.

A molecular dynamics study on pH response of protein adsorbed on peptide-modified polyvinyl alcohol hydrogel

The interactions between proteins and functional biomaterials under different physical and environmental conditions need to be understood when designing biomedical devices. Herein, we present a molecular dynamics simulation study of the fragment antigen-binding (Fab) of trastuzumab (a monoclonal antibody) and its complex with a peptide-modified polyvinyl alcohol (PVA) hydrogel at different pH values. Consistent with experiments, PVA when modified by charged ligands does shrink as a direct response to a drop in the pH. The protein maintains a stable conformation when adsorbed on the hydrogel matrix with a varied pH, showing no signs of denaturation in all simulated systems, suggesting that peptide-grafted PVA is a good biocompatible material. Under neutral conditions, the hydrogel alone stabilizes the interactions between the protein and the peptide ligands. Strikingly under acidic conditions the protein-ligand interactions are disrupted by a collective protonation of ligands. A sharp decrease in the interaction energies, accompanied by the sudden increase of the protein-ligand distance, indicates a rapid pH response in the protein-hydrogel complex. This will be important in protein delivery and purification. The effect of pH on the interactions and the dynamics of the protein and the sudden pH response of the hydrogel at the atomic level present a new functional perspective in developing new hydrogels with desirable properties.

On the loading mechanism of ssDNA into carbon nanotubes

Understanding of the mechanism and dynamics of DNA loading into carbon nanotubes (CNTs) is very important for the promising applications of CNTs in DNA sequencing, drug delivery and gene delivery systems etc. In this work, the loading mechanism and dynamics of different ssDNA oligomers into single-walled carbon nanotubes (SWNTs) was investigated through molecular dynamics simulations, steered molecular dynamics simulation and binding free energy calculations. Our simulation results showed that the loading of different ssDNA oligomers into the zigzag SWNT is much easier than for the armchair SWNT. Confined in both zigzag and armchair type SWNTs, ssDNA oligomers have a helical structure and their bases adapt the orientation parallel to the interior wall. From detailed analysis of the interaction energy, potential of mean force (PMF) of the unloading process and nucleotide binding free energy, our results show that the chirality of SWNTs has a large effect on the binding strength of nucleotides, and hence affects the loading dynamics of ssDNA into SWNTs.

Molecular dynamics simulations indicate that DNA bases using graphene nanopores can be identified by their translocation times

The improvement of the resolution of DNA sequencing by nanopore technology is very important for its real-life application. In this paper, we report our work on using molecular dynamics simulation to study the dependence of DNA sequencing on the translocation time of DNA through a graphene nanopore, using the single-strand DNA fragment translocation through graphene nanopores with diameters down to ∼2 nm as examples. We found that A, T, C, and G could be identified by the difference in the translocation time between different types of nucleotides through 2 nm graphene nanopores. In particular, the recognition of the graphene nanopore for different nucleotides can be greatly enhanced in a low electric field. Our study suggests that the recognition of a graphene nanopore by different nucleotides is the key factor for sequencing DNA by translocation time. Our study also indicates that the surface of a graphene nanopore can be modified to increase the recognition of nucleotides and to improve the resolution of DNA sequencing based on the DNA translocation time with a suitable electric field.

Contribution of Water Molecules in the Spontaneous Release of Protein by Graphene Sheets

Applications of graphene sheets in the fields of biosensors and biomedical devices are limited by the aqueous solubility of graphene. Consequently, understanding the role of water molecules in the aggregation or dispersion of graphene in aqueous solution with a biomolecule is of vital importance to its application. Herein, protein is spontaneously released by the layer-to-layer aggregation of two single-layer graphene sheets due to van der Waals force between the sheets. The properties of water molecules, including density and dynamics, are discussed in detail. The dynamic behavior of aggregation of graphene sheets is triggered by the dynamics of water molecules. To stabilize dispersed graphene sheets in aqueous solution, the density of water molecules between the graphene sheets should be larger than 0.83 g cm(-3), and graphene modified by hydroxyl groups could be a good choice. The stability of a model protein on the graphene sheet is studied to investigate the biological compatibility of graphene sheets. To be a material with good biocompatibility, graphene should be functionalized by hydrophilic groups. The results presented herein could be helpful in the research and application of graphene sheets in the fields of biomaterials, biosensors, and biomedical devices.

Trastuzumab-Peptide Interactions: Mechanism and Application in Structure-Based Ligand Design

Understanding of protein-ligand interactions and its influences on protein stability is necessary in the research on all biological processes and correlative applications, for instance, the appropriate affinity ligand design for the purification of bio-drugs. In this study, computational methods were applied to identify binding site interaction details between trastuzumab and its natural receptor. Trastuzumab is an approved antibody used in the treatment of human breast cancer for patients whose tumors overexpress the HER2 (human epidermal growth factor receptor 2) protein. However, rational design of affinity ligands to keep the stability of protein during the binding process is still a challenge. Herein, molecular simulations and quantum mechanics were used on protein-ligand interaction analysis and protein ligand design. We analyzed the structure of the HER2-trastuzumab complex by molecular dynamics (MD) simulations. The interaction energies of the mutated peptides indicate that trastuzumab binds to ligand through electrostatic and hydrophobic interactions. Quantitative investigation of interactions shows that electrostatic interactions play the most important role in the binding of the peptide ligand. Prime/MM-GBSA calculations were carried out to predict the binding affinity of the designed peptide ligands. A high binding affinity and specificity peptide ligand is designed rationally with equivalent interaction energy to the wild-type octadecapeptide. The results offer new insights into affinity ligand design.

Interaction between IGFBP7 and insulin: a theoretical and experimental study

Insulin-like growth factor binding protein 7 (IGFBP7) can bind to insulin with high affinity which inhibits the early steps of insulin action. Lack of recognition mechanism impairs our understanding of insulin regulation before it binds to insulin receptor. Here we combine computational simulations with experimental methods to investigate the interaction between IGFBP7 and insulin. Molecular dynamics simulations indicated that His200 and Arg198 in IGFBP7 were key residues. Verified by experimental data, the interaction remained strong in single mutation systems R198E and H200F but became weak in double mutation system R198E-H200F relative to that in wild-type IGFBP7. The results and methods in present study could be adopted in future research of discovery of drugs by disrupting protein–protein interactions in insulin signaling. Nevertheless, the accuracy, reproducibility, and costs of free-energy calculation are still problems that need to be addressed before computational methods can become standard binding prediction tools in discovery pipelines.