Accurate prediction of protonation state as a prerequisite for reliable MM-PB(GB)SA binding free energy calculations of HIV-1 protease inhibitors

Binding free energies were calculated for the inhibitors lopinavir, ritonavir, saquinavir, indinavir, amprenavir, and nelfinavir bound to HIV-1 protease. An MMPB/SA-type analysis was applied to conformational samples from 3 ns explicit solvent molecular dynamics simulations of the enzyme-inhibitor complexes. Binding affinities and the sampled conformations of the inhibitor and enzyme were compared between different HIV-1 protease protonation states to find the most likely protonation state of the enzyme in the complex with each of the inhibitors. The resulting set of protonation states leads to good agreement between calculated and experimental binding affinities. Results from the MMPB/SA analysis are compared with an explicit/implicit hybrid scheme and with MMGB/SA methods. It is found that the inclusion of explicit water molecules may offer a slight advantage in reproducing absolute binding free energies while the use of the Generalized Born approximation significantly affects the accuracy of the calculated binding affinities.     

(Thermo)dynamic Role of Receptor Flexibility,Entropy, and Motional Correlation in Protein–Ligand Binding

The binding of 2‐amino‐5‐methylthiazole to the W191G cavity mutant of cytochrome c peroxidase is an ideal test case to investigate the entropic contribution to the binding free energy due to changes in receptor flexibility. The dynamic and thermodynamic role of receptor flexibility are studied by 50 ns‐long explicit‐solvent molecular dynamics simulations of three separate receptor ensembles: W191G binding a K⁺ ion, W191G–2a5mt complex with a closed 190–195 gating loop, and apo with an open loop. We employ a method recently proposed to estimate accurate absolute single‐molecule configurational entropies and their differences for systems undergoing conformational transitions. We find that receptor flexibility plays a generally underestimated role in protein–ligand binding (thermo)dynamics and that changes of receptor motional correlation determine such large entropy contributions.     

新型二氟甲基磷酸类酪氨酸蛋白磷酸酯酶1B抑制剂的分子动力学模拟和结合自由能计算

通过分子对接建立了一系列含二氟甲基磷酸基团(DFMP)或二氟甲基硫酸基团(DFMS)的抑制剂与酪氨酸蛋白磷酸酯酶1B(PTP1B)的相互作用模式, 并通过1 ns的分子动力学模拟和molecular mechanics/generalized Born surface area (MM/GBSA)方法计算了其结合自由能. 计算获得的结合自由能排序和抑制剂与靶酶间结合能力排序一致; 通过基于主方程的自由能计算方法, 获得了抑制剂与靶酶残基间相互作用的信息, 这些信息显示DFMP/DFMS基团的负电荷中心与PTP1B的221位精氨酸正电荷中心之间的静电相互作用强弱决定了此类抑制剂的活性, 进一步的分析还显示位于DFMP/DFMS基团中的氟原子或其他具有适当原子半径的氢键供体原子会增进此类抑制剂与PTP1B活性位点的结合能力.     

Rapid and accurate prediction of binding free energies for saquinavir-bound HIV-1 proteases

To explain drug resistance by computer simulations at the molecular level, we first have to assess the accuracy of theoretical predictions. Herein we report an application of the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) technique to the ranking of binding affinities of the inhibitor saquinavir with the wild type (WT) and three resistant mutants of HIV-1 protease: L90M, G48V, and G48V/L90M. For each ligand-protein complex we report 10 ns of fully unrestrained molecular dynamics (MD) simulations with explicit solvent. We investigate convergence, internal consistency, and model dependency of MM/PBSA ligand binding energies. Converged enthalpy and entropy estimates produce ligand binding affinities within 1.5 kcal/mol of experimental values, with a remarkable level of correlation to the experimentally observed ranking of resistance levels. A detailed analysis of the enthalpic/entropic balance of drug-protease interactions explains resistance in L90M in terms of a higher vibrational entropy than in the WT complex, while G48V disrupts critical hydrogen bonds at the inhibitor's binding site and produces an altered, more unfavorable balance of Coulomb and polar desolvation energies.     

Actinomycin D binds to single stranded DNA oligomers which contain double GTC triplets

