Hot spot computational identification: Application to the complex formed between the hen egg white lysozyme (HEL) and the antibody HyHEL‐10

The definition and comprehension of the hot spots in an interface is a subject of primary interest for a variety of fields, including structure‐based drug design. Therefore, to achieve an alanine mutagenesis computational approach that is at the same time accurate and predictive, capable of reproducing the experimental mutagenesis values is a major challenge in the computational biochemistry field. Antibody/protein antigen complexes provide one of the greatest models to study protein–protein recognition process because they have three fundamentally features: specificity, high complementary association and a small epitope restricted to the diminutive complementary determining regions (CDR) region, while the remainder of the antibody is largely invariant. Thus, we apply a computational mutational methodological approach to the study of the antigen–antibody complex formed between the hen egg white lysozyme (HEL) and the antibody HyHEL‐10. A critical evaluation that focuses essentially on the limitations and advantages between different computational methods for hot spot determination, as well as between experimental and computational methodological approaches, is presented. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Fast and accurate predictions of binding free energies using MM‐PBSA and MM‐GBSA

In the drug discovery process, accurate methods of computing the affinity of small molecules with a biological target are strongly needed. This is particularly true for molecular docking and virtual screening methods, which use approximated scoring functions and struggle in estimating binding energies in correlation with experimental values. Among the various methods, MM‐PBSA and MM‐GBSA are emerging as useful and effective approaches. Although these methods are typically applied to large collections of equilibrated structures of protein‐ligand complexes sampled during molecular dynamics in water, the possibility to reliably estimate ligand affinity using a single energy‐minimized structure and implicit solvation models has not been explored in sufficient detail. Herein, we thoroughly investigate this hypothesis by comparing different methods for the generation of protein‐ligand complexes and diverse methods for free energy prediction for their ability to correlate with experimental values. The methods were tested on a series of structurally diverse inhibitors of Plasmodium falciparum DHFR with known binding mode and measured affinities. The results showed that correlations between MM‐PBSA or MM‐GBSA binding free energies with experimental affinities were in most cases excellent. Importantly, we found that correlations obtained with the use of a single protein‐ligand minimized structure and with implicit solvation models were similar to those obtained after averaging over multiple MD snapshots with explicit water molecules, with consequent save of computing time without loss of accuracy. When applied to a virtual screening experiment, such an approach proved to discriminate between true binders and decoy molecules and yielded significantly better enrichment curves. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Effect of sampling on BACE‐1 ligands binding free energy predictions via MM‐PBSA calculations

The BACE‐1 enzyme is a prime target to find a cure to Alzheimer's disease. In this article, we used the MM‐PBSA approach to compute the binding free energies of 46 reported ligands to this enzyme. After showing that the most probable protonation state of the catalytic dyad is mono‐protonated (on ASP32), we performed a thorough analysis of the parameters influencing the sampling of the conformational space (in total, more than 35 μs of simulations were performed). We show that ten simulations of 2 ns gives better results than one of 50 ns. We also investigated the influence of the protein force field, the water model, the periodic boundary conditions artifacts (box size), as well as the ionic strength. Amber03 with TIP3P, a minimal distance of 1.0 nm between the protein and the box edges and a ionic strength of I = 0.2 M provides the optimal correlation with experiments. Overall, when using these parameters, a Pearson correlation coefficient of R = 0.84 (R ² = 0.71) is obtained for the 46 ligands, spanning eight orders of magnitude of K _d (from 0.017 nm to 2000 μM, i.e., from −14.7 to −3.7 kcal/mol), with a ligand size from 22 to 136 atoms (from 138 to 937 g/mol). After a two‐parameter fit of the binding affinities for 12 of the ligands, an error of RMSD = 1.7 kcal/mol was obtained for the remaining ligands. © 2017 Wiley Periodicals, Inc.     

Comparison of thermodynamic integration and Bennett acceptance ratio for calculating relative protein‐ligand binding free energies

The performances of Bennett's acceptance ratio method and thermodynamic integration (TI) for the calculation of free energy differences in protein simulations are compared. For the latter, the standard trapezoidal rule, Simpson's rule, and Clenshaw‐Curtis integration are used as numerical integration methods. We evaluate the influence of the number and definition of intermediate states on the precision, accuracy, and efficiency of the free energy calculations. Our results show that non‐equidistantly spaced intermediate states are in some cases beneficial for the TI methods. Using several combinations of softness parameters and the λ power dependence, it is shown that these benefits are strongly dependent on the shape of the integrand. Although TI is more user‐friendly due to its simplicity, it was found that Bennett's acceptance ratio method is the more efficient method. It is also the least dependent on the choice of the intermediate states, making it more robust than TI. © 2013 Wiley Periodicals, Inc.     

