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

Structure based computational RNA design towards MafA transcriptional repressor implicated in multiple myeloma

Multiple myeloma is recognized as the second most common hematological cancer. MafA transcriptional repressor is an established mediator of myelomagenesis. While there are multitude of drugs available for targeting various effectors in multiple myeloma, current literature lacks a candidate RNA based MafA modulator. Thus, using the structure of MafA homodimer-consensus target DNA, a computational effort was implemented to design a novel RNA based chemical modulator against MafA. First, available MafA-consensus DNA structure was employed to generate an RNA library. This library was further subjected to global docking to select the most plausible RNA candidates, preferring to bind DNA binding region of MafA. Following global docking, MD-ready complexes that were prepared via local docking program, were subjected to 500 ns of MD simulations. First, each of these MD simulations were analyzed for relative binding free energy through MM-PBSA method, which pointed towards a strong RNA based MafA binder, RNA1. Second, through a detailed MD analysis, RNA1 was shown to prefer binding to a single monomer of the dimeric DNA binding domain of MafA using higher number of hydrophobic interactions compared with positive control MafA-DNA complex. At the final phase, a principal component analyses was conducted, which led us to identify the actual interaction region of RNA1 and MafA monomer. Overall, to our knowledge, this is the first computational study that presents an RNA molecule capable of potentially targeting MafA protein. Furthermore, limitations of our study together with possible future implications of RNA1 in multiple myeloma were also discussed.

Evaluation of surfactant-aided polycyclic aromatic hydrocarbon biodegradation by molecular docking and molecular dynamic simulation in the marine environment

Marine oil spills directly cause polycyclic aromatic hydrocarbons (PAHs) pollution and affect marine organisms due to their toxic property. Chemical and bio-based dispersants composed of surfactants and solvents are considered effective oil spill-treating agents. Dispersants enhance oil biodegradation in the marine environment by rapidly increasing their solubility in the water column. However, the effect of dispersants, especially surfactants, on PAHs degradation by enzymes produced by microorganisms has not been studied at the molecular level. The role of the cytochrome P450 (CYP) enzyme in converting contaminants into reactive metabolites during the biodegradation process has been evidenced, but the activity in the presence of surfactants is still ambiguous. Thus, this study focused on the evaluation of the impact of chemical and bio-surfactants (i.e., Tween 80 (TWE) and Surfactin (SUC)) on the biodegradation of naphthalene (NAP), chrysene (CHR), and pyrene (PYR), the representative components of PAHs, with CYP enzyme from microalgae Parachlorella kessleri using molecular docking and molecular dynamics (MD) simulation. The molecular docking analysis revealed that PAHs bound to residues at the CYP active site through hydrophobic interactions for biodegradation. The MD simulation showed that the surfactant addition changed the enzyme conformation in the CYP-PAH complexes to provide more interactions between the enzyme and PAHs. This led to an increase in the enzyme's capability to degrade PAHs. Binding free energy (ΔG‌‌Bind) calculations confirmed that surfactant treatment could enhance PAHs degradation by the enzyme. The SUC gave a better result on NAP and PYR biodegradation based on ΔG‌‌Bind, while TWE facilitated the biodegradation of CHR. The research outputs could greatly facilitate evaluating the behaviors of oil spill-treating agents and oil spill response operations in the marine environment.

In silico studies of the interaction of the minor groove binder Hoechst 33258 with B-DNA

Interaction of the minor groove binder, Hoechst 33258, with the Dickerson-Drew DNA dodecamer sequence has been investigated using docking, MM/QM, MM/GBSA and molecular dynamics computations to study the modes of binding and the interactions responsible for the binding. Besides the original Hoechst 33258 ligand (HT), a total of 12 ionization and stereochemical states for the ligand are obtained at the physiological pH and have been docked into B-DNA. These states have one or the other or both benzimidazole rings in protonated states, apart from the piperazine nitrogen, which has a quaternary nitrogen in all the states. Most of these states are found to exhibit good docking scores and free energy of binding with B-DNA. The best docked state has been taken further for molecular dynamics simulations and compared with the original HT. This state is protonated at both benzimidazole rings besides the piperazine ring and hence has very highly negative coulombic interaction energy. In both cases, there are strong coulombic interactions, but these are offset by the almost equally unfavorable solvation energies. Thus, the nonpolar forces, particularly van der Waals contacts, dominate the interaction, and the polar interactions highlight subtle changes in the binding energies, leading to more highly protonated states having more negative binding energies.Communicated by Ramaswamy H. Sarma.

Unusual similarity of DNA solvation dynamics in high-salinity crowding with divalent cations of varying concentrations

Double-stranded DNA bears the highest linear negative charge density (2e- per base-pair) among all biopolymers, leading to strong interactions with cations and dipolar water, resulting in the formation of a dense 'condensation layer' around DNA. Interactions involving proteins and ligands binding to DNA are primarily governed by strong electrostatic forces. Increased salt concentrations impede such electrostatic interactions - a situation that prevails in oceanic species due to their cytoplasm being enriched with salts. Nevertheless, how these interactions' dynamics are affected in crowded hypersaline environments remains largely unexplored. Here, we employ steady-state and time-resolved fluorescence Stokes shifts (TRFSS) of a DNA-bound ligand (DAPI) to investigate the static and dynamic solvation properties of DNA in the presence of two divalent cations, magnesium (Mg2+), and calcium (Ca2+) at varying high to very-high concentrations of 0.15 M, 1 M and 2 M. We compare the results to those obtained in physiological concentrations (0.15 M) of monovalent Na+ ions. Combining data from fluorescence femtosecond optical gating (FOG) and time-correlated single photon counting (TCSPC) techniques, dynamic fluorescence Stokes shifts in DNA are analysed over a broad range of time-scales, from 100 fs to 10 ns. We find that while divalent cation crowding strongly influences the DNA stability and ligand binding affinity to DNA, the dynamics of DNA solvation remain remarkably similar across a broad range of five decades in time, even in a high-salinity crowded environment with divalent cations, as compared to the physiological concentration of the Na+ ion. Steady-state and time-resolved data of the DNA-groove-bound ligand are seemingly unaffected by ion-crowding in hypersaline solution, possibly due to ions being mostly displaced by the DNA-bound ligand. Furthermore, the dynamic coupling of cations with nearby water may possibly contribute to a net-neutral effect on the overall collective solvation dynamics in DNA, owing to the strong anti-correlation of their electrostatic interaction energy fluctuations. Such dynamic scenarios may persist within the cellular environment of marine life and other biological cells that experience hypersaline conditions.

In vitro and in silico binding studies of phytochemical isochroman with calf thymus DNA using multi-spectroscopic and computational modelling techniques

