Nonpolar Contributions Research Articles

Inclusion of Water Multipoles into the Implicit Solvation Framework Leads to Accuracy Gains.

The current practical "workhorses" of the atomistic implicit solvation─the Poisson-Boltzmann (PB) and generalized Born (GB) models─face fundamental accuracy limitations. Here, we propose a computationally efficient implicit solvation framework, the Implicit Water Multipole GB (IWM-GB) model, that systematically incorporates the effects of multipole moments of water molecules in the first hydration shell of a solute, beyond the dipole water polarization already present at the PB/GB level. The framework explicitly accounts for coupling between polar and nonpolar contributions to the total solvation energy, which is missing from many implicit solvation models. An implementation of the framework, utilizing the GAFF force field and AM1-BCC atomic partial charges model, is parametrized and tested against the experimental hydration free energies of small molecules from the FreeSolv database. The resulting accuracy on the test set (RMSE ∼ 0.9 kcal/mol) is 12% better than that of the explicit solvation (TIP3P) treatment, which is orders of magnitude slower. We also find that the coupling between polar and nonpolar parts of the solvation free energy is essential to ensuring that several features of the IWM-GB model are physically meaningful, including the sign of the nonpolar contributions.

Effect of Composition and Size on Surface Properties of Anti-Cancer Nanoparticles.

Liposomal formulations offer significant advantages as anticancer drug carriers for targeted drug delivery; however, due to their complexity, clinical translation has been challenging. In addition, liposomal product manufacturing has been interrupted in the past, as was the case for Doxil® (doxorubicin hydrochloride liposome injection). Here, interfacial tension (IFT) measurements were investigated as a potential physicochemical characterization tool to aid in liposomal product characterization during development and manufacturing. A pendant drop method using an optical tensiometer was used to measure the interfacial tension of various analogues of Doxil® liposomal suspensions in air and in dodecane. The effect of liposome concentration, formulation (PEG and cholesterol content), presence of encapsulated drug, as well as average particle size was analyzed. It was observed that Doxil® analog liposomes demonstrate surfactant-like behavior with a sigmoidal-shape interfacial tension vs. concentration curve. This behavior was heavily dependent on PEG content, with a complete loss of surfactant-like behavior when PEG was removed from the formulation. In addition to interfacial tension, three data analyses were identified as able to distinguish between formulations with variations in PEG, cholesterol, and particle size: (i) polar and non-polar contribution to interfacial tension, (ii) liposomal concentration at which the polar and non-polar components were equal, and (iii) rate of interfacial tension decay after droplet formation, which is indicative of how quickly liposomes migrate from the bulk of the solution to the surface. We demonstrate for the first time that interfacial tension can be used to detect certain liposomal formulation changes, such as PEG content, encapsulated drug presence, and size variability, and may make a useful addition to physicochemical characterization during development and manufacturing of liposomal products.

Dissection of the Polar and Non-Polar Contributions to Aromatic Stacking Interactions in Solution.

Aromatic stacking interactions have been a matter of study and debate due to their crucial role in chemical and biological systems. The strong dependence on orientation and solvent together with the relatively small interaction energies have made evaluation and rationalization a challenge for experimental and theoretical chemists. We have used a supramolecular cage formed by two tris(pyridylmethyl)amines units to build chemical Double Mutant Cycles (DMC) for the experimental measurement of the free energies of π‐stacking interactions. Extrapolating the substituent effects to remove the contribution due to electrostatic interactions reveals that there is a substantial contribution to the measured stacking interaction energies which is due to non‐polar interactions (−3 to −6 kJ mol−1). The perfectly flat nature of the surface of an aromatic ring gives π‐stacking an inherent advantage over non‐polar interactions with alkyl groups and accounts for the wide‐spread prevalence of stacking interactions in Nature.

Elucidating specificity of an allosteric inhibitor WNK476 among With-No-Lysine kinase isoforms using molecular dynamic simulations.

