Structural Communication between the E. coli Chaperones DnaK and Hsp90

The 70 kDa and 90 kDa heat shock proteins Hsp70 and Hsp90 are two abundant and highly conserved ATP-dependent molecular chaperones that participate in the maintenance of cellular homeostasis. In Escherichia coli, Hsp90 (Hsp90Ec) and Hsp70 (DnaK) directly interact and collaborate in protein remodeling. Previous work has produced a model of the direct interaction of both chaperones. The locations of the residues involved have been confirmed and the model has been validated. In this study, we investigate the allosteric communication between Hsp90Ec and DnaK and how the chaperones couple their conformational cycles. Using elastic network models (ENM), normal mode analysis (NMA), and a structural perturbation method (SPM) of asymmetric and symmetric DnaK-Hsp90Ec, we extract biologically relevant vibrations and identify residues involved in allosteric signaling. When one DnaK is bound, the dominant normal modes favor biological motions that orient a substrate protein bound to DnaK within the substrate/client binding site of Hsp90Ec and release the substrate from the DnaK substrate binding domain. The presence of one DnaK molecule stabilizes the entire Hsp90Ec protomer to which it is bound. Conversely, the symmetric model of DnaK binding results in steric clashes of DnaK molecules and suggests that the Hsp90Ec and DnaK chaperone cycles operate independently. Together, this data supports an asymmetric binding of DnaK to Hsp90Ec.     

Structural Determination of Three Different Series of Compounds as Hsp90 Inhibitors Using 3D-QSAR Modeling, Molecular Docking and Molecular Dynamics Methods

Hsp90 is involved in correcting, folding, maturation and activation of a diverse array of client proteins; it has also been implicated in the treatment of cancer in recent years. In this work, comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), molecular docking and molecular dynamics were performed on three different series of Hsp90 inhibitors to build 3D-QSAR models, which were based on the ligand-based or receptor-based methods. The optimum 3D-QSAR models exhibited reasonable statistical characteristics with averaging internal q(2) > 0.60 and external r(2) (pred) > 0.66 for Benzamide tetrahydro-4H-carbazol-4-one analogs (BT), AT13387 derivatives (AT) and Dihydroxylphenyl amides (DA). The results revealed that steric effects contributed the most to the BT model, whereas H-bonding was more important to AT, and electrostatic, hydrophobic, H-bond donor almost contributed equally to the DA model. The docking analysis showed that Asp93, Tyr139 and Thr184 in Hsp90 are important for the three series of inhibitors. Molecular dynamics simulation (MD) further indicated that the conformation derived from docking is basically consistent with the average structure extracted from MD simulation. These results not only lead to a better understanding of interactions between these inhibitors and Hsp90 receptor but also provide useful information for the design of new inhibitors with a specific activity.     

The Chaperone System in Breast Cancer: Roles and Therapeutic Prospects of the Molecular Chaperones Hsp27, Hsp60, Hsp70, and Hsp90

Breast cancer (BC) is a major public health problem, with key pieces of information needed for developing preventive and curative measures still missing. For example, the participation of the chaperone system (CS) in carcinogenesis and anti-cancer responses is poorly understood, although it can be predicted to be a crucial factor in these mechanisms. The chief components of the CS are the molecular chaperones, and here we discuss four of them, Hsp27, Hsp60, Hsp70, and Hsp90, focusing on their pro-carcinogenic roles in BC and potential for developing anti-BC therapies. These chaperones can be targets of negative chaperonotherapy, namely the elimination/blocking/inhibition of the chaperone(s) functioning in favor of BC, using, for instance, Hsp inhibitors. The chaperones can also be employed in immunotherapy against BC as adjuvants, together with BC antigens. Extracellular vesicles (EVs) in BC diagnosis and management are also briefly discussed, considering their potential as easily accessible carriers of biomarkers and as shippers of anti-cancer agents amenable to manipulation and controlled delivery. The data surveyed from many laboratories reveal that, to enhance the understanding of the role of the CS in BS pathogenesis, one must consider the CS as a physiological system, encompassing diverse members throughout the body and interacting with the ubiquitin–proteasome system, the chaperone-mediated autophagy machinery, and the immune system (IS). An integrated view of the CS, including its functional partners and considering its highly dynamic nature with EVs transporting CS components to reach all the cell compartments in which they are needed, opens as yet unexplored pathways leading to carcinogenesis that are amenable to interference by anti-cancer treatments centered on CS components, such as the molecular chaperones.     

Bacterial J-Domains with C-Terminal Tags Contact the Substrate Binding Domain of DnaK and Sequester Chaperone Activity