A wide range of therapeutic molecules uses deoxyribonucleic acid (DNA) as an intracellular target. The interaction of small molecules to DNA is a key feature in pharmacology and plays a vital role in the development of novel and more efficient drugs with increased selective activity and enhanced therapeutic effectiveness. Isochroman (IC) is a constituent of Olea europea plant, which has been shown to exhibit several beneficial pharmacological activities. At present, its interaction studies using calf thymus DNA (ct-DNA) have not been explained. A set of multi-spectroscopic techniques has been performed to determine the interaction mechanism of isochroman with ct-DNA. Absorption spectra and quenching in fluorescence studies show that isochroman and ct-DNA form a complex. The static mode of quenching was determined by the Stern-Volmer plot. The value of binding constant, Kb = 4.0 × 103 M-1 revealed moderate type of binding. Effects of single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA) and ionic strength were studied to examine the isochroman binding to ct-DNA. Potassium iodide (KI) quenching effects and competitive binding studies clearly showed that isochroman binds in the minor groove of ct-DNA. Circular dichroic and DNA melting experiments also confirmed these results. The experimental outputs were further corroborated via in silico computational modelling studies. Lipinski's rule of 5 and SwissADME showed drug-likeness and oral bioavailability scores. Protox ІІ online software predicts oral and organ toxicity.Communicated by Ramaswamy H. Sarma.

X-ray Structure Characterization of the Selective Recognition of AT Base Pair Sequences

The rational design of small molecules that target specific DNA sequences is a promising strategy to modulate gene expression. This report focuses on a diamidinobenzimidazole compound, whose selective binding to the minor groove of AT DNA sequences holds broad significance in the molecular recognition of AT-rich human promoter sequences. The objective of this study is to provide a more detailed and systematized understanding, at an atomic level, of the molecular recognition mechanism of different AT-specific sequences by a rationally designed minor groove binder. The specialized method of X-ray crystallography was utilized to investigate how the sequence-dependent recognition properties in general, A-tract, and alternating AT sequences affect the binding of diamidinobenzimidazole in the DNA minor groove. While general and A-tract AT sequences give a narrower minor groove, the alternating AT sequences intrinsically have a wider minor groove which typically constricts upon binding. A strong and direct hydrogen bond between the N-H of the benzimidazole and an H-bond acceptor atom in the minor groove is essential for DNA recognition in all sequences described. In addition, the diamidine compound specifically utilizes an interfacial water molecule for its DNA binding. DNA complexes of AATT and AAAAAA recognition sites show that the diamidine compound can bind in two possible orientations with a preference for water-assisted hydrogen bonding at either cationic end. The complex structures of AAATTT, ATAT, ATATAT, and AAAA are bound in a singular orientation. Analysis of the helical parameters shows a minor groove expansion of about 1 Å across all the nonalternating DNA complexes. The results from this systematic approach will convey a greater understanding of the specific recognition of a diverse array of AT-rich sequences by small molecules and more insight into the design of small molecules with enhanced specificity to AT and mixed DNA sequences.

Assessment of Different Parameters on the Accuracy of Computational Alanine Scanning of Protein-Protein Complexes with the Molecular Mechanics/Generalized Born Surface Area Method

Computational alanine scanning with the molecular mechanics generalized Born surface area (MM/GBSA) method constitutes a widely used approach for identifying critical residues at protein-protein interfaces. Despite its popularity, the MM/GBSA method still has certain drawbacks due to its dependence on many factors. Here, we performed a systematical study on the impact of four different parameters, namely, the internal dielectric constant, the generalized Born model, the entropic term, and the inclusion of structural waters on the accuracy of computational alanine scanning calculations with the MM/GBSA method. Our results show that the internal dielectric constant is the most critical parameter for getting accurate predictions. The introduction of entropy and interfacial water molecules decreased the quality of the predictions, while the generalized Born model had little to no effect. Considering the significance of the internal dielectric value, we proposed a methodology based on the energetic predominance of a particular set of amino acids at the protein-protein interface for selecting an appropriate value for this variable. We hope that these results serve as a guideline for future studies of protein-protein complexes using the MM/GBSA method.

Encapsulation Dynamics of Neuromuscular Blocking Drugs by Sugammadex

BACKGROUND

The clinical actions of sugammadex have been well studied, but the detailed molecular mechanism of the drug encapsulation process has not been systematically documented. The hypothesis was that sugammadex would attract rocuronium and vecuronium via interaction with the sugammadex side-chain "tentacles," as previously suggested.

METHODS

Computational molecular dynamics simulations were done to investigate docking of sugammadex with rocuronium and vecuronium. To validate these methods, strength of binding was assessed between sugammadex and a heterogeneous group of nine other drugs, the binding affinities of which have been experimentally determined. These observations hinted that high concentrations of unbound sugammadex could bind to propofol, potentially altering its pharmacokinetic profile. This was tested experimentally in in vitro cortical slices.

RESULTS

Sugammadex encapsulation of rocuronium involved a sequential progression down a series of metastable states. After initially binding beside the sugammadex molecule (mean ± SD center-of-mass distance = 1.17 ± 0.13 nm), rocuronium then moved to the opposite side to that hypothesized, where it optimally aligned with the 16 hydroxyl groups (distance, 0.82 ± 0.04 nm) before entering the sugammadex cavity to achieve energetically stable encapsulation by approximately 120 ns (distance, 0.35 ± 0.12 nm). Vecuronium formed fewer hydrogen bonds with sugammadex than did rocuronium; hence, it was less avidly bound. For the other molecules, the computational results showed good agreement with the available experimental data, showing a clear bilogarithmic relation between the relative binding free energy and the association constant (R2 = 0.98). Weaker binding was manifest by periodic unbinding. The brain slice results confirmed the presence of a weak propofol-sugammadex interaction.

CONCLUSIONS

Computational simulations demonstrate the dynamics of neuromuscular blocking drug encapsulation by sugammadex occurring from the opposite direction to that hypothesized and also how high concentrations of unbound sugammadex can potentially weakly bind to other drugs given during general anesthesia.

EDITOR’S PERSPECTIVE

Engineered modular heterocyclic-diamidines for sequence-specific recognition of mixed AT/GC base pairs at the DNA minor groove