The present work carried out a study on the interactions between Actinomycin D (ActD) and some single-stranded DNA oligomers, which contain double GTC triplets separated by TTT sequence. The interactions of drugs with DNA oligomers were investigated by UV, circular dichroism (CD) spectroscopy and isothermal titration calorimetry (ITC). The results indicate that ActD binds to the single stranded DNA oligomers in the fold back binding model as supported by added A/T base at DNA strand terminal which facilitates the formation of hairpin. The apparent binding constant K_b, the apparent binding molar enthalpy ΔH⁰ and other thermodynamic data were obtained. The binding affinities are sequence dependent and related to the length of DNA strand. And the higher molar binding enthalpy indicates that the binding process is enthalpy driven.     

Molecular dynamics simulations and thermodynamics analysis of DNA-drug complexes. Minor groove binding between 4',6-diamidino-2-phenylindole and DNA duplexes in solution

An extended set of nanosecond-scale molecular dynamics simulations of DNA duplex sequences in explicit solvent interacting with the minor groove binding drug 4',6-diamidino-2-phenylindole (DAPI) are investigated for four different and sequence specific binding modes. Force fields for DAPI have been parametrized to properly reflect its internal nonplanarity. Sequences investigated include the binding modes observed experimentally, that is, AATT in d(CGCGAATTCGCG)(2) and ATTG in d(GGCCAATTGG)(2) and alternative shifted binding modes ATTC and AATT, respectively. In each case, stable MD simulations are obtained, well reproducing specific hydration patterns seen in the experiments. In contrast to the 2.4 A d(CGCGAATTCGCG)(2) crystal structure, the DAPI is nonplanar, consistent with its gas-phase geometry and the higher resolution crystal structure. The simulations also suggest that the DAPI molecule is able to adopt different conformational substates accompanied by specific hydration patterns that include long-residing waters. The MM_PBSA technology for estimating relative free energies was utilized. The most consistent free energy results were obtained with an approach that uses a single trajectory of the DNA-DAPI complex to estimate all free energy terms. It is demonstrated that explicit inclusion of a subset of bound water molecules shifts the calculated relative binding free energies in favor of both crystallographically observed binding modes, underlining the importance of structured hydration.     

Molecular dynamics of cryptophane and its complexes with tetramethylammonium and neopentane using a continuum solvent model

Time scales currently obtainable in explicit–solvent molecular dynamics simulations are inadequate for the study of many biologically important processes. This has led to increased interest in the use of continuum solvent models. For such models to be used effectively, it is important that their behavior relative to explicit simulation be clearly understood. Accordingly, 5 ns stochastic dynamics simulations of a derivative of cryptophane-E alone, and complexed with tetramethylammonium and neopentane were carried out. Solvation electrostatics were accounted for via solutions to the Poisson equation. Nonelectrostatic aspects of solvation were incorporated using a surface area-dependent energy term. Comparison of the trajectories to those from previously reported 25 ns explicit–solvent simulations shows that use of a continuum solvent model results in enhanced sampling. Use of the continuum solvent model also results in a considerable increase in computational efficiency. The continuum solvent model is found to predict qualitative structural characteristics that are similar to those observed in explicit solvent. However, some differences are significant, and optimization of the continuum parameterization will be required for this method to become an efficient alternative to explicit–solvent simulation. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 956–970, 1999     

Insights into the Ligand Binding to Bromodomain-Containing Protein 9 (BRD9): A Guide to the Selection of Potential Binders by Computational Methods

The estimation of the binding of a set of molecules against BRD9 protein was carried out through an in silico molecular dynamics-driven exhaustive analysis to guide the identification of potential novel ligands. Starting from eight crystal structures of this protein co-complexed with known binders and one apo form, we conducted an exhaustive molecular docking/molecular dynamics (MD) investigation. To balance accuracy and an affordable calculation time, the systems were simulated for 100 ns in explicit solvent. Moreover, one complex was simulated for 1 µs to assess the influence of simulation time on the results. A set of MD-derived parameters was computed and compared with molecular docking-derived and experimental data. MM-GBSA and the per-residue interaction energy emerged as the main indicators for the good interaction between the specific binder and the protein counterpart. To assess the performance of the proposed analysis workflow, we tested six molecules featuring different binding affinities for BRD9, obtaining promising outcomes. Further insights were reported to highlight the influence of the starting structure on the molecular dynamics simulations evolution. The data confirmed that a ranking of BRD9 binders using key parameters arising from molecular dynamics is advisable to discard poor ligands before moving on with the synthesis and the biological tests.     