Functional interactions between the molecular chaperone DnaK and cofactor J-proteins (DnaJs), as well as their homologs, are crucial to the maintenance of proteostasis across cell types. In the bacterial pathogen Mycobacterium tuberculosis, DnaK–DnaJ interactions are essential for cell growth and represent potential targets for antibiotic or adjuvant development. While the N-terminal J-domains of J-proteins are known to form important contacts with DnaK, C-terminal domains have varied roles. Here, we have studied the effect of adding C-terminal tags to N-terminal J-domain truncations of mycobacterial DnaJ1 and DnaJ2 to promote additional interactions with DnaK. We found that His₆ tags uniquely promote binding to additional sites in the substrate binding domain at the C-terminus of DnaK. Other C-terminal tags attached to J-domains, even peptides known to interact with DnaK, do not produce the same effects. Expression of C-terminally modified DnaJ1 or DnaJ2 J-domains in mycobacterial cells suppresses chaperone activity following proteotoxic stress, which is exaggerated in the presence of a small-molecule DnaK inhibitor. Hence, this work uncovers genetically encodable J-protein variants that may be used to study chaperone–cofactor interactions in other organisms.     

Molecular Chaperones and Thyroid Cancer

Thyroid cancers are the most common of the endocrine system malignancies and progress must be made in the areas of differential diagnosis and treatment to improve patient management. Advances in the understanding of carcinogenic mechanisms have occurred in various fronts, including studies of the chaperone system (CS). Components of the CS are found to be quantitatively increased or decreased, and some correlations have been established between the quantitative changes and tumor type, prognosis, and response to treatment. These correlations provide the basis for identifying distinctive patterns useful in differential diagnosis and for planning experiments aiming at elucidating the role of the CS in tumorigenesis. Here, we discuss studies of the CS components in various thyroid cancers (TC). The chaperones belonging to the families of the small heat-shock proteins Hsp70 and Hsp90 and the chaperonin of Group I, Hsp60, have been quantified mostly by immunohistochemistry and Western blot in tumor and normal control tissues and in extracellular vesicles. Distinctive differences were revealed between the various thyroid tumor types. The most frequent finding was an increase in the chaperones, which can be attributed to the augmented need for chaperones the tumor cells have because of their accelerated metabolism, growth, and division rate. Thus, chaperones help the tumor cell rather than protect the patient, exemplifying chaperonopathies by mistake or collaborationism. This highlights the need for research on chaperonotherapy, namely the development of means to eliminate/inhibit pathogenic chaperones.     

Molecular Dynamics Simulations of Acylpeptide Hydrolase Bound to Chlorpyrifosmethyl Oxon and Dichlorvos

Acylpeptide hydrolases (APHs) catalyze the removal of N-acylated amino acids from blocked peptides. Like other prolyloligopeptidase (POP) family members, APHs are believed to be important targets for drug design. To date, the binding pose of organophosphorus (OP) compounds of APH, as well as the different OP compounds binding and inducing conformational changes in two domains, namely, α/β hydrolase and β-propeller, remain poorly understood. We report a computational study of APH bound to chlorpyrifosmethyl oxon and dichlorvos. In our docking study, Val471 and Gly368 are important residues for chlorpyrifosmethyl oxon and dichlorvos binding. Molecular dynamics simulations were also performed to explore the conformational changes between the chlorpyrifosmethyl oxon and dichlorvos bound to APH, which indicated that the structural feature of chlorpyrifosmethyl oxon binding in APH permitted partial opening of the β-propeller fold and allowed the chlorpyrifosmethyl oxon to easily enter the catalytic site. These results may facilitate the design of APH-targeting drugs with improved efficacy.     

Exploring Large Domain Motions in Proteins Using Atomistic Molecular Dynamics with Enhanced Conformational Sampling

Conformational transitions in multidomain proteins are essential for biological functions. The Apo conformations are typically open and flexible, while the Holo states form more compact conformations stabilized by protein-ligand interactions. Unfortunately, the atomically detailed mechanisms for such open-closed conformational changes are difficult to be accessed experimentally as well as computationally. To simulate the transitions using atomistic molecular dynamics (MD) simulations, efficient conformational sampling algorithms are required. In this work, we propose a new approach based on generalized replica-exchange with solute tempering (gREST) for exploring the open-closed conformational changes in multidomain proteins. Wherein, selected surface charged residues in a target protein are defined as the solute region in gREST simulation and the solute temperatures are different in replicas and exchanged between them to enhance the domain motions. This approach is called gREST selected surface charged residues (gREST_SSCR) and is applied to the Apo and Holo states of ribose binding protein (RBP) in solution. The conformational spaces sampled with gREST_SSCR are much wider than those with the conventional MD, sampling open-closed conformational changes while maintaining RBP domains’ stability. The free-energy landscapes of RBP in the Apo and Holo states are drawn along with twist and hinge angles of the two moving domains. The inter-domain salt-bridges that are not observed in the experimental structures are also important in the intermediate states during the conformational changes.     