This report describes a breakthrough in a project to design minor groove binders to recognize any sequence of DNA. A key goal is to invent synthetic chemistry for compound preparation to recognize an adjacent GG sequence that has been difficult to target. After trying several unsuccessful compound designs, an N-alkyl-benzodiimidazole structure was selected to provide two H-bond acceptors for the adjacent GG-NH groups. Flanking thiophenes provide a preorganized structure with strong affinity, DB2831, and the structure is terminated by phenyl-amidines. The binding experimental results for DB2831 with a target AAAGGTTT sequence were successful and include a high ΔT m, biosensor SPR with a K D of 4 nM, a similar K D from fluorescence titrations and supporting competition mass spectrometry. MD analysis of DB2831 bound to an AAAGGTTT site reveals that the two unprotonated N of the benzodiimidazole group form strong H-bonds (based on distance) with the two central G-NH while the central -CH of the benzodiimidazole is close to the -C[double bond, length as m-dash]O of a C base. These three interactions account for the strong preference of DB2831 for a -GG- sequence. Surprisingly, a complex with one dynamic, interfacial water is favored with 75% occupancy.

gmx_MMPBSA: A New Tool to Perform End-State Free Energy Calculations with GROMACS

Molecular mechanics/Poisson-Boltzmann (Generalized-Born) surface area is one of the most popular methods to estimate binding free energies. This method has been proven to balance accuracy and computational efficiency, especially when dealing with large systems. As a result of its popularity, several programs have been developed for performing MM/PB(GB)SA calculations within the GROMACS community. These programs, however, present several limitations. Here we present gmx_MMPBSA, a new tool to perform end-state free energy calculations from GROMACS molecular dynamics trajectories. gmx_MMPBSA provides the user with several options, including binding free energy calculations with different solvation models (PB, GB, or 3D-RISM), stability calculations, computational alanine scanning, entropy corrections, and binding free energy decomposition. Noteworthy, several promising methodologies to calculate relative binding free energies such as alanine scanning with variable dielectric constant and interaction entropy have also been implemented in gmx_MMPBSA. Two additional tools-gmx_MMPBSA_test and gmx_MMPBSA_ana-have been integrated within gmx_MMPBSA to improve its usability. Multiple illustrating examples can be accessed through gmx_MMPBSA_test, while gmx_MMPBSA_ana provides fast, easy, and efficient access to different graphics plotted from gmx_MMPBSA output files. The latest version (v1.4.3, 26/05/2021) is available free of charge (documentation, test files, and tutorials included) at https://github.com/Valdes-Tresanco-MS/gmx_MMPBSA.

Evaluation of binding of potential ADMET/tox screened saquinavir analogues for inhibition of HIV-protease via molecular dynamics and binding free energy calculations

Developing novel drug molecules against HIV is a scientific quest necessitated by development of drug resistance against used drugs. We report comparative results of molecular dynamics simulation studies on 11 structural analogues of Saquinavir (SQV) - against HIV-protease that were earlier examined for pharmacodynamic and pharmacokinetic properties. We reported analogues S1, S5 and S8 to qualify the ADMET criterion and may be considered as potential lead molecules. In this study the designed molecules were successively docked with native HIV-protease at AutoDock. Docking scores established relative goodness of the 11 analogues against the benchmark for Saquinavir. The docked complexes were subjected to molecular dynamics simulation studies using GROMACS 4.6.2. Four parameters viz. H-bonding, RMSD, Binding energy and Protein-Ligand Distance were used for comparative analyses of the analogues relative to Saquinavir. The comparison and analysis of the results are indicative that analogues S8, S9 and S1 are promising candidates among all the analogues studied. From our earlier work and present study it is evident that among the three S8 and S1 qualify the ADMET criterion and between S1 and S8, the analogue S8 shows more target efficacy and specificity over S1 and have better molecular dynamics simulation results. Thus, of the 11 de novo Saquinavir analogues, the S8 appears to be the most promising candidate as lead molecule for HIV-protease inhibitor and is best suited for testing under biological system. Further validation of the proposed lead molecules through wet lab studies involving antiviral assays however is required.Communicated by Ramaswamy H. Sarma.

Gold-Etched Silver Nanowire Endoscopy: Toward a Widely Accessible Platform for Surface-Enhanced Raman Scattering-Based Analysis in Living Cells

Recently, our group introduced the use of silver nanowires (AgNWs) as novel non-invasive endoscopic probes for detecting intracellular Raman signals. This method, although innovative and promising, relies exclusively on the plasmonic waveguiding effect for signal enhancement. It, therefore, requires sophisticated operational tools and protocols, drastically limiting its applicability. Herein, an advanced strategy is offered to significantly enhance the performance of these endoscopic probes, making this approach widely accessible and versatile for cellular studies. By uniformly forming gold structures on the smooth AgNW surface via a galvanic replacement reaction, the density of the light coupling points along the whole probe surface is drastically increased, enabling high surface-enhanced Raman scattering (SERS) efficiency upon solely focusing the excitation light on the gold-etched AgNW. The applicability of these gold-etched AgNW probes for molecular sensing in cells is demonstrated by detecting site-specific and high-resolved SERS spectra of cell compartment-labeling dyes, namely, 4',6-diamidino-2-phenylindole in the nucleus and 3,3'-dioctadecyloxacarbocyanine on the membrane. The remarkable spectral sensitivity achieved provides essential structural information of the analytes, indicating the overall potential of the proposed approach for cellular studies of drug interactions with biomolecular items.

Tracing the driving forces responsible for the remarkable infectivity of 2019-nCoV: 1. Receptor binding domain in its bound and unbound states

The pandemic of COVID-19 has posed an urgent need to learn the dynamics of the virus and the mechanism of its contagion of host cells. By means of molecular dynamics simulations, this work addressed the behavior of 2019-nCoV in two aspects: the binding affinity of its receptor binding domain (RBD) with ACE2, and its potential conformation preferences in its unbound state. The results showed that the RBD of 2019-nCov bound much stronger with ACE2 than that of SARS-CoV due to a better organized hydrogen bond network between the former pair with most of the residues at the contact interface sharing the responsibility to hold the pair tightly. This is in contrast to the case of SARS-CoV, which strongly relied on the residues at the ends of the cleft. In its unbound state, the RBD of 2019-nCoV was found to fold part of its receptor binding motif (RBM) into a helical conformation and flip into a concave to minimize its contact with the external environment. This has the biological implication that the virus may achieve higher translational motion in the condensed phase and have a higher chance of survival by avoiding capture by the immune system before reaching its target receptor.

Structure based peptide design, molecular dynamics and MM-PBSA studies for targeting C terminal dimerization of NFAT5 DNA binding domain

