Lulu Ning

Exploring the Molecular Mechanism of Cross-Resistance to HIV-1 Integrase Strand Transfer Inhibitors by Molecular Dynamics Simulation and Residue Interaction Network Analysis

The rapid emergence of cross-resistance to the integrase strand transfer inhibitors (INSTIs) has become a serious problem in the therapy of human immunodeficiency virus type 1 (HIV-1) infection. Understanding the detailed molecular mechanism of INSTIs cross-resistance is therefore critical for the development of new effective therapy against cross-resistance. On the basis of the homology modeling constructed structure of tetrameric HIV-1 intasome, the detailed molecular mechanism of the cross-resistance mutation E138K/Q148K to three important INSTIs (Raltegravir (RAL, FDA approved in 2007), Elvitegravir (EVG, FDA approved in 2012), and Dolutegravir (DTG, phase III clinical trials)) was investigated by using molecular dynamics (MD) simulation and residue interaction network (RIN) analysis. The results from conformation analysis and binding free energy calculation can provide some useful information about the detailed binding mode and cross-resistance mechanism for the three INSTIs to HIV-1 intasome. Binding free energy decomposition analysis revealed that Pro145 residue in the 140s 1oop (Gly140 to Gly149) of the HIV-1 intasome had strong hydrophobic interactions with INSTIs and played an important role in the binding of INSTIs to HIV-1 intasome active site. A systematic comparison and analysis of the RIN proves that the communications between the residues in the resistance mutant is increased when compared with that of the wild-type HIV-1 intasome. Further analysis indicates that residue Pro145 may play an important role and is relevant to the structure rearrangement in HIV-1 intasome active site. In addition, the chelating ability of the oxygen atoms in INSTIs (e.g., RAL and EVG) to Mg(2+) in the active site of the mutated intasome was reduced due to this conformational change and is also responsible for the cross-resistance mechanism. Notably, the cross-resistance mechanism we proposed could give some important information for the future rational design of novel INSTIs overcoming cross-resistance. Furthermore, the combination use of molecular dynamics simulation and residue interaction network analysis can be generally applicable to investigate drug resistance mechanism for other biomolecular systems.

Influence of the pathogenic mutations T188K/R/A on the structural stability and misfolding of human prion protein: Insight from molecular dynamics simulations

BACKGROUND Prion diseases are associated with a conformational switch for PrP from PrP(C) to PrP(Sc). Many genetic mutations are linked with prion diseases, such as mutations T188K/R/A with fCJD. SCOPE OF REVIEW MD simulations for the WT PrP and its mutants were performed to explore the underlying dynamic effects of T188 mutations on human PrP. Although the globular domains are fairly conserved, the three mutations have diverse effects on the dynamics properties of PrP, including the shift of H1, the elongation of native β-sheet and the conversion of S2-H2 loop to a 3(10) helix. MAJOR CONCLUSIONS Our present study indicates that the three mutants for PrP may undergo different pathogenic mechanisms and the realistic atomistic simulations can provide insights into the effects of disease-associated mutations on PrP dynamics and stability, which can enhance our understanding of how mutations induce the conversion from PrP(C) to PrP(Sc). General significance Our present study helps to understand the effects of T188K/R/A mutations on human PrP: despite the three pathogenic mutations almost do not alter the native structure of PrP, but perturb its stability. This instability may further modulate the oligomerization pathways and determine the features of the PrP(Sc) assemblies.

Molecular mechanism of the inhibition and remodeling of human islet amyloid polypeptide (hIAPP(1-37)) oligomer by resveratrol from molecular dynamics simulation.

Natural polyphenols are one of the most actively investigated categories of amyloid inhibitors, and resveratrol has recently been reported to inhibit and remodel the human islet amyloid polypeptide (hIAPP) oligomers and fibrils. However, the exact mechanism of its action is still unknown, especially for the full-length hIAPP1-37. To this end, we performed all-atom molecular dynamics simulations for hIAPP1-37 pentamer with and without resveratrol. The obtained results show that the binding of resveratrol is able to cause remarkable conformational changes of hIAPP1-37 pentamer, in terms of secondary structures, order degree, and morphology. By clustering analysis, two possible binding sites of resveratrol on the hIAPP1-37 pentamer were found, located at the grooves of the top and bottom surfaces of β-sheet layer, respectively. After the binding free energy calculation and residue energy decomposition, it can be concluded that the bottom site is the more possible one, and that the nonpolar interactions act as the driving force for the binding of hIAPP1-37 to resveratrol. In addition, Arg11 is the most important residue for the binding of resveratrol. The full understanding of inhibitory mechanism of resveratrol on the hIAPP1-37 oligomer, and the identification of its binding sites on this protein are helpful for the future design and discovery of new amyloid inhibitors.

