1
|
Sun Z, He Q, Gong Z, Kalhor P, Huai Z, Liu Z. A General Picture of Cucurbit[8]uril Host–Guest Binding: Recalibrating Bonded Interactions. Molecules 2023; 28:molecules28073124. [PMID: 37049887 PMCID: PMC10095826 DOI: 10.3390/molecules28073124] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2023] [Revised: 03/15/2023] [Accepted: 03/27/2023] [Indexed: 04/03/2023] Open
Abstract
Atomic-level understanding of the dynamic feature of host–guest interactions remains a central challenge in supramolecular chemistry. The remarkable guest binding behavior of the Cucurbiturils family of supramolecular containers makes them promising drug carriers. Among Cucurbit[n]urils, Cucurbit[8]uril (CB8) has an intermediate portal size and cavity volume. It can exploit almost all host–guest recognition motifs formed by this host family. In our previous work, an extensive computational investigation of the binding of seven commonly abused and structurally diverse drugs to the CB8 host was performed, and a general dynamic binding picture of CB8-guest interactions was obtained. Further, two widely used fixed-charge models for drug-like molecules were investigated and compared in great detail, aiming at providing guidelines in choosing an appropriate charge scheme in host-guest modelling. Iterative refitting of atomic charges leads to improved binding thermodynamics and the best root-mean-squared deviation from the experimental reference is 2.6 kcal/mol. In this work, we focus on a thorough evaluation of the remaining parts of classical force fields, i.e., the bonded interactions. The widely used general Amber force fields are assessed and refitted with generalized force-matching to improve the intra-molecular conformational preference, and thus the description of inter-molecular host–guest interactions. The interaction pattern and binding thermodynamics show a significant dependence on the modelling parameters. The refitted system-specific parameter set improves the consistency of the modelling results and the experimental reference significantly. Finally, combining the previous charge-scheme comparison and the current force-field refitting, we provide general guidelines for the theoretical modelling of host–guest binding.
Collapse
|
2
|
Negi I, Singh B, Singh Mahmi A, Sharma P. Structural Properties of Hachimoji Nucleic Acids and Their Building Blocks: Comparison of Genetic Systems with Four, Six and Eight Alphabets. Chemphyschem 2023; 24:e202200714. [PMID: 36315394 DOI: 10.1002/cphc.202200714] [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/23/2022] [Revised: 10/29/2022] [Indexed: 11/07/2022]
Abstract
Expansion of the genetic alphabet is an ambitious goal. A recent breakthrough has led to the eight-base (hachimoji) genetics having canonical and unnatural bases. However, very little is known on the molecular-level features that facilitate the candidature of unnatural bases as genetic alphabets. Here we amalgamated DFT calculations and MD simulations to analyse the properties of the constituents of hachimoji DNA and RNA. DFT reveals the dominant syn conformation for isolated unnatural deoxyribonucleosides and at the 5'-end of oligonucleotides, although an anti/syn mixture is predicted at the nonterminal and 3'-terminal positions. However, isolated ribonucleotides prefer an anti/syn mixture, but mostly prefer anti conformation at the nonterminal positions. Further, the canonical base pairing combinations reveals significant strength, which may facilitate replication of hachimoji DNA. We also identify noncanonical base pairs that can better tolerate the substitution of unnatural pairs in RNA. Stacking strengths of 51 dimers reveals higher average stacking stabilization of dimers of hachimoji bases than canonical bases, which provides clues for choosing energetically stable sequences. A total of 14.4 μs MD simulations reveal the influence of solvent on the properties of hachimoji oligonucleotides and point to the likely fidelity of replication of hachimoji DNA. Our results pinpoint the features that explain the experimentally observed stability of hachimoji DNA.
Collapse
Affiliation(s)
- Indu Negi
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh, 160014, India
| | - Bimaldeep Singh
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh, 160014, India
| | - Amanpreet Singh Mahmi
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh, 160014, India
| | - Purshotam Sharma
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh, 160014, India.,Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Ave., Windsor, ON, N9B 3P4, Canada
| |
Collapse
|
3
|
Biaggne A, Kim YC, Melinger JS, Knowlton WB, Yurke B, Li L. Molecular dynamics simulations of cyanine dimers attached to DNA Holliday junctions. RSC Adv 2022; 12:28063-28078. [PMID: 36320263 PMCID: PMC9530999 DOI: 10.1039/d2ra05045e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2022] [Accepted: 09/20/2022] [Indexed: 11/07/2022] Open
Abstract
Dye aggregates and their excitonic properties are of interest for their applications to organic photovoltaics, non-linear optics, and quantum information systems. DNA scaffolding has been shown to be effective at promoting the aggregation of dyes in a controllable manner. Specifically, isolated DNA Holliday junctions have been used to achieve strongly coupled cyanine dye dimers. However, the structural properties of the dimers and the DNA, as well as the role of Holliday junction isomerization are not fully understood. To study the dynamics of cyanine dimers in DNA, molecular dynamics simulations were carried out for adjacent and transverse dimers attached to Holliday junctions in two different isomers. It was found that dyes attached to adjacent strands in the junction exhibit stronger dye-DNA interactions and larger inter-dye separations compared to transversely attached dimers, as well as end-to-end arrangements. Transverse dimers exhibit lower inter-dye separations and more stacked configurations. Furthermore, differences in Holliday junction isomer are analyzed and compared to dye orientations. For transverse dyes exhibiting the smaller inter-dye separations, excitonic couplings were calculated and shown to be in agreement with experiment. Our results suggested that dye attachment locations on DNA Holliday junctions affect dye-DNA interactions, dye dynamics, and resultant dye orientations which can guide the design of DNA-templated cyanine dimers with desired properties. Molecular dynamics simulations reveal dye attachment and DNA Holliday junction isomer effects on dye dimer orientations and excitonic couplings. These simulations can guide synthesis and experiments of dye-DNA structures for excitonic applications.![]()
Collapse
Affiliation(s)
- Austin Biaggne
- Micron School of Materials Science and Engineering, Boise State UniversityBoiseID 83725USA
| | - Young C. Kim
- Materials Science and Technology Division, U.S. Naval Research LaboratoryWashingtonDC20375USA
| | - Joseph. S. Melinger
- Electronics Science and Technology Division, U.S. Naval Research LaboratoryWashingtonDC20375USA
| | - William B. Knowlton
- Micron School of Materials Science and Engineering, Boise State UniversityBoiseID 83725USA,Department of Electrical and Computer Engineering, Boise State UniversityBoiseID 83725USA
| | - Bernard Yurke
- Micron School of Materials Science and Engineering, Boise State UniversityBoiseID 83725USA,Department of Electrical and Computer Engineering, Boise State UniversityBoiseID 83725USA
| | - Lan Li
- Micron School of Materials Science and Engineering, Boise State UniversityBoiseID 83725USA,Center for Advanced Energy StudiesIdaho FallsID 83401USA
| |
Collapse
|
4
|
Sun Z, Zheng L, Wang K, Huai Z, Liu Z. Primary vs secondary: Directionalized guest coordination in β-cyclodextrin derivatives. Carbohydr Polym 2022; 297:120050. [DOI: 10.1016/j.carbpol.2022.120050] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2022] [Revised: 08/17/2022] [Accepted: 08/25/2022] [Indexed: 02/01/2023]
|
5
|
Sun Z, Wang M, He Q, Liu Z. Molecular Modeling of Ionic Liquids: Force‐Field Validation and Thermodynamic Perspective from Large‐Scale Fast‐Growth Solvation Free Energy Calculations. ADVANCED THEORY AND SIMULATIONS 2022. [DOI: 10.1002/adts.202200274] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Zhaoxi Sun
- Beijing National Laboratory for Molecular Sciences College of Chemistry and Molecular Engineering Peking University Beijing 100871 China
| | - Mao Wang
- NCS Testing Technology Co., Ltd. No. 13, Gaoliangqiao Xiejie Beijing 100081 China
| | - Qiaole He
- AI Department of Enzymaster (Ningbo) Bio‐Engineering Co., Ltd. North Century Avenue 333 Ningbo 315100 China
| | - Zhirong Liu
- Beijing National Laboratory for Molecular Sciences College of Chemistry and Molecular Engineering Peking University Beijing 100871 China
| |
Collapse
|
6
|
Sun Z, He Q. Seeding the multi-dimensional nonequilibrium pulling for Hamiltonian variation: indirect nonequilibrium free energy simulations at QM levels. Phys Chem Chem Phys 2022; 24:8800-8819. [PMID: 35352744 DOI: 10.1039/d2cp00355d] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
The combination of free energy simulations in the alchemical and configurational spaces provides a feasible route to access the thermodynamic profiles under a computationally demanding target Hamiltonian. Normally, due to the significant differences between the computational cost of ab initio quantum mechanics (QM) calculations and those of semi-empirical quantum mechanics (SQM) and molecular mechanics (MM), this indirect method could be used to obtain the QM thermodynamics by combining the SQM or MM results and the SQM-to-QM or MM-to-QM corrections. In our previous work, a multi-dimensional nonequilibrium pulling framework for Hamiltonian variations was introduced based on bidirectional pulling and bidirectional reweighting. The method performs nonequilibrium free energy simulations in the configurational space to obtain the thermodynamic profile along the conformational change pathway under a selected computationally efficient Hamiltonian, and uses the nonequilibrium alchemical method to correct or perturb the thermodynamic profile to that under the target Hamiltonian. The BAR-based method is designed to achieve the best generality and transferability and thus leads to modest (∼20 fold) speedup. In this work, we explore the possibility of further accelerating the nonequilibrium free energy simulation by employing unidirectional pulling and using the selection criterion to obtain the initial configurations used to initiate nonequilibrium trajectories following the idea of adaptive steered molecular dynamics (ASMD). A single initial condition is used to seed the whole multi-dimensional nonequilibrium free energy simulation and the sampling is performed fully in the nonequilibrium ensemble. Introducing very short ps-length equilibrium sampling to grab more initial seeds could also be helpful. The ASMD scheme estimates the free energy difference with the unidirectional exponential average (EXP), but it does not follow exactly the requirements of the EXP estimator. Another deficiency of the seeding simulation is the inherently sequential or serial pulling due to the inter-segment dependency, which triggers some problems in the parallelizability of the simulation. Numerical tests are performed to grasp some insights and guidelines for using this selection-criterion-based ASMD scheme. The presented selection-criterion-based multi-dimensional ASMD scheme follows the same perturbation network of the BAR-based method, and thus could be used in various Hamiltonian-variation cases.
Collapse
Affiliation(s)
- Zhaoxi Sun
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Institute of Theoretical and Computational Chemistry, Peking University, Beijing 100871, China.
| | - Qiaole He
- AI Department of Enzymaster (Ningbo) Bio-Engineering Co., Ltd, North Century Avenue 333, 315100 Ningbo, China
| |
Collapse
|
7
|
Sun Z, Huai Z, He Q, Liu Z. A General Picture of Cucurbit[8]uril Host-Guest Binding. J Chem Inf Model 2021; 61:6107-6134. [PMID: 34818004 DOI: 10.1021/acs.jcim.1c01208] [Citation(s) in RCA: 32] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Describing, understanding, and designing complex interaction networks within macromolecular systems remain challenging in modern chemical research. Host-guest systems, despite their relative simplicity in both the structural feature and interaction patterns, still pose problems in theoretical modeling. The barrel-shaped supramolecular container cucurbit[8]uril (CB8) shows promising functionalities in various areas, e.g., catalysis and molecular recognition. It can stably coordinate a series of structurally diverse guests with high affinities. In this work, we examine the binding of seven commonly abused drugs to the CB8 host, aiming at providing a general picture of CB8-guest binding. Extensive sampling of the configurational space of these host-guest systems is performed, and the binding pathway and interaction patterns of CB8-guest complexes are investigated. A thorough comparison of widely used fixed-charge models for drug-like molecules is presented. Iterative refitting of the atomic charges suggests significant conformation dependence of charge generation. The initial model generated at the original conformation could be inaccurate for new conformations explored during conformational search, and the newly fitted charge set improves the prediction-experiment correlation significantly. Our investigations of the configurational space of CB8-drug complexes suggest that the host-guest interactions are more complex than expected. Despite the structural simplicities of these molecules, the conformational fluctuations of the host and the guest molecules and orientations of functional groups lead to the existence of an ensemble of binding modes. The insights of the binding thermodynamics, performance of fixed-charge models, and binding patterns of the CB8-guest systems are useful for studying and elucidating the binding mechanism of other host-guest complexes.
Collapse
Affiliation(s)
- Zhaoxi Sun
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Zhe Huai
- XtalPi-AI Research Center (XARC), 9F, Tower A, Dongsheng Building, No. 8, Zhongguancun East Road, Haidian District, Beijing 100083, P.R. China
| | - Qiaole He
- AI Department of Enzymaster (Ningbo) Bio-Engineering Co., Ltd., North Century Avenue 333, Ningbo 315100, China
| | - Zhirong Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
8
|
Sun Z, Kalhor P, Xu Y, Liu J. Extensive numerical tests of leapfrog integrator in middle thermostat scheme in molecular simulations. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2111242] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Zhaoxi Sun
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Institute of Theoretical and Computational Chemistry, Peking University, Beijing 100871, China
| | - Payam Kalhor
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Institute of Theoretical and Computational Chemistry, Peking University, Beijing 100871, China
| | - Yang Xu
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Institute of Theoretical and Computational Chemistry, Peking University, Beijing 100871, China
| | - Jian Liu
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Institute of Theoretical and Computational Chemistry, Peking University, Beijing 100871, China
| |
Collapse
|
9
|
Sun Z, Liu Z. BAR‐Based Multi‐Dimensional Nonequilibrium Pulling for Indirect Construction of QM/MM Free Energy Landscapes: Varying the QM Region. ADVANCED THEORY AND SIMULATIONS 2021. [DOI: 10.1002/adts.202100185] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Zhaoxi Sun
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering Peking University Beijing 100871 China
| | - Zhirong Liu
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering Peking University Beijing 100871 China
| |
Collapse
|
10
|
|
11
|
Huai Z, Yang H, Sun Z. Binding thermodynamics and interaction patterns of human purine nucleoside phosphorylase-inhibitor complexes from extensive free energy calculations. J Comput Aided Mol Des 2021; 35:643-656. [PMID: 33759016 DOI: 10.1007/s10822-021-00382-w] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Accepted: 03/13/2021] [Indexed: 11/29/2022]
Abstract
Human purine nucleoside phosphorylase (hPNP) plays a significant role in the catabolism of deoxyguanosine. The trimeric protein is an important target in the treatment of T-cell cancers and autoimmune disorders. Experimental studies on the inhibition of the hPNP observe that the first ligand bound to one of three subunits effectively inhibits the protein, while the binding of more ligands to the subsequent sites shows negative cooperativities. In this work, we performed extensive end-point and alchemical free energy calculations to determine the binding thermodynamics of the trimeric protein-ligand system. 13 Immucillin inhibitors with experimental results are under calculation. Two widely accepted charge schemes for small molecules including AM1-BCC and RESP are adopted for ligands. The results of RESP are in better agreement with the experimental reference. Further investigations of the interaction networks in the protein-ligand complexes reveal that several residues play significant roles in stabilizing the complex structure. The most commonly observed ones include PHE200, GLU201, MET219, and ASN243. The conformations of the protein in different protein-ligand complexes are observed to be similar. We expect these insights to aid the development of potent drugs targeting hPNP.
Collapse
Affiliation(s)
- Zhe Huai
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai, 200062, China
| | - Huaiyu Yang
- College of Engineering, Hebei Normal University, Shijiazhuang, 050024, China
| | - Zhaoxi Sun
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai, 200062, China.
| |
Collapse
|
12
|
Huai Z, Shen Z, Sun Z. Binding Thermodynamics and Interaction Patterns of Inhibitor-Major Urinary Protein-I Binding from Extensive Free-Energy Calculations: Benchmarking AMBER Force Fields. J Chem Inf Model 2020; 61:284-297. [PMID: 33307679 DOI: 10.1021/acs.jcim.0c01217] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Mouse major urinary protein (MUP) plays a key role in the pheromone communication system. The one-end-closed β-barrel of MUP-I forms a small, deep, and hydrophobic central cavity, which could accommodate structurally diverse ligands. Previous computational studies employed old protein force fields and short simulation times to determine the binding thermodynamics or investigated only a small number of structurally similar ligands, which resulted in sampled regions far from the experimental structure, nonconverged sampling outcomes, and limited understanding of the possible interaction patterns that the cavity could produce. In this work, extensive end-point and alchemical free-energy calculations with advanced protein force fields were performed to determine the binding thermodynamics of a series of MUP-inhibitor systems and investigate the inter- and intramolecular interaction patterns. Three series of inhibitors with a total of 14 ligands were simulated. We independently simulated the MUP-inhibitor complexes under two advanced AMBER force fields. Our benchmark test showed that the advanced AMBER force fields including AMBER19SB and AMBER14SB provided better descriptions of the system, and the backbone root-mean-square deviation (RMSD) was significantly lowered compared with previous computational studies with old protein force fields. Surprisingly, although the latest AMBER force field AMBER19SB provided better descriptions of various observables, it neither improved the binding thermodynamics nor lowered the backbone RMSD compared with the previously proposed and widely used AMBER14SB. The older but widely used AMBER14SB actually achieved better performance in the prediction of binding affinities from the alchemical and end-point free-energy calculations. We further analyzed the protein-ligand interaction networks to identify important residues stabilizing the bound structure. Six residues including PHE38, LEU40, PHE90, ALA103, LEU105, and TYR120 were found to contribute the most significant part of protein-ligand interactions, and 10 residues were found to provide favorable interactions stabilizing the bound state. The two AMBER force fields gave extremely similar interaction networks, and the secondary structures also showed similar behavior. Thus, the intra- and intermolecular interaction networks described with the two AMBER force fields are similar. Therefore, AMBER14SB could still be the default option in free-energy calculations to achieve highly accurate binding thermodynamics and interaction patterns.
Collapse
Affiliation(s)
- Zhe Huai
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
| | - Zhaoxi Shen
- Institute of Applied Physics and Materials Engineering, University of Macau, Avenida da Universidade, Taipa, Macau 999078, China
| | - Zhaoxi Sun
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
| |
Collapse
|
13
|
Sharma KD, Kathuria P, Wetmore SD, Sharma P. Can modified DNA base pairs with chalcogen bonding expand the genetic alphabet? A combined quantum chemical and molecular dynamics simulation study. Phys Chem Chem Phys 2020; 22:23754-23765. [PMID: 33063082 DOI: 10.1039/d0cp04921b] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A comprehensive (DFT and MD) computational study is presented with the goal to design and analyze model chalcogen-bonded modified nucleobase pairs that replace one (i.e., AXY:T, G:CXY, GXY:C) or two (GXY:CX'Y', X/X' = S, Se and Y/Y' = F, Cl, Br) Watson-Crick (WC) hydrogen bonds of the canonical A:T or G:C pair with chalcogen bond(s). DFT calculations on 18 base pair combinations that replace one WC hydrogen bond with a chalcogen bond reveal that the bases favorably interact in the gas phase (binding strengths up to -140 kJ mol-1) and water (up to -85 kJ mol-1). Although the remaining hydrogen bond(s) exhibits similar characteristics to those in the canonical base pairs, the structural features of the (Y-XO) chalcogen bond(s) change significantly with the identity of X and Y. The 36 doubly-substituted (GXY:CX'Y') base pairs have structural deviations from canonical G:C similar to those of the singly-substituted modifications (G:CXY or GXY:C). Furthermore, despite the replacement of two strong hydrogen bonds with chalcogen bonds, some GXY:CX'Y' pairs possess comparable binding energies (up to -132 kJ mol-1 in the gas phase and up to -92 kJ mol-1 in water) to the most stable G:CXY or GXY:C pairs, as well as canonical G:C. More importantly, G:C-modified pairs containing X = Se (high polarizability) and Y = F (high electronegativity) are the most stable, with comparable or slightly larger (by up to 13 kJ mol-1) binding energies than G:C. Further characterization of the chalcogen bonding in all modified base pairs (AIM, NBO and NCI analyses) reveals that the differences in the binding energies of modified base pairs are mainly dictated by the differences in the strengths of their chalcogen bonds. Finally, MD simulations on DNA oligonucleotides containing the most stable chalcogen-bonded base pair from each of the four classifications (AXY:T, G:CXY, GXY:C and GXY:CX'Y') reveal that the singly-modified G:C pairs best retain the local helical structure and pairing stability to a greater extent than the modified A:T pair. Overall, our study identifies two (G:CSeF and GSeF:C) promising pairs that retain chalcogen bonding in DNA and should be synthesized and further explored in terms of their potential to expand the genetic alphabet.
Collapse
Affiliation(s)
- Karan Deep Sharma
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh, 160014, India. and Department of Chemistry, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - Preetleen Kathuria
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh, 160014, India.
| | - Stacey D Wetmore
- Department of Chemistry and Biochemistry, University of Lethbridge, Alberta T1K 3M4, Canada.
| | - Purshotam Sharma
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh, 160014, India.
| |
Collapse
|
14
|
SAMPL7 TrimerTrip host-guest binding affinities from extensive alchemical and end-point free energy calculations. J Comput Aided Mol Des 2020; 35:117-129. [PMID: 33037549 DOI: 10.1007/s10822-020-00351-9] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2020] [Accepted: 09/30/2020] [Indexed: 12/12/2022]
Abstract
The prediction of host-guest binding affinities with computational modelling is still a challenging task. In the 7th statistical assessment of the modeling of proteins and ligands (SAMPL) challenge, a new host named TrimerTrip was synthesized and the thermodynamic parameters of 16 structurally diverse guests binding to the host were characterized. In the TrimerTrip-guest challenge, only structures of the host and the guests are provided, which indicates that the predictions of both the binding poses and the binding affinities are under assessment. In this work, starting from the binding poses obtained from our previous enhanced sampling simulations in the configurational space, we perform extensive alchemical and end-point free energy calculations to calculate the host-guest binding affinities retrospectively. The alchemical predictions with two widely accepted charge schemes (i.e. AM1-BCC and RESP) are in good agreement with the experimental reference, while the end-point estimates perform poorly in reproducing the experimental binding affinities. Aside from the absolute value of the binding affinity, the rank of binding free energies is also crucial in drug design. Surprisingly, the end-point MM/PBSA method seems very powerful in reproducing the experimental rank of binding affinities. Although the length of our simulations is long and the intermediate spacing is dense, the convergence behavior is not very good, which may arise from the flexibility of the host molecule. Enhanced sampling techniques in the configurational space may be required to obtain fully converged sampling. Further, as the length of sampling in alchemical free energy calculations already achieves several hundred ns, performing direct simulations of the binding/unbinding event in the physical space could be more useful and insightful. More details about the binding pathway and mechanism could be obtained in this way. The nonequilibrium method could also be a nice choice if one insists to use the alchemical method, as the intermediate sampling is avoided to some extent.
Collapse
|
15
|
Sun Z. SAMPL7 TrimerTrip host-guest binding poses and binding affinities from spherical-coordinates-biased simulations. J Comput Aided Mol Des 2020; 35:105-115. [PMID: 32776199 DOI: 10.1007/s10822-020-00335-9] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2020] [Accepted: 08/04/2020] [Indexed: 12/21/2022]
Abstract
Host-guest binding remains a major challenge in modern computational modelling. The newest 7th statistical assessment of the modeling of proteins and ligands (SAMPL) challenge contains a new series of host-guest systems. The TrimerTrip host binds to 16 structurally diverse guests. Previously, we have successfully employed the spherical coordinates as the collective variables coupled with the enhanced sampling technique metadynamics to enhance the sampling of the binding/unbinding event, search for possible binding poses and calculate the binding affinities in all three host-guest binding cases of the 6th SAMPL challenge. In this work, we report a retrospective study on the TrimerTrip host-guest systems by employing the same protocol to investigate the TrimerTrip host in the SAMPL7 challenge. As no binding pose is provided by the SAMPL7 host, our simulations initiate from randomly selected configurations and are proceeded long enough to obtain converged free energy estimates and search for possible binding poses. The calculated binding affinities are in good agreement with the experimental reference, and the obtained binding poses serve as a nice starting point for end-point or alchemical free energy calculations. Note that as the work is performed after the close of the SAMPL7 challenge, we do not participate in the challenge and the results are not formally submitted to the SAMPL7 challenge.
Collapse
Affiliation(s)
- Zhaoxi Sun
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China.
| |
Collapse
|
16
|
SAMPL6 host-guest binding affinities and binding poses from spherical-coordinates-biased simulations. J Comput Aided Mol Des 2020; 34:589-600. [PMID: 31974852 DOI: 10.1007/s10822-020-00294-1] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2019] [Accepted: 01/17/2020] [Indexed: 10/25/2022]
Abstract
Host-guest binding is a challenging problem in computer simulation. The prediction of binding affinities between hosts and guests is an important part of the statistical assessment of the modeling of proteins and ligands (SAMPL) challenges. In this work, the volume-based variant of well-tempered metadynamics is employed to calculate the binding affinities of the host-guest systems in the SAMPL6 challenge. By biasing the spherical coordinates describing the relative position of the host and the guest, the initial-configuration-induced bias vanishes and all possible binding poses are explored. The agreement between the predictions and the experimental results and the observation of new binding poses indicate that the volume-based technique serves as a nice candidate for the calculation of binding free energies and the search of the binding poses.
Collapse
|
17
|
Sun Z, Wang X, Zhang JZH. Theoretical understanding of the thermodynamics and interactions in transcriptional regulator TtgR-ligand binding. Phys Chem Chem Phys 2019; 22:1511-1524. [PMID: 31872826 DOI: 10.1039/c9cp05980f] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The transcriptional regulator TtgR belongs to the TetR family of transcriptional repressors. It depresses the transcription of the TtgABC operon and itself and thus regulates the extrusion of noxious chemicals with efflux pumps in bacterial cells. As the ligand-binding domain of TtgR is rather flexible, it can bind with a number of structurally diverse ligands, such as antibiotics, flavonoids and aromatic solvents. In the current work, we perform equilibrium and nonequilibrium alchemical free energy simulation to predict the binding affinities of a series of ligands targeting the TtgR protein and an agreement between the theoretical prediction and the experimental result is observed. End-point methods MM/PBSA and MM/GBSA are also employed for comparison. We further study the interaction maps and contacts between the protein and the ligand and identify important interactions in the protein-ligand binding cases. The dynamics fluctuation and secondary structures are also investigated. The current work sheds light on atomic and thermodynamic understanding of the TtgR-ligand interactions.
Collapse
Affiliation(s)
- Zhaoxi Sun
- Computational Biomedicine (IAS-5/INM-9), Forschungszentrum Jülich, Jülich 52425, Germany. and State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Xiaohui Wang
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China and Institute of Computational Science, Università della Svizzera italiana (USI), Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
| | - John Z H Zhang
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China and NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China and Department of Chemistry, New York University, NY, NY 10003, USA.
| |
Collapse
|
18
|
Sun Z. BAR-based multi-dimensional nonequilibrium pulling for indirect construction of QM/MM free energy landscapes: from semi-empirical to ab initio. Phys Chem Chem Phys 2019; 21:21942-21959. [PMID: 31552953 DOI: 10.1039/c9cp04113c] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
The indirect method for the construction of quantum mechanics (QM)/molecular mechanics (MM) free energy landscapes provides a cheaper alternative for free energy simulations at the QM level. The indirect method features a direct calculation of the free energy profile with a computationally efficient but less accurate Hamiltonian (i.e. low-level Hamiltonian) and a low-level-to-high-level correction. In the thermodynamic cycle, the direct low-level calculation along the physically meaningful reaction coordinate is corrected via the alchemical method, which is often achieved with perturbation-based techniques. In our previous work, a multi-dimensional nonequilibrium pulling framework is proposed for the indirect construction of QM/MM free energy landscapes. Previously, we focused on obtaining semi-empirical QM (SQM) results indirectly from direct MM simulations and MM to SQM corrections. In this work, we apply this method to obtain results under ab initio QM Hamiltonians by combining direct SQM results and SQM to QM corrections. A series of SQM and QM Hamiltonians are benchmarked. It is observed that PM6 achieves the best performance among the low-level Hamiltonians. Therefore, we recommend using PM6 as the low-level theory in the indirect free energy simulation. Considering its higher similarity to the high-level Hamiltonians, PM6 corrected with the bond charge correction could be more accurate than the existing AM1-BCC model. Another central result in the current work is a basic protocol of choosing the strength of restraints and an appropriate time step in nonequilibrium free energy simulation at the stiff spring limit. We provide theoretical derivations to emphasize the importance of using a sufficiently large force constant and choosing an appropriate time step. It is worth noting that a general rule of thumb for choosing the time step, according to our derivation, is that a time step of 1 fs or smaller should be used, as long as the stiff spring approximation is employed, even in simulations with constraints on bonds involving hydrogen atoms.
Collapse
Affiliation(s)
- Zhaoxi Sun
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| |
Collapse
|