NFAT5 as a transcription factor with an established role in osmotic stress response, has also been revealed to be active under numerous settings, including pathological conditions such as diabetic microvascular complications, chronic arthritis and cancer. Despite these links, current strategies for downregulating NFAT5 activity only relies on indirect modulators, not directly targeting NFAT5, itself. With this study, through using a computational approach, an original peptide was explored to directly target C terminal dimerization of NFAT5 RHR, located in its DNA binding domain. At first, homodimeric NFAT5 RHR bound to its consensus DNA was used for prediction of a preliminary peptide sequence. Possible amino acid replacements for this preliminary peptide were predicted for optimization, which was followed by addition of a cell penetrating peptide sequence. These attempts yielded a small peptide library, which was further investigated for peptide affinities towards C terminal of NFAT5 RHR through molecular docking, 50 ns and 250 ns molecular dynamics simulations, followed by estimation of MM-PBSA based relative binding free energies. Results indicated that after receiving mutations on the preliminary peptide sequence for optimization, a unique peptide could target C terminal dimerization region of NFAT5 RHR through using its cell penetrating peptide sequence. In conclusion, this is the first study presenting computational evidence on identification of a novel peptide capable of directly targeting NFAT5 dimerization. Besides, future implications of these observations were also discussed in terms of methodology and possible applications.

Extending the σ-Hole Motif for Sequence-Specific Recognition of the DNA Minor Groove

The majority of current drugs against diseases, such as cancer, can bind to one or more sites in a protein and inhibit its activity. There are, however, well-known limits on the number of druggable proteins, and complementary current drugs with compounds that could selectively target DNA or RNA would greatly enhance the availability of cellular probes and therapeutic progress. We are focusing on the design of sequence-specific DNA minor groove binders that, for example, target the promoter sites of transcription factors involved in a disease. We have started with AT-specific minor groove binders that are known to enter human cells and have entered clinical trials. To broaden the sequence-specific recognition of these compounds, several modules that have H-bond acceptors that strongly and specifically recognize G·C base pairs were identified. A lead module is a thiophene-N-alkyl-benzimidazole σ-hole-based system with terminal phenyl-amidines that have excellent affinity and selectivity for a G·C base pair in the minor groove. Efforts are now focused on optimizing this module. In this work, we are evaluating modifications to the compound aromatic system with the goal of improving GC selectivity and affinity. The lead compounds retain the thiophene-N-alkyl-BI module but have halogen substituents adjacent to an amidine group on the terminal phenyl-amidine. The optimum compounds must have strong affinity and specificity with a residence time of at least 100 s.

Small Sequence-Sensitive Compounds for Specific Recognition of the G⋅C Base Pair in DNA Minor Groove

A series of small diamidines with thiophene and modified N-alkylbenzimidazole σ-hole module represent specific binding to single G⋅C base pair (bp) DNA sequence. The variation of N-alkyl or aromatic rings were sensitive to microstructures of the DNA minor groove. Thirteen new compounds were synthesized to test their binding affinity and selectivity. The dicyanobenzimidazoles needed to synthesize the target diamidines were made via condensation/cyclization reactions of different aldehydes with different 3-amino-4-(alkyl- or phenyl-amino) benzonitriles. The final diamidines were synthesized using lithium bis-trimethylsilylamide (LiN[Si(CH3 )3 ]2 ) or Pinner methods. The newly synthesized compounds showed strong binding and selectivity to AAAGTTT compared to similar sequences AAATTT and AAAGCTTT investigated by several biophysical methods including biosensor-SPR, fluorescence spectroscopy, DNA thermal melting, ESI-MS spectrometry, circular dichroism, and molecular dynamics. The binding affinity results determined by fluorescence spectroscopy are in accordance with those obtained by biosensor-SPR. These small size single G⋅C bp highly specific binders extend the compound database for future biological applications.

Investigation of thermal properties of DNA structure with precise atomic arrangement via equilibrium and non-equilibrium molecular dynamics approaches

BACKGROUND AND OBJECTIVE

Thermal conductivity of Deoxyribonucleic acid molecules is important for nanotechnology applications. Theoretical simulations based on simple models predict thermal conductivity for these molecular structures.

METHODS

In this work, we calculate the thermal properties of Deoxyribonucleic acid with precise atomic arrangement via equilibrium and non-equilibrium molecular dynamics approaches. In these methods, each Deoxyribonucleic acid molecule is represented by C, N, O, and P atoms and implemented dreidng potential to describe their atomic interactions.

RESULTS

Our calculated rate for thermal conductivity via equilibrium and non-equilibrium molecular dynamics methods is 0.381 W/m K and 0.373 W/m K, respectively. By comparing results from these two methods, it was found that the results from equilibrium and non-equilibrium molecular dynamics methods are identical, approximately. On the other hand, the number of DNA molecules and the equilibrium temperature of the simulated structures were important factors in their thermal conductivity rates, and their thermal conductivity was calculated at 0.323 W/m K-0.381 W/m K intervals for equilibrium and 0.303 W/m K-0.373 W/m K interval for non-equilibrium calculations.

CONCLUSIONS

These results are in good agreement with thermal conductivity calculation with other research groups.

Origin of Slow Solvation Dynamics in DNA: DAPI in Minor Groove of Dickerson-Drew DNA

The measurement and understanding of collective solvation dynamics in DNA have vital biological implications, as protein and ligand binding to DNA can be directly controlled by complex electrostatic interactions of anionic DNA and surrounding dipolar water, and ions. Time-resolved fluorescence Stokes shift (TRFSS) experiments revealed anomalously slow solvation dynamics in DNA much beyond 100 ps that follow either power-law or slow multiexponential decay over several nanoseconds. The origin of such dispersed dynamics remains difficult to understand. Here we compare results of TRFSS experiments to molecular dynamics (MD) simulations of well-known 4',6-diamidino-2-phenylindole (DAPI)/Dickerson-Drew DNA complex over five decades of time from 100 fs to 10 ns to understand the origin of such dispersed dynamics. We show that the solvation time-correlation function (TCF) calculated from 200 ns simulation trajectory (total 800 ns) captures most features of slow dynamics as measured in TRFSS experiments. Decomposition of TCF into individual components unravels that slow dynamics originating from dynamically coupled DNA-water motion, although contribution from coupled water-Na⁺ motion is non-negligible. The analysis of residence time of water molecules around the probe (DAPI) reveals broad distribution from ∼6 ps to ∼3.5 ns: Several (49 nos.) water molecules show residences time greater than 500 ps, of which at least 14 water molecules show residence times of more than 1 ns in the first solvation shell of DAPI. Most of these slow water molecules are found to occupy two hydration sites in the minor groove near DAPI binding site. The residence time of Na⁺, however, is found to vary within ∼17-120 ps. Remarkably, we find that freezing the DNA fluctuations in simulation eliminates slower dynamics beyond ∼100 ps, where water and Na⁺ dynamics become faster, although strong anticorrelation exists between them. These results indicate that primary origin of slow dynamics lies within the slow fluctuations of DNA parts that couple with nearby slow water and ions to control the dispersed collective solvation dynamics in DNA minor groove.

