1
|
Ji X, Liu H, Zhang Y, Chen J, Chen HF. Personal Precise Force Field for Intrinsically Disordered and Ordered Proteins Based on Deep Learning. J Chem Inf Model 2023; 63:362-374. [PMID: 36533639 DOI: 10.1021/acs.jcim.2c01501] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Intrinsically disordered proteins (IDPs) are proteins without a fixed three-dimensional (3D) structure under physiological conditions and are associated with Parkinson's disease, Alzheimer's disease, cancer, cardiovascular disease, amyloidosis, diabetes, and other diseases. Experimental methods can hardly capture the ensemble of diverse conformations for IDPs. Molecular dynamics (MD) simulations can sample continuous conformations that might provide a valuable complement to experimental data. However, the accuracy of MD simulations depends on the quality of force field. In particular, the evolutionary conservation and coevolution of IDPs introduce that current force fields could not precisely reproduce the conformation of IDPs. In order to improve the performance of force field, deep learning and reweighting methods were used to automatically generate personal force field parameters for intrinsically disordered and ordered proteins. At first, the deep learning method predicted more accuracy φ/ψ dihedral of residue than the previous method. Then, reweighting optimized the personal force field parameters for each residue. Finally, typical representative systems such as IDPs, structure protein, and fast-folding protein were used to evaluate this force field. The results indicate that two personal force field parameters (named PPFF1 and PPFF1_af2) could better reproduce the experimental observables than ff03CMAP force field. In summary, this strategy will provide feasibility for the development of precise personal force fields.
Collapse
Affiliation(s)
- Xiaoyue Ji
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai200240, China
| | - Hao Liu
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai200240, China
| | - Yangpeng Zhang
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai200240, China
| | - Jun Chen
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai200240, China
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai200240, China.,Shanghai Center for Bioinformation Technology, Shanghai200235, China
| |
Collapse
|
2
|
Jiang Y, Chen HF. Performance evaluation of the balanced force field ff03CMAP for intrinsically disordered and ordered proteins. Phys Chem Chem Phys 2022; 24:29870-29881. [PMID: 36468450 DOI: 10.1039/d2cp04501j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
Intrinsically disordered proteins (IDPs) have been found to be closely associated with various human diseases. Because IDPs have no fixed tertiary structure under physiological conditions, current experimental methods, such as X-ray spectroscopy, NMR, and CryoEM, cannot capture all the dynamic conformations. Molecular dynamics simulation is an useful tool that is widely used to study the conformer distributions of IDPs and has become an important complementary tool for experimental methods. However, the accuracy of MD simulations directly depends on utilizing a precise force field. Recently a CMAP optimized force field based on the Amber ff03 force field (termed ff03CMAP herein) was developed for a balanced sampling of IDPs and folded proteins. In order to further evaluate the performance, more types of disordered and ordered proteins were used to test the ability for conformer sampling. The results showed that simulated chemical shifts, J-coupling, and Rg distribution with the ff03CMAP force field were in better agreement with NMR measurements and were more accurate than those with the ff03 force field. The sampling conformations by ff03CMAP were more diverse than those of ff03. At the same time, ff03CMAP could stabilize the conformers of the ordered proteins. These findings indicate that ff03CMAP can be widely used to sample diverse conformers for proteins, including the intrinsically disordered regions.
Collapse
Affiliation(s)
- Yuxin Jiang
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China.
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China. .,Shanghai Center for Bioinformation Technology, 200240, Shanghai, China
| |
Collapse
|
3
|
Tian X, Liu H, Chen HF. Catalytic mechanism of butane anaerobic oxidation for alkyl-coenzyme M reductase. Chem Biol Drug Des 2021; 98:701-712. [PMID: 34328701 DOI: 10.1111/cbdd.13931] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2021] [Revised: 07/02/2021] [Accepted: 07/24/2021] [Indexed: 12/18/2022]
Abstract
Methane is among the most potent of the greenhouse gases, which plays a key role in global climate change. As an excellent carbon and energy source, methane can be utilized by anaerobic methane oxidizing archaea and aerobic methane oxidizing bacteria. The previous work shows that an anaerobic thermophilic enrichment culture composed of dense consortia of archaea and bacteria apparently uses partly similar pathways to oxidize the C4 hydrocarbon butane. However, the catalytic mechanism of butane anaerobic oxidation for alkyl-coenzyme M reductase is still unknown. Therefore, molecular dynamics (MD) simulation was used to investigate the dynamics differences of catalytic mechanism between methane coenzyme M reductase (MCR) and alkyl-coenzyme M reductase (ACR). At first, the binding pocket of ACR is larger than that of MCR. Then, the complex of butane and ACR is more stable than that of methane and ACR. Protein conformation cloud suggests that the position of methane is dynamics and methane escapes from the binding pocket of ACR during most of the simulation time, while butane tightly binds in the pocket of ACR. The hydrophobic interactions between butane and ACR are more and stronger than those between methane and ACR. At the same time, the binding free energy between butane and ACR is significantly lower than that between methane and ACR. The dynamics correlation network indicates that the transformation of information flow for ACR-butane is smoother than that for ACR-methane. The shortest pathway for ACR-butane is from Gln144, Ala141, Hie135, Ile133, Ala160, Arg206, Asp97, Met94, Tyr347 to Phe345 with synergistic effect for two butane molecules. This study can insight into the catalytic mechanism for butane/ACR complex.
Collapse
Affiliation(s)
- Xiaopian Tian
- State Key Laboratory of Microbial Metabolism, Department of Bioinformatics and Biostatistics, SJTU-Yale Joint Center for Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, China
| | - Hao Liu
- Institute of Natural Sciences, Shanghai Jiao Tong University, Shanghai, China
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Department of Bioinformatics and Biostatistics, SJTU-Yale Joint Center for Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, China.,Shanghai Center for Bioinformation Technology, Shanghai, China
| |
Collapse
|
4
|
Cui X, Liu H, Rehman AU, Chen HF. Extensive evaluation of environment-specific force field for ordered and disordered proteins. Phys Chem Chem Phys 2021; 23:12127-12136. [PMID: 34032235 DOI: 10.1039/d1cp01385h] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Intrinsically disordered proteins (IDPs) have no fixed tertiary structure under physiological conditions and are associated with many human diseases. Because IDPs have the characteristic of possessing diverse conformations, current experimental methods cannot capture all the conformations of IDPs. However, molecular dynamics simulation can sample these atomistically diverse conformations as a valuable complement to experimental data. To accurately describe the properties of IDPs, the environment-specific precise force field (ESFF1) was successfully released to reproduce the conformer character of ordered and disordered proteins. Here, three typical IDPs and thirteen folded proteins were used to further evaluate the performance of this force field. The results indicate that the NMR observables of ESFF1 better approach experimental data than do those of ff14SB for IDPs. The sampling conformations by ESFF1 are more diverse than those of ff14SB. For folded proteins, these force fields have comparable performances for reproducing conformers. Therefore, ESFF1 can be used to reveal the model of sequence-disorder-function for IDPs.
Collapse
Affiliation(s)
- Xiaochen Cui
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China.
| | - Hao Liu
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China.
| | - Ashfaq Ur Rehman
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China.
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China. and Shanghai Center for Bioinformation Technology, Shanghai, 200235, China
| |
Collapse
|
5
|
Mondal M, Yang L, Cai Z, Patra P, Gao YQ. A perspective on the molecular simulation of DNA from structural and functional aspects. Chem Sci 2021; 12:5390-5409. [PMID: 34168783 PMCID: PMC8179617 DOI: 10.1039/d0sc05329e] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
As genetic material, DNA not only carries genetic information by sequence, but also affects biological functions ranging from base modification to replication, transcription and gene regulation through its structural and dynamic properties and variations. The motion and structural properties of DNA involved in related biological processes are also multi-scale, ranging from single base flipping to local DNA deformation, TF binding, G-quadruplex and i-motif formation, TAD establishment, compartmentalization and even chromosome territory formation, just to name a few. The sequence-dependent physical properties of DNA play vital role in all these events, and thus it is interesting to examine how simple sequence information affects DNA and the formation of the chromatin structure in these different hierarchical orders. Accordingly, molecular simulations can provide atomistic details of interactions and conformational dynamics involved in different biological processes of DNA, including those inaccessible by current experimental methods. In this perspective, which is mainly based on our recent studies, we provide a brief overview of the atomistic simulations on how the hierarchical structure and dynamics of DNA can be influenced by its sequences, base modifications, environmental factors and protein binding in the context of the protein-DNA interactions, gene regulation and structural organization of chromatin. We try to connect the DNA sequence, the hierarchical structures of DNA and gene regulation.
Collapse
Affiliation(s)
- Manas Mondal
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory 518055 Shenzhen China
| | - Lijiang Yang
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University 100871 Beijing China
| | - Zhicheng Cai
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University 100871 Beijing China.,Biomedical Pioneering Innovation Center, Peking University 100871 Beijing China
| | - Piya Patra
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory 518055 Shenzhen China .,Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University 100871 Beijing China
| | - Yi Qin Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory 518055 Shenzhen China .,Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University 100871 Beijing China.,Biomedical Pioneering Innovation Center, Peking University 100871 Beijing China.,Beijing Advanced Innovation Center for Genomics, Peking University 100871 Beijing China
| |
Collapse
|
6
|
Rahman MU, Rehman AU, Arshad T, Chen HF. Disaggregation mechanism of prion amyloid for tweezer inhibitor. Int J Biol Macromol 2021; 176:510-519. [PMID: 33607137 DOI: 10.1016/j.ijbiomac.2021.02.094] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2020] [Revised: 02/08/2021] [Accepted: 02/13/2021] [Indexed: 02/07/2023]
Abstract
The aggregation of amyloid has been an important event in the pathology of amyloidogenicity. A number of small molecules have been designed for Amyloidosis treatment. Molecular tweezer CLR01, a potential drug for misfolded β-amyloids inhibition, was reportedly bind directly to Lysine residues and interrupt oligomerization. However, the disaggregation mechanism of amyloid for this inhibitor is unclear. Here we used long timescale of molecular dynamic simulation to reveal the mechanism of disaggregation for pentamer prion amyloid. Molecular docking and molecular dynamics simulation demonstrate that CLR01 is attached with Lysine222 nitrogen by π-cation interaction of its nine aromatic rings and formation of salt bridge/hydrogen bond of one of the two rotatable peripheral anionic phosphate groups. Upon CLR01 binding, we found a major shifting occurs in initial conformation of the oligomer and stretch out the N-terminal chain A from the rest of the amyloid which seems to be the first stage of disaggregated the fibrils slowly yet efficiently. Moreover, the CLR01 remodelled the pentamer Prion220-272 into a compact structure which might be the resistant conformation for further oligomerization. Our work will contribute to better understand the interaction and deterioration mechanism of molecular tweezer for prions and similar amyloids, and offer significant insights into therapeutic development for Amyloidosis treatment.
Collapse
Affiliation(s)
- Mueed Ur Rahman
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Ashfaq Ur Rehman
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Taaha Arshad
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China; Shanghai Center for Bioinformation Technology, Shanghai 200235, China.
| |
Collapse
|
7
|
Rahman MU, Rehman AU, Liu H, Chen HF. Comparison and Evaluation of Force Fields for Intrinsically Disordered Proteins. J Chem Inf Model 2020; 60:4912-4923. [PMID: 32816485 DOI: 10.1021/acs.jcim.0c00762] [Citation(s) in RCA: 42] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Molecular dynamics (MD) simulations of six upgraded empirical force fields were compared and evaluated with short peptides, intrinsically disordered proteins, and folded proteins using trajectories of 1, 1.5, 5, or 10 μs (five replicates of 200 ns, 300 ns, 1 μs, or 2 μs) for each system. Previous studies have shown that different force fields, water models, simulation methods, and parameters can affect simulation outcomes. Here, the MD simulations were done in an explicit solvent with RS-peptide, HEWL19, HIV-rev, β amyloid (Aβ)-40, Aβ-42, phosphodiesterase-γ, CspTm, and ubiquitin using ff99IDPs, ff14IDPs, ff14IDPSFF, ff03w, CHARMM36m, and CHARMM22* force fields. The IDP ensembles generated by six all-atom empirical force fields were compared against NMR data. Despite using identical starting structures and simulation parameters, ensembles obtained with different force fields exhibit significant differences in NMR RMDs, secondary structure contents, and global properties such as the radius of gyration. The intrinsically disordered protein (IDP)-specific force fields could substantially reproduce the experimental observables in force field comparison: they have the lowest error in chemical shifts and J-couplings for short peptides/proteins, reasonably well for large IDPs and reasonably well with the radius of gyration. A high population of disorderness was observed in the IDP-specific force field for the IDP ensemble with a fraction of β sheets for β-amyloids. CHARMM22* performs better for many observables; however, it still has a preference toward the helicity for short peptides. The results of β-amyloid 42 starting from two different initial structures (Aβ421Z0Q and Aβ42model) were also compared with DSSP and NMR data. The results obtained with IDP-specific force fields within 2 μs simulation time are similar, even though starting from different structures. The current force fields perform equally well for folded proteins. The results of currently developed or modified force fields for IDPs are capable of enlightening the overall performance of the force field for disordered as well as folded proteins, thereby contributing to force field development.
Collapse
Affiliation(s)
- Mueed Ur Rahman
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Ashfaq Ur Rehman
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Hao Liu
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China.,Shanghai Center for Bioinformation Technology, Shanghai 200235, China
| |
Collapse
|
8
|
Nagarajan B, Sankaranarayanan NV, Desai UR. Rigorous analysis of free solution glycosaminoglycan dynamics using simple, new tools. Glycobiology 2020; 30:516-527. [PMID: 32080710 PMCID: PMC8179626 DOI: 10.1093/glycob/cwaa015] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2019] [Revised: 02/03/2020] [Accepted: 02/09/2020] [Indexed: 11/15/2022] Open
Abstract
Heparin/heparan sulfates (H/HS) are ubiquitous biopolymers that interact with many proteins to induce a range of biological functions. Unfortunately, how these biopolymers recognize their preferred protein targets remain poorly understood. It is suggested that computational simulations offer attractive avenues but a number of challenges, e.g., difficulty of selecting a comprehensive force field, few simple tools to interpret data, among others, remain. This work addresses several such challenges so as to help ease the implementation and analysis of computational experiments. First, this work presents a rigorous comparison of two different recent force fields, CHARMM36 and GLYCAM06, for H/HS studies. Second, it introduces two new straightforward parameters, i.e., end-to-end distance and minimum volume enclosing ellipsoid, to understand the myriad conformational forms of oligosaccharides that evolve over time in water. Third, it presents an application to elucidate the number and nature of inter and intramolecular, nondirect bridging water molecules, which help stabilize unique forms of H/HS. The results show that nonspecialists can use either CHARMM36 or GLYCAM06 force fields because both gave comparable results, albeit with small differences. The comparative study shows that the HS hexasaccharide samples a range of conformations with nearly equivalent energies, which could be the reason for its recognition by different proteins. Finally, analysis of the nondirect water bridges across the dynamics trajectory shows their importance in stabilization of certain conformational forms, which may become important for protein recognition. Overall, the work aids nonspecialists employ computational studies for understanding the solution behavior of H/HS.
Collapse
Affiliation(s)
- Balaji Nagarajan
- Institute for Structural Biology, Drug Discovery and Development, 800 E. Leigh Street, Suite 212, Richmond, VA 23219, USA
- Department of Medicinal Chemistry, 800 E. Leigh Street, Suite 205, Richmond, VA 23298, USA
| | - Nehru Viji Sankaranarayanan
- Institute for Structural Biology, Drug Discovery and Development, 800 E. Leigh Street, Suite 212, Richmond, VA 23219, USA
- Department of Medicinal Chemistry, 800 E. Leigh Street, Suite 205, Richmond, VA 23298, USA
| | - Umesh R Desai
- Institute for Structural Biology, Drug Discovery and Development, 800 E. Leigh Street, Suite 212, Richmond, VA 23219, USA
- Department of Medicinal Chemistry, 800 E. Leigh Street, Suite 205, Richmond, VA 23298, USA
| |
Collapse
|
9
|
Li J, Wu G, Fu Q, Ge H, Liu S, Li X, Cheng B. Exploring the influence of conserved lysine69 on the catalytic activity of the helicobacter pylori shikimate dehydrogenase: A combined QM/MM and MD simulations. Comput Biol Chem 2019; 83:107098. [DOI: 10.1016/j.compbiolchem.2019.107098] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2019] [Revised: 07/22/2019] [Accepted: 08/04/2019] [Indexed: 12/29/2022]
|
10
|
Hu G, Yu X, Bian Y, Cao Z, Xu S, Zhao L, Ji B, Wang W, Wang J. Atomistic Analysis of ToxN and ToxI Complex Unbinding Mechanism. Int J Mol Sci 2018; 19:E3524. [PMID: 30423909 PMCID: PMC6275071 DOI: 10.3390/ijms19113524] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2018] [Revised: 10/14/2018] [Accepted: 11/02/2018] [Indexed: 12/14/2022] Open
Abstract
ToxIN is a triangular structure formed by three protein toxins (ToxNs) and three specific noncoding RNA antitoxins (ToxIs). To respond to stimuli, ToxI is preferentially degraded, releasing the ToxN. Thus, the dynamic character is essential in the normal function interactions between ToxN and ToxI. Here, equilibrated molecular dynamics (MD) simulations were performed to study the stability of ToxN and ToxI. The results indicate that ToxI adjusts the conformation of 3' and 5' termini to bind to ToxN. Steered molecular dynamics (SMD) simulations combined with the recently developed thermodynamic integration in 3nD (TI3nD) method were carried out to investigate ToxN unbinding from the ToxIN complex. The potentials of mean force (PMFs) and atomistic pictures suggest the unbinding mechanism as follows: (1) dissociation of the 5' terminus from ToxN, (2) missing the interactions involved in the 3' terminus of ToxI without three nucleotides (G31, A32, and A33), (3) starting to unfold for ToxI, (4) leaving the binding package of ToxN for three nucleotides of ToxI, (5) unfolding of ToxI. This work provides information on the structure-function relationship at the atomistic level, which is helpful for designing new potent antibacterial drugs in the future.
Collapse
Affiliation(s)
- Guodong Hu
- Shandong Key Laboratory of Biophysics and Institutes of Biophysics, Dezhou University, Dezhou 253023, China.
| | - Xiu Yu
- Shandong Key Laboratory of Biophysics and Institutes of Biophysics, Dezhou University, Dezhou 253023, China.
| | - Yunqiang Bian
- Shandong Key Laboratory of Biophysics and Institutes of Biophysics, Dezhou University, Dezhou 253023, China.
| | - Zanxia Cao
- Shandong Key Laboratory of Biophysics and Institutes of Biophysics, Dezhou University, Dezhou 253023, China.
| | - Shicai Xu
- Shandong Key Laboratory of Biophysics and Institutes of Biophysics, Dezhou University, Dezhou 253023, China.
| | - Liling Zhao
- Shandong Key Laboratory of Biophysics and Institutes of Biophysics, Dezhou University, Dezhou 253023, China.
| | - Baohua Ji
- Shandong Key Laboratory of Biophysics and Institutes of Biophysics, Dezhou University, Dezhou 253023, China.
| | - Wei Wang
- National Laboratory of Solid State Microstructure and Department of Physics, Nanjing University, Nanjing 210093, China.
| | - Jihua Wang
- Shandong Key Laboratory of Biophysics and Institutes of Biophysics, Dezhou University, Dezhou 253023, China.
| |
Collapse
|
11
|
Collins AP, Anderson PC. Complete Coupled Binding-Folding Pathway of the Intrinsically Disordered Transcription Factor Protein Brinker Revealed by Molecular Dynamics Simulations and Markov State Modeling. Biochemistry 2018; 57:4404-4420. [PMID: 29990433 DOI: 10.1021/acs.biochem.8b00441] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Intrinsically disordered proteins (IDPs) make up a large class of proteins that lack stable structures in solution, existing instead as dynamic conformational ensembles. To perform their biological functions, many IDPs bind to other proteins or nucleic acids. Although IDPs are unstructured in solution, when they interact with binding partners, they fold into defined three-dimensional structures via coupled binding-folding processes. Because they frequently underlie IDP function, the mechanisms of this coupled binding-folding process are of great interest. However, given the flexibility inherent to IDPs and the sparse populations of intermediate states, it is difficult to reveal binding-folding pathways at atomic resolution using experimental methods. Computer simulations are another tool for studying these pathways at high resolution. Accordingly, we have applied 40 μs of unbiased molecular dynamics simulations and Markov state modeling to map the complete binding-folding pathway of a model IDP, the 59-residue C-terminal portion of the DNA binding domain of Drosophila melanogaster transcription factor Brinker (BrkDBD). Our modeling indicates that BrkDBD binds to its cognate DNA and folds in ∼50 μs by an induced fit mechanism, acquiring most of its stable secondary and tertiary structure only after it reaches the final binding site on the DNA. The protein follows numerous pathways en route to its bound and folded conformation, occasionally becoming stuck in kinetic traps. Each binding-folding pathway involves weakly bound, increasingly folded intermediate states located at different sites on the DNA surface. These findings agree with experimental data and provide additional insight into the BrkDBD folding mechanism and kinetics.
Collapse
Affiliation(s)
- Andrew P Collins
- Physical Sciences Division , University of Washington Bothell , Bothell , Washington 98011-8246 , United States
| | - Peter C Anderson
- Physical Sciences Division , University of Washington Bothell , Bothell , Washington 98011-8246 , United States
| |
Collapse
|
12
|
Nakatani K, Natsuhara N, Mori Y, Mukherjee S, Das B, Murata A. Synthesis of Naphthyridine Dimers with Conformational Restriction and Binding to DNA and RNA. Chem Asian J 2017; 12:3077-3087. [DOI: 10.1002/asia.201701293] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2017] [Revised: 09/30/2017] [Indexed: 12/11/2022]
Affiliation(s)
- Kazuhiko Nakatani
- Department of Regulatory Bioorganic Chemistry; The Institute of Scientific and Industrial Research; Osaka University; 8-1 Mihogaoka Ibaraki 567-0047 Japan
| | - Nozomi Natsuhara
- Department of Regulatory Bioorganic Chemistry; The Institute of Scientific and Industrial Research; Osaka University; 8-1 Mihogaoka Ibaraki 567-0047 Japan
| | - Yuki Mori
- Department of Regulatory Bioorganic Chemistry; The Institute of Scientific and Industrial Research; Osaka University; 8-1 Mihogaoka Ibaraki 567-0047 Japan
| | - Sanjukta Mukherjee
- Department of Regulatory Bioorganic Chemistry; The Institute of Scientific and Industrial Research; Osaka University; 8-1 Mihogaoka Ibaraki 567-0047 Japan
| | - Bimolendu Das
- Department of Regulatory Bioorganic Chemistry; The Institute of Scientific and Industrial Research; Osaka University; 8-1 Mihogaoka Ibaraki 567-0047 Japan
| | - Asako Murata
- Department of Regulatory Bioorganic Chemistry; The Institute of Scientific and Industrial Research; Osaka University; 8-1 Mihogaoka Ibaraki 567-0047 Japan
| |
Collapse
|
13
|
Chen J. Clarifying binding difference of ATP and ADP to extracellular signal-regulated kinase 2 by using molecular dynamics simulations. Chem Biol Drug Des 2016; 89:548-558. [PMID: 27696674 DOI: 10.1111/cbdd.12877] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2016] [Revised: 09/03/2016] [Accepted: 09/26/2016] [Indexed: 12/21/2022]
Abstract
Extracellular signal-regulated kinase 2 is a promising target for designs and development of anticancer drugs. Molecular dynamics simulations and molecular mechanics Poisson-Boltzmann method were applied to study binding difference of ADP and ATP to extracellular signal-regulated kinase 2. The results prove that the binding ability of ATP to extracellular signal-regulated kinase 2 is stronger than that of ADP. Principal component analysis performed by using molecular dynamics trajectories suggests that binding of ADP and ATP to extracellular signal-regulated kinase 2 change motion directions of two helices α1 and α2. Residue-based free energy decomposition method was adopted to calculate contributions of separate residues to associations of ADP and ATP with extracellular signal-regulated kinase 2. The results show that ADP and ATP produce strong CH-π interactions with five residues Ile29, Val37, Ala50, Leu105, and Leu154. In addition, five hydrogen bonding interactions of ADP and ATP with residues Lys52, Gln103, Asp104, and Met106 also stabilize bindings of ADP and ATP to extracellular signal-regulated kinase 2. Overall, the CH-π interactions of ATP with five residues Ile29, Val37, Ala50, Leu105, and Leu154 are stronger than ADP. This study is expected to contribute a significant theoretical hint for designs of anticancer drugs targeting extracellular signal-regulated kinase 2.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
14
|
Abstract
It is well-established that dynamics are central to protein function; their importance is implicitly acknowledged in the principles of the Monod, Wyman and Changeux model of binding cooperativity, which was originally proposed in 1965. Nowadays the concept of protein dynamics is formulated in terms of the energy landscape theory, which can be used to understand protein folding and conformational changes in proteins. Because protein dynamics are so important, a key to understanding protein function at the molecular level is to design experiments that allow their quantitative analysis. Nuclear magnetic resonance (NMR) spectroscopy is uniquely suited for this purpose because major advances in theory, hardware, and experimental methods have made it possible to characterize protein dynamics at an unprecedented level of detail. Unique features of NMR include the ability to quantify dynamics (i) under equilibrium conditions without external perturbations, (ii) using many probes simultaneously, and (iii) over large time intervals. Here we review NMR techniques for quantifying protein dynamics on fast (ps-ns), slow (μs-ms), and very slow (s-min) time scales. These techniques are discussed with reference to some major discoveries in protein science that have been made possible by NMR spectroscopy.
Collapse
|
15
|
Iacovelli F, Falconi M. Decoding the conformation-linked functional properties of nucleic acids by the use of computational tools. FEBS J 2015; 282:3298-310. [DOI: 10.1111/febs.13315] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2015] [Revised: 04/16/2015] [Accepted: 04/29/2015] [Indexed: 12/25/2022]
Affiliation(s)
| | - Mattia Falconi
- Department of Biology; University of Rome “Tor Vergata”; Italy
| |
Collapse
|
16
|
Binding Induced Intrinsically Disordered Protein Folding with Molecular Dynamics Simulation. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2015; 827:111-21. [DOI: 10.1007/978-94-017-9245-5_9] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
|
17
|
van der Vaart A. Coupled binding-bending-folding: The complex conformational dynamics of protein-DNA binding studied by atomistic molecular dynamics simulations. Biochim Biophys Acta Gen Subj 2014; 1850:1091-1098. [PMID: 25161164 DOI: 10.1016/j.bbagen.2014.08.009] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2014] [Revised: 08/14/2014] [Accepted: 08/18/2014] [Indexed: 12/21/2022]
Abstract
BACKGROUND Protein-DNA binding often involves dramatic conformational changes such as protein folding and DNA bending. While thermodynamic aspects of this behavior are understood, and its biological function is often known, the mechanism by which the conformational changes occur is generally unclear. By providing detailed structural and energetic data, molecular dynamics simulations have been helpful in elucidating and rationalizing protein-DNA binding. SCOPE OF REVIEW This review will summarize recent atomistic molecular dynamics simulations of the conformational dynamics of DNA and protein-DNA binding. A brief overview of recent developments in DNA force fields is given as well. MAJOR CONCLUSIONS Simulations have been crucial in rationalizing the intrinsic flexibility of DNA, and have been instrumental in identifying the sequence of binding events, the triggers for the conformational motion, and the mechanism of binding for a number of important DNA-binding proteins. GENERAL SIGNIFICANCE Molecular dynamics simulations are an important tool for understanding the complex binding behavior of DNA-binding proteins. With recent advances in force fields and rapid increases in simulation time scales, simulations will become even more important for future studies. This article is part of a Special Issue entitled Recent developments of molecular dynamics.
Collapse
Affiliation(s)
- Arjan van der Vaart
- Department of Chemistry, University of South Florida, 4202 East Fowler Avenue CHE 205, Tampa, FL 33620, USA.
| |
Collapse
|
18
|
Ye W, Yang J, Yu Q, Wang W, Hancy J, Luo R, Chen HF. Kink turn sRNA folding upon L7Ae binding using molecular dynamics simulations. Phys Chem Chem Phys 2014; 15:18510-22. [PMID: 24072031 DOI: 10.1039/c3cp53145g] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The kink-turn sRNA motif in archaea, whose combination with protein L7Ae initializes the assembly of small ribonucleoprotein particles (sRNPs), plays a key role in ribosome maturation and the translation process. Although many studies have been reported on this motif, the mechanism of sRNA folding coupled with protein binding is still poorly understood. Here, room and high temperature molecular dynamics (MD) simulations were performed on the complex of 25-nt kink-turn sRNA and L7Ae. The average RMSD values between the bound and corresponding apo structures and Kolmogorov-Smirnov P test analysis indicate that sRNA may follow an induced fit mechanism upon binding with L7Ae, both locally and globally. These conclusions are further supported by high-temperature unfolding kinetic analysis. Principal component analysis (PCA) found both closing and opening motions of the kink-turn sRNA. This might play a key role in the sRNP assembly and methylation catalysis. These combined computational methods can be used to study the specific recognition of other sRNAs and proteins.
Collapse
Affiliation(s)
- Wei Ye
- State Key Laboratory of Microbial Metabolism, Department of Bioinformatics and Biostatistics, College of Life Sciences and Biotechnology, Shanghai Jiaotong University, 800 Dongchuan Road, Shanghai, 200240, China
| | | | | | | | | | | | | |
Collapse
|
19
|
Zhang X, Liu Q, Xia T, Li N, He K, Wang C, Tan W, Fang X. Atomic force microscopy study of the effects of water-soluble fullerenes on the elasticity of living plant cells. Chem Asian J 2013; 8:2388-94. [PMID: 23929723 DOI: 10.1002/asia.201300522] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2013] [Revised: 05/13/2013] [Indexed: 11/11/2022]
Abstract
In this work, atomic force microscopy (AFM) was employed to characterize the elastic properties of a living suspension of Nicotiana tabacum L. cv. Bright Yellow (BY-2) cells and to investigate the changes in plant-cell elasticity that were induced by water-soluble C70 fullerene derivatives. The results revealed different effects of the three fullerene derivatives that had different numbers of carboxylic groups on the cell elasticity. BY-2 cells that were repressed by dimalonic-acid-modified C70 fullerenes (DiF70) and trimalonic-acid-modified C70 fullerenes (TriF70) showed a clear decrease in their Young's modulus. However, the Young's modulus of cells that were treated with tetramalonic-acid-modified C70 fullerenes (TetraF70) increased. Disruption of the actin cytoskeleton arrangement was observed following treatment with DiF70 and TriF70, but not with TetraF70. Moreover, the fullerene-induced cell-elasticity change was consistent with the change in cell-proliferation rate. This work provides a new approach and valuable information for the study of the biological effect of nanomaterials on plant cells.
Collapse
Affiliation(s)
- Xuejie Zhang
- Beijing National Laboratory for Molecular Sciences, Key Laboratory of Molecular Nanostructures and Nanotechnology, Institute of Chemistry, Chinese Academy of Sciences, 2 Zhongguancun Beiyijie, Beijing 100190 (P. R. China), Fax: (+86) 10-62653083
| | | | | | | | | | | | | | | |
Collapse
|
20
|
Baker CM, Best RB. Insights into the Binding of Intrinsically Disordered Proteins from Molecular Dynamics Simulation. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2013; 4:182-198. [PMID: 34354764 DOI: 10.1002/wcms.1167] [Citation(s) in RCA: 46] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
Intrinsically disordered proteins (IDPs) are a class of protein that, in the native state, possess no well-defined secondary or tertiary structure, existing instead as dynamic ensembles of conformations. They are biologically important, with approximately 20% of all eukaryotic proteins disordered, and found at the heart of many biochemical networks. To fulfil their biological roles, many IDPs need to bind to proteins and/or nucleic acids. And while unstructured in solution, IDPs typically fold into a well-defined three-dimensional structure upon interaction with a binding partner. The flexibility and structural diversity inherent to IDPs makes this coupled folding and binding difficult to study at atomic resolution by experiment alone, and computer simulation currently offers perhaps the best opportunity to understand this process. But simulation of coupled folding and binding is itself extremely challenging; these molecules are large and highly flexible, and their binding partners, such as DNA or cyclins, are also often large. Therefore, their study requires either or both simplified representations and advanced enhanced sampling schemes. It is not always clear that existing simulation techniques, optimized for studying folded proteins, are well-suited to IDPs. In this article, we examine the progress that has been made in the study of coupled folding and binding using molecular dynamics simulation. We summarise what has been learnt, and examine the state of the art in terms of both methodologies and models. We also consider the lessons to be learnt from advances in other areas of simulation and highlight the issues that remain of be addressed.
Collapse
Affiliation(s)
- Christopher M Baker
- University of Cambridge, Department of Chemistry, Lensfield Road, Cambridge, CB2 1EW, UK
| | - Robert B Best
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892-0520, USA
| |
Collapse
|
21
|
Wang W, Ye W, Yu Q, Jiang C, Zhang J, Luo R, Chen HF. Conformational Selection and Induced Fit in Specific Antibody and Antigen Recognition: SPE7 as a Case Study. J Phys Chem B 2013; 117:4912-23. [DOI: 10.1021/jp4010967] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Wei Wang
- State Key Laboratory of Microbial
metabolism, Department of Bioinformatics and Biostatistics, College
of Life Sciences and Biotechnology, Shanghai Jiaotong University, 800 Dongchuan Road, Shanghai, 200240, China
| | - Wei Ye
- State Key Laboratory of Microbial
metabolism, Department of Bioinformatics and Biostatistics, College
of Life Sciences and Biotechnology, Shanghai Jiaotong University, 800 Dongchuan Road, Shanghai, 200240, China
| | - Qingfen Yu
- State Key Laboratory of Microbial
metabolism, Department of Bioinformatics and Biostatistics, College
of Life Sciences and Biotechnology, Shanghai Jiaotong University, 800 Dongchuan Road, Shanghai, 200240, China
| | - Cheng Jiang
- State Key Laboratory of Microbial
metabolism, Department of Bioinformatics and Biostatistics, College
of Life Sciences and Biotechnology, Shanghai Jiaotong University, 800 Dongchuan Road, Shanghai, 200240, China
| | - Jian Zhang
- Department of Pathophysiology,
Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry
of Education, School of Medicine, Shanghai Jiaotong University, 280 Chongqing Road, Shanghai, 200025, China
| | - Ray Luo
- Departments of Molecular Biology
and Biochemistry and Biomedical Engineering, University of California, Irvine, California 92697-3900, United
States
| | - Hai-Feng Chen
- State Key Laboratory of Microbial
metabolism, Department of Bioinformatics and Biostatistics, College
of Life Sciences and Biotechnology, Shanghai Jiaotong University, 800 Dongchuan Road, Shanghai, 200240, China
- Shanghai Center for Bioinformation Technology, 100 Qinzhou Road, Shanghai,
200235, China
| |
Collapse
|
22
|
Chen L, Zhang JL, Yu LY, Zheng QC, Chu WT, Xue Q, Zhang HX, Sun CC. Influence of hyperthermophilic protein Cren7 on the stability and conformation of DNA: insights from molecular dynamics simulation and free energy analysis. J Phys Chem B 2012; 116:12415-25. [PMID: 23013198 DOI: 10.1021/jp305860h] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Cren7, a novel chromatin protein highly conserved among crenarchaea, plays an important role in genome packaging and gene regulation. However, the detail dynamical structural characteristic of the Cren7-DNA complex and the detail study of the DNA in the complex have not been done. Focused on two specific Cren7-DNA complexes (PDB codes 3LWH and 3LWI ), we applied molecular dynamics (MD) simulations at four different temperatures (300, 350, 400, and 450 K) and the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) free energy calculation at 300 and 350 K to examine the role of Cren7 protein in enhancing the stability of DNA duplexes via protein-DNA interactions, and to study the structural transition in DNA. The simulation results indicate that Cren7 stabilizes DNA duplex in a certain temperature range in the binary complex compared with the unbound DNA molecules. At the same time, DNA molecules were found to undergo B-like to A-like form transitions with increased temperature. The results of statistical analyses of the H-bond and hydrophobic contacts show that some residues have significant influence on the structure of DNA molecules. Our work can give important information to understand the interactions of proteins with nucleic acids and other ligands.
Collapse
Affiliation(s)
- Lin Chen
- State Key Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, Changchun 130023, People's Republic of China
| | | | | | | | | | | | | | | |
Collapse
|
23
|
Wang YJ, Wang JF, Ping J, Yu Y, Wang Y, Lian P, Li X, Li YX, Hao P. Computational studies on the substrate interactions of influenza A virus PB2 subunit. PLoS One 2012; 7:e44079. [PMID: 22957044 PMCID: PMC3434214 DOI: 10.1371/journal.pone.0044079] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2011] [Accepted: 07/30/2012] [Indexed: 12/28/2022] Open
Abstract
Influenza virus, which spreads around the world in seasonal epidemics and leads to large numbers of deaths every year, has several ribonucleoproteins in the central core of the viral particle. These viral ribonucleoproteins can specifically bind the conserved 3' and 5' caps of the viral RNAs with responsibility for replication and transcription of the viral RNA in the nucleus of infected cells. A fundamental question of most importance is that how the cap-binding proteins in the influenza virus discriminates between capped RNAs and non-capped ones. To get an answer, we performed molecular dynamics simulations and free energy calculations on the influenza A virus PB2 subunit, an important component of the RNP complexes, with a cap analog m7GTP. Our calculations showed that some key residues in the active site, such as Arg355, His357, Glu361 as well as Gln406, could offer significant hydrogen bonding and hydrophobic interactions with the guanine ring of the cap analog m7GTP to form an aromatic sandwich mechanism for the cap recognition and positioning in the active site. Subsequently, we applied this idea to a virtual screening procedure and identified 5 potential candidates that might be inhibitors against the PB2 subunit. Interestingly, 2 candidates Cpd1 and Cpd2 have been already reported to have inhibitory activities to the influenza virus cap-binding proteins. Further calculation also showed that they had comparatively higher binding affinities to the PB2 subunit than that of m7GTP. We believed that our findings could give an atomic insight into the deeper understanding of the cap recognition and binding mechanism, providing useful information for searching or designing novel drugs against influenza viruses.
Collapse
Affiliation(s)
- Ya-Jun Wang
- Key Laboratory of Systems Biomedicine (Ministry of Education), Shanghai Jiao Tong University, Shanghai, China
| | - Jing-Fang Wang
- Key Laboratory of Systems Biomedicine (Ministry of Education), Shanghai Jiao Tong University, Shanghai, China
- Shanghai Center for Bioinformation and Technology, Shanghai, China
| | - Jie Ping
- Key Laboratory of Systems Biomedicine (Ministry of Education), Shanghai Jiao Tong University, Shanghai, China
| | - Yao Yu
- Shanghai Center for Bioinformation and Technology, Shanghai, China
- Institute of Plant Physiology and Ecology, Chinese Academy of Sciences, Shanghai, China
| | - Ying Wang
- Key Laboratory of Systems Biomedicine (Ministry of Education), Shanghai Jiao Tong University, Shanghai, China
| | - Peng Lian
- Key Laboratory of Systems Biomedicine (Ministry of Education), Shanghai Jiao Tong University, Shanghai, China
| | - Xuan Li
- Shanghai Center for Bioinformation and Technology, Shanghai, China
- Institute of Plant Physiology and Ecology, Chinese Academy of Sciences, Shanghai, China
| | - Yi-Xue Li
- Shanghai Center for Bioinformation and Technology, Shanghai, China
- Bioinformatics Center, Key Laboratory of Systems Biology, Chinese Academy of Sciences, Shanghai, China
| | - Pei Hao
- Shanghai Center for Bioinformation and Technology, Shanghai, China
- Institute of Pasteur, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai, China
| |
Collapse
|
24
|
Atomistic mechanism of microRNA translation upregulation via molecular dynamics simulations. PLoS One 2012; 7:e43788. [PMID: 22952765 PMCID: PMC3428290 DOI: 10.1371/journal.pone.0043788] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2012] [Accepted: 07/24/2012] [Indexed: 01/16/2023] Open
Abstract
MicroRNAs are endogenous 23–25 nt RNAs that play important gene-regulatory roles in animals and plants. Recently, miR369-3 was found to upregulate translation of TNFα mRNA in quiescent (G0) mammalian cell lines. Knock down and immunofluorescence experiments suggest that microRNA-protein complexes (with FXR1 and AGO2) are necessary for the translation upregulation. However the molecular mechanism of microRNA translation activation is poorly understood. In this study we constructed the microRNA-mRNA-AGO2-FXR1 quadruple complex by bioinformatics and molecular modeling, followed with all atom molecular dynamics simulations in explicit solvent to investigate the interaction mechanisms for the complex. A combined analysis of experimental and computational data suggests that AGO2-FXR1 complex relocalize microRNA:mRNA duplex to polysomes in G0. The two strands of dsRNA are then separated upon binding of AGO2 and FXR1. Finally, polysomes may improve the translation efficiency of mRNA. The mutation research confirms the stability of microRNA-mRNA-FXR1 and illustrates importance of key residue of Ile304. This possible mechanism can shed more light on the microRNA-dependent upregulation of translation.
Collapse
|
25
|
Ye W, Chen Y, Wang W, Yu Q, Li Y, Zhang J, Chen HF. Insight into the stability of cross-β amyloid fibril from VEALYL short peptide with molecular dynamics simulation. PLoS One 2012; 7:e36382. [PMID: 22590535 PMCID: PMC3349666 DOI: 10.1371/journal.pone.0036382] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2012] [Accepted: 03/30/2012] [Indexed: 12/13/2022] Open
Abstract
Amyloid fibrils are found in many fatal neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, type II diabetes, and prion disease. The VEALYL short peptide from insulin has been confirmed to aggregate amyloid-like fibrils. However, the aggregation mechanism of amyloid fibril is poorly understood. Here, we utilized molecular dynamics simulation to analyse the stability of VEALYL hexamer. The statistical results indicate that hydrophobic residues play key roles in stabilizing VEALYL hexamer. Single point and two linkage mutants confirmed that Val1, Leu4, and Tyr5 of VEALYL are key residues. The consistency of the results for the VEALYL oligomer suggests that the intermediate states might be trimer (3-0) and pentamer(3-2). These results can help us to obtain an insight into the aggregation mechanism of amyloid fibril. These methods can be used to study the stability of amyloid fibril from other short peptides.
Collapse
Affiliation(s)
- Wei Ye
- State Key Laboratory of Microbial Metabolism, Department of Bioinformatics and Biostatistics, College of Life Sciences and Biotechnology, Shanghai Jiaotong University, Shanghai, China
| | - Yue Chen
- State Key Laboratory of Microbial Metabolism, Department of Bioinformatics and Biostatistics, College of Life Sciences and Biotechnology, Shanghai Jiaotong University, Shanghai, China
| | - Wei Wang
- State Key Laboratory of Microbial Metabolism, Department of Bioinformatics and Biostatistics, College of Life Sciences and Biotechnology, Shanghai Jiaotong University, Shanghai, China
| | - Qingfen Yu
- State Key Laboratory of Microbial Metabolism, Department of Bioinformatics and Biostatistics, College of Life Sciences and Biotechnology, Shanghai Jiaotong University, Shanghai, China
| | - Yixue Li
- Shanghai Center for Bioinformation Technology, Shanghai, China
- * E-mail: (HC); (YL); (JZ)
| | - Jian Zhang
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, School of Medicine, Shanghai Jiaotong University, Shanghai, China
- * E-mail: (HC); (YL); (JZ)
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Department of Bioinformatics and Biostatistics, College of Life Sciences and Biotechnology, Shanghai Jiaotong University, Shanghai, China
- Shanghai Center for Bioinformation Technology, Shanghai, China
- * E-mail: (HC); (YL); (JZ)
| |
Collapse
|
26
|
Binolfi A, Fernández CO, Sica MP, Delfino JM, Santos J. Recognition between a short unstructured peptide and a partially folded fragment leads to the thioredoxin fold sharing native-like dynamics. Proteins 2012; 80:1448-64. [DOI: 10.1002/prot.24043] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2011] [Revised: 12/27/2011] [Accepted: 01/11/2012] [Indexed: 11/09/2022]
|
27
|
Abstract
It has been known for decades that DNA is extremely flexible and polymorphic, but our knowledge of its accessible conformational space remains limited. Structural data, primarily from X-ray diffraction studies, is sparse in comparison to the manifold configurations possible, and direct experimental examinations of DNA's flexibility still suffer from many limitations. In the face of these shortcomings, molecular dynamics (MD) is now an essential tool in the study of DNA. It affords detailed structural and dynamical insights, which explains its recent transition from a small number of highly specialized laboratories to a large variety of groups dealing with challenging biological problems. MD is now making an irreversible journey to the mainstream of research in biology, with the attendant opportunities and challenges. But given the speed with which MD studies of DNA have spread, the roots remain somewhat shallow: in many cases, there is a lack of deep knowledge about the foundations, strengths, and limits of the technique. In this Account, we discuss how MD has become the most important source of structural and flexibility data on DNA, focusing on advances since 2007 of atomistic MD in the description of DNA under near-physiological conditions and highlighting the possibilities and shortcomings of the technique. The evolution in the field over the past four years is a prelude to the ongoing revolution. The technique has gained in robustness and predictive power, which when coupled with the spectacular improvements in software and hardware has enabled the tackling of systems of increasing complexity. Simulation times of microseconds have now been achieved, with even longer times when specialized hardware is used. As a result, we have seen the first real-time simulation of large conformational transitions, including folding and unfolding of short DNA duplexes. Noteworthy advances have also been made in the study of DNA-ligand interactions, and we predict that a global thermodynamic and kinetic picture of the binding landscape of DNA will become available in a few years. MD will become a crucial tool in areas such as biomolecular engineering and synthetic biology. MD has also been shown to be an excellent source of parameters for mesoscopic models of DNA flexibility. Such models can be refined through atomistic MD simulations on small duplexes and then applied to the study of entire chromosomes. Recent evidence suggests that MD-derived elastic models can successfully predict the position of regulatory regions in DNA and can help advance our understanding of nucleosome positioning and chromatin plasticity. If these results are confirmed, MD simulations can become the ultimate tool to decipher a physical code that can contribute to gene regulation. We are entering the golden age of MD simulations of DNA. Undoubtedly, the expectations are high, but the challenges are also enormous. These include the need for more accurate potential energy functionals and for longer and more complex simulations in more realistic systems. The joint research effort of several groups will be crucial for adapting the technique to the requirements of the coming decade.
Collapse
Affiliation(s)
- Alberto Pérez
- Joint IRB-BSC Program in Computational Biology, Institute of Research in Biomedicine Barcelona, Baldiri i Reixac 10, Barcelona 08028, Spain
| | - F. Javier Luque
- Department de Fisicoquímica and Institut de Biomedicina (IBUB), Facultat de Farmàcia, Universitat de Barcelona, Avgda Diagonal 643, Barcelona 08028, Spain
| | - Modesto Orozco
- Joint IRB-BSC Program in Computational Biology, Institute of Research in Biomedicine Barcelona, Baldiri i Reixac 10, Barcelona 08028, Spain
- Departament de Bioquímica, Universitat de Barcelona, Avgda Diagonal 647, Barcelona 08028, Spain, and Instituto Nacional de Bioinformàtica, Parc Científic de Barcelona, Baldiri i Reixac 10, Barcelona 08028, Spain
| |
Collapse
|
28
|
Qin F, Ye W, Chen Y, Chen X, Li Y, Zhang J, Chen HF. Specific recognition between intrinsically disordered LEF and DNA. Phys Chem Chem Phys 2012; 14:538-45. [DOI: 10.1039/c1cp22610j] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
29
|
Boehr DD. Promiscuity in protein-RNA interactions: conformational ensembles facilitate molecular recognition in the spliceosome: conformational diversity in U2AF⁶⁵ facilitates binding to diverse RNA sequences. Bioessays 2011; 34:174-80. [PMID: 22144099 DOI: 10.1002/bies.201100152] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Here I discuss findings that suggest a universal mechanism for proteins (and RNA) to recognize and interact with various binding partners by selectively binding to different conformations that pre-exist in the free protein's conformational ensemble. The tandem RNA recognition motif domains of splicing factor U2AF⁶⁵ fluctuate in solution between a predominately closed conformation in which the RNA binding site of one of the domains is blocked, and a lowly populated open conformation in which both RNA binding pockets are accessible. RNA binding to U2AF⁶⁵ may thus occur through the weakly populated open conformation, and the binding interaction stabilizes the open conformation. The conformational diversity observed in U2AF⁶⁵ might also facilitate binding to diverse RNA sequences as found in the polypyrimidine tracts that help define 3' splice sites. Similar binding pathways in other systems have important consequences in biological regulation, molecular evolution, and information storage.
Collapse
Affiliation(s)
- David D Boehr
- Department of Chemistry, The Pennsylvania State University, 240 Chemistry Building, University Park, PA, USA.
| |
Collapse
|