Specifically targeting the With-No-Lysine (WNK1) kinase, which is implicated in hypertension, renders a significant challenge in discovering competitive inhibitors due to the highly conserved ATP-binding pocket. However, an allosteric inhibitor may impart high specificity against the WNK kinase isoforms since it targets the less conserved site and can provide greater efficacy even under high physiological ATP concentration. In the current study, we have investigated the structural and energetic basis of the specificity of the allosteric inhibitor WNK476 against WNK kinase isoforms by combining molecular dynamics simulations and free energy calculations using molecular mechanics Poisson-Boltzmann surface area. Our study reveals that the conformational stabilization of αC-helix near the allosteric binding site, including conformational changes in activation and glycine-rich loop regions, favors the specificity of WNK476 toward WNK1. The MM/PBSA calculations suggest that the non-polar contribution from hydrophobic residues and polar solvation energy influences WNK/WNK476 complexation. Despite more favorable electrostatic and van der Waals interactions in WNK2/WNK476, WNK476 is more potent against WNK1 due to the lower contribution of disfavoring components-polar solvation and entropy. Further, we have identified that the hydrophobic residues of DLG, αC-helix, β4 , and β5 regions, and H-bond network near the β4 strand play a critical role in the specificity of WNK476 against WNK1. Finally, our study reveals that residues Leu272 , Val281 , Phe283 , and Leu369 of WNK1 actively contribute to the overall hydrophobic interactions for WNK1/WNK476. Overall, our study might help in the rational design of novel allosteric inhibitors against hypertension.

Efficient In Silico Saturation Mutagenesis of a Member of the Caspase Protease Family.

Rational-design methods have proven to be a valuable toolkit in the field of protein design. Numerical approaches such as free-energy calculations or QM/MM methods are fit to widen the understanding of a protein-sequence space but require large amounts of computational time and power. Here, we apply an efficient method for free-energy calculations that combines the one-step perturbation (OSP) with the third-power-fitting (TPF) approach. It is fit to calculate full free energies of binding from three different end states only. The nonpolar contribution to the free energies are calculated for a set of chosen amino acids from a single simulation of a judiciously chosen reference state. The electrostatic contributions, on the other hand, are predicted from simulations of the neutral and charged end states of the individual amino acids. We used this method to perform in silico saturation mutagenesis of two sites in human Caspase-2. We calculated relative binding free energies toward two different substrates that differ in their P1′ site and in their affinity toward the unmutated protease. Although being approximate, our approach showed very good agreement upon validation against experimental data. 76% of the predicted relative free energies of amino acid mutations was found to be true positives or true negatives. We observed that this method is fit to discriminate amino acid mutations because the rate of false negatives is very low (<1.5%). The approach works better for a substrate with medium/low affinity with a Matthews correlation coefficient (MCC) of 0.63, whereas for a substrate with very low affinity, the MCC was 0.38. In all cases, the combined TPF + OSP approach outperformed the linear interaction energy method.

Virtual Screening of the Indonesian Medicinal Plant and Zinc Databases for Potential Inhibitors of the RNA-Dependent RNA Polymerase (RdRp) of 2019 Novel Coronavirus

The novel coronavirus disease 19 (Covid-19) which is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been a pandemic across the world, which necessitate the need for the antiviral drug discovery. One of the potential protein targets for coronavirus treatment is RNA-dependent RNA polymerase. It is the key enzyme in the viral replication machinery, and it does not exist in human beings, therefore its targeting has been considered as a strategic approach. Here we describe the identification of potential hits from Indonesian Herbal and ZINC databases. The pharmacophore modeling was employed followed by molecular docking and dynamics simulation for 40 ns. 151 and 14480 hit molecules were retrieved from Indonesian herbal and ZINC databases, respectively. Three hits that were selected based on the structural analysis were stable during 40 ns, while binding energy prediction further implied that ZINC1529045114, ZINC169730811, and 9-Ribosyl-trans-zeatin had tighter binding affinities compared to Remdesivir. The ZINC169730811 had the strongest affinity toward RdRp compared to the other two hits including Remdesivir and its binding was corroborated by electrostatic, van der Waals, and nonpolar contribution for solvation energies. The present study offers three hits showing tighter binding to RdRp based on MM-PBSA binding energy prediction for further experimental verification.