A Comparative Linear Interaction Energy and MM/PBSA Study on SIRT1-Ligand Binding Free Energy Calculation

Binding free energy (ΔGbind) computation can play an important role in prioritizing compounds to be evaluated experimentally on their affinity for target proteins, yet fast and accurate ΔGbind calculation remains an elusive task. In this study, we compare the performance of two popular end-point methods, i.e., linear interaction energy (LIE) and molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA), with respect to their ability to correlate calculated binding affinities of 27 thieno[3,2-d]pyrimidine-6-carboxamide-derived sirtuin 1 (SIRT1) inhibitors with experimental data. Compared with the standard single-trajectory setup of MM/PBSA, our study elucidates that LIE allows to obtain direct ("absolute") values for SIRT1 binding free energies with lower compute requirements, while the accuracy in calculating relative values for ΔGbind is comparable (Pearson's r = 0.72 and 0.64 for LIE and MM/PBSA, respectively). We also investigate the potential of combining multiple docking poses in iterative LIE models and find that Boltzmann-like weighting of outcomes of simulations starting from different poses can retrieve appropriate binding orientations. In addition, we find that in this particular case study the LIE and MM/PBSA models can be optimized by neglecting the contributions from electrostatic and polar interactions to the ΔGbind calculations.

Structural approaches for the DNA binding motifs prediction in Bacillus thuringiensis sigma-E transcription factor (σETF)

The sigma-E transcription factor (σ^ETF) can be found in most of the bacteria cells including Bacillus thuringiensis. However, the cellular regulatory mechanisms of these transcription factors in the mass production of δ-endotoxins during sporulation stage are yet to be revealed. In addition, the recognition of DNA towards σ^ETF DNA binding motifs that led to the transcription activities is also being poorly studied. Therefore, this work studied the possible DNA binding motifs of σ^ETF by utilising in silico approaches. The structure of σ^ETF was first built via three different computational methods. A cognate DNA sequence was then docked to the predicted σ^ETF DNA-binding motifs. The binding free energy calculated using molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) for triplicate 50 ns simulation of σ^ETF-DNA complex revealed favourable binding energy of DNA to σ^ETF (average ∆G_bind = -34.57 kcal/mol) mainly driven by non-polar interactions. This study revealed that σ^ETF LYS131, ARG133, PHE138, TRP146, ARG222, LYS225 and ARG226 are most likely the key residues upon the binding and recognition of DNA prior to transcription actives. Since determination of genome-regulating protein which recognises specific DNA sequence is important to discriminate between the proteins preferences for different genes, this study might provide some understanding on the possible σ^ETF-DNA recognition prior to transcription initiated for the δ-endotoxins production.

Stereochemistry and innate immune recognition: (+)-norbinaltorphimine targets myeloid differentiation protein 2 and inhibits toll-like receptor 4 signaling

Deregulation of innate immune TLR4 signaling contributes to various diseases including neuropathic pain and drug addiction. Naltrexone is one of the rare TLR4 antagonists with good blood-brain barrier permeability and showing no stereoselectivity for TLR4. By linking 2 naltrexone units through a rigid pyrrole spacer, the bivalent ligand norbinaltorphimine was formed. Interestingly, (+)-norbinaltorphimine [(+)-1] showed ∼25 times better TLR4 antagonist activity than naltrexone in microglial BV-2 cell line, whereas (-)-norbinaltorphimine [(-)-1] lost TLR4 activity. The enantioselectivity of norbinaltorphimine was further confirmed in primary microglia, astrocytes, and macrophages. The activities of meso isomer of norbinaltorphimine and the molecular dynamic simulation results demonstrate that the stereochemistry of (+)-1 is derived from the (+)-naltrexone pharmacophore. Moreover, (+)-1 significantly increased and prolonged morphine analgesia in vivo. The efficacy of (+)-1 is long lasting. This is the first report showing enantioselective modulation of the innate immune TLR signaling.-Zhang, X., Peng, Y., Grace, P. M., Metcalf, M. D., Kwilasz, A. J., Wang, Y., Zhang, T., Wu, S., Selfridge, B. R., Portoghese, P. S., Rice, K. C., Watkins, L. R., Hutchinson, M. R., Wang, X. Stereochemistry and innate immune recognition: (+)-norbinaltorphimine targets myeloid differentiation protein 2 and inhibits toll-like receptor 4 signaling.

Small molecule binders recognize DNA microstructural variations via an induced fit mechanism

Regulation of gene-expression by specific targeting of protein-nucleic acid interactions has been a long-standing goal in medicinal chemistry. Transcription factors are considered "undruggable" because they lack binding sites well suited for binding small-molecules. In order to overcome this obstacle, we are interested in designing small molecules that bind to the corresponding promoter sequences and either prevent or modulate transcription factor association via an allosteric mechanism. To achieve this, we must design small molecules that are both sequence-specific and able to target G/C base pair sites. A thorough understanding of the relationship between binding affinity and the structural aspects of the small molecule-DNA complex would greatly aid in rational design of such compounds. Here we present a comprehensive analysis of sequence-specific DNA association of a synthetic minor groove binder using long timescale molecular dynamics. We show how binding selectivity arises from a combination of structural factors. Our results provide a framework for the rational design and optimization of synthetic small molecules in order to improve site-specific targeting of DNA for therapeutic uses in the design of selective DNA binders targeting transcription regulation.

Phosphorylation promotes binding affinity of Rap-Raf complex by allosteric modulation of switch loop dynamics

The effects of phosphorylation of a serine residue on the structural and dynamic properties of Ras-like protein, Rap, and its interactions with effector protein Ras binding domain (RBD) of Raf kinase, in the presence of GTP, are investigated via molecular dynamics simulations. The simulations show that phosphorylation significantly effects the dynamics of functional loops of Rap which participate in the stability of the complex with effector proteins. The effects of phosphorylation on Rap are significant and detailed conformational analysis suggest that the Rap protein, when phosphorylated and with GTP ligand, samples different conformational space as compared to non-phosphorylated protein. In addition, phosphorylation of SER11 opens up a new cavity in the Rap protein which can be further explored for possible drug interactions. Residue network analysis shows that the phosphorylation of Rap results in a community spanning both Rap and RBD and strongly suggests transmission of allosteric effects of local alterations in Rap to distal regions of RBD, potentially affecting the downstream signalling. Binding free energy calculations suggest that phosphorylation of SER11 residue increases the binding between Rap and Raf corroborating the network analysis results. The increased binding of the Rap-Raf complex can have cascading effects along the signalling pathways where availability of Raf can influence the oncogenic effects of Ras proteins. These simulations underscore the importance of post translational modifications like phosphorylation on the functional dynamics in proteins and can be an alternative to drug-targeting, especially in notoriously undruggable oncoproteins belonging to Ras-like GTPase family.