The adsorption mechanism and induced conformational changes of three typical proteins with different secondary structural features on graphene

Nanomaterials (NMs) have been widely used in the biomedical field. To explore the biological effects of graphene as one of the most widely used NMs, we studied the adsorption behavior and induced conformational changes of proteins representing different secondary structures on graphene: β-strands (WW domain), mixed α/β structure (BBA protein), and α-helices (λ-repressor). Our results indicate these model proteins were adsorbed onto the graphene surface quickly and tightly, however, varied degrees of conformational changes were observed. During the adsorption process, we found the β motif is a stiffer structural unit than the α-helix. Moreover, the level of conformational changes of the proteins is related not only to their sequence and structural properties but also to their orientation. Overall, from the different levels of intermolecular interaction, the protein adsorption was driven by van der Waals, hydrophobic and π–π stacking interactions. Our work suggests that classical molecular dynamics simulations and MM-GBSA calculations can provide useful information about the dynamics and energetics of the adsorption of proteins onto graphene. We believe that these findings will help us to further understand the adsorption of proteins on hydrophobic carbon nanomaterials at the atomic level.

In Silico Identification of Protein S-Palmitoylation Sites and Their Involvement in Human Inherited Disease.

S-Palmitoylation is a key regulatory mechanism controlling protein targeting, localization, stability, and activity. Since increasing evidence shows that its disruption is implicated in many human diseases, the identification of palmitoylation sites is attracting more attention. However, the computational methods that are published so far for this purpose have suffered from a poor balance of sensitivity and specificity; hence, it is difficult to get a good generalized prediction ability on an external validation set, which holds back the further analysis of associations between disruption of palmitoylation and human inherited diseases. In this work, we present a reliable identification method for protein S-palmitoylation sites, called SeqPalm, based on a series of newly composed features from protein sequences and the synthetic minority oversampling technique. With only 16 extracted key features, this approach achieves the most favorable prediction performance up to now with sensitivity, specificity, and Matthews correlation coefficient values of 95.4%, 96.3%, and 0.917, respectively. Then, all known disease-associated variations are studied by SeqPalm. It is found that 243 potential loss or gain of palmitoylation sites are highly associated with human inherited disease. The analysis presents several potential therapeutic targets for inherited diseases associated with loss or gain of palmitoylation function. There are even biological evidence that are coordinate with our prediction results. Therefore, this work presents a novel approach to discover the molecular basis of pathogenesis associated with abnormal palmitoylation. SeqPalm is now available online at http://lishuyan.lzu.edu.cn/seqpalm , which can not only annotate the palmitoylation sites of proteins but also distinguish loss or gain of palmitoylation sites by protein variations.

Stabilities and structures of islet amyloid polypeptide (IAPP22-28) oligomers: from dimer to 16-mer.

BACKGROUND The formation of amyloid fibrils is associated with many age-related degenerative diseases. Nevertheless, the molecular mechanism that directs the nucleation of these fibrils is not fully understood. METHODS Here, we performed MD simulations for the NFGAILS motif of hIAPP associated with the type II diabetes to estimate the stabilities of hIAPP22-28 protofibrils with different sizes: from 2 to 16 chains. In addition, to study the initial self-assembly stage, 4 and 8 IAPP22-28 chains in explicit solvent were also simulated. RESULTS Our results indicate that the ordered protofibrils with no more than 16 hIAPP22-28 chains will be structurally stable in two layers, while one-layer or three-layer models are not stable as expected. Furthermore, the oligomerization simulations show that the initial coil structures of peptides can quickly aggregate and convert to partially ordered β-sheet-rich oligomers. CONCLUSIONS Based on the obtained results, we found that the stability of an IAPP22-28 oligomer was not only related with its size but also with its morphology. The driving forces to form and stabilize an oligomer are the hydrophobic effects and backbone H-bond interaction. Our simulations also indicate that IAPP22-28 peptides tend to form an antiparallel strand orientation within the sheet. GENERAL SIGNIFICANCE Our finding can not only enhance the understanding about potential mechanisms of hIAPP nuclei formation and the extensive structural polymorphisms of oligomers, but also provide valuable information to develop potential β-sheet formation inhibitors against type II diabetes.

Exploring structural and thermodynamic stabilities of human prion protein pathogenic mutants D202N, E211Q and Q217R