Solvent reaction field potential inside an uncharged globular protein: a bridge between implicit and explicit solvent models?

The solvent reaction field potential of an uncharged protein immersed in simple point charge/extended explicit solvent was computed over a series of molecular dynamics trajectories, in total 1560 ns of simulation time. A finite, positive potential of 13-24 kbTec(-1) (where T=300 K), dependent on the geometry of the solvent-accessible surface, was observed inside the biomolecule. The primary contribution to this potential arose from a layer of positive charge density 1.0 A from the solute surface, on average 0.008 ec/A3, which we found to be the product of a highly ordered first solvation shell. Significant second solvation shell effects, including additional layers of charge density and a slight decrease in the short-range solvent-solvent interaction strength, were also observed. The impact of these findings on implicit solvent models was assessed by running similar explicit solvent simulations on the fully charged protein system. When the energy due to the solvent reaction field in the uncharged system is accounted for, correlation between per-atom electrostatic energies for the explicit solvent model and a simple implicit (Poisson) calculation is 0.97, and correlation between per-atom energies for the explicit solvent model and a previously published, optimized Poisson model is 0.99.     

Multiple protonation equilibria in electrostatics of protein-protein binding

All proteins contain groups capable of exchanging protons with their environment. We present here an approach, based on a rigorous thermodynamic cycle and the partition functions for energy levels characterizing protonation states of the associating proteins and their complex, to compute the electrostatic pH-dependent contribution to the free energy of protein-protein binding. The computed electrostatic binding free energies include the pH of the solution as the variable of state, mutual "polarization" of associating proteins reflected as changes in the distribution of their protonation states upon binding and fluctuations between available protonation states. The only fixed property of both proteins is the conformation; the structure of the monomers is kept in the same conformation as they have in the complex structure. As a reference, we use the electrostatic binding free energies obtained from the traditional Poisson-Boltzmann model, computed for a single macromolecular conformation fixed in a given protonation state, appropriate for given solution conditions. The new approach was tested for 12 protein-protein complexes. It is shown that explicit inclusion of protonation degrees of freedom might lead to a substantially different estimation of the electrostatic contribution to the binding free energy than that based on the traditional Poisson-Boltzmann model. This has important implications for the balancing of different contributions to the energetics of protein-protein binding and other related problems, for example, the choice of protein models for Brownian dynamics simulations of their association. Our procedure can be generalized to include conformational degrees of freedom by combining it with molecular dynamics simulations at constant pH. Unfortunately, in practice, a prohibitive factor is an enormous requirement for computer time and power. However, there may be some hope for solving this problem by combining existing constant pH molecular dynamics algorithms with so-called accelerated molecular dynamics algorithms.     

Secondary structure bias in generalized Born solvent models: comparison of conformational ensembles and free energy of solvent polarization from explicit and implicit solvation

The effects of the use of three generalized Born (GB) implicit solvent models on the thermodynamics of a simple polyalanine peptide are studied via comparing several hundred nanoseconds of well-converged replica exchange molecular dynamics (REMD) simulations using explicit TIP3P solvent to REMD simulations with the GB solvent models. It is found that when compared to REMD simulations using TIP3P the GB REMD simulations contain significant differences in secondary structure populations, most notably an overabundance of alpha-helical secondary structure. This discrepancy is explored via comparison of the differences in the electrostatic component of the free energy of solvation (DeltaDeltaG(pol)) between TIP3P (via thermodynamic Integration calculations), the GB models, and an implicit solvent model based on the Poisson equation (PE). The electrostatic components of the solvation free energies are calculated using each solvent model for four representative conformations of Ala10. Since the PE model is found to have the best performance with respect to reproducing TIP3P DeltaDeltaG(pol) values, effective Born radii from the GB models are compared to effective Born radii calculated with PE (so-called perfect radii), and significant and numerous deviations in GB radii from perfect radii are found in all GB models. The effect of these deviations on the solvation free energy is discussed, and it is shown that even when perfect radii are used the agreement of GB with TIP3P DeltaDeltaG(pol) values does not improve. This suggests a limit to the optimization of the effective Born radius calculation and that future efforts to improve the accuracy of GB models must extend beyond such optimizations.     