Evaluating the Solvent Stark Effect from Temperature-Dependent Solvatochromic Shifts of Anthracene.

The solvent Stark effect on the spectral shifts of anthracene is studied with temperature‐dependent solvatochromic measurements. The Stark contribution Δv Stark to the absorption shift Δv p in polar solvents is measured to be Δv Stark=(53±35) cm−1, in reasonable agreement with dielectric continuum theory estimate of 28 cm−1, whereas the major shift Δv p∼300 cm−1 presumably originates from the solute quadrupole. We pay attention to the accurate correction of Δv p for the nonpolar contribution that is crucial when the shifts are modest in magnitude.

Hydration free energies and solvation structures with molecular density functional theory in the hypernetted chain approximation.

The capability of molecular density functional theory in its lowest, second-order approximation, equivalent to the hypernetted chain approximation in integral equations, to predict accurately the hydration free-energies and microscopic structure of molecular solutes is explored for a variety of systems: spherical hydrophobic solutes, ions, water as a solute, and the Mobley's dataset of organic molecules. The successes and the caveats of the approach are carefully pinpointed. Compared to molecular simulations with the same force field and the same fixed solute geometries, the theory describes accurately the solvation of cations, less so that of anions or generally H-bond acceptors. Overall, the electrostatic contribution to solvation free-energies of neutral molecules is correctly reproduced. On the other hand, the cavity contribution is poorly described but can be corrected using scaled-particle theory ideas. Addition of a physically motivated, one-parameter cavity correction accounting for both pressure and surface effects in the nonpolar solvation contribution yields a precision of 0.8 kcal/mol for the overall hydration free energies of the whole Mobley's dataset. Inclusion of another one-parameter cavity correction for the electrostatics brings it to 0.6 kcal/mol, that is, kBT. This is accomplished with a three-orders of magnitude numerical speed-up with respect to molecular simulations.

Investigating specificity of the anti-hypertensive inhibitor WNK463 against With-No-Lysine kinase family isoforms via multiscale simulations

The With-No-Lysine (WNK) kinase family plays a significant role in regulating cation-chloride cotransporters, blood pressure and body fluid homeostasis. Mutations in the gene of WNK family, especially in WNK1 and WNK4 are responsible for pseudohypoaldosteronism type II (PHAII), characterized by hypertension. The selective inhibition of WNK1 over other isoforms has created an immense challenge in the design of an ATP competitive inhibitor due to their high conservatism. In this work, we have compared the selectivity of the inhibitor WNK463, which was designed for WNK1 with other WNK family isoforms by comprehensive molecular modeling, docking and molecular dynamics simulations in conjunction with the Molecular Mechanics Poisson-Boltzmann Surface Area method. Our calculations show that the affinity of the inhibitor decreases in the order WNK2 > WNK1 > WNK3 > WNK4, in agreement with the experiment. Our study reveals that the inhibitor is most selective to WNK2 due to decreased polar solvation and configurational entropy compared to other isoforms. Furthermore, our analyses indicated that the nonpolar contribution from the hydrophobic residues and hydrogen bonds in the hinge region gatekeeper residue Met304 of WNK1 and its equivalent residue from other kinases played a critical role in stabilizing the inhibitor against WNK kinases. Residues Lys233, Met304, Phe356 and Leu369 of WNK1 were the essential residue differences compared to other isoforms that led to specific interactions thereby forming the basis of molecular binding pattern of binding interactions. Overall, we have identified conserved WNK-inhibitor interactions and elucidated isoform-specific interactions that could be exploited in the design of more potent and selective WNK inhibitors.Communicated by Ramaswamy H. Sarma

