Binding free energy, energy and entropy calculations using simple model systems

Molecular docking and MD simulations reveal protease inhibitors block the catalytic residues in Prp8 intein of Aspergillus fumigatus: a potential target for antimycotics

Resistance to azoles and amphotericin B especially in Aspergillus fumigatus is a growing concern towards the treatment of invasive fungal infection. At this critical juncture, intein splicing would be a productive, and innovative target to establish therapies against resistant strains. Intein splicing is the central event for the activation of host protein, essential for the growth and survival of various microorganisms including A. fumigatus. The splicing process is a four-step protease-like nucleophilic cascade. Thus, we hypothesise that protease inhibitors would successfully halt intein splicing and potentially restrict the growth of the aforementioned pathogen. Using Rosetta Fold and molecular dynamics simulations, we modelled Prp8 intein structure; resembling classic intein fold with horse shoe shaped splicing domain. To fully comprehend the active site of Afu Prp8 intein, C1, T62, H65, H818, N819 from intein sequences and S820, the first C-extein residue are selected. Molecular docking shows that two FDA-approved drugs, i.e. Lufotrelvir and Remdesivir triphosphate efficiently interact with Prp8 intein from the assortment of 212 protease inhibitors. MD simulation portrayed that Prp8 undergoes conformational change upon ligand binding, and inferred the molecular recognition and stability of the docked complexes. Per-residue decomposition analysis confirms the importance of F: block R802, V803, and Q807 binding pocket in intein splicing domain towards recognition of inhibitors, along with active site residues through strong hydrogen bonds and hydrophobic contacts. However, in vitro and in vivo assays are required to confirm the inhibitory action on Prp8 intein splicing; which may pave the way for the development of new antifungals for A. fumigatus.Communicated by Ramaswamy H. Sarma.

Exploiting the potential of natural polyphenols as antivirals against monkeypox envelope protein F13 using machine learning and all-atoms MD simulations

The re-emergence of monkeypox (MPX), in the era of COVID-19 pandemic is a new global menace. Regardless of its leniency, there are chances of MPX expediting severe health deterioration. The role of envelope protein, F13 as a critical component for production of extracellular viral particles makes it a crucial drug target. Polyphenols, exhibiting antiviral properties have been acclaimed as an effective alternative to the traditional treatment methods for management of viral diseases. To facilitate the development of potent MPX specific therapeutics, herein, we have employed state-of-the-art machine learning techniques to predict a highly accurate 3-dimensional structure of F13 as well as identify binding hotspots on the protein surface. Additionally, we have effectuated high-throughput virtual screening methodology on 57 potent natural polyphenols having antiviral activities followed by all-atoms molecular dynamics (MD) simulations, to substantiate the mode of interaction of F13 protein and polyphenol complexes. The structure-based virtual screening based on Glide SP, XP and MM/GBSA scores enables the selection of six potent polyphenols having higher binding affinity towards F13. Non-bonded contact analysis, of pre- and post- MD complexes propound the critical role of Glu143, Asp134, Asn345, Ser321 and Tyr320 residues in polyphenol recognition, which is well supported by per-residue decomposition analysis. Close-observation of the structural ensembles from MD suggests that the binding groove of F13 is mostly hydrophobic in nature. Taken together, this structure-based analysis from our study provides a lead on Myricetin, and Demethoxycurcumin, which may act as potent inhibitors of F13. In conclusion, our study provides new insights into the molecular recognition and dynamics of F13-polyphenol bound states, offering new promises for development of antivirals to combat monkeypox. However, further in vitro and in vivo experiments are necessary to validate these results.

Solvent quality and solvent polarity in polypeptides