A semiempirical approach to ligand‐binding affinities: Dependence on the Hamiltonian and corrections

We present a combination of semiempirical quantum‐mechanical (SQM) calculations in the conductor‐like screening model with the MM/GBSA (molecular‐mechanics with generalized Born and surface‐area solvation) method for ligand‐binding affinity calculations. We test three SQM Hamiltonians, AM1, RM1, and PM6, as well as hydrogen‐bond corrections and two different dispersion corrections. As test cases, we use the binding of seven biotin analogues to avidin, nine inhibitors to factor Xa, and nine phenol‐derivatives to ferritin. The results vary somewhat for the three test cases, but a dispersion correction is mandatory to reproduce experimental estimates. On average, AM1 with the DH2 hydrogen‐bond and dispersion corrections gives the best results, which are similar to those of standard MM/GBSA calculations for the same systems. The total time consumption is only 1.3–1.6 times larger than for MM/GBSA. © 2012 Wiley Periodicals, Inc.     

Minimum MD simulation length required to achieve reliable results in free energy perturbation calculations: Case study of relative binding free energies of fructose‐1,6‐bisphosphatase inhibitors

In an attempt to establish the criteria for the length of simulation to achieve the desired convergence of free energy calculations, two studies were carried out on chosen complexes of FBPase‐AMP mimics. Calculations were performed for varied length of simulations and for different starting configurations using both conventional‐ and QM/MM‐FEP methods. The results demonstrate that for small perturbations, 1248 ps simulation time could be regarded a reasonable yardstick to achieve convergence of the results. As the simulation time is extended, the errors associated with free energy calculations also gradually tapers off. Moreover, when starting the simulation from different initial configurations of the systems, the results are not changed significantly, when performed for 1248 ps. This study carried on FBPase‐AMP mimics corroborates well with our previous successful demonstration of requirement of simulation time for solvation studies, both by conventional and ab initio FEP. The establishment of aforementioned criteria of simulation length serves a useful benchmark in drug design efforts using FEP methodologies, to draw a meaningful and unequivocal conclusion. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Interaction of Thymine DNA Glycosylase with Oxidised 5-Methyl-cytosines in Their Amino- and Imino-Forms

Thymine DNA Glycosylase (TDG) is an enzyme of the base excision repair mechanism and removes damaged or mispaired bases from DNA via hydrolysis of the glycosidic bond. Specificity is of high importance for such a glycosylase, so as to avoid the damage of intact DNA. Among the substrates reported for TDG are mispaired uracil and thymine but also formyl-cytosine and carboxyl-cytosine. Methyl-cytosine and hydroxylmethyl-cytosine are, in contrast, not processed by the TDG enzyme. We have in this work employed molecular dynamics simulations to explore the conformational dynamics of DNA carrying a formyl-cytosine or carboxyl-cytosine and compared those to DNA with the non-cognate bases methyl-cytosine and hydroxylmethyl-cytosine, as amino and imino tautomers. Whereas for the mispairs a wobble conformation is likely decisive for recognition, all amino tautomers of formyl-cytosine and carboxyl-cytosine exhibit the same Watson–Crick conformation, but all imino tautomers indeed form wobble pairs. The conformational dynamics of the amino tautomers in free DNA do not exhibit differences that could be exploited for recognition, and also complexation to the TDG enzyme does not induce any alteration that would indicate preferable binding to one or the other oxidised methyl-cytosine. The imino tautomers, in contrast, undergo a shift in the equilibrium between a closed and a more open, partially flipped state, towards the more open form upon complexation to the TDG enzyme. This stabilisation of the more open conformation is most pronounced for the non-cognate bases methyl-cytosine and hydroxyl-cytosine and is thus not a likely mode for recognition. Moreover, calculated binding affinities for the different forms indicate the imino forms to be less likely in the complexed DNA. These findings, together with the low probability of imino tautomers in free DNA and the indifference of the complexed amino tautomers, suggest that discrimination of the oxidised methyl-cytosines does not take place in the initial complex formation.     

Atomistic detailed free‐energy landscape of intrinsically disordered protein studied by multi‐scale divide‐and‐conquer molecular dynamics simulation