Bound Compound, Interfacial Water, and Phenyl Ring Rotation Dynamics of a Compound in the DNA Minor Groove

DB2277, a heterocyclic diamidine, is a successful design for mixed base pair (bp) DNA sequence recognition. The compound has a central aza-benzimidazole group that forms two H-bonds with a GC bp that has flanking AT bps. The nuclear magnetic resonance structure of the DB2277-DNA complex with an AAGATA recognition site sequence was determined, and here we report extended molecular dynamics (MD) simulations of the structure. DB2277 has two terminal phenyl-amidine groups, one of which is directly linked to the DB2277 heterocyclic core and the other through a flexible -OCH₂- group. The flexibly linked phenyl is too far from the minor groove floor to make direct H-bonds but is linked to an AT bp through water-mediated H-bonds. The flexibly linked phenyl-amidine with water-mediated H-bonds to the bases at the floor of the minor groove suggested that it might rotate in time spans accessible in MD. To test this idea, we conducted multimicrosecond MD simulations to determine if these phenyl rotations could be observed for a bound compound. In a 3 μs simulation, highly dynamic torsional motions were observed for the -OCH₂-linked phenyl but not for the other phenyl. The dynamics periodically reached a level to allow 180° rotation of the phenyl while it was still bound in the minor groove. This is the first observation of rotation of a phenyl bound to DNA, and the results provide mechanistic details about how a rotation can occur as well as how mixed bp recognition can occur for monomer compounds bound to the minor groove.

A Probabilistic Framework for Constructing Temporal Relations in Replica Exchange Molecular Trajectories

Knowledge of the structure and dynamics of biomolecules is essential for elucidating the underlying mechanisms of biological processes. Given the stochastic nature of many biological processes, like protein unfolding, it is almost impossible that two independent simulations will generate the exact same sequence of events, which makes direct analysis of simulations difficult. Statistical models like Markov chains, transition networks, etc. help in shedding some light on the mechanistic nature of such processes by predicting long-time dynamics of these systems from short simulations. However, such methods fall short in analyzing trajectories with partial or no temporal information, for example, replica exchange molecular dynamics or Monte Carlo simulations. In this work, we propose a probabilistic algorithm, borrowing concepts from graph theory and machine learning, to extract reactive pathways from molecular trajectories in the absence of temporal data. A suitable vector representation was chosen to represent each frame in the macromolecular trajectory (as a series of interaction and conformational energies), and dimensionality reduction was performed using principal component analysis (PCA). The trajectory was then clustered using a density-based clustering algorithm, where each cluster represents a metastable state on the potential energy surface (PES) of the biomolecule under study. A graph was created with these clusters as nodes with the edges learned using an iterative expectation maximization algorithm. The most reactive path is conceived as the widest path along this graph. We have tested our method on RNA hairpin unfolding trajectory in aqueous urea solution. Our method makes the understanding of the mechanism of unfolding in the RNA hairpin molecule more tractable. As this method does not rely on temporal data, it can be used to analyze trajectories from Monte Carlo sampling techniques and replica exchange molecular dynamics (REMD).

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Dissecting the Innate Immune Recognition of Opioid Inactive Isomer (+)-Naltrexone Derived Toll-like Receptor 4 (TLR4) Antagonists

The opioid inactive isomer (+)-naltrexone is one of the rare Toll-like receptor 4 (TLR4) antagonists with good blood-brain barrier (BBB) permeability, which is a lead with promising potential for treating neuropathic pain and drug addiction. (+)-Naltrexone targets the lipopolysaccharides (LPS) binding pocket of myeloid differentiation protein 2 (MD-2) and blocks innate immune TLR4 signaling. However, the details of the molecular interactions of (+)-naltrexone and its derivatives with MD-2 are not fully understood, which hinders the ligand-based drug discovery. Herein, in silico and in vitro assays were performed to elucidate the innate immune recognition of the opioid inactive (+)-isomers. The results showed that the conserved LPS binding pocket of MD-2 accommodated these opioid inactive (+)-isomers. The calculated binding free energies of (+)-naltrexone and its derivatives in complex with MD-2 correlated well with their experimental binding affinities and TLR4 antagonistic activities. Hydrophobic residues in the MD-2 cavity interacted directly with these (+)-naltrexone based TLR4 antagonists and principally participated in ligand binding. Increasing the hydrophobicity of substituted group at N-17 improved its TLR4 antagonistic activity, while charged groups disfavored the binding with MD-2. Molecular dynamics (MD) simulations showed the binding of (+)-naltrexone or its derivatives to MD-2 stabilized the "collapsed" conformation of MD-2, consequently blocking the binding and signaling of TLR4. Thermodynamics and dynamic analysis showed the topology of substituted group at N-17 of (+)-naltrexone affected the binding with MD-2 and TLR4 antagonistic activity. This study provides a molecular insight into the innate immune recognition of opioid inactive (+)-isomers, which would be of great help for the development of next-generation of (+)-opioid based TLR4 antagonists.

In Silico Study, Synthesis, and Cytotoxic Activities of Porphyrin Derivatives

Five known porphyrins, 5,10,15,20-tetrakis(p-tolyl)porphyrin (TTP), 5,10,15,20-tetrakis(p-bromophenyl)porphyrin (TBrPP), 5,10,15,20-tetrakis(p-aminophenyl)porphyrin (TAPP), 5,10,15-tris(tolyl)-20-mono(p-nitrophenyl)porphyrin (TrTMNP), 5,10,15-tris(tolyl)-20-mono(p-aminophenyl)porphyrin (TrTMAP), and three novel porphyrin derivatives, 5,15-di-[bis(3,4-ethylcarboxymethylenoxy)phenyl]-10,20-di(p-tolyl)porphyrin (DBECPDTP), 5,10-di-[bis(3,4-ethylcarboxymethylenoxy)phenyl]-15,20-di-(methylpyrazole-4-yl)porphyrin (cDBECPDPzP), 5,15-di-[bis(3,4-ethylcarboxymethylenoxy)phenyl]-10,20-di-(methylpyrazole-4-yl)porphyrin (DBECPDPzP), were used to study their interaction with protein targets (in silico study), and were synthesized. Their cytotoxic activities against cancer cell lines were tested using 3-(4,5-dimetiltiazol-2-il)-2,5-difeniltetrazolium bromide (MTT) assay. The interaction of porphyrin derivatives with carbonic anhydrase IX (CAIX) and REV-ERBβ proteins were studied by molecular docking and molecular dynamic simulation. In silico study results reveal that DBECPDPzP and TrTMNP showed the highest binding interaction with REV- ERBβ and CAIX, respectively, and both complexes of DBECPDPzP-REV-ERBβ and TrTMNP-CAIX showed good and comparable stability during molecular dynamic simulation. The studied porphyrins have selective growth inhibition activities against tested cancer cells and are categorized as marginally active compounds based on their IC₅₀.