Genes Identification,Molecular Docking and Dynamics Simulation Analysis of Laccases from Amylostereum areolatum Provides Molecular Basis of Laccase Bound to Lignin

An obligate mutualistic relationship exists between the fungus Amylostereum areolatum and woodwasp Sirex noctilio. The fungus digests lignin in the host pine, providing essential nutrients for the growing woodwasp larvae. However, the functional properties of this symbiosis are poorly described. In this study, we identified, cloned, and characterized 14 laccase genes from A. areolatum. These genes encoded proteins of 508 to 529 amino acids and contained three typical copper-oxidase domains, necessary to confer laccase activity. Besides, we performed molecular docking and dynamics simulation of the laccase proteins in complex with lignin compounds (monomers, dimers, trimers, and tetramers). AaLac2, AaLac3, AaLac6, AaLac8, and AaLac10 were found that had low binding energies with all lignin model compounds tested and three of them could maintain stability when binding to these compounds. Among these complexes, amino acid residues ALA, GLN, LEU, PHE, PRO, and SER were commonly present. Our study reveals the molecular basis of A. areolatum laccases interacting with lignin, which is essential for understanding how the fungus provides nutrients to S. noctilio. These findings might also provide guidance for the control of S. noctilio by informing the design of enzyme mutants that could reduce the efficiency of lignin degradation.     

Membrane Binding and Insertion of a pHLIP Peptide Studied by All-Atom Molecular Dynamics Simulations

Recent experiments in function mechanism study reported that a pH low-insertion peptide (pHLIP) can insert into a zwitterionic palmitoyloleoylphosphatidylcholine (POPC) lipid bilayer at acidic pH while binding to the bilayer surface at basic pH. However, the atomic details of the pH-dependent interaction of pHLIP with a POPC bilayer are not well understood. In this study, we investigate the detailed interactions of pHLIP with a POPC bilayer at acidic and basic pH conditions as those used in function mechanism study, using all-atom molecular dynamics (MD) simulations. Simulations have been performed by employing the initial configurations, where pHLIP is placed in aqueous solution, parallel to bilayer surface (system S), partially-inserted (system P), or fully-inserted (system F) in POPC bilayers. On the basis of multiple 200-ns MD simulations, we found (1) pHLIP in system S can spontaneously insert into a POPC bilayer at acidic pH, while binding to the membrane surface at basic pH; (2) pHLIP in system P can insert deep into a POPC bilayer at acidic pH, while it has a tendency to exit, and stays at bilayer surface at basic pH; (3) pHLIP in system F keeps in an α-helical structure at acidic pH while partially unfolding at basic pH. This study provides at atomic-level the pH-induced insertion of pHLIP into POPC bilayer.     

基于反应分子动力学的1-丁基-3-甲基咪唑硝酸盐与煤反应机理研究

使用FTT锥形量热仪对烟煤和烟煤-[Bmim][NO₃](1-丁基-3-甲基咪唑硝酸盐)混合物进行燃烧性能测试,并结合分子反应动力学方法(ReaxFF)分析1200K,1600K,1800K和2200K四个温度条件下[Bmim][NO₃]与煤的化学反应过程。实验结果表明:[Bmim][NO₃]的添加会使体系达到热释放速率(HRR)峰值的时间提前,即促进烟煤的氧化燃烧反应。模拟结果表明:当模拟温度为1200K时,仅产生了少量H₂O和NO₂;当模拟温度升高至1600K时,开始产生小分子烷烃、CO和CO₂,但数目较少且产生时间延后;当模拟温度继续升高至1800K时,开始产生少量的H₂等气体无机物,CO与CO₂的数目较1600K时增加一倍,且初次出现时间提前约100ps;当模拟温度为2200K时,反应物中的N元素转化为NH₃和HCN,CO数目继续增加,CO₂数目逐渐降低。总体来看,温度较低时反应过程主要受[NO₃]-影响,[NO₃]-通过与煤分子中羧基反应使煤分子脱氢形成得电子结构,从而易于发生分解反应;当温度继续升高,[Bmim]+开始显著参与反应过程,[Bmim]+热解产生的自由基与煤分子发生反应,生成大量NH₃和HCN,从而会对人体和环境造成危害。     