Calcineurin (CaN) is a eukaryotic serine/threonine protein phosphatase activated by both Ca²⁺ and calmodulin (CaM), including intrinsically disordered region (IDR). The region undergoes folding into an α‐helix form in the presence Ca²⁺‐loaded CaM. To sample the ordered structure of the IDR by conventional all atom model (AAM) molecular dynamics (MD) simulation, the IDR and Ca²⁺‐loaded CaM must be simultaneously treated. However, it is time‐consuming task because the coupled folding and binding should include repeated binding and dissociation. Then, in this study, we propose novel multi‐scale divide‐and‐conquer MD (MSDC‐MD), which combines AAM‐MD and coarse‐grained model MD (CGM‐MD). To speed up the conformation sampling, MSDC‐MD simulation first treats the IDR by CGM to sample conformations from wide conformation space; then, multiple AAM‐MD in a limited area is initiated using the resultant CGM conformation, which is reconstructed by homology modeling method. To investigate performance, we sampled the ordered conformation of the IDR using MSDC‐MD; the root‐mean‐square distance (RMSD) with respect to the experimental structure was 2.23 Å.     

Computational Characterization of Binding of Small Molecule Inhibitors to HIV‐1 gp41

Developing orally available small molecule inhibitors of HIV‐1 fusion has attracted significant interest over many years. Frey had recently reported several synthetic compounds which are experimentally shown to inhibit cell‐cell fusion in the low micromolar range. We carried out computational study to help identify possible binding modes by docking these compounds onto the hydrophobic pocket on gp41 and to characterize structures of binding complexes. The detailed gp41‐molecule binding interactions and free energies of binding are obtained through molecular dynamics simulation and MM‐PBSA calculation. Specific molecular interactions in the gp41‐inhibitor complexes are identified. The present computational study complements the corresponding experimental investigation and helps establish a good starting point for further refinement of small molecular gp41 inhibitors.     

Net charge changes in the calculation of relative ligand‐binding free energies via classical atomistic molecular dynamics simulation

The calculation of binding free energies of charged species to a target molecule is a frequently encountered problem in molecular dynamics studies of (bio‐)chemical thermodynamics. Many important endogenous receptor‐binding molecules, enzyme substrates, or drug molecules have a nonzero net charge. Absolute binding free energies, as well as binding free energies relative to another molecule with a different net charge will be affected by artifacts due to the used effective electrostatic interaction function and associated parameters (e.g., size of the computational box). In the present study, charging contributions to binding free energies of small oligoatomic ions to a series of model host cavities functionalized with different chemical groups are calculated with classical atomistic molecular dynamics simulation. Electrostatic interactions are treated using a lattice‐summation scheme or a cutoff‐truncation scheme with Barker–Watts reaction‐field correction, and the simulations are conducted in boxes of different edge lengths. It is illustrated that the charging free energies of the guest molecules in water and in the host strongly depend on the applied methodology and that neglect of correction terms for the artifacts introduced by the finite size of the simulated system and the use of an effective electrostatic interaction function considerably impairs the thermodynamic interpretation of guest‐host interactions. Application of correction terms for the various artifacts yields consistent results for the charging contribution to binding free energies and is thus a prerequisite for the valid interpretation or prediction of experimental data via molecular dynamics simulation. Analysis and correction of electrostatic artifacts according to the scheme proposed in the present study should therefore be considered an integral part of careful free‐energy calculation studies if changes in the net charge are involved. © 2013 The Authors Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

Direct Detection of DNA and DNA‐Ligand Interaction by Impedance Spectroscopy

In this work we present an impedimetric detection system for DNA‐ligand interactions. The sensor system consists of thiol‐modified single‐stranded DNA chemisorbed to gold. Impedance measurements in the presence of the redox system ferri‐/ferrocyanide show an increase in charge transfer resistance (R_ct) after hybridisation of a complementary target. Different amounts of capture strands, used for gold electrode modification, result in surface coverages between 3 and 15 pmol/cm² ssDNA. The relative change in R_ct upon hybridisation increases with increasing amount of capture probe on the electrode from 1.5‐ to 4.5‐fold. Impedimetric detection of binding events of a metal‐intercalator ([Ru(phen)₃]²⁺) and a groove binder (spermine) to double‐stranded DNA is demonstrated. Binding of [Ru(phen)₃]²⁺ and spermine exhibits a decrease in charge transfer resistance. Here, the ligand’s interaction leads to electrostatic shielding of the negatively charged DNA backbone. The impedance changes have been evaluated in dependence on the concentration of both DNA binders. Furthermore, the association of a single‐stranded binding protein (SSBP) is found to cause an increase in charge transfer resistance only when incubated with single‐stranded DNA. The specific binding of an anti‐dsDNA antibody to the dsDNA‐modified electrode surface decreases in contrast the interfacial impedance.     