Insights into the structure-affinity relationships and solvation effects between OfHex1 and inhibitors using molecular dynamics simulations

OfHex1 is a potential target for the rational design of pesticides. TMG-chitotriomycin is one of the most highly specific known inhibitors of chitinolytic β-GlcNAcases from bacteria, fungi and insects. TMG-chitotriomycin and its analogues show different activities to OfHex1, dependent on the number of GlcNAc units. Subsequently, it is essential to explore how these GlcNAc unit number changes cause alterations in activity. In this study, we examined the free energy patterns and per residue decomposition of binding within the complexes of OfHex1 and a series of inhibitors, utilizing restricted molecular dynamics (MD) and water-mediated MM/GBSA calculations. The results indicated Glu328 could form a stronger polar interaction with OfHex1 inhibitors, while Trp448 and Trp490 had important non-polar contributions. Interestingly, the conformation of Trp448 was different in the open or closed state, when OfHex1 bound different inhibitors. Moreover, the water molecule that mediates the GlcNAc Ⅱ and Trp490 may be critical to stabilizing the hydrophobic interaction. Further study showed that isomerization of TMG-chitotriomycin analogs did not decrease binding affinity, however, there was a highly positive correlation between the calculated binding affinities and the experimental activity data (r2 = 0.92) when water molecules were explicitly taken into account. Moreover, the water molecules that mediated GlcNAc II and Trp490 might be critical to the stabilization of the hydrophobic interaction and cause the activity difference between TMG-(GlcNAc)2 and TMG-(GlcNAc).

Molecular dynamics investigation on the Asciminib resistance mechanism of I502L and V468F mutations in BCR-ABL

Asciminib, a highly selective non-ATP competitive inhibitor of BCR-ABL, has demonstrated to be a promising drug for patients with chronic myeloid leukemia. It is a pity that two resistant mutations (I502L and V468F) have been found during the clinical trial, which is a challenge for the curative effect of Asciminib. In this study, molecular dynamics simulations and molecular mechanics generalized Born surface area (MM-GB/SA) calculations were performed to investigate the molecular mechanism of Asciminib resistance induced by the two mutants. The obtained results indicate that the mutations have adversely influence on the binding of Asciminib to BCR-ABL, as the nonpolar contributions decline in the two mutants. In addition, I502L mutation causes α-helix I’ (αI′) to shift away from the helical bundle composed of αE, αF, and αH, making the distance between αI’ and Asciminib increased. For V468F mutant, the side chain of Phe468 occupies the bottom of the myristoyl pocket (MP), which drives Asciminib to shift toward the outside of MP. Our results provide the molecular insights of Asciminib resistance mechanism in BCR-ABL mutants, which may help the design of novel inhibitors.

Surface properties of the composite films based on poly(vinyl alcohol) and nanodiamonds as studied by wetting techniques and autoradiography

Surface concentration of oxidized nanodiamonds (NDs) and NDs modified with oleylamine in NDs–poly(vinyl alcohol) composite films has been determined by autoradiography using tritium-labeled NDs. It has been demonstrated that oxidized NDs are distributed uniformly in the composite films, and their surface concentration is the same at both sides of the film. After modification with oleylamine, ND concentration at the bottom side of the film exceeds that at the top side approximately tenfold. Zisman critical surface tension of the composite films as well as polar and non-polar contributions to the surface tension has been found. It has been shown that surface tension is determined by the polymer, and the addition of NDs does not affect its value.

Residue-level determinants of angiopoietin-2 interactions with its receptor Tie2.

