1
|
Tolokh IS, Folescu DE, Onufriev AV. Inclusion of Water Multipoles into the Implicit Solvation Framework Leads to Accuracy Gains. J Phys Chem B 2024; 128:5855-5873. [PMID: 38860842 PMCID: PMC11194828 DOI: 10.1021/acs.jpcb.4c00254] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2024] [Revised: 05/28/2024] [Accepted: 05/29/2024] [Indexed: 06/12/2024]
Abstract
The current practical "workhorses" of the atomistic implicit solvation─the Poisson-Boltzmann (PB) and generalized Born (GB) models─face fundamental accuracy limitations. Here, we propose a computationally efficient implicit solvation framework, the Implicit Water Multipole GB (IWM-GB) model, that systematically incorporates the effects of multipole moments of water molecules in the first hydration shell of a solute, beyond the dipole water polarization already present at the PB/GB level. The framework explicitly accounts for coupling between polar and nonpolar contributions to the total solvation energy, which is missing from many implicit solvation models. An implementation of the framework, utilizing the GAFF force field and AM1-BCC atomic partial charges model, is parametrized and tested against the experimental hydration free energies of small molecules from the FreeSolv database. The resulting accuracy on the test set (RMSE ∼ 0.9 kcal/mol) is 12% better than that of the explicit solvation (TIP3P) treatment, which is orders of magnitude slower. We also find that the coupling between polar and nonpolar parts of the solvation free energy is essential to ensuring that several features of the IWM-GB model are physically meaningful, including the sign of the nonpolar contributions.
Collapse
Affiliation(s)
- Igor S. Tolokh
- Department
of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Dan E. Folescu
- Department
of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
- Department
of Mathematics, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Alexey V. Onufriev
- Department
of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
- Department
of Physics, Virginia Tech, Blacksburg, Virginia 24061, United States
- Center
for Soft Matter and Biological Physics, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
2
|
Bass L, Elder LH, Folescu DE, Forouzesh N, Tolokh IS, Karpatne A, Onufriev AV. Improving the Accuracy of Physics-Based Hydration-Free Energy Predictions by Machine Learning the Remaining Error Relative to the Experiment. J Chem Theory Comput 2024; 20:396-410. [PMID: 38149593 PMCID: PMC10950260 DOI: 10.1021/acs.jctc.3c00981] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2023]
Abstract
The accuracy of computational models of water is key to atomistic simulations of biomolecules. We propose a computationally efficient way to improve the accuracy of the prediction of hydration-free energies (HFEs) of small molecules: the remaining errors of the physics-based models relative to the experiment are predicted and mitigated by machine learning (ML) as a postprocessing step. Specifically, the trained graph convolutional neural network attempts to identify the "blind spots" in the physics-based model predictions, where the complex physics of aqueous solvation is poorly accounted for, and partially corrects for them. The strategy is explored for five classical solvent models representing various accuracy/speed trade-offs, from the fast analytical generalized Born (GB) to the popular TIP3P explicit solvent model; experimental HFEs of small neutral molecules from the FreeSolv set are used for the training and testing. For all of the models, the ML correction reduces the resulting root-mean-square error relative to the experiment for HFEs of small molecules, without significant overfitting and with negligible computational overhead. For example, on the test set, the relative accuracy improvement is 47% for the fast analytical GB, making it, after the ML correction, almost as accurate as uncorrected TIP3P. For the TIP3P model, the accuracy improvement is about 39%, bringing the ML-corrected model's accuracy below the 1 kcal/mol threshold. In general, the relative benefit of the ML corrections is smaller for more accurate physics-based models, reaching the lower limit of about 20% relative accuracy gain compared with that of the physics-based treatment alone. The proposed strategy of using ML to learn the remaining error of physics-based models offers a distinct advantage over training ML alone directly on reference HFEs: it preserves the correct overall trend, even well outside of the training set.
Collapse
Affiliation(s)
- Lewis Bass
- Department of Computer Engineering, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Luke H Elder
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Dan E Folescu
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
- Department of Mathematics, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Negin Forouzesh
- Department of Computer Science, California State University, Los Angeles, California 90032, United States
| | - Igor S Tolokh
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Anuj Karpatne
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Alexey V Onufriev
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
- Department of Physics, Virginia Tech, Blacksburg, Virginia 24061, United States
- Center for Soft Matter and Biological Physics, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
3
|
Zhang J, Zhang M, Wang Q, Wen H, Liu Z, Wang F, Wang Y, Yao F, Song N, Kou Z, Li Y, Guo F, Zhu S. Distinct structure and gating mechanism in diverse NMDA receptors with GluN2C and GluN2D subunits. Nat Struct Mol Biol 2023; 30:629-639. [PMID: 36959261 DOI: 10.1038/s41594-023-00959-z] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2022] [Accepted: 03/01/2023] [Indexed: 03/25/2023]
Abstract
N-methyl-D-aspartate (NMDA) receptors are heterotetramers comprising two GluN1 and two alternate GluN2 (N2A-N2D) subunits. Here we report full-length cryo-EM structures of the human N1-N2D di-heterotetramer (di-receptor), rat N1-N2C di-receptor and N1-N2A-N2C tri-heterotetramer (tri-receptor) at a best resolution of 3.0 Å. The bilobate N-terminal domain (NTD) in N2D intrinsically adopts a closed conformation, leading to a compact NTD tetramer in the N1-N2D receptor. Additionally, crosslinking the ligand-binding domain (LBD) of two N1 protomers significantly elevated the channel open probability (Po) in N1-N2D di-receptors. Surprisingly, the N1-N2C di-receptor adopted both symmetric (minor) and asymmetric (major) conformations, the latter further locked by an allosteric potentiator, PYD-106, binding to a pocket between the NTD and LBD in only one N2C protomer. Finally, the N2A and N2C subunits in the N1-N2A-N2C tri-receptor display a conformation close to one protomer in the N1-N2A and N1-N2C di-receptors, respectively. These findings provide a comprehensive structural understanding of diverse function in major NMDA receptor subtypes.
Collapse
Affiliation(s)
- Jilin Zhang
- Institute of Neuroscience, State Key Laboratory of Neuroscience, CAS Center for Excellence in Brain Science and Intelligence Technology, Chinese Academy of Sciences, Shanghai, China
- University of Chinese Academy of Sciences, Beijing, China
| | - Ming Zhang
- University of Chinese Academy of Sciences, Beijing, China
- Center for Neurological and Psychiatric Research and Drug Discovery, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
| | | | - Han Wen
- DP Technology, Beijing, China
| | - Zheyi Liu
- CAS Key Laboratory of Separation Sciences for Analytical Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, China
| | - Fangjun Wang
- CAS Key Laboratory of Separation Sciences for Analytical Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, China
| | | | - Fenyong Yao
- Institute of Neuroscience, State Key Laboratory of Neuroscience, CAS Center for Excellence in Brain Science and Intelligence Technology, Chinese Academy of Sciences, Shanghai, China
| | - Nan Song
- Institute of Neuroscience, State Key Laboratory of Neuroscience, CAS Center for Excellence in Brain Science and Intelligence Technology, Chinese Academy of Sciences, Shanghai, China
| | - Zengwei Kou
- Institute of Neuroscience, State Key Laboratory of Neuroscience, CAS Center for Excellence in Brain Science and Intelligence Technology, Chinese Academy of Sciences, Shanghai, China
| | - Yang Li
- University of Chinese Academy of Sciences, Beijing, China
- Center for Neurological and Psychiatric Research and Drug Discovery, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
| | - Fei Guo
- University of Chinese Academy of Sciences, Beijing, China
- Center for Neurological and Psychiatric Research and Drug Discovery, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
| | - Shujia Zhu
- Institute of Neuroscience, State Key Laboratory of Neuroscience, CAS Center for Excellence in Brain Science and Intelligence Technology, Chinese Academy of Sciences, Shanghai, China.
- University of Chinese Academy of Sciences, Beijing, China.
| |
Collapse
|
4
|
Parkin D, Takano M. The Intrinsic Radius as a Key Parameter in the Generalized Born Model to Adjust Protein-Protein Electrostatic Interaction. Int J Mol Sci 2023; 24:ijms24054700. [PMID: 36902130 PMCID: PMC10003186 DOI: 10.3390/ijms24054700] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2023] [Revised: 01/30/2023] [Accepted: 02/19/2023] [Indexed: 03/05/2023] Open
Abstract
The generalized Born (GB) model is an extension of the continuum dielectric theory of Born solvation energy and is a powerful method for accelerating the molecular dynamic (MD) simulations of charged biological molecules in water. While the effective dielectric constant of water that varies as a function of the separation distance between solute molecules is incorporated into the GB model, adjustment of the parameters is indispensable for accurate calculation of the Coulomb (electrostatic) energy. One of the key parameters is the lower limit of the spatial integral of the energy density of the electric field around a charged atom, known as the intrinsic radius ρ. Although ad hoc adjustment of ρ has been conducted to improve the Coulombic (ionic) bond stability, the physical mechanism by which ρ affects the Coulomb energy remains unclear. Via energetic analysis of three differently sized systems, here, we clarify that the Coulomb bond stability increases with increasing ρ and that the increased stability is caused by the interaction energy term, not by the self-energy (desolvation energy) term, as was supposed previously. Our results suggest that the use of larger values for the intrinsic radii of hydrogen and oxygen atoms, together with the use of a relatively small value for the spatial integration cutoff in the GB model, can better reproduce the Coulombic attraction between protein molecules.
Collapse
Affiliation(s)
- Dan Parkin
- Research Institute for Science and Engineering, Waseda University, Okubo 3-4-1, Sinjuku-ku, Tokyo 169-8555, Japan
| | - Mitsunori Takano
- Research Institute for Science and Engineering, Waseda University, Okubo 3-4-1, Sinjuku-ku, Tokyo 169-8555, Japan
- Department of Pure and Applied Physics, Waseda University, Okubo 3-4-1, Sinjuku-ku, Tokyo 169-8555, Japan
- Correspondence: ; Tel.: +81-3-5286-3512
| |
Collapse
|
5
|
Molecular Insights into Substrate Binding of the Outer Membrane Enzyme OmpT. Catalysts 2023. [DOI: 10.3390/catal13020214] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
The enzyme OmpT of the outer membrane of Escherichia coli shows proteolytic activity and cleaves peptides and proteins. Using molecular dynamics simulations in a fully hydrated lipid bilayer on a time scale of hundreds of nanoseconds, we draw a detailed atomic picture of substrate recognition in the OmpT-holo enzyme complex. Hydrogen bonds and salt bridges are essential for maintaining the integrity of the active site and play a central role for OmpT in recognizing its substrate. Electrostatic interactions are critical at all stages from approaching the substrate to docking at the active site. Computational alanine scanning based on the Molecular Mechanics Generalized Born Surface Area (MM-GBSA) approach confirms the importance of multiple residues in the active site that form salt bridges. The substrate fluctuates along the axis of the β-barrel, which is associated with oscillations of the binding cleft formed by the residue pairs D210-H212 and D83-D85. Principal component analysis suggests that substrate and protein movements are correlated. We observe the transient presence of putative catalytic water molecules near the active site, which may be involved in the nucleophilic attack on the cleavable peptide bond of the substrate.
Collapse
|
6
|
Wang K, Li B, Ge L, Xie Y. Molecular insight into the systematic affinity and selectivity of partner recognition sites between the WW1 and WW2 domains of human KIBRA neuroprotein. J Mol Graph Model 2022; 116:108258. [PMID: 35810735 DOI: 10.1016/j.jmgm.2022.108258] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2021] [Revised: 05/22/2022] [Accepted: 06/21/2022] [Indexed: 12/14/2022]
Abstract
Human KIBRA, a member of the WWC family proteins, is an upstream regulator of the Salvador/Warts/Hippo (SWH) signaling pathway and predominately expressed in nervous system. The protein has two functionally regulatory domains WW1 and WW2 at N-terminal region, which recognize and bind to the PY-motif segments of their partner proteins to serve as a signaling scaffold role in the SWH pathway. The two domains are highly conserved, but their downstream ligands and biological functions may not be fully consistent. In this study, we attempted to systematically profile the PY-motif affinity to and selectivity between KIBRA WW1 and WW2 domains involved in partner recognition sites. Ontology mining was used to enrich the KIBRA-interacting proteins in literature libraries, from which a variety of PY-motif peptide segments were identified, and their binding behavior to each domain was then analyzed by integrating computational modeling and experimental assay. Most PY-motif peptides were found to interact potently with WW1 and WW2, but they generally only exhibit a moderate or modest selectivity between the two domains. Subsequently, several representative peptides were further examined in detail to elucidate the molecular mechanism underlying their affinity and selectivity. It is revealed that the middle motif region of PY-motif peptides is primarily responsible for the affinity and stability of peptide binding, but only contributes marginally to peptide selectivity. Instead, the N-terminal region and, particularly, C-terminal region of PY-motif peptides play a crucial role in the selectivity. Hydrophobic contacts and hydrogen bonds confer stability and specificity to the domain-peptide interaction, respectively.
Collapse
Affiliation(s)
- Kai Wang
- Department of Anesthesiology, Rizhao People's Hospital, Affiliated to Jining Medical University, Rizhao, 276827, China
| | - Baoqiang Li
- Department of Anesthesiology, Rizhao People's Hospital, Affiliated to Jining Medical University, Rizhao, 276827, China
| | - Lei Ge
- Department of Emergency, Rizhao People's Hospital, Affiliated to Jining Medical University, Rizhao, 276827, China
| | - Yi Xie
- Department of Anesthesiology, Zibo Central Hospital, Affiliated to Binzhou Medical University, Zibo, 255020, China.
| |
Collapse
|
7
|
Cong Y, Li M, Qi Y, Zhang JZH. A fast-slow method to treat solute dynamics in explicit solvent. Phys Chem Chem Phys 2022; 24:14498-14510. [PMID: 35665790 DOI: 10.1039/d2cp00732k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Aiming to reduce the computational cost in the current explicit solvent molecular dynamics (MD) simulation, this paper proposes a fast-slow method for the fast MD simulation of biomolecules in explicit solvent. This fast-slow method divides the entire system into two parts: a core layer (typically solute or biomolecule) and a peripheral layer (typically solvent molecules). The core layer is treated using standard MD method but the peripheral layer is treated by a slower dynamics method to reduce the computational cost. We compared four different simulation models in testing calculations for several small proteins. These include gas-phase, implicit solvent, fast-slow explicit solvent and standard explicit solvent MD simulations. Our study shows that gas-phase and implicit solvent models do not provide a realistic solvent environment and fail to correctly produce reliable dynamic structures of proteins. On the other hand, the fast-slow method can essentially reproduce the same solvent effect as the standard explicit solvent model while gaining an order of magnitude in efficiency. This fast-slow method thus provides an efficient approach for accelerating the MD simulation of biomolecules in explicit solvent.
Collapse
Affiliation(s)
- Yalong Cong
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Key Laboratory of Green Chemistry & Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University at Shanghai, 200062, China.
| | - Mengxin Li
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Key Laboratory of Green Chemistry & Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University at Shanghai, 200062, China.
| | - Yifei Qi
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Key Laboratory of Green Chemistry & Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University at Shanghai, 200062, China.
| | - John Z H Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Key Laboratory of Green Chemistry & Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University at Shanghai, 200062, China. .,CAS Key Laboratory of Quantitative Engineering Biology, Shenzhen Institute of Synthetic Biology, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, Guangdong, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China.,Department of Chemistry, New York University, NY, NY 10003, USA.,Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi, 030006, China
| |
Collapse
|
8
|
Zhou Y, Jiang Y, Chen SJ. RNA-ligand molecular docking: advances and challenges. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2022; 12:e1571. [PMID: 37293430 PMCID: PMC10250017 DOI: 10.1002/wcms.1571] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/26/2021] [Accepted: 07/20/2021] [Indexed: 12/16/2022]
Abstract
With rapid advances in computer algorithms and hardware, fast and accurate virtual screening has led to a drastic acceleration in selecting potent small molecules as drug candidates. Computational modeling of RNA-small molecule interactions has become an indispensable tool for RNA-targeted drug discovery. The current models for RNA-ligand binding have mainly focused on the docking-and-scoring method. Accurate docking and scoring should tackle four crucial problems: (1) conformational flexibility of ligand, (2) conformational flexibility of RNA, (3) efficient sampling of binding sites and binding poses, and (4) accurate scoring of different binding modes. Moreover, compared with the problem of protein-ligand docking, predicting ligand binding to RNA, a negatively charged polymer, is further complicated by additional effects such as metal ion effects. Thermodynamic models based on physics-based and knowledge-based scoring functions have shown highly encouraging success in predicting ligand binding poses and binding affinities. Recently, kinetic models for ligand binding have further suggested that including dissociation kinetics (residence time) in ligand docking would result in improved performance in estimating in vivo drug efficacy. More recently, the rise of deep-learning approaches has led to new tools for predicting RNA-small molecule binding. In this review, we present an overview of the recently developed computational methods for RNA-ligand docking and their advantages and disadvantages.
Collapse
Affiliation(s)
- Yuanzhe Zhou
- Department of Physics and Astronomy, Department of Biochemistry, Institute of Data Sciences and Informatics, University of Missouri, Columbia, MO 65211-7010, USA
| | - Yangwei Jiang
- Department of Physics and Astronomy, Department of Biochemistry, Institute of Data Sciences and Informatics, University of Missouri, Columbia, MO 65211-7010, USA
| | - Shi-Jie Chen
- Department of Physics and Astronomy, Department of Biochemistry, Institute of Data Sciences and Informatics, University of Missouri, Columbia, MO 65211-7010, USA
| |
Collapse
|
9
|
Golyshev VM, Pyshnyi DV, Lomzov AA. Calculation of Energy for RNA/RNA and DNA/RNA Duplex Formation by Molecular Dynamics Simulation. Mol Biol 2021. [DOI: 10.1134/s002689332105006x] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
Abstract
The development of approaches for predictive calculation of hybridization properties of various nucleic acid (NA) derivatives is the basis for the rational design of the NA-based constructs. Modern advances in computer modeling methods provide the feasibility of these calculations. We have analyzed the possibility of calculating the energy of DNA/RNA and RNA/RNA duplex formation using representative sets of complexes (65 and 75 complexes, respectively). We used the classical molecular dynamics (MD) method, the MMPBSA or MMGBSA approaches to calculate the enthalpy (ΔH°) component, and the quasi-harmonic approximation (Q-Harm) or the normal mode analysis (NMA) methods to calculate the entropy (ΔS°) contribution to the Gibbs energy ($$\Delta G_{{37}}^{^\circ }$$ ) of the NA complex formation. We have found that the MMGBSA method in the analysis of the MD trajectory of only the NA duplex and the empirical linear approximation allow calculation of the enthalpy of formation of the DNA, RNA, and hybrid duplexes of various lengths and GC content with an accuracy of 8.6%. Within each type of complex, the combination of rather efficient MMGBSA and Q-Harm approaches being applied to the trajectory of only the bimolecular complex makes it possible to calculate the $$\Delta G_{{37}}^{^\circ }$$ of the duplex formation with an error value of 10%. The high accuracy of predictive calculation for different types of natural complexes (DNA/RNA, DNA/RNA, and RNA/RNA) indicates the possibility of extending the considered approach to analogs and derivatives of nucleic acids, which gives a fundamental opportunity in the future to perform rational design of new types of NA-targeted sequence-specific compounds.
Collapse
|
10
|
Zannoni A, Pelliciari S, Musiani F, Chiappori F, Roncarati D, Scarlato V. Definition of the Binding Architecture to a Target Promoter of HP1043, the Essential Master Regulator of Helicobacter pylori. Int J Mol Sci 2021; 22:ijms22157848. [PMID: 34360614 PMCID: PMC8345958 DOI: 10.3390/ijms22157848] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2021] [Revised: 07/13/2021] [Accepted: 07/14/2021] [Indexed: 11/30/2022] Open
Abstract
HP1043 is an essential orphan response regulator of Helicobacter pylori orchestrating multiple crucial cellular processes. Classified as a member of the OmpR/PhoB family of two-component systems, HP1043 exhibits a highly degenerate receiver domain and evolved to function independently of phosphorylation. Here, we investigated the HP1043 binding mode to a target sequence in the hp1227 promoter (Php1227). Scanning mutagenesis of HP1043 DNA-binding domain and consensus sequence led to the identification of residues relevant for the interaction of the protein with a target DNA. These determinants were used as restraints to guide a data-driven protein-DNA docking. Results suggested that, differently from most other response regulators of the same family, HP1043 binds in a head-to-head conformation to the Php1227 target promoter. HP1043 interacts with DNA largely through charged residues and contacts with both major and minor grooves of the DNA are required for a stable binding. Computational alanine scanning on molecular dynamics trajectory was performed to corroborate our findings. Additionally, in vitro transcription assays confirmed that HP1043 positively stimulates the activity of RNA polymerase.
Collapse
Affiliation(s)
- Annamaria Zannoni
- Department of Pharmacy and Biotechnology (FaBiT), University of Bologna, 40126 Bologna, Italy; (A.Z.); (S.P.); (F.M.)
| | - Simone Pelliciari
- Department of Pharmacy and Biotechnology (FaBiT), University of Bologna, 40126 Bologna, Italy; (A.Z.); (S.P.); (F.M.)
| | - Francesco Musiani
- Department of Pharmacy and Biotechnology (FaBiT), University of Bologna, 40126 Bologna, Italy; (A.Z.); (S.P.); (F.M.)
| | - Federica Chiappori
- Istituto di Tecnologie Biomediche-Consiglio Nazionale delle Ricerche (ITB-CNR), 20054 Segrate, Italy;
| | - Davide Roncarati
- Department of Pharmacy and Biotechnology (FaBiT), University of Bologna, 40126 Bologna, Italy; (A.Z.); (S.P.); (F.M.)
- Correspondence: (D.R.); (V.S.)
| | - Vincenzo Scarlato
- Department of Pharmacy and Biotechnology (FaBiT), University of Bologna, 40126 Bologna, Italy; (A.Z.); (S.P.); (F.M.)
- Correspondence: (D.R.); (V.S.)
| |
Collapse
|
11
|
Nguyen TB, Lane DP, Verma CS. Can Glycosylation Mask the Detection of MHC Expressing p53 Peptides by T Cell Receptors? Biomolecules 2021; 11:biom11071056. [PMID: 34356680 PMCID: PMC8301869 DOI: 10.3390/biom11071056] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2021] [Revised: 06/28/2021] [Accepted: 07/07/2021] [Indexed: 12/02/2022] Open
Abstract
Proteins of the major histocompatibility complex (MHC) class I, or human leukocyte antigen (HLA) in humans interact with endogenous peptides and present them to T cell receptors (TCR), which in turn tune the immune system to recognize and discriminate between self and foreign (non-self) peptides. Of especial importance are peptides derived from tumor-associated antigens. T cells recognizing these peptides are found in cancer patients, but not in cancer-free individuals. What stimulates this recognition, which is vital for the success of checkpoint based therapy? A peptide derived from the protein p53 (residues 161–169 or p161) was reported to show this behavior. T cells recognizing this unmodified peptide could be further stimulated in vitro to create effective cancer killing CTLs (cytotoxic T lymphocytes). We hypothesize that the underlying difference may arise from post-translational glycosylation of p161 in normal individuals, likely masking it against recognition by TCR. Defects in glycosylation in cancer cells may allow the presentation of the native peptide. We investigate the structural consequences of such peptide glycosylation by investigating the associated structural dynamics.
Collapse
Affiliation(s)
- Thanh Binh Nguyen
- Division of Biomolecular Structure to Mechanism, Bioinformatics Institute, Agency for Science, Technology and Research (A*STAR), Singapore 138671, Singapore;
| | - David P. Lane
- p53 Laboratory, Agency for Science, Technology and Research (A*STAR), Singapore 138648, Singapore;
| | - Chandra S. Verma
- Division of Biomolecular Structure to Mechanism, Bioinformatics Institute, Agency for Science, Technology and Research (A*STAR), Singapore 138671, Singapore;
- School of Biological Sciences, College of Science, Nanyang Technological University, Singapore 637551, Singapore
- Department of Biological Sciences, Faculty of Science, National University of Singapore, Singapore 117543, Singapore
- Correspondence: ; Tel.: +65-6478-8273; Fax: +65-6478-9048
| |
Collapse
|
12
|
Zou T, Woodrum BW, Halloran N, Campitelli P, Bobkov AA, Ghirlanda G, Ozkan SB. Local Interactions That Contribute Minimal Frustration Determine Foldability. J Phys Chem B 2021; 125:2617-2626. [PMID: 33687216 DOI: 10.1021/acs.jpcb.1c00364] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Earlier experiments suggest that the evolutionary information (conservation and coevolution) encoded in protein sequences is necessary and sufficient to specify the fold of a protein family. However, there is no computational work to quantify the effect of such evolutionary information on the folding process. Here we explore the role of early folding steps for sequences designed using coevolution and conservation through a combination of computational and experimental methods. We simulated a repertoire of native and designed WW domain sequences to analyze early local contact formation and found that the N-terminal β-hairpin turn would not form correctly due to strong non-native local contacts in unfoldable sequences. Through a maximum likelihood approach, we identified five local contacts that play a critical role in folding, suggesting that a small subset of amino acid pairs can be used to solve the "needle in the haystack" problem to design foldable sequences. Thus, using the contact probability of those five local contacts that form during the early stage of folding, we built a classification model that predicts the foldability of a WW sequence with 81% accuracy. This classification model was used to redesign WW domain sequences that could not fold due to frustration and make them foldable by introducing a few mutations that led to the stabilization of these critical local contacts. The experimental analysis shows that a redesigned sequence folds and binds to polyproline peptides with a similar affinity as those observed for native WW domains. Overall, our analysis shows that evolutionary-designed sequences should not only satisfy the folding stability but also ensure a minimally frustrated folding landscape.
Collapse
Affiliation(s)
- Taisong Zou
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, Arizona 85287, United States
| | - Brian W Woodrum
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287, United States
| | - Nicholas Halloran
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287, United States
| | - Paul Campitelli
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, Arizona 85287, United States
| | - Andrey A Bobkov
- Conrad Prebys Center for Chemical Genomics, Sanford Burnham Prebys Medical Discovery Institute, La Jolla, California 92037, United States
| | - Giovanna Ghirlanda
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287, United States
| | - Sefika Banu Ozkan
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, Arizona 85287, United States
| |
Collapse
|
13
|
Blake S, Hemming I, Heng JIT, Agostino M. Structure-Based Approaches to Classify the Functional Impact of ZBTB18 Missense Variants in Health and Disease. ACS Chem Neurosci 2021; 12:979-989. [PMID: 33621064 DOI: 10.1021/acschemneuro.0c00758] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
The Cys2His2 type zinc finger is a motif found in many eukaryotic transcription factor proteins that facilitates binding to genomic DNA so as to influence cellular gene expression. One such transcription factor is ZBTB18, characterized as a repressor that orchestrates the development of mammalian tissues including skeletal muscle and brain during embryogenesis. In humans, it has been recognized that disease-associated ZBTB18 missense variants mapping to the coding sequence of the zinc finger domain influence sequence-specific DNA binding, disrupt transcriptional regulation, and impair neural circuit formation in the brain. Furthermore, general population ZBTB18 missense variants that influence DNA binding and transcriptional regulation have also been documented within this domain; however, the molecular traits that explain why some variants cause disease while others do not are poorly understood. Here, we have applied five structure-based approaches to evaluate their ability to discriminate between disease-associated and general population ZBTB18 missense variants. We found that thermodynamic integration and Residue Scanning in the Schrodinger Biologics Suite were the best approaches for distinguishing disease-associated variants from general population variants. Our results demonstrate the effectiveness of structure-based approaches for the functional characterization of missense alleles to DNA binding, zinc finger transcription factor protein-coding genes that underlie human health and disease.
Collapse
Affiliation(s)
- Steven Blake
- Curtin Health Innovation Research Institute, Curtin University, Bentley, Western Australia 6102, Australia
- Ralph and Patricia Sarich Neuroscience Research Institute, Nedlands, Western Australia 6009, Australia
- School of Pharmacy and Biomedical Sciences, Curtin University, Bentley, Western Australia 6845, Australia
| | - Isabel Hemming
- Curtin Health Innovation Research Institute, Curtin University, Bentley, Western Australia 6102, Australia
- Ralph and Patricia Sarich Neuroscience Research Institute, Nedlands, Western Australia 6009, Australia
- The Faculty of Health and Medical Sciences, Medical School, The University of Western Australia, Crawley, Western Australia 6009, Australia
| | - Julian Ik-Tsen Heng
- Curtin Health Innovation Research Institute, Curtin University, Bentley, Western Australia 6102, Australia
- Ralph and Patricia Sarich Neuroscience Research Institute, Nedlands, Western Australia 6009, Australia
| | - Mark Agostino
- Curtin Health Innovation Research Institute, Curtin University, Bentley, Western Australia 6102, Australia
- School of Pharmacy and Biomedical Sciences, Curtin University, Bentley, Western Australia 6845, Australia
- Curtin Institute for Computation, Curtin University, Bentley, Western Australia, Australia
| |
Collapse
|
14
|
MD-TSPC4: Computational Method for Predicting the Thermal Stability of I-Motif. Int J Mol Sci 2020; 22:ijms22010061. [PMID: 33374624 PMCID: PMC7793491 DOI: 10.3390/ijms22010061] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2020] [Revised: 12/16/2020] [Accepted: 12/17/2020] [Indexed: 11/23/2022] Open
Abstract
I-Motif is a tetrameric cytosine-rich DNA structure with hemi-protonated cytosine: cytosine base pairs. Recent evidence showed that i-motif structures in human cells play regulatory roles in the genome. Therefore, characterization of novel i-motifs and investigation of their functional implication are urgently needed for comprehensive understanding of their roles in gene regulation. However, considering the complications of experimental investigation of i-motifs and the large number of putative i-motifs in the genome, development of an in silico tool for the characterization of i-motifs in the high throughput scale is necessary. We developed a novel computation method, MD-TSPC4, to predict the thermal stability of i-motifs based on molecular modeling and molecular dynamic simulation. By assuming that the flexibility of loops in i-motifs correlated with thermal stability within certain temperature ranges, we evaluated the correlation between the root mean square deviations (RMSDs) of model structures and the thermal stability as the experimentally obtained melting temperature (Tm). Based on this correlation, we propose an equation for Tm prediction from RMSD. We expect this method can be useful for estimating the overall structure and stability of putative i-motifs in the genome, which can be a starting point of further structural and functional studies of i-motifs.
Collapse
|
15
|
Wang X, Zhang H, Li W. DNA-binding mechanisms of human and mouse cGAS: a comparative MD and MM/GBSA study. Phys Chem Chem Phys 2020; 22:26390-26401. [PMID: 33179635 DOI: 10.1039/d0cp04162a] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Cyclic GMP-AMP synthase (cGAS) can detect the presence of cytoplasmic DNA and activate the innate immune system via the cGAS-STING pathway. Although several structures of cGAS-DNA complexes were resolved recently, the molecular mechanism of cGAS in its recognition of DNA has not yet been fully understood. In order to reveal the subtle differences between human and mouse cGAS in terms of their DNA-binding mechanisms, four systems, both human and mouse cGAS in complex with two different DNA sequences of equal length, were studied by molecular dynamics simulations and molecular mechanics/generalized Born surface area analysis. Several residues, including ARG176/ARG161, ARG195/ARG180, ASN210/ASN196, LYS384/LYS372, CYM397/CYM385, LYS403/LYS391, LYS407/LYS395, and LYS411/LYS399, were identified to be the common key residues in the recognition of DNA for cGAS in both humans and mice. In addition, four residue pairs LYS173/ARG158, ASP177/LYS162, CYS199/LYS184, and GLU398/SER387 were suggested to be the major residues that make human cGAS and mouse cGAS different in terms of their binding to DNA. Besides the well-known zinc-thumb domain, two residues at the kink of the spine helix were also proposed for the first time to be the major binding motifs in cGAS-DNA interaction.
Collapse
Affiliation(s)
- Xiaowen Wang
- Institute for Advanced Study, Shenzhen University, Room 341, Administration Building, Shenzhen 518060, China.
| | | | | |
Collapse
|
16
|
Forouzesh N, Mukhopadhyay A, Watson LT, Onufriev AV. Multidimensional Global Optimization and Robustness Analysis in the Context of Protein-Ligand Binding. J Chem Theory Comput 2020; 16:4669-4684. [PMID: 32450041 PMCID: PMC8594251 DOI: 10.1021/acs.jctc.0c00142] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Accuracy of protein-ligand binding free energy calculations utilizing implicit solvent models is critically affected by parameters of the underlying dielectric boundary, specifically, the atomic and water probe radii. Here, a global multidimensional optimization pipeline is developed to find optimal atomic radii specifically for protein-ligand binding calculations in implicit solvent. The computational pipeline has these three key components: (1) a massively parallel implementation of a deterministic global optimization algorithm (VTDIRECT95), (2) an accurate yet reasonably fast generalized Born implicit solvent model (GBNSR6), and (3) a novel robustness metric that helps distinguish between nearly degenerate local minima via a postprocessing step of the optimization. A graph-based "kT-connectivity" approach to explore and visualize the multidimensional energy landscape is proposed: local minima that can be reached from the global minimum without exceeding a given energy threshold (kT) are considered to be connected. As an illustration of the capabilities of the optimization pipeline, we apply it to find a global optimum in the space of just five radii: four atomic (O, H, N, and C) radii and water probe radius. The optimized radii, ρW = 1.37 Å, ρC = 1.40 Å, ρH = 1.55 Å, ρN = 2.35 Å, and ρO = 1.28 Å, lead to a closer agreement of electrostatic binding free energies with the explicit solvent reference than two commonly used sets of radii previously optimized for small molecules. At the same time, the ability of the optimizer to find the global optimum reveals fundamental limits of the common two-dielectric implicit solvation model: the computed electrostatic binding free energies are still almost 4 kcal/mol away from the explicit solvent reference. The proposed computational approach opens the possibility to further improve the accuracy of practical computational protocols for binding free energy calculations.
Collapse
Affiliation(s)
- Negin Forouzesh
- Department of Computer Science, Virginia Polytechnic Institute & State University, Blacksburg, Virginia 24061, United States
| | - Abhishek Mukhopadhyay
- Department of Physics, Virginia Polytechnic Institute & State University, Blacksburg, Virginia 24061, United States
| | - Layne T Watson
- Department of Computer Science, Virginia Polytechnic Institute & State University, Blacksburg, Virginia 24061, United States
- Department of Mathematics, Virginia Polytechnic Institute & State University, Blacksburg, Virginia 24061, United States
- Department of Aerospace and Ocean Engineering, Virginia Polytechnic Institute & State University, Blacksburg, Virginia 24061, United States
- Center for Soft Matter and Biological Physics, Virginia Polytechnic Institute & State University, Blacksburg, Virginia 24061, United States
| | - Alexey V Onufriev
- Department of Computer Science, Virginia Polytechnic Institute & State University, Blacksburg, Virginia 24061, United States
- Department of Physics, Virginia Polytechnic Institute & State University, Blacksburg, Virginia 24061, United States
- Center for Soft Matter and Biological Physics, Virginia Polytechnic Institute & State University, Blacksburg, Virginia 24061, United States
| |
Collapse
|
17
|
Röder K, Stirnemann G, Dock-Bregeon AC, Wales DJ, Pasquali S. Structural transitions in the RNA 7SK 5' hairpin and their effect on HEXIM binding. Nucleic Acids Res 2020; 48:373-389. [PMID: 31732748 PMCID: PMC7145557 DOI: 10.1093/nar/gkz1071] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2019] [Revised: 10/23/2019] [Accepted: 10/31/2019] [Indexed: 12/25/2022] Open
Abstract
7SK RNA, as part of the 7SK ribonucleoprotein complex, is crucial to the regulation of transcription by RNA-polymerase II, via its interaction with the positive transcription elongation factor P-TEFb. The interaction is induced by binding of the protein HEXIM to the 5′ hairpin (HP1) of 7SK RNA. Four distinct structural models have been obtained experimentally for HP1. Here, we employ computational methods to investigate the relative stability of these structures, transitions between them, and the effects of mutations on the observed structural ensembles. We further analyse the results with respect to mutational binding assays, and hypothesize a mechanism for HEXIM binding. Our results indicate that the dominant structure in the wild type exhibits a triplet involving the unpaired nucleotide U40 and the base pair A43-U66 in the GAUC/GAUC repeat. This conformation leads to an open major groove with enough potential binding sites for peptide recognition. Sequence mutations of the RNA change the relative stability of the different structural ensembles. Binding affinity is consequently lost if these changes alter the dominant structure.
Collapse
Affiliation(s)
- Konstantin Röder
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK
| | - Guillaume Stirnemann
- CNRS Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique, PSL University, Université de Paris, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Anne-Catherine Dock-Bregeon
- Laboratoire de Biologie Intégrative des Modèles Marins, UMR CNRS 8227, Sorbonne Université, Station Biologique de Roscoff, 29680 Roscoff, France
| | - David J Wales
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK
| | - Samuela Pasquali
- Laboratoire CiTCoM, CNRS UMR 8038, Université de Paris, 4 Avenue de l'observatoire, 75270 Paris, France
| |
Collapse
|
18
|
Zakrevsky P, Kasprzak WK, Heinz WF, Wu W, Khant H, Bindewald E, Dorjsuren N, Fields EA, de Val N, Jaeger L, Shapiro BA. Truncated tetrahedral RNA nanostructures exhibit enhanced features for delivery of RNAi substrates. NANOSCALE 2020; 12:2555-2568. [PMID: 31932830 DOI: 10.1039/c9nr08197f] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Using RNA as a material for nanoparticle construction provides control over particle size and shape at the nano-scale. RNA nano-architectures have shown promise as delivery vehicles for RNA interference (RNAi) substrates, allowing multiple functional entities to be combined on a single particle in a programmable fashion. Rather than employing a completely bottom-up approach to scaffold design, here multiple copies of an existing synthetic supramolecular RNA nano-architecture serve as building blocks along with additional motifs for the design of a novel truncated tetrahedral RNA scaffold, demonstrating that rationally designed RNA assemblies can themselves serve as modular pieces in the construction of larger rationally designed structures. The resulting tetrahedral scaffold displays enhanced characteristics for RNAi-substrate delivery in comparison to similar RNA-based scaffolds, as evidenced by its increased functional capacity, increased cellular uptake and ultimately an increased RNAi efficacy of its adorned Dicer substrate siRNAs. The unique truncated tetrahedral shape of the nanoparticle core appears to contribute to this particle's enhanced function, indicating the physical characteristics of RNA scaffolds merit significant consideration when designing platforms for delivery of functional RNAs via RNA nanoparticles.
Collapse
Affiliation(s)
- Paul Zakrevsky
- RNA Biology Laboratory, National Cancer Institute, Frederick, MD 21702, USA.
| | - Wojciech K Kasprzak
- Basic Science Program, Frederick National Laboratory for Cancer Research, Frederick, MD 21702, USA
| | - William F Heinz
- Optical Microscopy and Analysis Laboratory, Cancer Research Technology Program, Frederick National Laboratory for Cancer Research, Frederick, MD 21702, USA
| | - Weimin Wu
- Center for Molecular Microscopy, Center for Cancer Research, National Cancer Institute, National Institutes of Health, Frederick, MD 21702, USA
| | - Htet Khant
- Center for Molecular Microscopy, Center for Cancer Research, National Cancer Institute, National Institutes of Health, Frederick, MD 21702, USA
| | - Eckart Bindewald
- Basic Science Program, Frederick National Laboratory for Cancer Research, Frederick, MD 21702, USA
| | - Nomongo Dorjsuren
- RNA Biology Laboratory, National Cancer Institute, Frederick, MD 21702, USA.
| | - Eric A Fields
- RNA Biology Laboratory, National Cancer Institute, Frederick, MD 21702, USA.
| | - Natalia de Val
- Center for Molecular Microscopy, Center for Cancer Research, National Cancer Institute, National Institutes of Health, Frederick, MD 21702, USA and Cancer Research Technology Program, Frederick National Laboratory for Cancer Research, Leidos Biomedical Research Inc, Frederick, MD 21701, USA
| | - Luc Jaeger
- Department of Chemistry and Biochemistry, Biomolecular Science and Engineering Program, University of California, Santa Barbara, CA 93106-9510, USA.
| | - Bruce A Shapiro
- RNA Biology Laboratory, National Cancer Institute, Frederick, MD 21702, USA.
| |
Collapse
|
19
|
Wei H, Luo A, Qiu T, Luo R, Qi R. Improved Poisson-Boltzmann Methods for High-Performance Computing. J Chem Theory Comput 2019; 15:6190-6202. [PMID: 31525962 DOI: 10.1021/acs.jctc.9b00602] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
Implicit solvent models based on the Poisson-Boltzmann equation (PBE) have been widely used to study electrostatic interactions in biophysical processes. These models often treat the solvent and solute regions as high and low dielectric continua, leading to a large jump in dielectrics across the molecular surface which is difficult to handle. Higher order interface schemes are often needed to seek higher accuracy for PBE applications. However, these methods are usually very liberal in the use of grid points nearby the molecular surface, making them difficult to use on high-performance computing platforms. Alternatively, the harmonic average (HA) method has been used to approximate dielectric interface conditions near the molecular surface with surprisingly good convergence and is well suited for high-performance computing. By adopting a 7-point stencil, the HA method is advantageous in generating simple 7-banded coefficient matrices, which greatly facilitate linear system solution with dense data parallelism, on high-performance computing platforms such as a graphics processing unit (GPU). However, the HA method is limited due to its lower accuracy. Therefore, it would be of great interest for high-performance applications to develop more accurate methods while retaining the simplicity and effectiveness of the 7-point stencil discretization scheme. In this study, we have developed two new algorithms based on the spirit of the HA method by introducing more physical interface relations and imposing the discretized Poisson's equation to the second order, respectively. Our testing shows that, for typical biomolecules, the new methods significantly improve the numerical accuracy to that comparable to the second-order solvers and with ∼65% overall efficiency gain on widely available high-performance GPU platforms.
Collapse
|
20
|
Sun Z, Wang X. Thermodynamics of Helix formation in small peptides of varying lengthin vacuo, implicit solvent and explicit solvent: Comparison between AMBER force fields. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2019. [DOI: 10.1142/s0219633619500159] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
Helix formation is of great significance in protein folding. The helix-forming tendencies of amino acids are accumulated along the sequence to determine the helix-forming tendency of peptides. Computer simulation can be used to model this process in atomic details and give structural insights. In the current work, we employ equilibrate-state free energy simulation to systematically study the folding/unfolding thermodynamics of a series of mutated peptides. Two AMBER force fields including AMBER99SB and AMBER14SB are compared. The new 14SB force field uses refitted torsion parameters compared with 99SB and they share the same atomic charge scheme. We find that in vacuo the helix formation is mutation dependent, which reflects the different helix propensities of different amino acids. In general, there are helix formers, helix indifferent groups and helix breakers. The helical structure becomes more favored when the length of the sequence becomes longer, which arises from the formation of additional backbone hydrogen bonds in the lengthened sequence. Therefore, the helix indifferent groups and helix breakers will become helix formers in long sequences. Also, protonation-dependent helix formation is observed for ionizable groups. In 14SB, the helical structures are more stable than in 99SB and differences can be observed in their grouping schemes, especially in the helix indifferent group. In solvents, all mutations are helix indifferent due to protein–solvent interactions. The decrease in the number of backbone hydrogen bonds is the same with the increase in the number of protein–water hydrogen bonds. The 14SB in explicit solvent is able to capture the free energy minima in the helical state while 14SB in implicit solvent, 99SB in explicit solvent and 99SB in implicit solvent cannot. The helix propensities calculated under 14SB agree with the corresponding experimental values, while the 99SB results obviously deviate from the references. Hence, implicit solvent models are unable to correctly describe the thermodynamics even for the simple helix formation in isolated peptides. Well-developed force fields and explicit solvents are needed to correctly describe the protein dynamics. Aside from the free energy, differences in conformational ensemble under different force fields in different solvent models are observed. The numbers of hydrogen bonds formed under different force fields agree and they are mostly determined by the solvent model.
Collapse
Affiliation(s)
- Zhaoxi Sun
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, P. R. China
- Computational Biomedicine (IAS-5/INM-9), Forschungszentrum Jülich, Jülich 52425, Germany
| | - Xiaohui Wang
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, P. R. China
- Institute of Computational Science, Universitàdella Svizzeraitaliana (USI), Via Giuseppe Buffi 13, CH-6900 Lugano, Ticino, Switzerland
| |
Collapse
|
21
|
Abstract
It would often be useful in computer simulations to use an implicit description of solvation effects, instead of explicitly representing the individual solvent molecules. Continuum dielectric models often work well in describing the thermodynamic aspects of aqueous solvation and can be very efficient compared to the explicit treatment of the solvent. Here, we review a particular class of so-called fast implicit solvent models, generalized Born (GB) models, which are widely used for molecular dynamics (MD) simulations of proteins and nucleic acids. These approaches model hydration effects and provide solvent-dependent forces with efficiencies comparable to molecular-mechanics calculations on the solute alone; as such, they can be incorporated into MD or other conformational searching strategies in a straightforward manner. The foundations of the GB model are reviewed, followed by examples of newer, emerging models and examples of important applications. We discuss their strengths and weaknesses, both for fidelity to the underlying continuum model and for the ability to replace explicit consideration of solvent molecules in macromolecular simulations.
Collapse
Affiliation(s)
- Alexey V Onufriev
- Departments of Computer Science and Physics, Center for Soft Matter and Biological Physics, Virginia Tech, Blacksburg, Virginia 24060, USA;
| | - David A Case
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA;
| |
Collapse
|
22
|
Luo A, Li X, Zhang X, Zhan H, Du H, Zhang Y, Peng X. Identification of AtHsp90.6 involved in early embryogenesis and its structure prediction by molecular dynamics simulations. ROYAL SOCIETY OPEN SCIENCE 2019; 6:190219. [PMID: 31218061 PMCID: PMC6550000 DOI: 10.1098/rsos.190219] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/17/2019] [Accepted: 04/02/2019] [Indexed: 05/29/2023]
Abstract
Heat-shock protein of 90 kDa (Hsp90) is a key molecular chaperone involved in folding the synthesized protein and controlling protein quality. Conformational dynamics coupled to ATPase activity in N-terminal domain is essential for Hsp90's function. However, the relevant process is still largely unknown in plant Hsp90s, especially those required for plant embryogenesis which is inextricably tied up with human survival. Here, AtHsp90.6, a member of Hsp90 family in Arabidopsis, was firstly identified as a protein essential for embryogenesis. Thus we modelled AtHsp90.6 in its functionally closed 'lid-down' and open 'lid-up' states, exploring the nucleotide binding mechanism in these two states. Free energy landscape and electrostatic potential analysis revealed the switching mechanism between these two states. Collectively, this study quantitatively analysed the conformational changes of AtHsp90.6 bound to ATP or ADP. This result may help us understand the mechanism of action of AtHsp90.6 in future.
Collapse
Affiliation(s)
- An Luo
- College of Life Science, Yangtze University, Jingzhou 434023, People's Republic of China
| | - Xinbo Li
- College of Life Science, State Key Laboratory of Hybrid Rice, Wuhan University, Wuhan 430072, People's Republic of China
- Center for Tissue Engineering and Regenerative Medicine, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan 430072, People's Republic of China
| | - Xuecheng Zhang
- College of Life Science, State Key Laboratory of Hybrid Rice, Wuhan University, Wuhan 430072, People's Republic of China
| | - Huadong Zhan
- College of Life and Environment Sciences, Shanghai Normal University, Shanghai, 200234, People's Republic of China
| | - Hewei Du
- College of Life Science, Yangtze University, Jingzhou 434023, People's Republic of China
| | - Yubo Zhang
- Department of Food Science, Foshan University, Foshan 528231, People's Republic of China
| | - Xiongbo Peng
- College of Life Science, State Key Laboratory of Hybrid Rice, Wuhan University, Wuhan 430072, People's Republic of China
| |
Collapse
|
23
|
Lu Y, Yuan L, Chen X, Zhang A, Zhang P, Zou D. Systematic analysis and identification of unexpected interactions from the neuroprotein drug interactome in hydrocephalus pharmacological intervention. J Bioinform Comput Biol 2019; 17:1950002. [PMID: 30866733 DOI: 10.1142/s0219720019500021] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Hydrocephalus is a neurological condition caused by an abnormal accumulation of cerebrospinal fluid; pharmacological intervention of the disease has been found to elicit a variety of adverse drug reactions (ADRs) in central nervous system (CNS) by unexpectedly targeting certain functional neuroproteins. Here, a systematic neuroprotein drug interactome (SNDI) is created for 11 hydrocephalus drugs/metabolites plus 20 control drugs across 518 druggable pockets on the surface of 472 CNS neuroproteins via a large-scale molecular docking approach. Heuristic clustering analysis of the SNDI profile divides the 31 investigated drug ligands into a distinct panel and a background panel; the former consists of two hydrocephalus drugs (Furosemide and Triamterene) and their respective metabolites (Furosemide glucuronide and Hydroxytriamterene) that are inferred to have generally high affinity towards the whole array of neuroprotein pockets. A total of 13 neuroproteins are enriched in gene ontology semantic mining as putative unexpected targets of the distinct panel, and their intermolecular interactions with hydrocephalus drugs/metabolites are investigated in detail using dynamics simulation and energetics analysis. We also perform kinase assay and viability test to substantiate the interactome analysis. It is found that the Furosemide and Triamterene have significant cytotoxic effects on normal human astrocytes, in which the Triamterene can inhibit the neurokinase ROCK2, a representative of putative unexpected targets, with a high activity, which is comparable with the sophisticated ROCK2 inhibitor Fasudil.
Collapse
Affiliation(s)
- Youming Lu
- 1 Department of Neurosurgery, Shanghai Jiao Tong University Affiliated Sixth People's Hospital, Shanghai 200233, P. R. China
| | - Lei Yuan
- 1 Department of Neurosurgery, Shanghai Jiao Tong University Affiliated Sixth People's Hospital, Shanghai 200233, P. R. China
| | - Xin Chen
- 1 Department of Neurosurgery, Shanghai Jiao Tong University Affiliated Sixth People's Hospital, Shanghai 200233, P. R. China
| | - Aijun Zhang
- 1 Department of Neurosurgery, Shanghai Jiao Tong University Affiliated Sixth People's Hospital, Shanghai 200233, P. R. China
| | - Pengqi Zhang
- 1 Department of Neurosurgery, Shanghai Jiao Tong University Affiliated Sixth People's Hospital, Shanghai 200233, P. R. China
| | - Dongdong Zou
- 1 Department of Neurosurgery, Shanghai Jiao Tong University Affiliated Sixth People's Hospital, Shanghai 200233, P. R. China
| |
Collapse
|
24
|
Yu Z, Liu Y, Zhu J, Han J, Tian X, Han W, Zhao L. Insights from molecular dynamics simulations and steered molecular dynamics simulations to exploit new trends of the interaction between HIF-1α and p300. J Biomol Struct Dyn 2019; 38:1-12. [DOI: 10.1080/07391102.2019.1580616] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- Zhengfei Yu
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, Changchun, China
| | - Ye Liu
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, Changchun, China
| | - Jingxuan Zhu
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, Changchun, China
| | - Jiarui Han
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, Changchun, China
| | - Xiaopian Tian
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, Changchun, China
| | - Weiwei Han
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, Changchun, China
| | - Li Zhao
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, Changchun, China
| |
Collapse
|
25
|
Wei H, Luo R, Qi R. An efficient second-order poisson-boltzmann method. J Comput Chem 2019; 40:1257-1269. [PMID: 30776135 PMCID: PMC6422926 DOI: 10.1002/jcc.25783] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2018] [Revised: 11/29/2018] [Accepted: 01/06/2019] [Indexed: 11/08/2022]
Abstract
Immersed interface method (IIM) is a promising high-accuracy numerical scheme for the Poisson-Boltzmann model that has been widely used to study electrostatic interactions in biomolecules. However, the IIM suffers from instability and slow convergence for typical applications. In this study, we introduced both analytical interface and surface regulation into IIM to address these issues. The analytical interface setup leads to better accuracy and its convergence closely follows a quadratic manner as predicted by theory. The surface regulation further speeds up the convergence for nontrivial biomolecules. In addition, uncertainties of the numerical energies for tested systems are also reduced by about half. More interestingly, the analytical setup significantly improves the linear solver efficiency and stability by generating more precise and better-conditioned linear systems. Finally, we implemented the bottleneck linear system solver on GPUs to further improve the efficiency of the method, so it can be widely used for practical biomolecular applications. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Haixin Wei
- Department of Chemical Engineering and Materials Science, University of California, Irvine, California, 92697
| | - Ray Luo
- Department of Chemical Engineering and Materials Science, University of California, Irvine, California, 92697.,Department of Molecular Biology and Biochemistry, University of California, Irvine, California, 92697.,Department of Biomedical Engineering, University of California, Irvine, California, 92697
| | - Ruxi Qi
- Department of Molecular Biology and Biochemistry, University of California, Irvine, California, 92697
| |
Collapse
|
26
|
Nishikawa N, Sakae Y, Gouda T, Tsujimura Y, Okamoto Y. Structural Analysis of a Trimer of β 2-Microgloblin Fragment by Molecular Dynamics Simulations. Biophys J 2019; 116:781-790. [PMID: 30771855 DOI: 10.1016/j.bpj.2018.11.3143] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2017] [Revised: 06/08/2018] [Accepted: 11/06/2018] [Indexed: 01/22/2023] Open
Abstract
A peptide β2-m21-31, which is a fragment from residue 21 to residue 31 of β2-microgloblin, is experimentally known to self-assemble and form amyloid fibrils. In order to understand the mechanism of amyloid fibril formations, we applied the replica-exchange molecular dynamics method to the system consisting of three fragments of β2-m21-31. From the analyses on the temperature dependence, we found that there is a clear phase transition temperature in which the peptides aggregate with each other. Moreover, we found by the free energy analyses that there are two major stable states: One of them is like amyloid fibrils and the other is amorphous aggregates.
Collapse
Affiliation(s)
- Naohiro Nishikawa
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Aichi, Japan; Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, Okazaki, Aichi, Japan
| | - Yoshitake Sakae
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Aichi, Japan
| | - Takuya Gouda
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Aichi, Japan
| | - Yuichiro Tsujimura
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Aichi, Japan
| | - Yuko Okamoto
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Aichi, Japan; Structural Biology Research Center, Graduate School of Science, Nagoya University, Nagoya, Aichi, Japan; Center for Computational Science, Graduate School of Engineering, Nagoya University, Nagoya, Aichi, Japan; Information Technology Center, Nagoya University, Nagoya, Aichi, Japan; JST-CREST, Nagoya, Aichi, Japan.
| |
Collapse
|
27
|
Combining Structural Data with Computational Methodologies to Investigate Structure-Function Relationships in TRP Channels. Methods Mol Biol 2019; 1987:65-82. [PMID: 31028674 DOI: 10.1007/978-1-4939-9446-5_5] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/08/2022]
Abstract
Since the emergence of high-resolution three-dimensional structures of membrane proteins, and the increasing availability of state-of-the-art algorithms and high-performance-computing facilities, classical molecular dynamics (MD) simulations have become a routine device to explore the molecular behavior of these proteins. The rise of cryo-electron microscopy (cryo-EM) as a credible experimental tool to resolve structures at an atomic level has revolutionized structural biology in recent years, culminating in the disclosure of the first high-resolution three-dimensional structure of a transient receptor potential (TRP) channel, the vanilloid receptor 1 (TRPV1). As a result, the number of research articles investigating the molecular behavior of TRP channels using macromolecular simulation techniques has proliferated. This review provides an overview of the current state of this field, including our understanding of TRP channel structure, the framework of classical MD simulations, and how to perform such simulations to investigate structure-function relationships in TRP channels.
Collapse
|
28
|
Glycan binding and specificity of viral influenza neuraminidases by classical molecular dynamics and replica exchange molecular dynamics simulations. J Biomol Struct Dyn 2018; 37:3354-3365. [DOI: 10.1080/07391102.2018.1514326] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
29
|
Şimşek Y, Brown A. Two-Photon Absorption Properties of Gold Fluorescent Protein: A Combined Molecular Dynamics and Quantum Chemistry Study. J Phys Chem B 2018; 122:5738-5748. [PMID: 29741903 DOI: 10.1021/acs.jpcb.8b00885] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
Molecular dynamics (MD) simulations were carried out to obtain the conformational changes of the chromophore in the gold fluorescent protein (PDB ID: 1OXF ). To obtain two-photon absorption (TPA) cross-sections, time dependent density functional theory (TD-DFT) computations were performed for chromophore geometries sampled along the trajectory. The TD-DFT computations used the CAM-B3LYP functional and 6-31+G(d) basis set. Results showed that two dihedral angles change remarkably over the simulation time. TPA cross-sections were found to average 13.82 GM for the excitation to S1 computed from the equilibrium geometries; however, extending the structures with a water molecule and GLU residue, which make H bonds with the chromophore molecule, increased excitation energies and TPA cross-sections significantly. Besides the effects of the surrounding residues and the dihedrals on the spectroscopic properties, some bond lengths affected the excitation energies and the TPA cross-sections significantly (up to ±25-30%), while the effects of the bond angles were smaller (±5%). Overall the present results provide insight into the effects of the conformational flexibility on TPA (with gold fluorescent protein as a specific example) and suggest that further experimental measurements of TPA for the gold fluorescent protein should be undertaken.
Collapse
Affiliation(s)
- Yusuf Şimşek
- Vocational School of Health Services , Gazi University , 06560 Ankara , Turkey.,Department of Chemistry , University of Alberta , Edmonton , AB T6G 2G2 , Canada
| | - Alex Brown
- Department of Chemistry , University of Alberta , Edmonton , AB T6G 2G2 , Canada
| |
Collapse
|
30
|
Meng L, Huang Z. In silico-in vitro discovery of untargeted kinase-inhibitor interactions from kinase-targeted therapies: A case study on the cancer MAPK signaling pathway. Comput Biol Chem 2018; 75:196-204. [PMID: 29803964 DOI: 10.1016/j.compbiolchem.2018.05.012] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2018] [Revised: 05/09/2018] [Accepted: 05/13/2018] [Indexed: 12/14/2022]
Abstract
Protein kinase inhibitors have been widely used as therapeutic agents to treat a variety of diseases, but many of them may cause off-target effects by unexpectedly targeting other noncognate kinases due to high conversion across the protein kinase family. The mitogen-activated protein kinase (MAPK) signaling pathway plays an essential role in tumorigenesis, which has been recognized as a high priority in the druggable target candidates of anticancer therapy. Here, we attempt to investigate the untargeted kinase-inhibitor interactions (UKIIs) of kinase-targeted therapies for the cancer MAPK signaling cascade via an integration of biomolecular modeling, cell viability assay and kinase inhibition analysis. A systematic kinase-inhibitor interaction profile is created for 28 FDA-approved kinase inhibitor drugs across 9 caner-related MAPK kinases. The created profile is analyzed at structural, energetic and dynamic levels and, consequently, totally 18 promising UKII pairs with high theoretical affinity are derived, from which the noncognate inhibitors Cabozantinib, Regorafenib and Crizotinib are selected to test their cytotoxic effects on human epithelial colorectal adenocarcinoma Caco-2 cell line and inhibition activity against the recombinant protein of human p38α kinase domain. The obtained results are compared with two cognate MAPK inhibitors JNK-IN-8 and BIRB796. As might be expected, the Regorafenib, Crizotinib and Cabozantinib exhibit high, moderate and low cytotoxicities, respectively. In addition, the Regorafenib is determined to have a potent p38α-inhibitory activity. This is basically in line with the test results of positive controls JNK-IN-8 and BIRB796 and can be well confirmed by computational modeling.
Collapse
Affiliation(s)
- Li Meng
- College of Pharmaceutical Sciences, Jiangsu Vocational College of Medicine, Yancheng 224008, China
| | - Zhijun Huang
- Department of General Surgery, Yancheng First People's Hospital, Yancheng 224005, China.
| |
Collapse
|
31
|
Kumbhar BV, Panda D, Kunwar A. Interaction of microtubule depolymerizing agent indanocine with different human αβ tubulin isotypes. PLoS One 2018; 13:e0194934. [PMID: 29584771 PMCID: PMC5870988 DOI: 10.1371/journal.pone.0194934] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2017] [Accepted: 03/13/2018] [Indexed: 12/02/2022] Open
Abstract
Tubulin isotypes are known to regulate the stability and dynamics of microtubules, and are also involved in the development of resistance against microtubule-targeted cancer drugs. Indanocine, a potent microtubule depolymerizing agent, is highly active against multidrug-resistant (MDR) cancer cells without affecting normal cells. It is known to disrupt microtubule dynamics in cells and induce apoptotic cell death. Indanocine is reported to bind to tubulin at the colchicine site i.e. at the interface of αβ tubulin heterodimer. However, it’s precise binding mode, involved molecular interactions and the binding affinities with different αβ-tubulin isotypes present in MDR cells are not well understood. Here, the binding affinities of human αβ-tubulin isotypes with indanocine were examined, employing the molecular modeling approach i.e. docking, molecular dynamics simulation and binding energy calculations. Multiple sequence analysis suggests that the amino acid sequences are different in the indanocine binding pockets of βI, βIIa, βIII and βVI isotypes. However, such differences are not observed in the amino acid sequences of βIVa, βIVb, and βV tubulin isotypes at indanocine binding pockets. Docking and molecular dynamics simulation results show that indanocine prefers the interface binding pocket of αβIIa, αβIII, αβIVb, αβV, and αβVI tubulin isotypes; whereas it is expelled from the interface binding pocket of αβIVa and αβI-tubulin isotypes. Further, binding free energy calculations show that αβVI has the highest binding affinity and αβI has the lowest binding affinity for indanocine among all β-tubulin isotypes. The binding free energy decreases in the order of αβVI > αβIVb > αβIIa > αβIII > αβV > αβIVa > αβI. Thus, our study provides a significant understanding of involved molecular interactions of indanocine with tubulin isotypes, which may help to design potent indanocine analogues for specific tubulin isotypes in MDR cells in future.
Collapse
Affiliation(s)
- Bajarang Vasant Kumbhar
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai, Maharashtra, India
| | - Dulal Panda
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai, Maharashtra, India
| | - Ambarish Kunwar
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai, Maharashtra, India
- * E-mail:
| |
Collapse
|
32
|
Devi M, Chingbiaknem E, Lyngdoh RHD. A molecular mechanics study on GA codon box translation. J Theor Biol 2018; 441:28-43. [PMID: 29305181 DOI: 10.1016/j.jtbi.2018.01.002] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2017] [Revised: 01/01/2018] [Accepted: 01/02/2018] [Indexed: 11/28/2022]
Abstract
The GA codon box incorporates the two-fold degeneracy of aspartic acid and of glutamic acid. Using the molecular mechanics approach of the AMBER suite, the four codons of the GA box are paired via H-bonding with two aspartic acid anticodons and two glutamic acid anticodons to yield 8 cognate and 11 non-cognate codon-anticodon duplexes. In addition four select non-cognate duplexes between the GA box codons and three alanine anticodons are also studied. These 23 duplexes display a variety of base-pairing possibilities at the wobble position. Cognate duplexes are differentiated from non-cognate duplexes on the grounds of structure and stability (chiefly the former). The results are in line with Crick's wobble hypothesis, and corroborate the observed reading properties of the aspartic acid anticodons GUC and QUC and of the glutamic acid anticodons CUC and SmnUC.
Collapse
Affiliation(s)
- Martina Devi
- Department of Chemistry, North-Eastern Hill University, Shillong 793022, India
| | - Esther Chingbiaknem
- Department of Chemistry, North-Eastern Hill University, Shillong 793022, India
| | - R H Duncan Lyngdoh
- Department of Chemistry, North-Eastern Hill University, Shillong 793022, India.
| |
Collapse
|
33
|
Zhuang RJ, Jin WD, Wang XY, Wu XM. Identification and characterization of the druggable kinase targets of olmesartan and its analogues from a systematic kinase-chemical interaction profile in atherosclerosis. J Mol Graph Model 2018; 80:211-216. [PMID: 29414040 DOI: 10.1016/j.jmgm.2018.01.016] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2017] [Revised: 01/12/2018] [Accepted: 01/22/2018] [Indexed: 10/17/2022]
Abstract
Olmesartan (OL) is the pharmacologically active metabolite of Olmesartan medoxomil (OM), an FDA-approved angiotensin II receptor antagonist for administrating cardiovascular diseases. The drug has been found to have potential effects on diverse protein kinase signaling involved in the pathogenesis of atherosclerosis, either by directly inhibiting the hub kinases or by indirectly modulating marginal members in the signaling pathways. In the present study, we computationally model the kinase-chemical Interaction Profile between six OL-related chemicals (i.e. OL, OM, Valsartan [VL], Losartan [LS], Candesartan [CD] and Telmisartan [TL]) and 23 human protein kinases in atherosclerosis. The profile is analyzed systematically at molecular level to identify unexpected kinase targets for OL. There is a good consistence between co-citation frequency and affinity scoring for the chemical association with kinase candidates; the OL and its analogs VL and LS exhibit a similar binding profile to the atherosclerosis kinase spectrum. It is suggested that the Ser/Thr-specific kinases PI3Kα and ROCK1 are potential druggable targets of OL for atherosclerosis therapy. As a paradigm, kinase assays reveal that the inhibitory potency of OL and Y-27632 (positive control) on ROCK1 is determined at micromolar level, while the OM (negative control) possesses no detectable activity for the kinase.
Collapse
Affiliation(s)
- Rui-Juan Zhuang
- Department of Cardioloy, Affiliated Third Hospital of Nantong University, Wuxi 214041, China
| | - Wei-Dong Jin
- Department of Cardioloy, Affiliated Third Hospital of Nantong University, Wuxi 214041, China
| | - Xiao-Yan Wang
- Department of Cardioloy, Affiliated Third Hospital of Nantong University, Wuxi 214041, China
| | - Xue-Ming Wu
- Department of Cardioloy, Affiliated Third Hospital of Nantong University, Wuxi 214041, China.
| |
Collapse
|
34
|
Yang HC, Ge YC, Yang CH, Chao WC. Substrate Channeling of Prostaglandin H2 on the Stereochemical Control of a Cascade Cyclization Route. ACS Catal 2018. [DOI: 10.1021/acscatal.7b03687] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Hsiao-Ching Yang
- Department of Chemistry, Fu Jen Catholic University, New Taipei City 24205, Taiwan
| | - Yung-Chi Ge
- Department of Chemistry, Fu Jen Catholic University, New Taipei City 24205, Taiwan
| | - Cheng-Han Yang
- Department of Chemistry, Fu Jen Catholic University, New Taipei City 24205, Taiwan
| | - Wei-Chih Chao
- Department of Chemistry, Fu Jen Catholic University, New Taipei City 24205, Taiwan
| |
Collapse
|
35
|
Izadi S, Harris RC, Fenley MO, Onufriev AV. Accuracy Comparison of Generalized Born Models in the Calculation of Electrostatic Binding Free Energies. J Chem Theory Comput 2018; 14:1656-1670. [PMID: 29378399 DOI: 10.1021/acs.jctc.7b00886] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
The need for accurate yet efficient representation of the aqueous environment in biomolecular modeling has led to the development of a variety of generalized Born (GB) implicit solvent models. While many studies have focused on the accuracy of available GB models in predicting solvation free energies, a systematic assessment of the quality of these models in binding free energy calculations, crucial for rational drug design, has not been undertaken. Here, we evaluate the accuracies of eight common GB flavors (GB-HCT, GB-OBC, GB-neck2, GBNSR6, GBSW, GBMV1, GBMV2, and GBMV3), available in major molecular dynamics packages, in predicting the electrostatic binding free energies ( ΔΔ Gel) for a diverse set of 60 biomolecular complexes belonging to four main classes: protein-protein, protein-drug, RNA-peptide, and small complexes. The GB flavors are examined in terms of their ability to reproduce the results from the Poisson-Boltzmann (PB) model, commonly used as accuracy reference in this context. We show that the agreement with the PB of ΔΔ Gel estimates varies widely between different GB models and also across different types of biomolecular complexes, with R2 correlations ranging from 0.3772 to 0.9986. A surface-based "R6" GB model recently implemented in AMBER shows the closest overall agreement with reference PB ( R2 = 0.9949, RMSD = 8.75 kcal/mol). The RNA-peptide and protein-drug complex sets appear to be most challenging for all but one model, as indicated by the large deviations from the PB in ΔΔ Gel. Small neutral complexes present the least challenge for most of the GB models tested. The quantitative demonstration of the strengths and weaknesses of the GB models across the diverse complex types provided here can be used as a guide for practical computations and future development efforts.
Collapse
Affiliation(s)
- Saeed Izadi
- Early Stage Pharmaceutical Development , Genentech Inc. , 1 DNA Way , South San Francisco , California 94080 , United States
| | - Robert C Harris
- Department of Pharmaceutical Sciences , University of Maryland School of Pharmacy , Baltimore , Maryland 21201 , United States
| | - Marcia O Fenley
- Institute of Molecular Biophysics , Florida State University , Tallahassee , Florida 32306-3408 , United States
| | | |
Collapse
|
36
|
Zhou HX, Pang X. Electrostatic Interactions in Protein Structure, Folding, Binding, and Condensation. Chem Rev 2018; 118:1691-1741. [PMID: 29319301 DOI: 10.1021/acs.chemrev.7b00305] [Citation(s) in RCA: 454] [Impact Index Per Article: 75.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Charged and polar groups, through forming ion pairs, hydrogen bonds, and other less specific electrostatic interactions, impart important properties to proteins. Modulation of the charges on the amino acids, e.g., by pH and by phosphorylation and dephosphorylation, have significant effects such as protein denaturation and switch-like response of signal transduction networks. This review aims to present a unifying theme among the various effects of protein charges and polar groups. Simple models will be used to illustrate basic ideas about electrostatic interactions in proteins, and these ideas in turn will be used to elucidate the roles of electrostatic interactions in protein structure, folding, binding, condensation, and related biological functions. In particular, we will examine how charged side chains are spatially distributed in various types of proteins and how electrostatic interactions affect thermodynamic and kinetic properties of proteins. Our hope is to capture both important historical developments and recent experimental and theoretical advances in quantifying electrostatic contributions of proteins.
Collapse
Affiliation(s)
- Huan-Xiang Zhou
- Department of Chemistry and Department of Physics, University of Illinois at Chicago , Chicago, Illinois 60607, United States.,Department of Physics and Institute of Molecular Biophysics, Florida State University , Tallahassee, Florida 32306, United States
| | - Xiaodong Pang
- Department of Physics and Institute of Molecular Biophysics, Florida State University , Tallahassee, Florida 32306, United States
| |
Collapse
|
37
|
Zhong XF, Zhang YB, Huang GD, Ouyang YZ, Liao DJ, Peng JW, Huang WZ. Proteomic analysis of stachyose contribution to the growth of Lactobacillus acidophilus CICC22162. Food Funct 2018; 9:2979-2988. [DOI: 10.1039/c8fo00528a] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Stachyose is a functional oligosaccharide, acting as a potential prebiotic for colonic fermentation.
Collapse
Affiliation(s)
- Xian-feng Zhong
- Department of Food Science
- Foshan University
- Foshan 528231
- China
- Foshan Engineering Research Center for Brewing Technology
| | - Yu-bo Zhang
- Department of Food Science
- Foshan University
- Foshan 528231
- China
- Foshan Engineering Research Center for Brewing Technology
| | - Gui-dong Huang
- Department of Food Science
- Foshan University
- Foshan 528231
- China
- Foshan Engineering Research Center for Brewing Technology
| | - Yong-zhong Ouyang
- School of Environmental and Chemical Engineering
- Foshan University
- Foshan 528231
- China
| | | | - Jia-wei Peng
- Department of Food Science
- Foshan University
- Foshan 528231
- China
- Foshan Engineering Research Center for Brewing Technology
| | - Wei-zhi Huang
- Department of Food Science
- Foshan University
- Foshan 528231
- China
- Foshan Engineering Research Center for Brewing Technology
| |
Collapse
|
38
|
Xiao Y, McElheny D, Hoshi M, Ishii Y. Solid-State NMR Studies of Amyloid Materials: A Protocol to Define an Atomic Model of Aβ(1-42) in Amyloid Fibrils. Methods Mol Biol 2018; 1777:407-428. [PMID: 29744851 DOI: 10.1007/978-1-4939-7811-3_26] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Intense efforts have been made to understand the molecular structures of misfolded amyloid β (Aβ) in order to gain insight into the pathological mechanism of Alzheimer's disease. Solid-state NMR spectroscopy (SSNMR) is considered a primary tool for elucidating the structures of insoluble and noncrystalline amyloid fibrils and other amyloid assemblies. In this chapter, we describe a detailed protocol to obtain the first atomic model of the 42-residue human Aβ peptide Aβ(1-42) in structurally homogeneous amyloid fibrils from our recent SSNMR study (Nat Struct Mol Biol 22:499-505, 2015). Despite great biological and clinical interest in Aβ(1-42) fibrils, their structural details have been long-elusive until this study. The protocol is divided into four sections. First, the solid-phase peptide synthesis (SPPS) and purification of monomeric Aβ(1-42) is described. We illustrate a controlled incubation method to prompt misfolding of Aβ(1-42) into homogeneous amyloid fibrils in an aqueous solution with fragmented Aβ(1-42) fibrils as seeds. Next, we detail analysis of Aβ(1-42) fibrils by SSNMR to obtain structural restraints. Finally, we describe methods to construct atomic models of Aβ(1-42) fibrils based on SSNMR results through two-stage molecular dynamics calculations.
Collapse
Affiliation(s)
- Yiling Xiao
- Department of Chemistry, University of Illinois at Chicago, Chicago, IL, USA
| | - Dan McElheny
- Department of Chemistry, University of Illinois at Chicago, Chicago, IL, USA
| | - Minako Hoshi
- Institute of Biomedical Research and Innovation, FBRI, Kobe, Japan
- Department of Anatomy and Developmental Biology, Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Yoshitaka Ishii
- Department of Chemistry, University of Illinois at Chicago, Chicago, IL, USA.
- School of Life Science and Technology, Tokyo Institute of Technology, Yokohama, Japan.
| |
Collapse
|
39
|
Cacciotto P, Ramaswamy VK, Malloci G, Ruggerone P, Vargiu AV. Molecular Modeling of Multidrug Properties of Resistance Nodulation Division (RND) Transporters. Methods Mol Biol 2018; 1700:179-219. [PMID: 29177832 DOI: 10.1007/978-1-4939-7454-2_11] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Efflux pumps of the resistance nodulation division (RND) superfamily are among the major contributors to intrinsic and acquired multidrug resistance in Gram-negative bacteria. Structural information on AcrAB-TolC and MexAB-OprM, major efflux pumps of Escherichia coli and Pseudomonas aeruginosa respectively, boosted intensive research aimed at understanding the molecular mechanisms ruling the active extrusion processes. In particular, several studies were devoted to the understanding of the determinants behind the extraordinary broad specificity of the RND transporters AcrB and MexB. In this chapter, we discuss the ever-growing role computational methods have been playing in deciphering key structural and dynamical features of these transporters and of their interaction with substrates and inhibitors. We further discuss and illustrate examples from our lab of how molecular docking, homology modeling, all-atom molecular dynamics simulations and in silico free energy estimations can all together give precious insights into the processes of recognition and extrusion of substrates, as well as on the possible inhibition strategies.
Collapse
Affiliation(s)
- Pierpaolo Cacciotto
- Department of Physics, University of Cagliari, s.p. 8, 09042, Monserrato, CA, Italy
| | - Venkata K Ramaswamy
- Department of Physics, University of Cagliari, s.p. 8, 09042, Monserrato, CA, Italy
| | - Giuliano Malloci
- Department of Physics, University of Cagliari, s.p. 8, 09042, Monserrato, CA, Italy
| | - Paolo Ruggerone
- Department of Physics, University of Cagliari, s.p. 8, 09042, Monserrato, CA, Italy
| | - Attilio V Vargiu
- Department of Physics, University of Cagliari, s.p. 8, 09042, Monserrato, CA, Italy.
| |
Collapse
|
40
|
Zhao K, Zhou X, Ding M. Molecular insight into mutation-induced conformational change in metastasic bowel cancer BRAF kinase domain and its implications for selective inhibitor design. J Mol Graph Model 2017; 79:59-64. [PMID: 29145034 DOI: 10.1016/j.jmgm.2017.11.005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2017] [Revised: 10/31/2017] [Accepted: 11/06/2017] [Indexed: 11/30/2022]
Abstract
Oncogenic BRAF V600E mutation confers constitutive activation for the kinase and is closely related to the pathogenesis of metastasic bowel cancer (MBC). Here, the V600E-induced conformational change in MBC BRAF kinase domain is characterized systematically at structural, energetic and dynamic levels. The mutation is observed to cause a conformational conversion of the kinase's activation loop from DFG-out to DFG-in, thus activating the kinase. Electrostatic force is primarily responsible for the conformational conversion and stabilization of DFG-in associated with the mutation. Molecular docking calculations are employed to analyze the binding mode difference of mutant-selective inhibitors between the DFG-out and DFG-in conformations of BRAF kinase. It is revealed that the mutation can reshape inhibitor selectivity profile by altering kinase loop conformation. Several compounds are determined to have a high or moderate selectivity for mutant over wild-type kinase. The selectivity is primarily originated from hydrogen bond interactions of inhibitor ligands with mutant rather than wild type due to the conformational difference in kinase domain.
Collapse
Affiliation(s)
- Kai Zhao
- Department of Gastroenterology, People's Hospital of Jintan, Changzhou 213200, China
| | - Xin Zhou
- Department of Gastroenterology, People's Hospital of Jintan, Changzhou 213200, China
| | - Ming Ding
- Department of Respiration, The Affiliated Hospital of Jiangsu University, Zhenjiang 212001, China.
| |
Collapse
|
41
|
Onufriev AV, Izadi S. Water models for biomolecular simulations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1347] [Citation(s) in RCA: 94] [Impact Index Per Article: 13.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Alexey V. Onufriev
- Department of Physics; Virginia Tech; Blacksburg VA USA
- Department of Computer Science; Virginia Tech; Blacksburg VA USA
- Center for Soft Matter and Biological Physics; Virginia Tech; Blacksburg VA USA
| | - Saeed Izadi
- Early Stage Pharmaceutical Development; Genentech Inc.; South San Francisco, CA USA
| |
Collapse
|
42
|
Gantchev TG, Petkov PS, Hunting DJ. Conformational rearrangement of 1,2-d(GG) intrastrand cis-diammineplatinum crosslinked DNA is driven by counter-ion penetration within the minor groove of the modified site. J Mol Model 2017; 23:278. [PMID: 28913561 DOI: 10.1007/s00894-017-3445-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2016] [Accepted: 08/25/2017] [Indexed: 11/27/2022]
Abstract
The major structural aberrations of DNA induced by a cis-diammineplatinum (II) 1,2-d(GG) intrastrand cross-link (CPT) have been known for decades. To gain deeper insights into the structural dynamics of the sequence-dependent DNA distortions adjacent to the CPT adduct, we employed molecular modeling and molecular dynamics (MD) simulations. The structural dynamics of native (N-DNA) and cisPt 1,2-d(GG) crosslinked (CPT-DNA) in the form of symmetric 36 nt d(G2T15G*G*T15G2)●C2A15CCA15C2) oligonucleotide duplexes is compared. The selected sequence context enabled tracking of the origin of the DNA axis curvature at the YpR flexible points (N-DNA), the enhancement of axis bending, and further distortions due to steric/electrostatic perturbations arising from the CPT-crosslink. In addition to the known structural distortions of CPT-DNA: helix bend towards the major groove; local helix unwinding; high roll angle between cross-linked guanine bases; and adoption of A-form DNA on the 5'-side of the CPT-crosslink (TpG junction); our results show the existence of a singular irreversible and reproducible conformational rearrangement, not previously observed, resulting in two stable CPT-DNA1 and CPT-DNA2 conformers. The CPT-DNA2 conformation presents an enhanced DNA axis bend and a wider and shallower minor grove with increased solvent accessibility within the modified site. It is concluded that the polymorphous (unstable) DNA environment near the cisPt 1,2-d(GG) unit in synergy with specific dynamic events, such as prolonged minor groove retention of particular Na+ ions and water redistribution within the d(TG*G*T) site, together with the formation of extra and more stable H-bonds between Pt(NH3)2 amines and neighboring nucleotides, are cooperatively responsible for the initiation of the conformational rearrangement leading to the CPT-DNA2 conformer, which, surprisingly, closely resembles the HMGB1-bound CPT-DNA structure. Graphical abstract Superimposed averaged structures of normal (N-DNA, green) and cisplatin intrastrand cross-linked (CPT-DNA, orange). Global DNA axes: N-DNA (blue beads); CPT-DNA (red beads); PT (yellow sphere); crosslinked dGs viewed from the minor groove (bold).
Collapse
Affiliation(s)
- Tsvetan G Gantchev
- Department of Nuclear Medicine & Radiobiology, Faculty of Medicine, Université de Sherbrooke, Sherbrooke, QC, J1H 5N4, Canada. .,"Roumen Tsanev" Institute of Molecular Biology, Bulgarian Academy of Sciences, 1113, Sofia, Bulgaria.
| | - Peicho St Petkov
- Department of Atomic Physics, Faculty of Physics, Sofia University, 5 James Bourchier Blvd., 1164, Sofia, Bulgaria
| | - Darel J Hunting
- Department of Nuclear Medicine & Radiobiology, Faculty of Medicine, Université de Sherbrooke, Sherbrooke, QC, J1H 5N4, Canada
| |
Collapse
|
43
|
Nishikawa N, Sakae Y, Gouda T, Tsujimura Y, Okamoto Y. Two major stable structures of amyloid-forming peptides: amorphous aggregates and amyloid fibrils. MOLECULAR SIMULATION 2017. [DOI: 10.1080/08927022.2017.1359746] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Naohiro Nishikawa
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Japan
- Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, Okazaki, Japan
| | - Yoshitake Sakae
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Japan
| | - Takuya Gouda
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Japan
| | - Yuichiro Tsujimura
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Japan
| | - Yuko Okamoto
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Japan
- Structural Biology Research Center, Graduate School of Science, Nagoya University, Nagoya, Japan
- Center for Computational Science, Graduate School of Engineering, Nagoya University, Nagoya, Japan
- Information Technology Center, Nagoya University, Nagoya, Japan
- JST-CREST, Nagoya, Japan
| |
Collapse
|
44
|
Sun Y, Ai X, Hou J, Ye X, Liu R, Shen S, Li Z, Lu S. Integrated discovery of FOXO1-DNA stabilizers from marine natural products to restore chemosensitivity to anti-EGFR-based therapy for metastatic lung cancer. MOLECULAR BIOSYSTEMS 2017; 13:330-337. [PMID: 27966721 DOI: 10.1039/c6mb00678g] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
The transcription factor forkhead box O1 (FOXO1) negatively regulates activated EGFR signaling by turning on the gene expression of tumor suppressor Kruppel-like factor 6. Here, we propose that the chemosensitivity to anti-EGFR-based lung cancer therapy can be restored by stabilization of the FOXO1-DNA complex architecture using small-molecule marine natural medicines. A synthetic protocol that integrates computational ligand-protein-DNA binding analysis and an experimental fluorescence binding assay was applied against a large library of structurally diverse, drug-like marine natural products to discover novel stabilizers of DNA-bound FOXO1 conformation. The screening utilized chemical similarity analysis to exclude structurally redundant compounds, and then carried out high-throughput molecular docking and computational binding analysis to identify potential marine natural product candidates. Consequently, eight commercially available hits were selected and tested in vitro, from which four marine natural product compounds (tanzawaic acid D, hymenidin, cribrostatin 6 and barbamide) were found to have high or moderate potency to selectively bind to the FOXO1 DNA-binding domain (DBD) in the presence of its cognate DNA partner. Atomistic molecular dynamics (MD) simulations revealed that the identified stabilizers do not directly interact with DNA; instead, they can effectively stabilize the free FOXO1 DBD domain in the DNA-bound conformation and thus promote the binding of FOXO1 to DNA.
Collapse
Affiliation(s)
- Yingjia Sun
- Oncology Department, Shanghai Chest Hospital, Shanghai Jiaotong University, Shanghai 200030, P. R. China.
| | - Xinghao Ai
- Oncology Department, Shanghai Chest Hospital, Shanghai Jiaotong University, Shanghai 200030, P. R. China.
| | - Jingwen Hou
- Instrumental Analysis Center, Shanghai Jiaotong University, Shanghai 200240, P. R. China
| | - Xiangyun Ye
- Oncology Department, Shanghai Chest Hospital, Shanghai Jiaotong University, Shanghai 200030, P. R. China.
| | - Ruijun Liu
- Oncology Department, Shanghai Chest Hospital, Shanghai Jiaotong University, Shanghai 200030, P. R. China.
| | - Shengping Shen
- Oncology Department, Shanghai Chest Hospital, Shanghai Jiaotong University, Shanghai 200030, P. R. China.
| | - Ziming Li
- Oncology Department, Shanghai Chest Hospital, Shanghai Jiaotong University, Shanghai 200030, P. R. China.
| | - Shun Lu
- Oncology Department, Shanghai Chest Hospital, Shanghai Jiaotong University, Shanghai 200030, P. R. China.
| |
Collapse
|
45
|
Mizuhara Y, Parkin D, Umezawa K, Ohnuki J, Takano M. Over-Destabilization of Protein-Protein Interaction in Generalized Born Model and Utility of Energy Density Integration Cutoff. J Phys Chem B 2017; 121:4669-4677. [PMID: 28426223 DOI: 10.1021/acs.jpcb.7b01438] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The generalize Born (GB) model is frequently used in MD simulations of biomolecular systems in aqueous solution. The GB model is usually based on the so-called Coulomb field approximation (CFA) for the energy density integration. In this study, we report that the GB model with CFA overdestabilizes the long-range electrostatic attraction between oppositely charged molecules (ionic bond forming two-helix system and kinesin-tubulin system) when the energy density integration cutoff, rmax, which is used to calculate the Born energy, is set to a large value. We show that employing large rmax, which is usually expected to make simulation results more accurate, worsens the accuracy so that the attraction is changed into repulsion. It is demonstrated that the overdestabilization is caused by the overestimation of the desolvation penalty upon binding that originates from CFA. We point out that the overdestabilization can be corrected by employing a relatively small cutoff (rmax = 10-15 Å), affirming that the GB models, even with CFA, can be used as a powerful tool to theoretically study the protein-protein interaction, particularly on its dynamical aspect, such as binding and unbinding.
Collapse
Affiliation(s)
- Yukinobu Mizuhara
- Department of Pure and Applied Physics, Waseda University , Okubo 3-4-1, Shinjuku-Ku, Tokyo 169-8555, Japan
| | - Dan Parkin
- Department of Advanced Science and Engineering, Waseda University , Okubo 3-4-1, Shinjuku-Ku, Tokyo 169-8555, Japan
| | | | - Jun Ohnuki
- Department of Pure and Applied Physics, Waseda University , Okubo 3-4-1, Shinjuku-Ku, Tokyo 169-8555, Japan
| | - Mitsunori Takano
- Department of Pure and Applied Physics, Waseda University , Okubo 3-4-1, Shinjuku-Ku, Tokyo 169-8555, Japan
| |
Collapse
|
46
|
Devi M, Lyngdoh RD. Favored and less favored codon–anticodon duplexes arising from the GC codon family box encoding for alanine: some computational perspectives. J Biomol Struct Dyn 2017; 36:1029-1049. [DOI: 10.1080/07391102.2017.1308886] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Martina Devi
- Department of Chemistry, North-Eastern Hill University, Shillong 793022, India
| | - R.H. Duncan Lyngdoh
- Department of Chemistry, North-Eastern Hill University, Shillong 793022, India
| |
Collapse
|
47
|
Aytenfisu A, Spasic A, Grossfield A, Stern HA, Mathews DH. Revised RNA Dihedral Parameters for the Amber Force Field Improve RNA Molecular Dynamics. J Chem Theory Comput 2017; 13:900-915. [PMID: 28048939 PMCID: PMC5312698 DOI: 10.1021/acs.jctc.6b00870] [Citation(s) in RCA: 81] [Impact Index Per Article: 11.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2016] [Indexed: 12/27/2022]
Abstract
The backbone dihedral parameters of the Amber RNA force field were improved by fitting using multiple linear regression to potential energies determined by quantum chemistry calculations. Five backbone and four glycosidic dihedral parameters were fit simultaneously to reproduce the potential energies determined by a high-level density functional theory calculation (B97D3 functional with the AUG-CC-PVTZ basis set). Umbrella sampling was used to determine conformational free energies along the dihedral angles, and these better agree with the population of conformations observed in the protein data bank for the new parameters than for the conventional parameters. Molecular dynamics simulations performed on a set of hairpin loops, duplexes and tetramers with the new parameter set show improved modeling for the structures of tetramers CCCC, CAAU, and GACC, and an RNA internal loop of noncanonical pairs, as compared to the conventional parameters. For the tetramers, the new parameters largely avoid the incorrect intercalated structures that dominate the conformational samples from the conventional parameters. For the internal loop, the major conformation solved by NMR is stable with the new parameters, but not with the conventional parameters. The new force field performs similarly to the conventional parameters for the UUCG and GCAA hairpin loops and the [U(UA)6A]2 duplex.
Collapse
Affiliation(s)
- Asaminew
H. Aytenfisu
- Department
of Biochemistry & Biophysics, University
of Rochester Medical Center, Rochester, New York 14642, United States
- Center
for RNA Biology, University of Rochester
Medical Center, Rochester, New York 14642, United States
| | - Aleksandar Spasic
- Department
of Biochemistry & Biophysics, University
of Rochester Medical Center, Rochester, New York 14642, United States
- Center
for RNA Biology, University of Rochester
Medical Center, Rochester, New York 14642, United States
| | - Alan Grossfield
- Department
of Biochemistry & Biophysics, University
of Rochester Medical Center, Rochester, New York 14642, United States
| | - Harry A. Stern
- Center
for Integrated Research Computing, University
of Rochester, Rochester, New York 14642, United
States
| | - David H. Mathews
- Department
of Biochemistry & Biophysics, University
of Rochester Medical Center, Rochester, New York 14642, United States
- Center
for RNA Biology, University of Rochester
Medical Center, Rochester, New York 14642, United States
- Department
of Biostatistics & Computational Biology, University of Rochester Medical Center, Rochester, New York 14642, United States
| |
Collapse
|
48
|
Le MT, Kasprzak WK, Kim T, Gao F, Young MYL, Yuan X, Shapiro BA, Seog J, Simon AE. Folding behavior of a T-shaped, ribosome-binding translation enhancer implicated in a wide-spread conformational switch. eLife 2017; 6:e22883. [PMID: 28186489 PMCID: PMC5336357 DOI: 10.7554/elife.22883] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2016] [Accepted: 02/07/2017] [Indexed: 11/29/2022] Open
Abstract
Turnip crinkle virus contains a T-shaped, ribosome-binding, translation enhancer (TSS) in its 3'UTR that serves as a hub for interactions throughout the region. The viral RNA-dependent RNA polymerase (RdRp) causes the TSS/surrounding region to undergo a conformational shift postulated to inhibit translation. Using optical tweezers (OT) and steered molecular dynamic simulations (SMD), we found that the unusual stability of pseudoknotted element H4a/Ψ3 required five upstream adenylates, and H4a/Ψ3 was necessary for cooperative association of two other hairpins (H5/H4b) in Mg2+. SMD recapitulated the TSS unfolding order in the absence of Mg2+, showed dependence of the resistance to pulling on the 3D orientation and gave structural insights into the measured contour lengths of the TSS structure elements. Adenylate mutations eliminated one-site RdRp binding to the 3'UTR, suggesting that RdRp binding to the adenylates disrupts H4a/Ψ3, leading to loss of H5/H4b interaction and promoting a conformational switch interrupting translation and promoting replication.
Collapse
Affiliation(s)
- My-Tra Le
- Department of Cell Biology and Molecular Genetics, University of Maryland, College Park, United States
| | - Wojciech K Kasprzak
- Basic Science Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research, Frederick, United States
| | - Taejin Kim
- RNA Biology Laboratory, Center for Cancer Research, National Cancer Institute, Frederick, United States
| | - Feng Gao
- Department of Cell Biology and Molecular Genetics, University of Maryland, College Park, United States
| | - Megan YL Young
- Department of Cell Biology and Molecular Genetics, University of Maryland, College Park, United States
| | - Xuefeng Yuan
- Department of Cell Biology and Molecular Genetics, University of Maryland, College Park, United States
| | - Bruce A Shapiro
- RNA Biology Laboratory, Center for Cancer Research, National Cancer Institute, Frederick, United States
| | - Joonil Seog
- Department of Materials Science and Engineering, University of Maryland, College Park, United States
| | - Anne E Simon
- Department of Cell Biology and Molecular Genetics, University of Maryland, College Park, United States
| |
Collapse
|
49
|
Chen X, Xie B, Cao L, Zhu F, Chen B, Lv H, Fan X, Han L, Bie L, Cao X, Shen X, Cao F. Direct binding of microRNA-21 pre-element with Regorafenib: An alternative mechanism for anti-colorectal cancer chemotherapy? J Mol Graph Model 2017; 73:48-53. [PMID: 28236743 DOI: 10.1016/j.jmgm.2017.02.004] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/01/2017] [Revised: 02/05/2017] [Accepted: 02/08/2017] [Indexed: 02/07/2023]
Abstract
The Regorafenib is a broad-spectrum kinase inhibitor that has been approved to treat colorectal cancer (CRC). However, evidences have shown that the agent is also implicated in drug interaction with microRNA-21 (miR-21), an oncogenic miRNA which plays a key role in resisting programmed cell death in CRC cells. Here, we supposed that, instead of kinase inhibition, Regorafenib can directly bind to and then stabilize miR-21 pre-element, thus preventing RNase Dicer-meditated cleavage of the pre-element to mature miR-21. In order to verify the notion, an in silico-in vitro integrated investigation of the direct intermolecular interaction between Regorafenib and miR-21 pre-element was performed by using active pocket identification, RNA-ligand docking, molecular dynamics (MD) simulation, binding energetic analysis, and fluorescence-based assay. It was revealed that the Regorafenib can bind at the major groove-like stem region of miR-21 pre-element through three geometrically satisfactory hydrogen bonds (H-bonds) as well as a number of hydrophobic forces and π-π stacking, conferring strong specificity and high stability to the RNA-ligand complex system (Kd=0.73μM). Separate inversion mutation of two base pairs (G6C, C12G) and (A13U, U4A) that are involved in the H-bonding can considerably impair the affinity of Regorafenib to miR-21 pre-element, with Kd increase to 27 and 96μM, respectively. All these supported that Regorafenib can directly bind to miR-21 pre-element at molecular level and the binding mode can be properly modeled by using the proposed integrated strategy. This study would provide a potential, alternative mechanism for anti-colorectal cancer chemotherapy with Regorafenib.
Collapse
Affiliation(s)
- Xiaobing Chen
- Department of Internal Oncology, Hennan Cancer Hospital, The Affiliated Cancer Hospital of Zhengzhou University, Zhengzhou 450008, China.
| | - Bojian Xie
- Department of Surgical Oncology, Taizhou Hospital, Wenzhou Medical University, Taizhou 317000, China
| | - Liang Cao
- The Second Affiliated Hospital of Zhengzhou University, Zhengzhou 450014, China
| | - Feng Zhu
- Taizhou Hospital, Wenzhou Medical University, Taizhou 317000, China
| | - Beibei Chen
- Department of Internal Oncology, Hennan Cancer Hospital, The Affiliated Cancer Hospital of Zhengzhou University, Zhengzhou 450008, China
| | - Huifang Lv
- Department of Internal Oncology, Hennan Cancer Hospital, The Affiliated Cancer Hospital of Zhengzhou University, Zhengzhou 450008, China
| | - Xingxing Fan
- Department of Internal Oncology, Hennan Cancer Hospital, The Affiliated Cancer Hospital of Zhengzhou University, Zhengzhou 450008, China
| | - Lili Han
- Department of Internal Oncology, Hennan Cancer Hospital, The Affiliated Cancer Hospital of Zhengzhou University, Zhengzhou 450008, China
| | - Liangyu Bie
- Department of Internal Oncology, Hennan Cancer Hospital, The Affiliated Cancer Hospital of Zhengzhou University, Zhengzhou 450008, China
| | - Xinguang Cao
- Department of Internal Oncology, Hennan Cancer Hospital, The Affiliated Cancer Hospital of Zhengzhou University, Zhengzhou 450008, China
| | - Xiaokun Shen
- Department of Surgical Oncology, Taizhou Hospital, Wenzhou Medical University, Taizhou 317000, China.
| | - Feilin Cao
- Department of Surgical Oncology, Taizhou Hospital, Wenzhou Medical University, Taizhou 317000, China.
| |
Collapse
|
50
|
Zhang R, Erler J, Langowski J. Histone Acetylation Regulates Chromatin Accessibility: Role of H4K16 in Inter-nucleosome Interaction. Biophys J 2017; 112:450-459. [PMID: 27931745 PMCID: PMC5300776 DOI: 10.1016/j.bpj.2016.11.015] [Citation(s) in RCA: 92] [Impact Index Per Article: 13.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2016] [Revised: 09/30/2016] [Accepted: 11/14/2016] [Indexed: 12/17/2022] Open
Abstract
The N-terminal tail of histone H4 is an indispensable mediator for inter-nucleosome interaction, which is required for chromatin fiber condensation. H4K16 acetylation (H4K16Ac) activates gene transcription by influencing both chromatin structure and interplay with nonhistone proteins. To understand the influence of H4K16Ac on inter-nucleosome interaction, we performed a simulation study for the H4 tail in the context of two nucleosomes in neighboring unit cells in the crystal structure. The binding conformation of H4 tail with/without K16Ac was sampled by replica exchange with solute tempering, and the free energy landscape was explored by metadynamics. The results indicate two important features of H4K16: 1) it is the first button to anchor the H4 tail on the adjacent nucleosome; and 2) it is the only acetylation site interacting with the acidic patch. H4K16Ac disrupts the electrostatic interactions of K16, weakens H4 tail-acidic patch binding, and significantly increases H4 tail conformation diversity. Our study suggests that H4K16Ac directly reduces the inter-nucleosome interaction mediated by the H4 tail, which might further encourage the binding of nonhistone proteins on the acidic patch.
Collapse
Affiliation(s)
- Ruihan Zhang
- Division Biophysics of Macromolecules, German Cancer Research Center, Heidelberg, Germany
| | - Jochen Erler
- Division Biophysics of Macromolecules, German Cancer Research Center, Heidelberg, Germany
| | - Jörg Langowski
- Division Biophysics of Macromolecules, German Cancer Research Center, Heidelberg, Germany.
| |
Collapse
|