Analyzing the robustness of the MM/PBSA free energy calculation method: application to DNA conformational transitions

The ability to predict and characterize free energy differences associated with conformational equilibria or the binding of biomolecules is vital to understanding the molecular basis of many important biological functions. As biological studies focus on larger molecular complexes and properties of the genome, proteome, and interactome, the development and characterization of efficient methods for calculating free energy becomes increasingly essential. The aim of this study is to examine the robustness of the end-point free energy method termed the molecular mechanics Poisson-Boltzmann solvent accessible surface area (MM/PBSA) method. Specifically, applications of MM/PBSA to the conformational equilibria of nucleic acid (NA) systems are explored. This is achieved by comparing A to B form DNA conformational free energy differences calculated using MM/PBSA with corresponding free energy differences determined with a more rigorous and time-consuming umbrella sampling algorithm. In addition, the robustness of NA MM/PBSA calculations is also evaluated in terms of the sensitivity towards the choice of force field and the choice of solvent model used during conformational sampling. MM/PBSA calculations of the free energy difference between A-form and B-form DNA are shown to be in very close agreement with the PMF result determined using an umbrella sampling approach. Further, it is found that the MM/PBSA conformational free energy differences were also in agreement using either the CHARMM or AMBER force field. The influence of ionic strength on conformational stability was particularly insensitive to the choice of force field. Finally, it is also shown that the use of a generalized Born implicit solvent during conformational sampling results in free energy estimates that deviate slightly from those obtained using explicitly solvated MD simulations in these NA systems.     

Computation of the binding affinities of catechol‐O‐methyltransferase inhibitors: Multisubstate relative free energy calculations

Alchemical free energy simulations are amongst the most accurate techniques for the computation of the free energy changes associated with noncovalent protein–ligand interactions. A procedure is presented to estimate the relative binding free energies of several ligands to the same protein target where multiple, low‐energy configurational substates might coexist, as opposed to one unique structure. The contributions of all individual substates were estimated, explicitly, with the free energy perturbation method, and combined in a rigorous fashion to compute the overall relative binding free energies and dissociation constants. It is shown that, unless the most stable bound forms are known a priori, inaccurate results may be obtained if the contributions of multiple substates are ignored. The method was applied to study the complex formed between human catechol‐O‐methyltransferase and BIA 9‐1067, a newly developed tight‐binding inhibitor that is currently under clinical evaluation for the therapy of Parkinson's disease. Our results reveal an exceptionally high‐binding affinity (K_d in subpicomolar range) and provide insightful clues on the interactions and mechanism of inhibition. The inhibitor is, itself, a slowly reacting substrate of the target enzyme and is released from the complex in the form of O‐methylated product. By comparing the experimental catalytic rate (k_cat) and the estimated dissociation rate (k_off) constants of the enzyme‐inhibitor complex, one can conclude that the observed inhibition potency (K_i) is primarily dependent on the catalytic rate constant of the inhibitor's O‐methylation, rather than the rate constant of dissociation of the complex. © 2012 Wiley Periodicals, Inc.     

Molecular dynamics simulations and free energy calculations on the enzyme 4‐hydroxyphenylpyruvate dioxygenase

4‐Hydroxyphenylpyruvate dioxygenase is a relevant target in both pharmaceutical and agricultural research. We report on molecular dynamics simulations and free energy calculations on this enzyme, in complex with 12 inhibitors for which experimental affinities were determined. We applied the thermodynamic integration approach and the more efficient one‐step perturbation. Even though simulations seem well converged and both methods show excellent agreement between them, the correlation with the experimental values remains poor. We investigate the effect of slight modifications on the charge distribution of these highly conjugated systems and find that accurate models can be obtained when using improved force field parameters. This study gives insight into the applicability of free energy methods and current limitations in force field parameterization. © 2011 Wiley Periodicals, Inc. J Comput Chem 2011     

Binding modes identification through molecular dynamic simulations: A case study with carnosine enantiomers and the Teicoplanin A2‐2‐based chiral stationary phase