A modular design for minor groove binding and recognition of mixed base pair sequences of DNA

The design and synthesis of compounds that target mixed, AT/GC, DNA sequences is described. The design concept connects two N-methyl-benzimidazole-thiophene single GC recognition units with a flexible linker that lets the compound fit the shape and twist of the DNA minor groove while covering a full turn of the double helix.

Experimental Probing and Molecular Dynamics Simulation of the Molecular Recognition of DNA Duplexes by the Flavonoid Luteolin

Luteolin (C₁₅H₁₀O₆) is an important flavonoid found in many fruits, plants, medicinal herbs, and vegetables exhibiting many pharmacological properties. The anticancer, antitumor, antioxidant, and anti-inflammatory activities of luteolin have been reported. The pharmacological action of small molecules is dependent upon its interaction with biomacromolecules. The interactions of small molecules with DNA play a major role in the transcription and translation process. In this work, we explored the energetic profile of DNA-luteolin interaction by isothermal titration calorimetry (ITC). The effect of temperature and salt concentration on DNA binding was examined by UV-Vis method. The mode of interaction was further probed by UV melting temperature analysis and differential scanning calorimetry. An atomic level insight on the recognition of luteolin with DNA was achieved by employing molecular dynamics (MD) simulation on luteolin in complex with AT- and GC-rich DNA sequences. AMBER force field proves to be appropriate in providing an understanding on the binding mode and specificity of luteolin with duplex DNA. MD results suggest a minor groove binding of luteolin with DNA and the binding free energy obtained is in agreement with the experimental results.

Focus on PNA Flexibility and RNA Binding using Molecular Dynamics and Metadynamics

Peptide Nucleic Acids (PNAs) can efficiently target DNA or RNA acting as chemical tools for gene regulation. Their backbone modification and functionalization is often used to increase the affinity for a particular sequence improving selectivity. The understanding of the trading forces that lead the single strand PNA to bind the DNA or RNA sequence is preparatory for any further rational design, but a clear and unique description of this process is still not complete. In this paper we report further insights into this subject, by a computational investigation aiming at the characterization of the conformations of a single strand PNA and how these can be correlated to its capability in binding DNA/RNA. Employing Metadynamics we were able to better define conformational pre-organizations of the single strand PNA and γ-modified PNA otherwise unrevealed through classical molecular dynamics. Our simulations driven on backbone modified PNAs lead to the conclusion that this γ-functionalization affects the single strand preorganization and targeting properties to the DNA/RNA, in agreement with circular dichroism (CD) spectra obtained for this class of compounds. MD simulations on PNA:RNA dissociation and association mechanisms allowed to reveal the critical role of central bases and preorganization in the binding process.

Computational and Experimental Characterization of Patient Derived Mutations Reveal an Unusual Mode of Regulatory Spine Assembly and Drug Sensitivity in EGFR Kinase

The catalytic activation of protein kinases requires precise positioning of key conserved catalytic and regulatory motifs in the kinase core. The Regulatory Spine (RS) is one such structural motif that is dynamically assembled upon kinase activation. The RS is also a mutational hotspot in cancers; however, the mechanisms by which cancer mutations impact RS assembly and kinase activity are not fully understood. In this study, through mutational analysis of patient derived mutations in the RS of EGFR kinase, we identify an activating mutation, M766T, at the RS3 position. RS3 is located in the regulatory αC-helix, and a series of mutations at the RS3 position suggest a strong correlation between the amino acid type present at the RS3 position and ligand (EGF) independent EGFR activation. Small polar amino acids increase ligand independent activity, while large aromatic amino acids decrease kinase activity. M766T relies on the canonical asymmetric dimer for full activation. Molecular modeling and molecular dynamics simulations of WT and mutant EGFR suggest a model in which M766T activates the kinase domain by disrupting conserved autoinhibitory interactions between M766 and hydrophobic residues in the activation segment. In addition, a water mediated hydrogen bond network between T766, the conserved K745-E762 salt bridge, and the backbone amide of the DFG motif is identified as a key determinant of M766T-mediated activation. M766T is resistant to FDA approved EGFR inhibitors such as gefitinib and erlotinib, and computational estimation of ligand binding free energy identifies key residues associated with drug sensitivity. In sum, our studies suggest an unusual mode of RS assembly and oncogenic EGFR activation, and provide new clues for the design of allosteric protein kinase inhibitors.

Molecular dynamics simulations provide insights into the substrate specificity of FAOX family members

Enzymatic assays based on Fructosyl Amino Acid Oxidases (FAOX) represent a potential, rapid and economical strategy to measure glycated hemoglobin (HbA1c), which is in turn a reliable method to monitor the insurgence and the development of diabetes mellitus.

Molecular Dynamics: New Frontier in Personalized Medicine

The field of drug discovery has witnessed infinite development over the last decade with the demand for discovery of novel efficient lead compounds. Although the development of novel compounds in this field has seen large failure, a breakthrough in this area might be the establishment of personalized medicine. The trend of personalized medicine has shown stupendous growth being a hot topic after the successful completion of Human Genome Project and 1000 genomes pilot project. Genomic variant such as SNPs play a vital role with respect to inter individual's disease susceptibility and drug response. Hence, identification of such genetic variants has to be performed before administration of a drug. This process requires high-end techniques to understand the complexity of the molecules which might bring an insight to understand the compounds at their molecular level. To sustenance this, field of bioinformatics plays a crucial role in revealing the molecular mechanism of the mutation and thereby designing a drug for an individual in fast and affordable manner. High-end computational methods, such as molecular dynamics (MD) simulation has proved to be a constitutive approach to detecting the minor changes associated with an SNP for better understanding of the structural and functional relationship. The parameters used in molecular dynamic simulation elucidate different properties of a macromolecule, such as protein stability and flexibility. MD along with docking analysis can reveal the synergetic effect of an SNP in protein-ligand interaction and provides a foundation for designing a particular drug molecule for an individual. This compelling application of computational power and the advent of other technologies have paved a promising way toward personalized medicine. In this in-depth review, we tried to highlight the different wings of MD toward personalized medicine.