We combined computational and experimental methods to interrogate the binding determinants of angiopoietin-2 (Ang2) to its receptor tyrosine kinase (RTK) Tie2-a central signaling system in angiogenesis, inflammation, and tumorigenesis. We used physics-based electrostatic and surface-area calculations to identify the subset of interfacial Ang2 and Tie2 residues that can affect binding directly. Using random and site-directed mutagenesis and yeast surface display (YSD), we validated these predictions and identified additional Ang2 positions that affected receptor binding. We then used burial-based calculations to classify the larger set of Ang2 residues that are buried in the Ang2 core, whose mutations can perturb the Ang2 structure and thereby affect interactions with Tie2 indirectly. Our analysis showed that the Ang2-Tie2 interface is dominated by nonpolar contributions, with only three Ang2 and two Tie2 residues that contribute electrostatically to intermolecular interactions. Individual interfacial residues contributed only moderately to binding, suggesting that engineering of this interface will require multiple mutations to reach major effects. Conversely, substitutions in substantially buried Ang2 residues were more prevalent in our experimental screen, reduced binding substantially, and are therefore more likely to have a deleterious effect that might contribute to oncogenesis. Computational analysis of additional RTK-ligand complexes, c-Kit-SCF and M-CSF-c-FMS, and comparison to previous YSD results, further show the utility of our combined methodology.

Molecular mechanism study of several inhibitors binding to BRD9 bromodomain based on molecular simulations

Bromodomain-containing protein 9 (BRD9) has been employed as a potential target for anticancer drugs in recent years. In this work, molecular docking, molecular dynamics (MD) simulations, binding free energy calculations, and per residue energy decomposition approaches were performed to elucidate the different binding modes between four pyridinone-like scaffold inhibitors and BRD9 bromodomain. Analysis results indicate that non-polar contribution mainly deriving from van der Waals energy is a critical impact on binding affinity of inhibitors against BRD9. Some key residues Phe44, Phe47, Val49, and Ile53 (at ZA loop) enhance the binding energy of inhibitors in BRD9 by means of providing hydrophobic interactions. Moreover, it is observed that BRD9 is anchored by the formation of a stable hydrogen bond between the carbonyl of the inhibitors and the residue Asn100 (at BC loop), and a strong π–π stacking interaction formed between the residue Tyr106 (at BC loop) and the inhibitors. The existence of dimethoxyphenyl structure and the aromatic ring merged to pyridinone scaffold are useful to enhance the BRD9 binding affinity. These findings should guide the rational design of more prospective inhibitors targeting BRD9.Communicated by Ramaswamy H. Sarma

Probing the “fingers” domain binding pocket of Hepatitis C virus NS5B RdRp and D559G resistance mutation via molecular docking, molecular dynamics simulation and binding free energy calculations

The NS5B RdRp polymerase is a prominent enzyme for the replication of Hepatitis C virus (HCV). During the HCV replication, the template RNA binding takes place in the “fingers” sub-domain of NS5B. The “fingers” domain is a new emerging allosteric site for the HCV drug development. The inhibitors of the “fingers” sub-domain adopt a new antiviral mechanism called RNA intervention. The details of essential amino acid residues, binding mode of the ligand, and the active site intermolecular interactions of RNA intervention reflect that this mechanism is ambiguous in the experimental study. To elucidate these details, we performed molecular docking analysis of the fingers domain inhibitor quercetagetin (QGN) with NS5B polymerase. The detailed analysis of QGN-NS5B intermolecular interactions was carried out and found that QGN interacts with the binding pocket amino acid residues Ala97, Ala140, Ile160, Phe162, Gly283, Gly557, and Asp559; and also forms π⋯π stacking interaction with Phe162 and hydrogen bonding interaction with Gly283. These are found to be the essential interactions for the RNA intervention mechanism. Among the strong hydrogen bonding interactions, the QGN⋯Ala140 is a newly identified important hydrogen bonding interaction by the present work and this interaction was not resolved by the previously reported crystal structure. Since D559G mutation at the fingers domain was reported for reducing the inhibition percentage of QGN to sevenfold, we carried out molecular dynamics (MD) simulation for wild and D559G mutated complexes to study the stability of protein conformation and intermolecular interactions. At the end of 50 ns MD simulation, the π⋯π stacking interaction of Phe162 with QGN found in the wild-type complex is altered into T-shaped π stacking interaction, which reduces the inhibition strength. The origin of the D559G resistance mutation was studied using combined MD simulation, binding free energy calculations and principal component analysis. The results were compared with the wild-type complex. The mutation D559G reduces the binding affinity of the QGN molecule to the fingers domain. The free energy decomposition analysis of each residue of wild-type and mutated complexes revealed that the loss of non-polar energy contribution is the origin of the resistance.Communicated by Ramaswamy H. Sarma