Using molecular dynamics and thermodynamic integration, we report on the solvation process of seven polypeptides (GLY, ALA, ILE, ASN, LYS, ARG, GLU) in water and in cyclohexane. The polypeptides are selected to cover the full hydrophobic scale while varying their chain length from tri- to undeca-homopeptides, providing indications on possible non-additivity effects as well as the role of the peptide backbone in the overall stability of the polypeptides. The use of different solvents and different polypeptides allows us to investigate the relation between solvent quality - the capacity of a given solvent to fold/unfold a given biopolymer often described on a scale ranging from "good" to "poor"; and solvent polarity - related to the specific interactions of any solvent with respect to a reference solvent. Undeca-glycine is found to be the only polypeptide to have a stable collapse in water (polar solvent), with the other hydrophobic polypeptides displaying repeated folding and unfolding events in water, with polar polypeptides presenting even more complex behavior. By contrast, all polypeptides are found to keep an extended conformation in cyclohexane, irrespective of their polarity. All considered polypeptides are also found to have favorable solvation free energy independent of the solvent polarity and their intrinsic hydrophobicity, clearly highlighting the prominent stabilizing role of the peptide backbone - with the solvation process largely enthalpically dominated in polar polypeptides and partially entropically driven for hydrophobic polypeptides. Our study thus reveals the complexity of the solvation process of polypeptides defying the common view "like dissolves like", with the solute polarity playing the most prominent role. The absence of mirror symmetry upon the inversion of polarities of both the solvent and the polypeptides is confirmed.

Galantamine Based Novel Acetylcholinesterase Enzyme Inhibitors: A Molecular Modeling Design Approach

Acetylcholinesterase (AChE) enzymes play an essential role in the development of Alzheimer's disease (AD). Its excessive activity causes several neuronal problems, particularly psychopathies and neuronal cell death. A bioactive pose on the hAChE B site of the human acetylcholinesterase (hAChE) enzyme employed in this investigation, which was obtained from the Protein Data Bank (PDB ID 4EY6), allowed for the prediction of the binding affinity and free binding energy between the protein and the ligand. Virtual screening was performed to obtain structures similar to Galantamine (GNT) with potential hAChE activity. The top 200 hit compounds were prioritized through the use of filters in ZincPharmer, with special features related to the pharmacophore. Critical analyses were carried out, such as hierarchical clustering analysis (HCA), ADME/Tox predictions, molecular docking, molecular simulation studies, synthetic accessibility (SA), lipophilicity, water solubility, and hot spots to confirm the stable binding of the two promising molecules (ZINC16951574-LMQC2, and ZINC08342556-LMQC5). The metabolism prediction, with metabolites M3-2, which is formed by Glutathionation reaction (Phase II), M1-2, and M2-2 formed from the reaction of S-oxidation and Aliphatic hydroxylation (Phase I), were both reactive but with no side effects. Theoretical synthetic routes and prediction of synthetic accessibility for the most promising compounds are also proposed. In conclusion, this study shows that in silico modeling can be used to create new drug candidate inhibitors for hAChE. The compounds ZINC16951574-LMQC2, and ZINC08342556-LMQC5 are particularly promising for oral administration because they have a favorable drug-likeness profile, excellent lipid solubility, high bioavailability, and adequate pharmacokinetics.

Hierarchical Self-Assembly Pathways of Peptoid Helices and Sheets

Peptoids (N-substituted glycines) are a class of tailorable synthetic peptidomic polymers. Amphiphilic diblock peptoids have been engineered to assemble 2D crystalline lattices with applications in catalysis and molecular separations. Assembly is induced in an organic solvent/water mixture by evaporating the organic phase, but the assembly pathways remain uncharacterized. We conduct all-atom molecular dynamics simulations of Nbrpe6Nc6 as a prototypical amphiphilic diblock peptoid comprising an NH₂-capped block of six hydrophobic N-((4-bromophenyl)ethyl)glycine residues conjugated to a polar NH₃(CH₂)₅CO tail. We identify a thermodynamically controlled assembly mechanism by which monomers assemble into disordered aggregates that self-order into 1D chiral helical rods then 2D achiral crystalline sheets. We support our computational predictions with experimental observations of 1D rods using small-angle X-ray scattering, circular dichroism, and atomic force microscopy and 2D crystalline sheets using X-ray diffraction and atomic force microscopy. This work establishes a new understanding of hierarchical peptoid assembly and principles for the design of peptoid-based nanomaterials.