In the present study, an in silico methodology able to define the binding modes adopted by carnosine enantiomers in the setting of the chiral recognition process is described. The inter‐ and intramolecular forces involved in the enantioseparation process with the Teicoplanin A2‐2 chiral selector and carnosine as model compound are successfully identified. This approach fully rationalizes, at a molecular level, the (S) < (R) enantiomeric elution order obtained under reversed‐phase conditions. Consistent explanations were achieved by managing molecular dynamics results with advanced techniques of data analysis. As a result, the time‐dependent identification of all the interactions simultaneously occurring in the chiral selector‐enantiomeric analyte binding process was obtained. Accordingly, it was found that only (R)‐carnosine is able to engage a stabilizing charge–charge interaction through its ionized imidazole ring with the carboxylate counter‐part on the chiral selector. Instead, (S)‐carnosine establishes intramolecular contacts between its ionized functional groups, that limit its conformational freedom and impair the association with the chiral selector unit.     

DNA Photosensitization by an “Insider”: Photophysics and Triplet Energy Transfer of 5‐Methyl‐2‐pyrimidone Deoxyribonucleoside

The main chromophore of (6‐4) photoproducts, namely, 5‐methyl‐2‐pyrimidone (Pyo), is an artificial noncanonical nucleobase. This chromophore has recently been reported as a potential photosensitizer that induces triplet damage in thymine DNA. In this study, we investigate the spectroscopic properties of the Pyo unit embedded in DNA by means of explicit solvent molecular‐dynamics simulations coupled to time‐dependent DFT and quantum‐mechanics/molecular‐mechanics techniques. Triplet‐state transfer from the Pyo to the thymine unit was monitored in B‐DNA by probing the propensity of this photoactive pyrimidine analogue to induce a Dexter‐type triplet photosensitization and subsequent DNA damage.     

Molecular dynamics study of DNA binding by INT‐DBD under a polarized force field

The DNA binding domain of transposon Tn916 integrase (INT‐DBD) binds to DNA target site by positioning the face of a three‐stranded antiparallel β‐sheet within the major groove. As the negatively charged DNA directly interacts with the positively charged residues (such as Arg and Lys) of INT‐DBD, the electrostatic interaction is expected to play an important role in the dynamical stability of the protein–DNA binding complex. In the current work, the combined use of quantum‐based polarized protein‐specific charge (PPC) for protein and polarized nucleic acid‐specific charge (PNC) for DNA were employed in molecular dynamics simulation to study the interaction dynamics between INT‐DBD and DNA. Our study shows that the protein–DNA structure is stabilized by polarization and the calculated protein–DNA binding free energy is in good agreement with the experimental data. Furthermore, our study revealed a positive correlation between the measured binding energy difference in alanine mutation and the occupancy of the corresponding residue's hydrogen bond. This correlation relation directly relates the contribution of a specific residue to protein–DNA binding energy to the strength of the hydrogen bond formed between the specific residue and DNA. © 2013 Wiley Periodicals, Inc.     

Interaction of Metal Complexes with G‐Quadruplex DNA

Guanine‐rich sequences of DNA can assemble into tetrastranded structures known as G‐quadruplexes. It has been suggested that these secondary DNA structures could be involved in the regulation of several key biological processes. In the human genome, guanine‐rich sequences with the potential to form G‐quadruplexes exist in the telomere as well as in promoter regions of certain oncogenes. The identification of these sequences as novel targets for the development of anticancer drugs has sparked great interest in the design of molecules that can interact with quadruplex DNA. While most reported quadruplex DNA binders are based on purely organic templates, numerous metal complexes have more recently been shown to interact effectively with this DNA secondary structure. This Review provides an overview of the important roles that metal complexes can play as quadruplex DNA binding molecules, highlighting the unique properties metals can confer to these molecules.     

On the induced‐fit mechanism of substrate‐enzyme binding structures of nylon‐oligomer hydrolase

We present a detailed computational investigation of the induced‐fit motion in a nylon‐oligomer hydrolase (NylB) upon substrate binding. To this aim, we resort on the recently introduced parallel cascade selection molecular dynamics approach, allowing for an accelerated access to the set of conformational changes from an open‐ to a closed‐state structure to form the enzyme‐substrate complex in a specific induce‐fit mechanism. The structural investigation is quantitatively complemented by free energy analyses within the umbrella sampling algorithm accompanied by weighted histogram analysis. We find that the stabilization free energy is about 1.4 kcal/mol, whereas the highest free energy barrier to be overcome is about 2.3 kcal/mol. Conversely, the energetic contribution for the substrate binding is about 20 kcal/mol, as estimated from Generalized Born/Surface Area. This means that the open‐close induced‐fit motion could occur frequently once the substrate binds to the open state of NylB. © 2014 Wiley Periodicals, Inc.