Fragment Molecular Orbital Calculations with Implicit Solvent Based on the Poisson-Boltzmann Equation: Implementation and DNA Study.

In this study, an ab initio fragment molecular orbital (FMO) methodology was developed to evaluate the solvent effects on electrostatic interactions, which make a significant contribution to the physical and chemical processes occurring in biological systems. Here, a fully polarizable solute consisting of the FMO electron density was electrostatically coupled with an implicit solvent based on the Poisson-Boltzmann (PB) equation; in addition, the nonpolar contributions empirically obtained from the molecular surface area (SA) were added. Interaction analysis considering solvent-screening and dispersion effects is now available as a powerful tool to determine the local stabilities inside solvated biomolecules. This methodology is applied to a deoxyribonucleic acid (DNA) duplex known as the Dickerson dodecamer. We found that excessively large electrostatic interactions inside the duplex are effectively damped by the screening, and the frontier molecular orbital energies are also successfully lowered. These observations indicate the stability of highly charged DNA duplexes in solution. Moreover, the solvation free energies in the implicit model show fairly good agreement with those in the explicit model while avoiding the costly statistical sampling of the electrolyte distribution. Consequently, our FMO-PBSA approach could yield new insights into biological phenomena and pharmacological problems via this ab initio methodology.

Breaking the polar-nonpolar division in solvation free energy prediction.

Implicit solvent models divide solvation free energies into polar and nonpolar additive contributions, whereas polar and nonpolar interactions are inseparable and nonadditive. We present a feature functional theory (FFT) framework to break this ad hoc division. The essential ideas of FFT are as follows: (i) representability assumption: there exists a microscopic feature vector that can uniquely characterize and distinguish one molecule from another; (ii) feature-function relationship assumption: the macroscopic features, including solvation free energy, of a molecule is a functional of microscopic feature vectors; and (iii) similarity assumption: molecules with similar microscopic features have similar macroscopic properties, such as solvation free energies. Based on these assumptions, solvation free energy prediction is carried out in the following protocol. First, we construct a molecular microscopic feature vector that is efficient in characterizing the solvation process using quantum mechanics and Poisson-Boltzmann theory. Microscopic feature vectors are combined with macroscopic features, that is, physical observable, to form extended feature vectors. Additionally, we partition a solvation dataset into queries according to molecular compositions. Moreover, for each target molecule, we adopt a machine learning algorithm for its nearest neighbor search, based on the selected microscopic feature vectors. Finally, from the extended feature vectors of obtained nearest neighbors, we construct a functional of solvation free energy, which is employed to predict the solvation free energy of the target molecule. The proposed FFT model has been extensively validated via a large dataset of 668 molecules. The leave-one-out test gives an optimal root-mean-square error (RMSE) of 1.05 kcal/mol. FFT predictions of SAMPL0, SAMPL1, SAMPL2, SAMPL3, and SAMPL4 challenge sets deliver the RMSEs of 0.61, 1.86, 1.64, 0.86, and 1.14 kcal/mol, respectively. Using a test set of 94 molecules and its associated training set, the present approach was carefully compared with a classic solvation model based on weighted solvent accessible surface area. © 2017 Wiley Periodicals, Inc.

Discovery of a ROCK inhibitor, FPND, which prevents cerebral hemorrhage through maintaining vascular integrity by interference with VE-cadherin