Molecular Docking and Molecular Dynamics Simulation Studies of Triterpenes from Vernonia patula with the Cannabinoid Type 1 Receptor

A molecular docking approach was employed to evaluate the binding affinity of six triterpenes, namely epifriedelanol, friedelin, α-amyrin, α-amyrin acetate, β-amyrin acetate, and bauerenyl acetate, towards the cannabinoid type 1 receptor (CB1). Molecular docking studies showed that friedelin, α-amyrin, and epifriedelanol had the strongest binding affinity towards CB1. Molecular dynamics simulation studies revealed that friedelin and α-amyrin engaged in stable non-bonding interactions by binding to a pocket close to the active site on the surface of the CB1 target protein. The studied triterpenes showed a good capacity to penetrate the blood–brain barrier. These results help to provide some evidence to justify, at least in part, the previously reported antinociceptive and sedative properties of Vernonia patula.     

NPT系综分子动力学模拟水的密度和扩散系数

采用常温常压(NPT)系综分子动力学方法模拟了298.15K、0.1013MPa条件下水的密度和自扩散系数,模拟结果与实验值一致,相对误差分别为1.74%和3.83%。     

Promotion and Inhibition of Amyloid-β Peptide Aggregation: Molecular Dynamics Studies

Aggregates of amyloid-β (Aβ) peptides are known to be related to Alzheimer’s disease. Their aggregation is enhanced at hydrophilic–hydrophobic interfaces, such as a cell membrane surface and air-water interface, and is inhibited by polyphenols, such as myricetin and rosmarinic acid. We review molecular dynamics (MD) simulation approaches of a full-length Aβ peptide, Aβ40, and Aβ(16–22) fragments in these environments. Since these peptides have both hydrophilic and hydrophobic amino acid residues, they tend to exist at the interfaces. The high concentration of the peptides accelerates the aggregation there. In addition, Aβ40 forms a β-hairpin structure, and this structure accelerates the aggregation. We also describe the inhibition mechanism of the Aβ(16–22) aggregation by polyphenols. The aggregation of Aβ(16–22) fragments is caused mainly by the electrostatic attraction between charged amino acid residues known as Lys16 and Glu22. Since polyphenols form hydrogen bonds between their hydroxy and carboxyl groups and these charged amino acid residues, they inhibit the aggregation.     

Identification of Bacterial Membrane Selectivity of Romo1-Derived Antimicrobial Peptide AMPR-22 via Molecular Dynamics

The abuse or misuse of antibiotics has caused the emergence of extensively drug-resistant (XDR) bacteria, rendering most antibiotics ineffective and increasing the mortality rate of patients with bacteremia or sepsis. Antimicrobial peptides (AMPs) are proposed to overcome this problem; however, many AMPs have attenuated antimicrobial activities with hemolytic toxicity in blood. Recently, AMPR-11 and its optimized derivative, AMPR-22, were reported to be potential candidates for the treatment of sepsis with a broad spectrum of antimicrobial activity and low hemolytic toxicity. Here, we performed molecular dynamics (MD) simulations to clarify the mechanism of lower hemolytic toxicity and higher efficacy of AMPR-22 at an atomic level. We found four polar residues in AMPR-11 bound to a model mimicking the bacterial inner/outer membranes preferentially over eukaryotic plasma membrane. AMPR-22 whose polar residues were replaced by lysine showed a 2-fold enhanced binding affinity to the bacterial membrane by interacting with bacterial specific lipids (lipid A or cardiolipin) via hydrogen bonds. The MD simulations were confirmed experimentally in models that partially mimic bacteremia conditions in vitro and ex vivo. The present study demonstrates why AMPR-22 showed low hemolytic toxicity and this approach using an MD simulation would be helpful in the development of AMPs.     

Molecular Dynamics Simulation on the Conformational Transition of the Mad2 Protein from the Open to the Closed State

The Mad2 protein, with two distinct conformations of open- and closed-states, is a key player in the spindle checkpoint. The closed Mad2 state is more active than the open one. We carried out conventional and targeted molecular dynamics simulations for the two stable Mad2 states and their conformational transition to address the dynamical transition mechanism from the open to the closed state. The intermediate structure in the transition process shows exposure of the β6 strand and an increase of space around the binding sites of β6 strand due to the unfolding of the β7/8 sheet and movement of the β6/4/5 sheet close to the αC helix. Therefore, Mad2 binding to the Cdc20 protein in the spindle checkpoint is made possible. The interconversion between these two states might facilitate the functional activity of the Mad2 protein. Motion correlation analysis revealed the allosteric network between the β1 strand and β7/8 sheet via communication of the β5-αC loop and the β6/4/5 sheet in this transition process.     