Molecular dynamics simulations of large integral membrane proteins with an implicit membrane model

The heterogeneous dielectric generalized Born (HDGB) methodology is an the extension of the GBMV model for the simulation of integral membrane proteins with an implicit membrane environment. Three large integral membrane proteins, the bacteriorhodopsin monomer and trimer and the BtuCD protein, were simulated with the HDGB model in order to evaluate how well thermodynamic and dynamic properties are reproduced. Effects of the truncation of electrostatic interactions were examined. For all proteins, the HDGB model was able to generate stable trajectories that remained close to the starting experimental structures, in excellent agreement with explicit membrane simulations. Dynamic properties evaluated through a comparison of B-factors are also in good agreement with experiment and explicit membrane simulations. However, overall flexibility was slightly underestimated with the HDGB model unless a very large electrostatic cutoff is employed. Results with the HDGB model are further compared with equivalent simulations in implicit aqueous solvent, demonstrating that the membrane environment leads to more realistic simulations.     

Terahertz absorption of DNA decamer duplex

This work combines experimental and theoretical approaches to investigate terahertz absorption spectra of the DNA formed by the sequence oligomer 5'-CCGGCGCCGG-3'. The three-dimensional structure of this self-complimentary DNA decamer has been well-studied, permitting us to perform direct identification of the low-frequency phonon modes associated with specific conformation and to conduct comprehensive computer simulations. Two modeling techniques, normal-mode analysis and nanosecond molecular dynamics with explicit solvent molecules, were employed to extract the low-frequency vibrational modes based on which the absorption spectra were calculated. The absorption spectra of the DNA decamer in aqueous solution were measured in the frequency range 10-25 cm(-1) using the terahertz Fourier transform infrared spectroscopy. Multiple well-resolved and reproducible resonance modes were observed. When calculated and experimental spectra were compared, the spectrum based on molecular dynamics simulations showed a better correlation with the experimental spectra than the one based on normal-mode analysis. These results demonstrate that there exist a considerable number of active low-frequency phonon modes in this short DNA duplex.     

甘氨酸二肽分子酰胺-Ⅰ带光谱与结构相关性

系统探索了蛋白质二肽模型分子——甘氨酸二肽(GLYD)在气相与水溶液中的结构与光谱特性。从分子动力学轨迹中提取具有代表性结构的GLYD-D₂O聚集体的瞬态结构开展简正模式分析,获取了对蛋白质二级结构敏感的酰胺-Ⅰ带的振动光谱参数,建立起振动光谱与特征基团结构间的相关性。将溶剂作用以静电势场的形式投影至二肽分子骨架中,与酰胺-Ⅰ带在气/液相中的频率差相关联,并引入酰胺-Ⅰ带简正模式随二级结构变化的规律,将各个构象态可能存在的振动耦合包含在内,构建具有二级结构敏感性的静电频率转换图,实现溶液相中多肽骨架酰胺-Ⅰ带的快速准确预测。     

Double-stranded DNA dissociates into single strands when dragged into a poor solvent

DNA displays a richness of biologically relevant supramolecular structures, which depend on both sequence and ambient conditions. The effect of dragging double-stranded DNA (dsDNA) from water into poor solvent on the double-stranded structure is still unclear because of condensation. Here, we employed single molecule techniques based on atomic force microscopy and molecular dynamics (MD) simulations to investigate the change in structure and mechanics of DNA during the ambient change. We found that the two strands are split apart when the dsDNA is pulled at one strand from water into a poor solvent. The findings were corroborated by MD simulations where dsDNA was dragged from water into poor solvent, revealing details of the strand separation at the water/poor solvent interface. Because the structure of DNA is of high polarity, all poor solvents show a relatively low polarity. We speculate that the principle of spontaneous unwinding/splitting of dsDNA by providing a low-polarity (in other word, hydrophobic) micro-environment is exploited as one of the catalysis mechanisms of helicases.     

