RESUMEN
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has led to the global spread of the coronavirus disease (COVID-19), which has caused great loss of life and property worldwide. We investigated the regulatory mechanism with the antibody targeting the N-terminal domain (NTD) of the S protein by molecular dynamic simulation. It was found that the structure of the S1-4A8 complex experienced the largest change when the receptor binding domain (RBD) of S1 was in the Up state. By calculating the angle between domains of S1 in the Down and Up states, we found that the RBD angle changed more in the Up state. We further performed binding free energy calculations for S1-4A8 complexes in both Up and Down states, and the results showed that 4A8 has a stronger affinity with NTD in the Up state. These results indicate that 4A8 plays a stronger regulatory role in the RBD Up state. The N3 and N5 loops on the NTD are the main antigen-antibody binding sites, and residues on the antibody complementarity determining region 3 (CDR3) in the Up state can penetrate deeper into the hydrophobic pocket at the bottom of the N5 loop to form a tighter binding. Through the tICA method, we found that except the residues at the binding interface, distant residues including A609, V610, G652, and A653 at the linker region of subdomain 2, and residues S359 and N360 near the bottom of RBD are able to influence the regulatory effect in the long-range. This work provides new insights into the neutralization mechanism of targeting NTD antibodies in SARS-CoV-2.
Asunto(s)
COVID-19 , SARS-CoV-2 , Humanos , Unión Proteica , AnticuerposRESUMEN
Transthyretin amyloidosis (ATTR) is a progressive and systemic disease caused by the misfolding and amyloid aggregation of transthyretin (TTR). Stabilizing the TTR tetramers and disrupting the formed TTR aggregation are treated as a promising strategy for the treatment of ATTR. Previous studies have reported that epigallocatechin gallate (EGCG) can participate in the whole process of TTR aggregation to prevent ATTR. However, the interaction mechanism of EGCG in this process is still obscure. In this work, we performed molecular dynamics simulations to investigate the interactions between EGCG and TTR tetramers, and between EGCG and TTR aggregates formed by the V30M mutation. The obtained results suggest that EGCG at the binding site of the V30M TTR tetramer can form stable hydrogen bonds with residues in the flexible AB-loop and EF-helix-loop, which reduces the structural mobility of these regions significantly. Additionally, the polyaromatic property of EGCG contributes to the increasement of hydrophobicity at the binding site and thus makes the tetramer difficult to be solvated and dissociated. For V30M-TTR-generated aggregates, EGCG can promote the dissociation of boundary ß-strands by destroying key residue interactions of TTR aggregates. Moreover, EGCG is capable of inserting into the side-chain of residues of neighboring ß-strands and disrupting the highly structured aggregates. Taken together, this study elucidates the role of EGCG in preventing TTR amyloidosis, which can provide important theoretical support for the future of drug design for ATTR.
RESUMEN
The dissociation of the transthyretin (TTR) tetramer into a monomer is closely related to various TTR amyloidoses in humans. While the tetramer dissociation has been reported to be the rate-limiting step for TTR aggregation, few details are known about the mechanism. Here, molecular dynamics (MD) simulations were performed by combining conventional MD and biased metadynamics to investigate the mechanism for the wild-type (WT) and mutant (T119M) structures. Both were found to have a great deal in common. Conventional MD simulations reveal that interfacial hydrophobic interactions contribute significantly to stabilize the tetramer. Interfacial residues including L17, V20, L110, and V121 with close contacts form a hydrophobic channel. Metadynamics simulations indicate that the mouth opening of the hydrophobic channel is the first and the most difficult step for dissociation. Interactions of V20 between opposing dimers lock four monomers into the tetramer, and disruption of the interactions is found to be involved in the final step. During the dissociation, an increasing extent of solvation was observed by calculating the radial distribution functions of water around interfacial hydrophobic residues, suggesting that water plays a role in driving the tetramer dissociation. Moreover, compared to T119, residue M119 has a longer side chain that extends into the hydrophobic channel, making solvation more difficult, consistent with a higher energy barrier for dissociation of the T119M tetramer. This result provides a good explanation for the protective role of the T119M mutation. Overall, this study can provide atomic-level insights to better understand the pathogenesis of TTR amyloidosis and guide rational drug design in the future.
Asunto(s)
Simulación de Dinámica Molecular , Prealbúmina , Humanos , Prealbúmina/química , Prealbúmina/genética , Mutación , Diseño de FármacosRESUMEN
SARS-CoV-2 has led to a global pandemic of new crown pneumonia, which has had a tremendous impact on human society. Antibody drug therapy is one of the most effective way of combating SARS-CoV-2. In order to design potential antibody drugs with high affinity, we used antibody S309 from patients with SARS-CoV as the target antibody and RBD of S protein as the target antigen. Systems with RBD glycosylated and non-glycosylated were constructed to study the influence of glycosylation. From the results of molecular dynamics simulations, the steric effects of glycans on the surface of RBD plays a role of "wedge", which makes the L335-E340 region of RBD close to the CDR3 region of the heavy chain of antibody and increases the contact area between antigen and antibody. By mutating the key residues of antibody at the interaction interface, we found that the binding affinities of antibody mutants G103A, P28W and Y100W were all stronger than that of the wild-type, especially for the G103A mutant. G103A significantly reduces the distance between the binding region of L335-K356 in the antigen and P28-Y32 of heavy chain in the antibody through structural transition. Taken together, the antibody design method described in this work can provide theoretical guidance and a time-saving method for antibody drug design.
Asunto(s)
COVID-19 , SARS-CoV-2 , Humanos , Simulación de Dinámica Molecular , Anticuerpos , Diseño de Fármacos , Unión ProteicaRESUMEN
The development of novel metal organic framework (MOF) friction power generation materials with high stability is important. This paper reports the first example of a double-helix metal chain organic framework with a network structure (ZUT-8). ZUT-8 shows high chemical stability, functional adjustability, and excellent output performance of friction power generation, which is superior to traditional coordination polymer materials. The cathodic protection system with ZUT-8 can prevent metal corrosion significantly. The output performance can be improved effectively by enhancing the conjugate effect of the linker. The theoretical calculation results showed that an increase in the degree of conjugation could significantly reduce the band gap, thereby affecting the friction power output signal. This study opens the door to constructing MOF materials with a double-helix metal chain and will promote their potential applications in self-powered electrochemical cathodic protection.
RESUMEN
Traditional Chinese medicine treatment of diseases has been recognized, but the material basis and mechanisms are not clear. In this study, target prediction of the antigastric cancer (GC) effect of Guiqi Baizhu (GQBZP) and the analysis of potential key compounds, key targets, and key pathways for the therapeutic effects against GC were carried out based on the method of network analysis and Kyoto Encyclopedia of Genes and Genomes enrichment. There were 33 proteins shared between GQBZP and GC, and 131 compounds of GQBZP had a high correlation with these proteins, indicating that the PI3K-AKT signaling pathway might play a key role in GC. From these studies, we selected human epidermal growth factor receptor 2 (HER2) and programmed cell death 1-ligand 1 (PD-L1) for docking; the results showed that 385 and 189 compounds had high docking scores with HER2 and PD-L1, respectively. Six compounds were selected for microscale thermophoresis (MST). Daidzein/quercetin and isorhamnetin/formononetin had the highest binding affinity for HER2 and PD-L1, with Kd values of 3.7 µmol/L and 490, 667, and 355 nmol/L, respectively. Molecular dynamics simulation studies based on the docking complex structures as the initial conformation yielded the binding free energy between daidzein/quercetin with HER2 and isorhamnetin/formononetin with PD-L1, calculated by molecular mechanics Poisson-Boltzmann surface area, of -26.55, -14.18, -19.41, and -11.86 kcal/mol, respectively, and were consistent with the MST results. In vitro experiments showed that quercetin, daidzein, and isorhamnetin had potential antiproliferative effects in MKN-45 cells. Enzyme activity assays showed that quercetin could inhibit the activity of HER2 with an IC50 of 570.07 nmol/L. Our study provides a systematic investigation to explain the material basis and molecular mechanism of traditional Chinese medicine in treating diseases.
Asunto(s)
Antígeno B7-H1/metabolismo , Medicamentos Herbarios Chinos/metabolismo , Proteínas de Neoplasias/metabolismo , Receptor ErbB-2/metabolismo , Neoplasias Gástricas/metabolismo , Antígeno B7-H1/química , Línea Celular Tumoral , Proliferación Celular/efectos de los fármacos , Medicamentos Herbarios Chinos/uso terapéutico , Humanos , Isoflavonas/metabolismo , Isoflavonas/farmacología , Simulación del Acoplamiento Molecular/métodos , Proteínas de Neoplasias/química , Fosfatidilinositol 3-Quinasas/metabolismo , Estructura Cuaternaria de Proteína , Estructura Terciaria de Proteína , Proteínas Proto-Oncogénicas c-akt/metabolismo , Quercetina/análogos & derivados , Quercetina/metabolismo , Quercetina/farmacología , Receptor ErbB-2/antagonistas & inhibidores , Receptor ErbB-2/química , Transducción de Señal , Neoplasias Gástricas/tratamiento farmacológicoRESUMEN
The specific properties of carbon nanoparticles (NPs) have attracted great attention in applications in biotechnology and biomedicine, e.g., in the field of amyloidosis. To date, it is still indefinable whether carbon NPs would promote or inhibit the fibril formation of amyloid proteins. Here, to uncover the effects of carbon nanoparticles (NPs) including graphene and carbon nanotubes on the aggregation of prion proteins, whose misfolding and aggregation will lead to prion diseases, a ThT fluorescence assay and a molecular dynamics (MD) simulation were performed. The ThT fluorescence assay reveals that both graphene and carbon nanotubes can inhibit the fibril formation of prion proteins, especially graphene. Further MD simulation of the PrP127-147 tetramer with or without carbon NPs suggests that the interactions between prion proteins and carbon NPs reduce the aggregation tendency of PrP127-147 by decreasing the interpeptide interactions and thus inhibiting ß-sheet formation. Meanwhile, aromatic residues greatly contribute to the inhibition effects of carbon NPs by a π-π stacking interaction. The obtained results can increase our understanding on the interaction between nanoparticles and amyloid-related proteins.
Asunto(s)
Carbono/farmacología , Nanopartículas , Nanotubos de Carbono , Proteínas Priónicas/metabolismo , Agregado de Proteínas , Carbono/química , Grafito/química , Grafito/farmacología , Humanos , Simulación de Dinámica Molecular , Nanomedicina , Nanopartículas/química , Nanotubos de Carbono/química , Enfermedades por Prión/metabolismo , Enfermedades por Prión/terapia , Proteínas Priónicas/química , Agregado de Proteínas/efectos de los fármacos , Agregación Patológica de Proteínas/metabolismo , Agregación Patológica de Proteínas/terapiaRESUMEN
As co-chaperones of the 90-kDa heat shock protein(HSP90), FK506 binding protein 51 (FKBP51) and FK506 binding protein 52 (FKBP52) modulate the maturation of steroid hormone receptor through their specific FK1 domains (FKBP12-like domain 1). The inhibitors targeting FK1 domains are potential therapies for endocrine-related physiological disorders. However, the structural conservation of the FK1 domains between FKBP51 and FKBP52 make it difficult to obtain satisfactory selectivity in FK506-based drug design. Fortunately, a series of iFit ligands synthesized by Hausch et al exhibited excellent selectivity for FKBP51, providing new opportunity for design selective inhibitors. We performed molecular dynamics simulation, binding free energy calculation and unbinding pathway analysis to reveal selective mechanism for the inhibitor iFit4 binding with FKBP51 and FKBP52. The conformational stability evaluation of the "Phe67-in" and "Phe67-out" states implies that FKBP51 and FKBP52 have different preferences for "Phe67-in" and "Phe67-out" states, which we suggest as the determinant factor for the selectivity for FKBP51. The binding free energy calculations demonstrate that nonpolar interaction is favorable for the inhibitors binding, while the polar interaction and entropy contribution are adverse for the inhibitors binding. According to the results from binding free energy decomposition, the electrostatic difference of residue 85 causes the most significant thermodynamics effects on the binding of iFit4 to FKBP51 and FKBP52. Furthermore, the importance of substructure units on iFit4 were further evaluated by unbinding pathway analysis and residue-residue contact analysis between iFit4 and the proteins. The results will provide new clues for the design of selective inhibitors for FKBP51.
Asunto(s)
Simulación de Dinámica Molecular , Proteínas de Unión a Tacrolimus/química , Secuencia de Aminoácidos , Sitios de Unión , Humanos , Enlace de Hidrógeno , Ligandos , Estructura Molecular , Unión Proteica , Conformación Proteica , Isoformas de Proteínas/antagonistas & inhibidores , Isoformas de Proteínas/química , Proteínas de Unión a Tacrolimus/antagonistas & inhibidores , TermodinámicaRESUMEN
BACKGROUND: The inhibitors blocking the interaction between programmed cell death protein 1(PD-1) and programmed death-ligand 1(PD-L1) can activate the immune response of T cell and eliminate cancer cells. The crystallographic studies have provided structural insights of the interactive interfaces between PD-L1 and its protein ligands. However, the hotspot residues on PD-L1 as well as structural and energetic basis for different protein ligands still need to be further investigated. METHODS: Molecular modeling methods including molecular dynamics simulation, per-residue free energy decomposition, virtual alanine scanning mutagenesis and residue-residue contact analysis were used to qualitatively and quantitatively analyze the interactions between PD-L1 and different protein ligands. RESULTS: The results of virtual alanine scanning mutagenesis suggest that Y56, Q66, M115, D122, Y123, R125 are the hotspot residues on PD-L1. The residue-residue contact analysis further shows that PD-1 interacts with PD-L1 mainly by F and G strands while monoclonal antibodies like avelumab and BMS-936559 mainly interact with PD-L1 by CDR2 and CDR3 loops of the heavy chain. CONCLUSIONS: A structurally similar ß-hairpin peptide with 13 or 14 residues was extracted from each protein ligand and these ß-hairpin peptides were found tightly binding to the putative hotspot residues on PD-L1. GENERAL SIGNIFICANCE: This study recognizes the hotspot residues on PD-L1 and uncovers the common structural and energetic basis of different protein ligands binding to PD-L1. These results will be valuable for the design of small molecule or peptide inhibitors targeting on PD-L1.
Asunto(s)
Anticuerpos Monoclonales/farmacología , Antígeno B7-H1/química , Simulación de Dinámica Molecular , Secuencias de Aminoácidos , Aminoácidos/química , Anticuerpos Monoclonales/química , Anticuerpos Monoclonales/metabolismo , Anticuerpos Monoclonales Humanizados , Antineoplásicos Inmunológicos/química , Antineoplásicos Inmunológicos/metabolismo , Antígeno B7-H1/genética , Antígeno B7-H1/inmunología , Antígeno B7-H1/metabolismo , Cristalografía por Rayos X , Diseño de Fármacos , Humanos , Ligandos , Modelos Moleculares , Mutagénesis Sitio-Dirigida , Fragmentos de Péptidos/química , Fragmentos de Péptidos/metabolismo , Receptor de Muerte Celular Programada 1/química , Receptor de Muerte Celular Programada 1/metabolismo , Conformación Proteica , Mapeo de Interacción de ProteínasRESUMEN
The structural transition of prion proteins from a native α-helix (PrPC) to a misfolded ß-sheet-rich conformation (PrPSc) is believed to be the main cause of a number of prion diseases in humans and animals. Understanding the molecular basis of misfolding and aggregation of prion proteins will be valuable for unveiling the etiology of prion diseases. However, due to the limitation of conventional experimental techniques and the heterogeneous property of oligomers, little is known about the molecular architecture of misfolded PrPSc and the mechanism of structural transition from PrPC to PrPSc. The prion fragment 127-147 (PrP127-147) has been reported to be a critical region for PrPSc formation in Gerstmann-Straussler-Scheinker (GSS) syndrome and thus has been used as a model for the study of prion aggregation. In the present study, we employ molecular dynamics (MD) simulation techniques to study the conformational change of this fragment that could be relevant to the PrPC-PrPSc transition. Employing extensive replica exchange molecular dynamics (REMD) and conventional MD simulations, we sample a huge number of conformations of PrP127-147. Using the Markov state model (MSM), we identify the metastable conformational states of this fragment and the kinetic network of transitions between the states. The resulting MSM reveals that disordered random-coiled conformations are the dominant structures. A key metastable folded state with typical extended ß-sheet structures is identified with Pro137 being located in a turn region, consistent with a previous experimental report. Conformational analysis reveals that intrapeptide hydrophobic interaction and two key residue interactions, including Arg136-His140 and Pro137-His140, contribute a lot to the formation of ordered extended ß-sheet states. However, network pathway analysis from the most populated disordered state indicates that the formation of extended ß-sheet states is quite slow (at the millisecond level), as large structural rearrangement is needed from disordered states. We speculate that the formation process of the extended ß-sheet folded states may represent an important event during the early formation of prion oligomers and the results of our study provide insights into the molecular details of the early stage of prion aggregation.
RESUMEN
Eleven new abietane type (1â11), and one new kaurane (12), diterpenes, together with eleven known compounds (13-23), were isolated and identified from the stems of Tripterygium regelii, which has been used as a traditional folk Chinese medicine for the treatment of rheumatoid arthritis in China. The structures of new compounds were characterized by means of the interpretation of high-resolution electrospray ionization mass spectrometry (HRESIMS), extensive nuclear magnetic resonance (NMR) spectroscopic data and comparisons of their experimental CD spectra with calculated electronic circular dichroism (ECD) spectra. Compound 1 is the first abietane type diterpene with an 18â1 lactone ring. Compound 19 was isolated from the plants of the Tripterygium genus for the first time, and compounds 14-17 were isolated from T. regelii for the first time. Triregelin I (9) showed significant cytotoxicity against A2780 and HepG2 with IC50 values of 5.88 and 11.74 µM, respectively. It was found that this compound was inactive against MCF-7 cells. The discovery of these twelve new diterpenes not only provided information on chemical substances of T. regelii, but also contributed to the chemical diversity of natural terpenoids.
Asunto(s)
Abietanos/aislamiento & purificación , Diterpenos de Tipo Kaurano/aislamiento & purificación , Tallos de la Planta/química , Tripterygium/química , Abietanos/química , Abietanos/farmacología , Espectroscopía de Resonancia Magnética con Carbono-13 , Muerte Celular/efectos de los fármacos , Línea Celular Tumoral , Dicroismo Circular , Diterpenos de Tipo Kaurano/química , Diterpenos de Tipo Kaurano/farmacología , Humanos , Espectroscopía de Protones por Resonancia MagnéticaRESUMEN
The aggregation of human islet amyloid polypeptide (hIAPP) is closely related with the occurrence of type 2 diabetes (T2D). Natural flavonoid morin was confirmed to not only inhibit the amyloid formation of hIAPP, but disaggregate its preformed amyloid fibrils. In this study, with the goal of elucidating the molecular mechanism of inhibition and destabilization of morin on the full-length hIAPP(1-37) oligomer, molecular dynamics simulations were performed for hIAPP(1-37) pentamer in the presence and absence of morin. The obtained results show that during the protein-inhibitor interaction, morin can notably alter the structural properties of hIAPP(1-37) pentamer, such as morphology, solvent accessible surface area and secondary structure. Moreover, we identified three possible binding sites of morin on hIAPP, all of which located near the amyloidogenic region of this protein. From the binding free energy calculations, we found that Site II was the most possible one. Further conformational analysis together with energy decomposition showed that the residues His18, Phe23 and Ile26 play a key role in the binding with morin by hydrogen bond, π-π and hydrophobic interactions. The proposal of the theoretical mechanism of morin against hIAPP aggregation will provide valuable information for the development of new drugs to inhibit hIAPP aggregation.
Asunto(s)
Flavonoides/química , Polipéptido Amiloide de los Islotes Pancreáticos/química , Sitios de Unión , Humanos , Enlace de Hidrógeno , Interacciones Hidrofóbicas e Hidrofílicas , Polipéptido Amiloide de los Islotes Pancreáticos/metabolismo , Simulación de Dinámica Molecular , Agregado de Proteínas/fisiología , Estructura Secundaria de Proteína , Solventes/química , Electricidad Estática , TermodinámicaRESUMEN
Nine new monoterpenoid indole alkaloids, ervatamines A-I (1-9), and five known ones (10-14), were isolated from Ervatamia hainanensis. The new structures were elucidated by extensive spectroscopic analysis and comparison to known compounds. Their absolute configurations were determined by various methods including computational methods, X-ray diffraction analysis, and electronic circular dichroism spectroscopy, as well as chemical transformations. Ervatamine A (1) is a ring-C-contracted ibogan-type monoterpenoid indole alkaloid with an unusual 6/5/6/6/6 pentacyclic rearranged ring system. Ervatamines B-E (2-5) display a nitrogen-containing 9/6 ring system, which is rarely observed in nature. The epimeric ervatamines B (2) and C (3) possess a 22-nor-monoterpenoid indole alkaloid carbon skeleton, which was only found in deformylstemmadenine. Compounds 10 and 14 exhibited significant anti-inflammatory activities, with IC50 values of 25.5 and 41.5 µM, respectively, while the IC50 value of indomethacin as a positive control was found to be 42.6 µM. Additionally, compound 9 showed mild activity against 786-O and HL-60 cell lines.
Asunto(s)
Antiinflamatorios/aislamiento & purificación , Antiinflamatorios/farmacología , Apocynaceae/química , Medicamentos Herbarios Chinos/aislamiento & purificación , Medicamentos Herbarios Chinos/farmacología , Alcaloides de Triptamina Secologanina/aislamiento & purificación , Alcaloides de Triptamina Secologanina/farmacología , Antiinflamatorios/química , Medicamentos Herbarios Chinos/química , Células HL-60 , Humanos , Estructura Molecular , Alcaloides de Triptamina Secologanina/químicaRESUMEN
Channelrhodopsin-2 (ChR2) is a cationic channel protein that has been extensively studied in optogenetics. The ion channel is opened via a series of proton transfers and H-bond changes during the photocycle but the detailed mechanism is still unknown. Molecular dynamics (MD) simulations with enhanced sampling were performed on the dark-adapted state (i.e., D470) and two photocycle intermediates (P1 500 and P2 390) to study the proton transfer path of the Schiff base and the subsequent conformational changes. The results suggest there are two possible proton transfer pathways from the Schiff base to proton acceptors (i.e., E123 or D253), depending on the protonation of E90. If E90 is protonated in the P1 500 state, the proton on the Schiff base will transfer to E123. The polyene chain of 13-cis retinal tilts and opens the channel that detours the blocking central gate (CG) and forms a narrow channel through the transmembrane helices (TM) 2, 3, 6 and 7. In contrast, if E90 deprotonates after retinal isomerization, the primary proton acceptor is D253, and an almost-open channel through TM1, 2, 3 and 7 is generated. The channel diameter is very close to the experimental value. The potential mean force (PMF) suggests that the free energy is extremely low for ions passing through this channel.
RESUMEN
Conformational transition from the normal cellular form of prion protein (PrPC) to the pathogenic "scrapie" form (PrPSc) is considered to be a key event in the occurrence of prion disease. Additionally, the H2 C-terminus is widely considered to be a vital site for PrP conformational transition, which can be used as an important region to explore the potential mechanism of PrP misfolding. Therefore, to study the misfolding mechanism of PrP, 500 ns well-tempered metadynamics simulations were performed by focusing on the H2 C-terminus of PrP. For comparison, three systems were designed in total, including PrP in neutral and acidic conditions, as well as H187R mutant. The resulting free energy surfaces (FESs) obtained from metadynamics simulations reveal that acidic conditions and H187R mutation can facilitate PrP misfolding by decreasing free energy barriers for conformational transition and forming energy stable conformational states. Further analyses aimed at H2 C-terminus show that due to the increase of positive charge on residue 187 in both acidic and H187R systems, the electrostatic repulsion of residue 187 and R136/R156 increases greatly, which disrupts the electrostatic interaction network around H2 C-terminus and exposes the hydrophobic core to the solvent. Taken together, acidic conditions and H187R mutation can accelerate PrP misfolding mainly by forming more energetically stable metastable conformations with lower free energy barriers, and electrostatic network disruption involving residue 187 drives the initial misfolding of H2 C-terminus. This study provides quantitative insight into the related function of the H2 C-terminus in the PrP misfolding process, which may guide H2 C-terminus mediated drug design in the future.
Asunto(s)
Enfermedades por Prión , Priones , Humanos , Interacciones Hidrofóbicas e Hidrofílicas , Proteínas Priónicas/genética , Priones/genética , Conformación Proteica , Pliegue de ProteínaRESUMEN
Amyloid transthyretin (ATTR) amyloidosis is a widespread and fatal systemic amyloidosis characterized by the misfolding and amyloid aggregation of transthyretin (TTR). Studies suggest that dissociation of the TTR tetramer is the key step for its misfolding. Because of the importance of tetramer dissociation on ATTR amyloidosis, many TTR stabilizers have been discovered to stabilize the tetramer structure. This paper describes the application conventional molecular dynamics and metadynamics simulations to investigate the binding and unbinding mechanisms of two TTR stabilizers, including AG10 and tafamidis. AG10 has been granted an orphan drug designation by the U.S. Food and Drug Administration (FDA), and tafamidis was the first FDA-approved treatment for ATTR cardiomyopathy. The conventional molecular dynamics simulations reveal that both AG10 and tafamidis can stabilize the TTR tetramer through different mechanisms. AG10 stabilizes TTR tetramer by forming H-bonds with S117 to mimic the protective effect of T119M. Tafamidis stabilizes the tetramer by forming H-bond with S52 in the flexible CD loop to increase its structural stability. Despite the strong binding affinity of tafamidis, the free-energy surface constructed from metadynamics simulation suggests that tafamidis unbinds more readily than AG10 with lower free-energy barriers between the binding state and other intermediates. Finally, by performing pharmacophore analysis, we found two common important moieties of the studied compounds for their binding on the pockets, which can provide valuable guidance for future lead compounds' optimization in designing drugs for ATTR amyloidosis.
Asunto(s)
Neuropatías Amiloides Familiares , Prealbúmina , Benzoxazoles/farmacología , Humanos , Simulación de Dinámica MolecularRESUMEN
The microtubule-associated protein tau is critical for the development and maintenance of the nervous system. Tau dysfunction is associated with a variety of neurodegenerative diseases called tauopathies, which are characterized by neurofibrillary tangles formed by abnormally aggregated tau protein. Studying the aggregation mechanism of tau protein is of great significance for elucidating the etiology of tauopathies. The hexapeptide 306VQIVYK311 (PHF6) of R3 has been shown to play a vital role in promoting tau aggregation. In this study, long-term all-atom molecular dynamics simulations in explicit solvent were performed to investigate the mechanisms of spontaneous aggregation and template-induced misfolding of PHF6, and the dimerization at the early stage of nucleation was further specifically analyzed by the Markov state model (MSM). Our results show that PHF6 can spontaneously aggregate to form multimers enriched with ß-sheet structure and the ß-sheets in multimers prefer to exist in a parallel way. It is observed that PHF6 monomer can be induced to form a ß-sheet structure on either side of the template but in a different way. In detail, the ß-sheet structure is easier to form on the left side but does not extend well, but on the right side, the monomer can form the extended ß-sheet structure. Furthermore, MSM analysis shows that the formation of dimer mainly occurs in three steps. First, the separated monomers collide with each other at random orientations, and then a dimer with short ß-sheet structure at the N-terminal forms; finally, ß-sheets elongate to form an extended parallel ß-sheet dimer. During these processes, multiple intermediate states are identified and multiple paths can form a parallel ß-sheet dimer from the disordered coil structure. Moreover, the residues I308, V309, and Y310 play an essential role in the dimerization. In a word, our results uncover the aggregation and misfolding mechanism of PHF6 from the atomic level, which can provide useful theoretical guidance for rational design of effective therapeutic drugs against tauopathies.
Asunto(s)
Agregado de Proteínas , Agregación Patológica de Proteínas/metabolismo , Proteínas tau/química , Secuencia de Aminoácidos , Sitios de Unión , Dimerización , Humanos , Enlace de Hidrógeno , Interacciones Hidrofóbicas e Hidrofílicas , Cadenas de Markov , Microtúbulos/metabolismo , Modelos Moleculares , Simulación de Dinámica Molecular , Ovillos Neurofibrilares/metabolismo , Conformación Proteica , Dominios Proteicos , Pliegue de Proteína , Estructura Secundaria de ProteínaRESUMEN
Recently, small-molecule compounds have been reported to block the PD-1/PD-L1 interaction by inducing the dimerization of PD-L1. All these inhibitors had a common scaffold and interacted with the cavity formed by two PD-L1 monomers. This special interactive mode provided clues for the structure-based drug design, however, also showed limitations for the discovery of small-molecule inhibitors with new scaffolds. In this study, we revealed the structure-activity relationship of the current small-molecule inhibitors targeting dimerization of PD-L1 by predicting their binding and unbinding mechanism via conventional molecular dynamics and metadynamics simulation. During the binding process, the representative inhibitors (BMS-8 and BMS-1166) tended to have a more stable binding mode with one PD-L1 monomer than the other and the small-molecule inducing PD-L1 dimerization was further stabilized by the non-polar interaction of Ile54, Tyr56, Met115, Ala121, and Tyr123 on both monomers and the water bridges involved in ALys124. The unbinding process prediction showed that the PD-L1 dimerization kept stable upon the dissociation of ligands. It's indicated that the formation and stability of the small-molecule inducing PD-L1 dimerization was the key factor for the inhibitory activities of these ligands. The contact analysis, R-group based quantitative structure-activity relationship (QSAR) analysis and molecular docking further suggested that each attachment point on the core scaffold of ligands had a specific preference for pharmacophore elements when improving the inhibitory activities by structural modifications. Taken together, the results in this study could guide the structural optimization and the further discovery of novel small-molecule inhibitors targeting PD-L1.
RESUMEN
The conformational transition of prion protein (PrP) from a native form PrPC to a pathological isoform PrPSc is the main cause of a number of prion diseases in human and animals. Thus, understanding the molecular basis of conformational transition of PrP will be valuable for unveiling the etiology of PrP-related diseases. Here, to explore the potential misfolding mechanism of PrP under the acidic condition, which is known to promote PrP misfolding and trigger its aggregation, the conventional and accelerated molecular dynamics (MD) simulations combined with the Markov state model (MSM) analysis were performed. The conventional MD simulations reveal that, at an acidic pH, the globular domain of PrP is partially unfolded, particularly for the α2 C-terminus. Structural analysis of the key macrostates obtained by MSM indicates that the α2 C-terminus and the ß2-α2 loop may serve as important sites for the pH-induced PrP misfolding. Meanwhile, the α1 may also participate in the pH-induced structural conversion by moving away from the α2-α3 subdomain. Notably, dynamical network analysis of the key metastable states indicates that the protonated H187 weakens the interactions between the α2 C-terminus, α1-ß2 loop, and α2-α3 loop, leading these domains, especially the α2 C-terminus, to become unstable and to begin to misfold. Therefore, the α2 C-terminus plays a key role in the PrP misfolding process and serves as a potential site for drug targeting. Overall, our findings can deepen the understanding of the pathogenesis related to PrP and provide useful guidance for the future drug discovery.
Asunto(s)
Simulación de Dinámica Molecular , Proteínas PrPC/química , Proteínas PrPSc/química , Pliegue de Proteína , Humanos , Concentración de Iones de Hidrógeno , Cadenas de MarkovRESUMEN
The level of tau aggregation into neurofibrillary tangles, including paired helical filament (PHF) and straight filament (SF), is closely associated with Alzheimer's disease. Despite the pathological importance of misfolding and aggregation of tau, the corresponding mechanism remains unclear. Therefore, to uncover the misfolding mechanism of the tau monomer upon induction of formed PHF and SF, in this study, a conventional molecular dynamics simulation combined with a steered molecular dynamics simulation was performed to study the dissociation of the boundary chain. Interestingly, our results show that the dissociation mechanisms of the boundary chain in PHF and SF are different. In PHF, the boundary chain begins to dissociate from regions ß2 and ß3 and ends at ß8. However, in SF, it is simultaneously dissociated from ß1 and ß8 and ends at ß5. The dissociation of the boundary chain is the reverse of template-induced misfolding of the monomer. Therefore, we can deduce the misfolding mechanism of the monomer upon induction of the template. For PHF, ß8 first interacts with the template by hydrophobic interaction. Then ß7, ß6, ß5, ß4, and ß1 sequentially bind to the template by electrostatic and hydrophobic interactions. After ß1 binds to the template, ß2 and ß3 very quickly bind to the template through hydrophobic interaction. For SF, ß5 of the monomer first interacts with the template by electrostatic attraction. Then ß4 and ß6, ß3 and ß7, and ß2 and ß8 bind to the template in turn. Finally, ß1 and ß8 are fully bound to the template by hydrophobic interaction. The obtained results will be vital for understanding the earlier events during misfolding and aggregation of tau.