A Benchmark of Electrostatic Method Performance in Relative Binding Free Energy Calculations

Relative free energy calculations are fast becoming a critical part of early stage pharmaceutical design, making it important to know how to obtain the best performance with these calculations in applications that could span hundreds of calculations and molecules. In this work, we compared two different treatments of long-range electrostatics, Particle Mesh Ewald (PME) and Reaction Field (RF), in relative binding free energy calculations using a nonequilibrium switching protocol. We found simulations using RF achieve comparable results to those using PME but gain more efficiency when using CPU and similar performance using GPU. The results from this work encourage more use of RF in molecular simulations.

Ligand design by targeting a binding site water

Solvent reorganization is a major driving force of protein–ligand association, but the contribution of binding site waters to ligand affinity is poorly understood. We investigated how altered interactions with a water network can influence ligand binding to a receptor. A series of ligands of the A_2A adenosine receptor, which either interacted with or displaced an ordered binding site water, were studied experimentally and by molecular dynamics simulations. An analog of the endogenous ligand that was unable to hydrogen bond to the ordered water lost affinity and this activity cliff was captured by molecular dynamics simulations. Two compounds designed to displace the ordered water from the binding site were then synthesized and evaluated experimentally, leading to the discovery of an A_2A agonist with nanomolar activity. Calculation of the thermodynamic profiles resulting from introducing substituents that interacted with or displaced the ordered water showed that the gain of binding affinity was enthalpy driven. Detailed analysis of the energetics and binding site hydration networks revealed that the enthalpy change was governed by contributions that are commonly neglected in structure-based drug optimization. In particular, simulations suggested that displacement of water from a binding site to the bulk solvent can lead to large energy contributions. Our findings provide insights into the molecular driving forces of protein–ligand binding and strategies for rational drug design.

Solvent reorganization is a major driving force of protein–ligand association, but the contribution of binding site waters to ligand affinity is poorly understood.

Experimental characterization of the association of β-cyclodextrin and eight novel cyclodextrin derivatives with two guest compounds

We investigate the binding of native β-cyclodextrin (β-CD) and eight novel β-CD derivatives with two different guest compounds, using isothermal calorimetry and 2D NOESY NMR. In all cases, the stoichiometry is 1:1 and binding is exothermic. Overall, modifications at the 3' position of β-CD, which is at the secondary face, weaken binding by several kJ/mol relative to native β-CD, while modifications at the 6' position (primary face) maintain or somewhat reduce the binding affinity. The variations in binding enthalpy are larger than the variations in binding free energy, so entropy-enthalpy compensation is observed. Characterization of the bound conformations with NOESY NMR shows that the polar groups of the guests may be situated at either face, depending on the host molecule, and, in some cases, both orientations are populated. The present results were used in the SAMPL7 blinded prediction challenge whose results are detailed in the same special issue of JCAMD.

Rigid Scaffolds: Synthesis of 2,6-Bridged Piperazines with Functional Groups in all three Bridges