多柔比星稀土金属离子配合物与DNA相互作用的研究

用紫外分光光度法和荧光光谱法研究了多柔比星( Adriamycin,ADM)稀土金属离子配合物(ADM-M)与DNA的相互作用.结果发现,在pH=7.0时,ADM与Eu3+、yb3+能形成稳定配合物,该配合物可使DNA的最大吸收产生明显的减色效应及红移,并能够竞争置换溴化乙锭(EB)与DNA的结合点.KI猝灭试验发现D...     

Modeling and molecular dynamics simulation of the human gonadotropin-releasing hormone receptor in a lipid bilayer

In the present study, a model for the human gonadotropin-releasing hormone receptor embedded in an explicit lipid bilayer was developed. The final conformation was obtained by extensive molecular dynamics simulations of a homology model based on the bovine rhodopsin crystal structure. The analysis of the receptor structure allowed us to detect a number of specific contacts between different amino acid residues, as well as water- and lipid-mediated interactions. These interactions were stable in six additional independent 35 ns long simulations at 310 and 323 K, which used the refined model as the starting structure. All loops, particularly the extracellular loop 2 and the intracellular loop 3, exhibited high fluctuations, whereas the transmembrane helices were more static. Although other models of this receptor have been previously developed, none of them have been subjected to extensive molecular dynamics simulations, and no other three-dimensional structure is publicly available. Our results suggest that the presence of ions as well as explicit solvent and lipid molecules are critical for the structure of membrane protein models, and that molecular dynamics simulations are certainly useful for their refinement.     

Allosteric and Chelate Cooperativity in Divalent Crown Ether/Ammonium Complexes with Strong Binding Enhancement

A thorough thermodynamic analysis by isothermal titration calorimetry of allosteric and chelate cooperativity effects in divalent crown ether/ammonium complexes is combined with DFT calculations including implicit solvent on the one hand and large‐scale molecular dynamics simulations with explicit solvent molecules on the other. The complexes studied exhibit binding constants up to 2×10⁶ m ^?1 with large multivalent binding enhancements and thus strong chelate cooperativity effects. Slight structural changes in the spacers, that is, the exchange of two ether oxygen atoms by two isoelectronic methylene groups, cause significantly stronger binding and substantially increased chelate cooperativity. The analysis is complemented by the examination of solvent effects and allosteric cooperativity. Such a detailed understanding of the binding processes will help to efficiently design and construct larger supramolecular architectures with multiple multivalent building blocks.     

Specific ion binding to macromolecules: effects of hydrophobicity and ion pairing

Using molecular dynamics simulations in an explicit aqueous solvent, we examine the binding of fluoride versus iodide to a spherical macromolecule with both hydrophobic and positively charged patches. Rationalizing our observations, we divide the ion association interaction into two mechanisms: (1) poorly solvated iodide ions are attracted to hydrophobic surface patches, while (2) the strongly solvated fluoride and to a minor extent also iodide bind via cation-anion interactions. Quantitatively, the binding affinities vary significantly with the accessibility of the charged groups as well as the surface potential; therefore, we expect the ion-macromolecule association to be modulated by the local surface characteristics of the (bio-)macromolecule. The observed cation-anion pairing preference is in excellent agreement with experimental data.     

Numerical simulation of the effect of solvent viscosity on the motions of a beta-peptide heptamer

This report examines the effect of a decrease in solvent viscosity on the simulated folding behaviour of a beta-peptide heptamer in methanol. Simulations of the molecular dynamics of the heptamer H-beta3-HVal-beta3-HAla-beta3-HLeu-(S,S)-beta3-HAla(alphaMe)-beta3-HVal-beta3-HAla-beta3-HLeu-OH in methanol, with an explicit representation of the methanol molecules, were performed for 80 ns at various solvent viscosities. The simulations indicate that at a solvent viscosity of one third of that of methanol, only the dynamic aspects of the folding process are altered, and that the rate of folding is increased. At a viscosity of one tenth of that of methanol, insufficient statistics are obtained within the 80 ns period. We suggest that 80 ns is an insufficient time to reach conformational equilibrium at very low viscosity because the dependence of the folding rate of a beta-peptide on solvent viscosity has two regimes; a result that was observed in another computational study for alpha-peptides.