Power-Law Solvation Dynamics in G-Quadruplex DNA: Role of Hydration Dynamics on Ligand Solvation inside DNA

G-quadruplex DNA (GqDNA) structures act as promising anticancer targets for small-molecules (ligands). Solvation dynamics of a ligand (DAPI: 4',6-diamidino-2-phenylindole) inside antiparallel-GqDNA is studied through direct comparison of time-resolved experiments to molecular dynamics (MD) simulation. Dynamic Stokes shifts of DAPI in GqDNA prepared in H2O buffer and D2O are compared to find the effect of water on ligand solvation. Experimental dynamics (in H2O) is then directly compared with the dynamics computed from 65 ns simulation on the same DAPI-GqDNA complex. Ligand solvation follows power-law relaxation (summed with fast exponential relaxation) from ~100 fs to 10 ns. Simulation results show relaxation below ~5 ps is dominated by water motion, while both water and DNA contribute comparably to dictate long-time power-law dynamics. Ion contribution is, however, found to be negligible. Simulation results also suggest that anomalous solvation dynamics may have origin in subdiffusive motion of perturbed water near GqDNA.

Molecular docking and dynamics simulations on the interaction of cationic porphyrin-anthraquinone hybrids with DNA G-quadruplexes

A series of cationic porphyrin-anthraquinone hybrids bearing either pyridine, imidazole, or pyrazole rings at the meso-positions have been investigated for their interaction with DNA G-quadruplexes by employing molecular docking and molecular dynamics simulations. Three types of DNA G-quadruplexes were utilized, which comprise parallel, antiparallel, and mixed hybrid topologies. The porphyrin hybrids have a preference to bind with parallel and mixed hybrid structures compared to the antiparallel structure. This preference arises from the end stacking of porphyrin moiety following G-stem and loop binding of anthraquinone tail, which is not found in the antiparallel due to the presence of diagonal and lateral loops that crowd the G-quartet. The binding to the antiparallel, instead, occurred with poorer affinity through both the loop and wide groove. All sites of porphyrin binding were confirmed by 6 ns molecular dynamics simulation, as well as by the negative value of the total binding free energies that were calculated using the MMPBSA method. Free energy analysis shows that the favorable contribution came from the electrostatic term, which supposedly originated from the interaction of either cationic pyridinium, pyrazole, or imidazole groups and the anionic phosphate backbone, and also from the van der Waals energy, which primarily contributed through end stacking interaction.

The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities

INTRODUCTION

The molecular mechanics energies combined with the Poisson-Boltzmann or generalized Born and surface area continuum solvation (MM/PBSA and MM/GBSA) methods are popular approaches to estimate the free energy of the binding of small ligands to biological macromolecules. They are typically based on molecular dynamics simulations of the receptor-ligand complex and are therefore intermediate in both accuracy and computational effort between empirical scoring and strict alchemical perturbation methods. They have been applied to a large number of systems with varying success.

AREAS COVERED

The authors review the use of MM/PBSA and MM/GBSA methods to calculate ligand-binding affinities, with an emphasis on calibration, testing and validation, as well as attempts to improve the methods, rather than on specific applications.

EXPERT OPINION

MM/PBSA and MM/GBSA are attractive approaches owing to their modular nature and that they do not require calculations on a training set. They have been used successfully to reproduce and rationalize experimental findings and to improve the results of virtual screening and docking. However, they contain several crude and questionable approximations, for example, the lack of conformational entropy and information about the number and free energy of water molecules in the binding site. Moreover, there are many variants of the method and their performance varies strongly with the tested system. Likewise, most attempts to ameliorate the methods with more accurate approaches, for example, quantum-mechanical calculations, polarizable force fields or improved solvation have deteriorated the results.

Probing the relationship between anti-Pneumocystis carinii activity and DNA binding of bisamidines by molecular dynamics simulations

The anti-Pneumocystis carinii activity of 13 synthetic pentamidine analogs was analyzed. The experimental differences in melting points of DNA dodecamer 5'-(CGCGAATTCGCG)₂-3' complexes (ΔT_m), and in the biological activity measured using ATP bioluminescence assay (IC₅₀) together with the theoretical free energy of DNA-ligand binding estimated by the proposed computational protocol, showed that the experimental activity of the tested pentamidines appeared to be due to the binding to the DNA minor groove with extended AT sequences. The effect of heteroatoms in the aliphatic linker, and the sulfonamide or methoxy substituents on the compound inducing changes in the interactions with the DNA minor groove was examined and was correlated with biological activity. In computational analysis, the explicit solvent approximation with the discrete water molecules was taken into account, and the role of water molecules in the DNA-ligand complexes was defined.

Physics-based potentials for the coupling between backbone- and side-chain-local conformational states in the UNited RESidue (UNRES) force field for protein simulations

The UNited RESidue (UNRES) model of polypeptide chains is a coarse-grained model in which each amino-acid residue is reduced to two interaction sites, namely, a united peptide group (p) located halfway between the two neighboring α-carbon atoms (Cαs), which serve only as geometrical points, and a united side chain (SC) attached to the respective Cα. Owing to this simplification, millisecond molecular dynamics simulations of large systems can be performed. While UNRES predicts overall folds well, it reproduces the details of local chain conformation with lower accuracy. Recently, we implemented new knowledge-based torsional potentials (Krupa et al. J. Chem. Theory Comput. 2013, 9, 4620–4632) that depend on the virtual-bond dihedral angles involving side chains: Cα···Cα···Cα···SC (τ(1)), SC···Cα···Cα···Cα (τ(2)), and SC···Cα···Cα···SC (τ(3)) in the UNRES force field. These potentials resulted in significant improvement of the simulated structures, especially in the loop regions. In this work, we introduce the physics-based counterparts of these potentials, which we derived from the all-atom energy surfaces of terminally blocked amino-acid residues by Boltzmann integration over the angles λ(1) and λ(2) for rotation about the Cα···Cα virtual-bond angles and over the side-chain angles χ. The energy surfaces were, in turn, calculated by using the semiempirical AM1 method of molecular quantum mechanics. Entropy contribution was evaluated with use of the harmonic approximation from Hessian matrices. One-dimensional Fourier series in the respective virtual-bond-dihedral angles were fitted to the calculated potentials, and these expressions have been implemented in the UNRES force field. Basic calibration of the UNRES force field with the new potentials was carried out with eight training proteins, by selecting the optimal weight of the new energy terms and reducing the weight of the regular torsional terms. The force field was subsequently benchmarked with a set of 22 proteins not used in the calibration. The new potentials result in a decrease of the root-mean-square deviation of the average conformation from the respective experimental structure by 0.86 Å on average; however, improvement of up to 5 Å was observed for some proteins.