The activity of pharmacologically active compounds can be increased by presenting a drug in a defined conformation, which fits exactly into the binding pocket of its target. Herein, the piperazine scaffold was conformationally restricted by substituted C2- or C3-bridges across the 2- and 6-position. At first, a three-step, one-pot procedure was developed to obtain reproducibly piperazine-2,6-diones with various substituents at the N-atoms in high yields. Three strategies for bridging of piperazine-2,6-diones were pursued: 1. The bicyclic mixed ketals 8-benzyl-6-ethoxy-3-(4-methoxybenzyl)-6-(trimethylsilyloxy)-3,8-diazabicyclo[3.2.1]octane-2,4-diones were prepared by Dieckmann analogous cyclization of 2-(3,5-dioxopiperazin-2-yl)acetates. 2. Stepwise allylation, hydroboration and oxidation of piperazine-2,6-diones led to 3-(3,5-dioxopiperazin-2-yl)propionaldehydes. Whereas reaction of such an aldehyde with base provided the bicyclic alcohol 9-benzyl-6-hydroxy-3-(4-methoxybenzyl)-3,9-diazabicyclo[3.3.1]nonane-2,4-dione in only 10 % yield, the corresponding sulfinylimines reacted with base to give N-(2,4-dioxo-3,9-diazabicyclo[3.3.1]nonan-6-yl)-2-methylpropane-2-sulfinamides in >66 % yield. 3. Transformation of a piperazine-2,6-dione with 1,4-dibromobut-2-ene and 3-halo-2-halomethylprop-1-enes provided 3,8-diazabicyclo[3.2.1]octane-2,4-dione and 3,9-diazabicyclo[3.3.1]nonane-2,4-dione with a vinyl group at the C2- or a methylene group at the C3-bridge, respectively. Since bridging via sulfinylimines and the one-pot bridging with 3-bromo-2-bromomethylprop-1-ene gave promising yields, these strategies will be exploited for the synthesis of novel receptor ligands bearing various substituents in a defined orientation at the carbon bridge.

Entropy deepens loading chemical potentials of small alcohols by narrow carbon nanotubes

Small alcohol confinement within narrow carbon nanotubes has been extensively and systematically studied via rigorous free-energy calculations.

Water in protein hydration and ligand recognition

This review describes selected basics of water in biomolecular recognition. We focus on a qualitative understanding of the most important physical aspects, how these change in magnitude between bulk water and protein environment, and how the roles that water plays for proteins arise from them. These roles include mechanical support, thermal coupling, dielectric screening, mass and charge transport, and the competition with a ligand for the occupation of a binding site. The presence or absence of water has ramifications that range from the thermodynamic binding signature of a single ligand up to cellular survival. The large inhomogeneity in water density, polarity and mobility around a solute is hard to assess in experiment. This is a source of many difficulties in the solvation of protein models and computational studies that attempt to elucidate or predict ligand recognition. The influence of water in a protein binding site on the experimental enthalpic and entropic signature of ligand binding is still a point of much debate. The strong water‐water interaction in enthalpic terms is counteracted by a water molecule's high mobility in entropic terms. The complete arrest of a water molecule's mobility sets a limit on the entropic contribution of a water displacement process, while the solvent environment sets limits on ligand reactivity.

Free-Energy Calculations for Bioisosteric Modifications of A3 Adenosine Receptor Antagonists

Adenosine receptors are a family of G protein-coupled receptors with increased attention as drug targets on different indications. We investigate the thermodynamics of ligand binding to the A₃ adenosine receptor subtype, focusing on a recently reported series of diarylacetamidopyridine inhibitors via molecular dynamics simulations. With a combined approach of thermodynamic integration and one-step perturbation, we characterize the impact of the charge distribution in a central heteroaromatic ring on the binding affinity prediction. Standard charge distributions according to the GROMOS force field yield values in good agreement with the experimental data and previous free energy calculations. Subsequently, we examine the thermodynamics of inhibitor binding in terms of the energetic and entropic contributions. The highest entropy penalties are found for inhibitors with methoxy substituents in meta position of the aryl groups. This bulky group restricts rotation of aromatic rings attached to the pyrimidine core which leads to two distinct poses of the ligand. Our predictions support the previously proposed binding pose for the o-methoxy ligand, yielding in this case a very good correlation with the experimentally measured affinities with deviations below 4 kJ/mol.

Redox thermodynamics of B-class dye-decolorizing peroxidases