The central event in the pathogenesis of prion protein (PrP) is a profound conformational change from its α-helical (PrP(C)) to its β-sheet-rich isoform (PrP(Sc)). Many single amino acid mutations of PrP are associated with familial prion diseases, such as D202N, E211Q, and Q217R mutations located at the third native α-helix of human PrP. In order to explore the underlying structural and dynamic effects of these mutations, we performed all-atom molecular dynamics (MD) simulations for the wild-type (WT) PrP and its mutants. The obtained results indicate that these amino acid substitutions have subtle effects on the protein structures, but show large changes of the overall electrostatic potential distributions. We can infer that the changes of PrP electrostatic surface due to the studied mutations may influence the intermolecular interactions during the aggregation process. In addition, the mutations also affect the thermodynamic stabilities of PrP.

Structural Diversity and Initial Oligomerization of PrP106–126 Studied by Replica-Exchange and Conventional Molecular Dynamics Simulations

Prion diseases are marked by cerebral accumulation of the abnormal isoform of the prion protein. A fragment of prion protein composed of residues 106–126 (PrP106–126) exhibits similar properties to full length prion and plays a key role in the conformational conversion from cellular prion to its pathogenic pattern. Soluble oligomers of PrP106–126 have been proposed to be responsible for neurotoxicity. However, the monomeric conformational space and initial oligomerization of PrP106–126 are still obscure, which are very important for understanding the conformational conversion of PrP106–126. In this study, replica exchange molecular dynamics simulations were performed to investigate monomeric and dimeric states of PrP106–126 in implicit solvent. The structural diversity of PrP106–126 was observed and this peptide did not acquire stable structure. The dimeric PrP106–126 also displayed structural diversity and hydrophobic interaction drove the dimerization. To further study initial oligomerization of PrP106–126, 1 µs conventional molecular dynamics simulations of trimer and tetramer formation were carried out in implicit solvent. We have observed the spontaneous formation of several basic oligomers and stable oligomers with high β-sheet contents were sampled in the simulations of trimer and tetramer formation. The β-hairpin formed in hydrophobic tail of PrP106–126 with residues 118–120 in turn may stabilize these oligomers and seed the formation oligomers. This study can provide insight into the detailed information about the structure of PrP106–126 and the dynamics of aggregation of monomeric PrP106–126 into oligomers in atomic level.

The role of Cys179–Cys214 disulfide bond in the stability and folding of prion protein: insights from molecular dynamics simulations

AbstractPrion diseases are associated with misfolding and aggregation of prion protein (PrP). Cellular prion protein contains a disulfide bond linking Cys residues at positions 179 and 214. It has been proposed that this disulfide bond plays an important role in the conversion between cellular (PrPC) and the scrapie form of prion protein (PrPSc). To probe the role of this disulfide bond in the stability and folding of prion protein, we employed molecular dynamics simulations to study the reduced prion protein and a variant of PrP in which the two cysteines were replaced by alanines residues. The simulations highlighted the changes that occurred upon breakage of the disulfide bond. Breakage of the disulfide bond resulted in a shift of H1, elongation of the native β-sheet and perturbation of the hydrophobic core of huPrP. The changes are similar to the conformational transitions of prion protein in low pH, in denaturing conditions or with pathogenic mutations, which indicate that rupture of the disulfide bond may lead to the misfolding of prion protein. FigureThe breakage of Cys179-Cys214 disulfide bond resulted in conformational transitions of prion protein

Effects of the Pathogenic Mutation A117V and the Protective Mutation H111S on the Folding and Aggregation of PrP106-126: Insights from Replica Exchange Molecular Dynamics Simulations

The fragment 106-126 of prion protein exhibits similar properties to full-length prion. Experiments have shown that the A117V mutation enhances the aggregation of PrP106-126, while the H111S mutation abolishes the assembly. However, the mechanism of the change in the aggregation behavior of PrP106-126 upon the two mutations is not fully understood. In this study, replica exchange molecular dynamics simulations were performed to investigate the conformational ensemble of the WT PrP106-126 and its two mutants A117V and H111S. The obtained results indicate that the three species are all intrinsically disordered but they have distinct morphological differences. The A117V mutant has a higher propensity to form β-hairpin structures than the WT, while the H111S mutant has a higher population of helical structures. Furthermore, the A117V mutation increases the hydrophobic solvent accessible surface areas of PrP106-126 and the H111S mutation reduces the exposure of hydrophobic residues. It can be concluded that the difference in populations of β-hairpin structures and the change of hydrophobic solvent accessible areas may induce the different aggregation behaviors of the A117V and the H111S mutated PrP106-126. Understanding why the two mutations have contrary effects on the aggregation of PrP106-126 is very meaningful for further elucidation of the mechanism underlying aggregation and design of inhibitor against aggregation process.