Hemorrhagic stroke occurs when a weakened vessel ruptures and bleeds into the surrounding brain, leading to high rates of death and disability worldwide. A series of complex pathophysiological cascades contribute to the risk of hemorrhagic stroke, and no therapies have proven effective to prevent hemorrhagic stroke. Stabilization of vascular integrity has been considered as a potential therapeutic target for hemorrhagic stroke. ROCKs, which belong to the serine/threonine protein kinase family and participate in the organization of actin cytoskeleton, have become attractive targets for the treatment of strokes. In this study, in vitro enzyme-based assays revealed that a new compound (FPND) with a novel scaffold identified by docking-based virtual screening could inhibit ROCK1 specifically at low micromolar concentration. Molecular modeling showed that FPND preferentially interacted with ROCK1, and the difference between the binding affinity of FPND toward ROCK1 and ROCK2 primarily resulted from non-polar contributions. Furthermore, FPND significantly prevented statin-induced cerebral hemorrhage in a zebrafish model. In addition, in vitro studies using the xCELLigence RTCA system, immunofluorescence and western blotting revealed that FPND prevented statin-induced cerebral hemorrhage by enhancing endothelial cell–cell junctions through inhibiting the ROCK-mediated VE-cadherin signaling pathway. As indicated by the extremely low toxicity of FPND against mice, it is safe and can potentially prevent vascular integrity loss-related diseases, such as hemorrhagic stroke.

Explicit Solvent Hydration Benchmark for Proteins with Application to the PBSA Method.

Explicit and implicit solvent models have a proven record of delivering hydration free energies of small, druglike solutes in reasonable agreement with experiment. Hydration of macromolecules, such as proteins, is to a large extent uncharted territory, with few results shedding light on quantitative consistency between different solvent models, let alone their ability to reproduce real water. In this work, based on extensive explicit solvent simulations employing TIP3P and SPC/E water models we analyze hydration free energy changes between fixed conformations of 5 diverse proteins, including large multidomain structures. For the two solvent models we find better agreement in electrostatic rather than nonpolar contributions (RMSE of 2.3 and 2.7 kcal/mol, respectively), even though absolute values of the latter are typically an order of magnitude smaller. We also highlight the importance of finite size corrections to relative protein hydration free energies, which turn out to be rather large, on the order of several kcal/mol, and are necessary for proper interpretation of results obtained under periodic boundary conditions. We further compare gathered data with predictions of the implicit solvent approach based on the Poisson equation and the surface or volume based nonpolar term. We find definitely lesser consistency than between the two explicit models (RMSE between implicit and TIP3 results of 11.3 and 8.4 kcal/mol for electrostatic and nonpolar contributions, respectively). In the process we determine the value of the protein dielectric constant and the geometric model for the dielectric boundary that provide for the best agreement. Finally, we evaluate the usefulness of surface and volume based models of nonpolar contributions to hydration free energy of large biomolecules.

Exploring the binding mechanism of Heteroaryldihydropyrimidines and Hepatitis B Virus capsid combined 3D-QSAR and molecular dynamics

The non-nucleoside drugs have been developed to treat HBV infection owing to their increased efficacy and lesser side effects, in which heteroaryldihydropyrimidines (HAPs) have been identified as effective inhibitors of HBV capsid. In this paper, the binding mechanism of HAPs targeting on HBV capsid protein was explored through three-dimensional quantitative structure-activity relationship, molecular dynamics and binding free energy decompositions. The obtained models of comparative molecular field analysis and comparative molecular similarity indices analysis enable the sufficient interpretation of structure-activity relationship of HAPs-HBV. The binding free energy analysis correlates with the experimental data. The computational results disclose that the non-polar contribution is the major driving force and Y132A mutation enhances the binding affinity for inhibitor 2 bound to HBV. The hydrogen bond interactions between the inhibitors and Trp102 help to stabilize the conformation of HAPs-HBV. The study provides insight into the binding mechanism of HAPs-HBV and would be useful for the rational design and modification of new lead compounds of HAP drugs.