With >5000 annotated genes dye-decolorizing peroxidases (DyPs) represent a heme b peroxidase family of broad functional diversity. Bacterial B-class DyPs are poor peroxidases of unknown physiological function. Hydrogen peroxide efficiently mediates the rapid formation of Compound I in B-class DyPs, which, however, is stable and shows modest reactivity towards organic and inorganic electron donors. To understand these characteristics, we have investigated the redox thermodynamics of the one-electron reduction of the ferric high-spin form of wild-type B-class DyP from the pathogenic bacterium Klebsiella pneumoniae (KpDyP) and the variants D143A, R232A and D143A/R232A. These distal amino acids are fully conserved in all DyPs and play important roles in Compound I formation and maintenance of the heme cavity architecture and substrate access route(s). The E°' values of the respective redox couples Fe(III)/Fe(II) varied from -350 mV (wild-type KpDyP) to -299 mV (D143A/R232A) at pH 7.0. Variable-temperature spectroelectrochemical experiments revealed that the reduction reaction of B-class DyPs is enthalpically unfavored but entropically favored with significant differences in enthalpic and entropic contributions to E°' between the four proteins. Molecular dynamics simulations demonstrated the impact of solvent reorganization on the entropy change during reduction reaction and revealed the dynamics and restriction of substrate access channels. Obtained data are discussed with respect to the poor peroxidase activities of B-class DyPs and compared with heme peroxidases from other (super)families as well as with chlorite dismutases, which do not react with hydrogen peroxide but share a similar fold and heme cavity architecture.

Consensus Conformations of Dinucleoside Monophosphates Described with Well-Converged Molecular Dynamics Simulations

Dinucleoside monophosphates (DNMPs) have been described using various experimental approaches as flexible molecules which generate ensembles populating at least a small set of different conformations in solution. However, due to limitations of each approach in its ability to delineate the ensemble of conformations, an accurate and quantitative description of certain conformational features has not been performed for all DNMPs. Here, we apply a temperature replica-exchange molecular dynamics approach to fully and quickly converge conformational distributions of all RNA DNMPs immersed in the TIP3P water model using the AMBER ff14 force field. For a selection of DNMPs, the conformational ensembles were also generated when immersed in the OPC water model using alternative AMBER and CHARMM force fields. The OPC water model and other force field choices did not introduce new conformational classes but shifted the populations among existing conformations. Except for pyrimidine-pyrimidine dinucleosides, all other DNMPs populated four major conformations (which are defined in the main text and labeled A-form, Ladder, Inverted, and Sheared), in addition to an Extended form. Pyrimidine-pyrimidines did not generate the Sheared conformation. Distinguishing features and stabilizing factors of each conformation were identified and assessed based on the known experimental interpretations. The configuration of the glycosidic bond and the nonbonding interactions of hydrogen bond acceptors with the 2'-hydroxyl group were found to play determining roles in stabilizing particular conformations which could serve as a guide for potential force field modifications to improve the accuracy. Additionally, we computed stacking free energies based on the DNMP conformational distributions and found significant discrepancies with a previous study. Our investigation determined that the AMBER force field was incorrectly implemented in the previous study. In the future, this simulation approach can be used to quickly analyze the effects of new force field modifications in shifting the conformational populations of DNMPs, and can can be further applied to foresee such effects in larger RNA motifs including tetranucleotides and tetraloops.

Validation of Molecular Simulation: An Overview of Issues

Computer simulation of molecular systems enables structure-energy-function relationships of molecular processes to be described at the sub-atomic, atomic, supra-atomic, or supra-molecular level. To interpret results of such simulations appropriately, the quality of the calculated properties must be evaluated. This depends on the way the simulations are performed and on the way they are validated by comparison to values Q^exp of experimentally observable quantities Q. One must consider 1) the accuracy of Q^exp , 2) the accuracy of the function Q(r^N ) used to calculate a Q-value based on a molecular configuration r^N of N particles, 3) the sensitivity of the function Q(r^N ) to the configuration r^N , 4) the relative time scales of the simulation and experiment, 5) the degree to which the calculated and experimental properties are equivalent, and 6) the degree to which the system simulated matches the experimental conditions. Experimental data is limited in scope and generally corresponds to averages over both time and space. A critical analysis of the various factors influencing the apparent degree of (dis)agreement between simulations and experiment is presented and illustrated using examples from the literature. What can be done to enhance the validation of molecular simulation is also discussed.