Folding of Trp-cage Mini Protein Using Temperature and Biasing Potential Replica��Exchange Molecular Dynamics Simulations

The folding process of the 20 residue Trp-cage mini-protein was investigated using standard temperature replica exchange molecular dynamics (T-RexMD) simulation and a biasing potential RexMD (BP-RexMD) method. In contrast to several conventional molecular dynamics simulations, both RexMD methods sampled conformations close to the native structure after 10–20 ns simulation time as the dominant conformational states. In contrast, to T-RexMD involving 16 replicas the BP-RexMD method achieved very similar sampling results with only five replicas. The result indicates that the BP-RexMD method is well suited to study folding processes of proteins at a significantly smaller computational cost, compared to T-RexMD. Both RexMD methods sampled not only similar final states but also agreed on the sampling of intermediate conformations during Trp-cage folding. The analysis of the sampled potential energy contributions indicated that Trp-cage folding is favored by both van der Waals and to a lesser degree electrostatic contributions. Folding does not introduce any significant sterical strain as reflected by similar energy distributions of bonded energy terms (bond length, bond angle and dihedral angle) of folded and unfolded Trp-cage structures.     

Characterization of the Interaction between Gallic Acid and Lysozyme by Molecular Dynamics Simulation and Optical Spectroscopy

The binding interaction between gallic acid (GA) and lysozyme (LYS) was investigated and compared by molecular dynamics (MD) simulation and spectral techniques. The results from spectroscopy indicate that GA binds to LYS to generate a static complex. The binding constants and thermodynamic parameters were calculated. MD simulation revealed that the main driving forces for GA binding to LYS are hydrogen bonding and hydrophobic interactions. The root-mean-square deviation verified that GA and LYS bind to form a stable complex, while the root-mean-square fluctuation results showed that the stability of the GA-LYS complex at 298 K was higher than that at 310 K. The calculated free binding energies from the molecular mechanics/Poisson-Boltzmann surface area method showed that van der Waals forces and electrostatic interactions are the predominant intermolecular forces. The MD simulation was consistent with the spectral experiments. This study provides a reference for future study of the pharmacological mechanism of GA.     

聚环氧乙烷、聚四氢呋喃及其共聚醚力学性能的MD模拟

为寻求预示复合材料力学性能的理论方法,对优选NEPE推进剂的黏合剂预聚物提供基础信息,在COMPASS力场下,取NVT系综,对聚环氧乙烷、聚四氢呋喃和它们的无规共聚醚进行了分子动力学(MD)模拟,并计算出弹性系数、模量、泊松比和柯西压等力学参数.经比较发现,聚环氧乙烷的刚性和延展性最优,共聚醚次之,聚四氢呋喃最差.     

A Combination of 3D-QSAR, Molecular Docking and Molecular Dynamics Simulation Studies of Benzimidazole-Quinolinone Derivatives as iNOS Inhibitors

Inducible Nitric Oxide Synthase (iNOS) has been involved in a variety of diseases, and thus it is interesting to discover and optimize new iNOS inhibitors. In previous studies, a series of benzimidazole-quinolinone derivatives with high inhibitory activity against human iNOS were discovered. In this work, three-dimensional quantitative structure-activity relationships (3D-QSAR), molecular docking and molecular dynamics (MD) simulation approaches were applied to investigate the functionalities of active molecular interaction between these active ligands and iNOS. A QSAR model with R² of 0.9356, Q² of 0.8373 and Pearson-R value of 0.9406 was constructed, which presents a good predictive ability in both internal and external validation. Furthermore, a combined analysis incorporating the obtained model and the MD results indicates: (1) compounds with the proper-size hydrophobic substituents at position 3 in ring-C (R₃ substituent), hydrophilic substituents near the X₆ of ring-D and hydrophilic or H-bond acceptor groups at position 2 in ring-B show enhanced biological activities; (2) Met368, Trp366, Gly365, Tyr367, Phe363, Pro344, Gln257, Val346, Asn364, Met349, Thr370, Glu371 and Tyr485 are key amino acids in the active pocket, and activities of iNOS inhibitors are consistent with their capability to alter the position of these important residues, especially Glu371 and Thr370. The results provide a set of useful guidelines for the rational design of novel iNOS inhibitors.