Computational Calorimetry: High-Precision Calculation of Host-Guest Binding Thermodynamics

We present a strategy for carrying out high-precision calculations of binding free energy and binding enthalpy values from molecular dynamics simulations with explicit solvent. The approach is used to calculate the thermodynamic profiles for binding of nine small molecule guests to either the cucurbit[7]uril (CB7) or β-cyclodextrin (βCD) host. For these systems, calculations using commodity hardware can yield binding free energy and binding enthalpy values with a precision of ∼0.5 kcal/mol (95% CI) in a matter of days. Crucially, the self-consistency of the approach is established by calculating the binding enthalpy directly, via end point potential energy calculations, and indirectly, via the temperature dependence of the binding free energy, i.e., by the van't Hoff equation. Excellent agreement between the direct and van't Hoff methods is demonstrated for both host-guest systems and an ion-pair model system for which particularly well-converged results are attainable. Additionally, we find that hydrogen mass repartitioning allows marked acceleration of the calculations with no discernible cost in precision or accuracy. Finally, we provide guidance for accurately assessing numerical uncertainty of the results in settings where complex correlations in the time series can pose challenges to statistical analysis. The routine nature and high precision of these binding calculations opens the possibility of including measured binding thermodynamics as target data in force field optimization so that simulations may be used to reliably interpret experimental data and guide molecular design.

Toward Improved Force-Field Accuracy through Sensitivity Analysis of Host-Guest Binding Thermodynamics

Improving the capability of atomistic computer models to predict the thermodynamics of noncovalent binding is critical for successful structure-based drug design, and the accuracy of such calculations remains limited by nonoptimal force field parameters. Ideally, one would incorporate protein-ligand affinity data into force field parametrization, but this would be inefficient and costly. We now demonstrate that sensitivity analysis can be used to efficiently tune Lennard-Jones parameters of aqueous host-guest systems for increasingly accurate calculations of binding enthalpy. These results highlight the promise of a comprehensive use of calorimetric host-guest binding data, along with existing validation data sets, to improve force field parameters for the simulation of noncovalent binding, with the ultimate goal of making protein-ligand modeling more accurate and hence speeding drug discovery.

Thermodynamic characterization of new positive allosteric modulators binding to the glutamate receptor A2 ligand-binding domain: combining experimental and computational methods unravels differences in driving forces

Positive allosteric modulation of the ionotropic glutamate receptor GluA2 presents a potential treatment of cognitive disorders, for example, Alzheimer's disease. In the present study, we describe the synthesis, pharmacology, and thermodynamic studies of a series of monofluoro-substituted 3,4-dihydro-2H-1,2,4-benzothiadiazine 1,1-dioxides. Measurements of ligand binding by isothermal titration calorimetry (ITC) showed similar binding affinities for the modulator series at the GluA2 LBD but differences in the thermodynamic driving forces. Binding of 5c (7-F) and 6 (no-F) is enthalpy driven, and 5a (5-F) and 5b (6-F) are entropy driven. For 5d (8-F), both quantities were equal in size. Thermodynamic integration (TI) and one-step perturbation (OSP) were used to calculate the relative binding affinity of the modulators. The OSP calculations had a higher predictive power than those from TI, and combined with the shorter total simulation time, we found the OSP method to be more effective for this setup. Furthermore, from the molecular dynamics simulations, we extracted the enthalpies and entropies, and along with the ITC data, this suggested that the differences in binding free energies are largely explained by the direct ligand-surrounding enthalpies. Furthermore, we used the OSP setup to predict binding affinities for a series of polysubstituted fluorine compounds and monosubstituted methyl compounds and used these predictions to characterize the modulator binding pocket for this scaffold of positive allosteric modulators.

On the thermodynamics of carbon nanotube single-file water loading: free energy, energy and entropy calculations

Single-file water chains confined in carbon nanotubes are studied using molecular dynamics simulations and free energy calculations to characterize the pore loading in terms of energy, entropy and free energy.