Radial Behavior of Electrostatic Potentials of Atoms and Ions Revisited: Isotropy and Anisotropy.

In this paper we revisit earlier work relating to monoatomic atoms and ions published by pioneers in the area of electrostatic potentials. We include plots of the radial distributions of the electrostatic potentials for spherically symmetric atoms and cations, and for singly, doubly and triply negative anions. For atoms with anisotropy in their densities and electrostatic potentials, such as the halonium cations, it is shown how the molecular surface approach for plotting electrostatic potentials complements that achieved by directional radial distributions.

Virtual screening, MMGBSA, and molecular dynamics approaches for identification of natural products from South African biodiversity as potential Onchocerca volvulus pi-class glutathione S-transferase inhibitors.

We investigated 1012 molecules from natural products previously isolated from the South African biodiversity (SANCDB, https://sancdb.rubi.ru.ac.za/), for putative inhibition of Onchocerca volvulus pi-class glutathione S-transferase (Ov-GST2) by virtual screening, MMGBSA, and molecular dynamics approaches. ADMET, docking, and MMGBSA shortlisted 12 selected homoisoflavanones-type hit molecules, among which two namely SANC00569, and SANC00689 displayed high binding affinities of -46.09 and -46.26 kcal mol-1, respectively towards π-class Ov-GST2, respectively. The molecular dynamics results of SANC00569 showed the presence of intermolecular H-bonding, hydrophobic interactions between the ligand and key amino acids of Ov-GST2, throughout the simulation period. This hit molecule had a stable binding pose and occupied the binding pockets throughout the 200 ns simulation. To the best of our knowledge, there is no report of any alleged anti-onchocerciasis activity referring to homoisoflavanones or flavonoids. Nevertheless, homoisoflavanones, which are a subclass of flavonoids, exhibit a plethora of biological activities. All these results led to the conclusion that SANC00569 is the most hypothetical Ov-GST2, which could lead the development of new drugs against Onchocerca volvulus pi-class glutathione S-transferase. Further validation of these findings through in vitro and in vivo studies is required.

Bovine serum albumin uptake and polypeptide disaggregation studies of hypoglycemic ruthenium(II) uracil Schiff-base complexes.

Our prior studies have illustrated that the uracil ruthenium(II) diimino complex, [Ru(H3ucp)Cl(PPh3)] (1) (H4ucp = 2,6-bis-((6-amino-1,3-dimethyluracilimino)methylene)pyridine) displayed high hypoglycemic effects in diet-induced diabetic rats. To rationalize the anti-diabetic effects of 1, three new derivatives have been prepared, cis-[Ru(bpy)2(urdp)]Cl2 (2) (urdp = 2,6-bis-((uracilimino)methylene)pyridine), trans-[RuCl2(PPh3)(urdp)] (3), and cis-[Ru(bpy)2(H4ucp)](PF6)2 (4). Various physicochemical techniques were utilized to characterize the structures of the novel ruthenium compounds. Prior to biomolecular interactions or in vitro studies, the stabilities of 1-4 were monitored in anhydrous DMSO, aqueous phosphate buffer containing 2% DMSO, and dichloromethane (DCM) via UV-Vis spectrophotometry. Time-dependent stability studies showed ligand exchange between DMSO nucleophiles and chloride co-ligands of 1 and 3, which was suppressed in the presence of an excess amount of chloride ions. In addition, the metal complexes 1 and 3 are stable in both DCM and an aqueous phosphate buffer containing 2% DMSO. In the case of compounds 2 and 4 with no chloride co-ligands within their coordination spheres, high stability in aqueous phosphate buffer containing 2% DMSO was observed. Fluorescence emission titrations of the individual ruthenium compounds with bovine serum albumin (BSA) showed that the metal compounds interact non-discriminately within the protein's hydrophobic cavities as moderate to strong binders. The metal complexes were capable of disintegrating mature amylin amyloid fibrils. In vivo glucose metabolism studies in liver (Chang) cell lines confirmed enhanced glucose metabolism as evidenced by the increased glucose utilization and glycogen synthesis in liver cell lines in the presence of complexes 2-4.

Anisotropies in electronic densities and electrostatic potentials of Halonium Ions: focus on Chlorine, Bromine and Iodine.

CONTEXT: Why are the halonium cations so effective in forming strongly-bound complexes? We directed our research to address this question and we present electrostatic potential data for the valence-state halogen atoms X and halonium cations X+, where X = Cl, Br, I. The electron densities and electrostatic potentials of the halonium cations show considerably greater anisotropy than do the valence state halogens. The distances from the electrostatic potential surface maxima to the halogen nuclei are about 0.5 Å smaller than the distances from the electrostatic potential surface minima to the nuclei, giving the halonium cations each a more disk-like shape than the corresponding neutral valence state halogens. Their surface electrostatic potentials are totally consistent with the directionalities of halonium cations in complexes and the strengths of their interactions. To add perspective to this brief report, we have included calculations of the isotropic cation K+ and noble gas Kr. METHODS: The calculations of the electrostatic potentials of the valence states of the halogen atoms Cl, Br and I and the halonium cations Cl+, Br+ and I+, as well as K+ and Kr, on 0.001 au contours of their electronic densities were carried out with Gaussian O9 and the Wave Function Analysis - Surface Analysis Suite (WFA-SAS) at the M06-2X/6-31 + G(d,p) and M06-2X/3-21G* levels.

Interaction between escitalopram and ibuprofen or paracetamol: DFT and molecular docking on the drug-drug interactions.

A large number of drugs are introduced each year to treat different diseases. Most of the time, patients suffer from more than one health problem which makes it necessary to take multiple drugs. When drugs are combined, the problem of drug-drug interaction becomes relevant. In this work, we studied the drug-drug interaction between escitalopram and ibuprofen or paracetamol using density functional theory and quantum theory of atoms in molecules. The results suggest that following the interactions, the activity of drugs changes according to site of interaction. Most reactive and most stable interactions would be preferable for the purpose of use. The in silico drug-likeness studies show that escitalopram and paracetamol couple is more bioavailable than escitalopram and ibuprofen couple. Moreover, in order to gain additional insights into the mentioned drugs' interactions, the drugs were docked separately and jointly against the potential targets for antidepressants and NSAIDs, namely 6HIS and 2PXX. The molecular docking results showed a potential improvement of the effectiveness of the drugs after combining by forming hydrogen bonds, hydrophobic contacts and π…π stacking.Communicated by Ramaswamy H. Sarma.

Synthesis, Structure, Hirshfeld Surface Analysis, Non-Covalent Interaction, and In Silico Studies of 4-Hydroxy-1-[(4-Nitrophenyl)Sulfonyl]Pyrrolidine-2-Carboxyllic Acid.

The new compound 4-hydroxy-1-[(4-nitrophenyl)sulfonyl]pyrrolidine-2-carboxyllic acid was obtained by the reaction of 4-hydroxyproline with 4-nitrobenzenesulfonyl chloride. The compound was characterized using single crystal X-ray diffraction studies. Spectroscopic methods including NMR, FTIR, ES-MS, and UV were employed for further structural analysis of the synthesized compound. The title compound was found to have crystallized in an orthorhombic crystal system with space group P212121. The S1-N1 bond length of 1.628 (2) Å was a strong indication of the formation of the title compound. The absence of characteristic downfield 1H NMR peak of pyrrolidine ring and the presence of S-N stretching vibration at 857.82 cm-1 on the FTIR are strong indications for the formation of the sulfonamide. The experimental study was complemented with computations at the B3LYP/6-311G + + (d,p) level of theory to gain more understanding of interactions in the compound at the molecular level. Noncovalent interaction, Hirsfeld surface analysis and interaction energy calculations were employed in the analysis of the supramolecular architecture of the compound. Predicted ADMET parameters, awarded suitable bioavailability credentials, while the molecular docking study indicated that the compound enchants promising inhibition prospects against dihydropteroate synthase, DNA topoisomerase, and SARS-CoV-2 spike. Graphical Abstract: Herein we present the solid state structure, noncovalent interaction and spectroscopic analysis of a prospective bioactive compound 4-hydroxy-1-[(4-nitrophenyl)sulphonyl]pyrrolidine-2-carboxyllic acid. Supplementary Information: The online version contains supplementary material available at 10.1007/s10870-023-00978-0.

Computational study of the interaction of the psychoactive amphetamine with 1,2-indanedione and 1,8-diazafluoren-9-one as fingerprinting reagents.

In this study, we used computational methods to investigate the interaction of amphetamine (AMP) with 1,2-indanedione (IND) and 1,8-diazafluoren-9-one (DFO) so as to understand whether AMP can be detected in latent fingerprints using either of these reagents. The results show that the binding energies of AMP with IND and DFO were enhanced by the presence of amino acid from -9.29 to -12.35 kcal mol-1 and -7.98 to -10.65 kcal mol-1, respectively. The physical origins of these interactions could be better understood by symmetry-adapted perturbation theory. The excited state properties of the binding structures with IND demonstrate distinguishable absorption peaks in the UV-vis spectra but zero fluorescence. Furthermore, the UV-vis spectra of the possible reaction products between AMP and the reagents reveal absorption peaks in the visible spectrum. Therefore, we could predict that reaction of AMP with IND would be observable by a reddish colour while with DFO, a colour change to violet is expected. To conclude, the reagents IND and DFO may be used to detect AMP by UV-vis spectroscopy and if their reactions are allowed, the reagents may then act as a potentially rapid, affordable and easy colorimetric test for AMP in latent fingerprints without destruction of the fingerprint sample.

A study of the thermodynamics and mechanisms of the atmospherically relevant reaction dimethyl sulphide (DMS) with atomic chlorine (Cl) in the absence and presence of water, using electronic structure methods.

The thermodynamics and mechanisms of the atmospherically relevant reaction dimethyl sulphide (DMS) + atomic chlorine (Cl) were investigated in the absence and presence of a single water molecule, using electronic structure methods. Stationary points on each reaction surface were located using density functional theory (DFT) with the M06-2X functional with aug-cc-pVDZ (aVDZ) and aug-cc-pVTZ (aVTZ) basis sets. Then fixed point calculations were carried out using the UM06-2X/aVTZ optimised stationary point geometries, with aug-cc-pVnZ basis sets (n = T and Q), using the coupled cluster method [CCSD(T)], as well as the domain-based local pair natural orbitals coupled cluster [DLPNO-UCCSD(T)] approach. Four reaction channels are possible, formation of (A) CH3SCH2 + HCl, (B) CH3S + CH3Cl, (C) CH3SCl + CH3, and (C') CH3S(Cl)CH3. The results show that, in the absence of water, channels A and C' are the dominant channels. In the presence of water, the calculations show that the reaction mechanisms for A and C formation change significantly. Channel A occurs via submerged TSs and is expected to be rapid. Channel B occurs via TSs which present significant energy barriers indicating that this channel is not significant in the presence of water relative to CH3SCH2 + HCl and DMS·Cl adduct formation, as is the case in the absence of water. Channel C was not considered as it is endothermic in the absence of water. In the presence of water, pathways which proceed via (a) DMS·H2O + Cl, (b) Cl·H2O + DMS and (c) DMS·Cl + H2O were considered. It was found that under tropospheric conditions, reactions via pathway (b) are of minor importance relative to those that proceed via pathways (a) and (c). This study has shown that water changes the mechanisms of the DMS + Cl reactions significantly but the presence of water is not expected to affect the overall reaction rate coefficient under atmospheric conditions as the DMS + Cl reaction has a rate coefficient at room temperature close to the collisional limit.

Theoretical insights of the drug-drug interaction between favipiravir and ibuprofen: a DFT, QTAIM and drug-likeness investigation.

One of the approaches adopted, as the first treatment method, in previous pandemic viruses before the production of new vaccines or drugs, is to evaluate the effectiveness of existing drugs. Sometimes, it is known that the combination of drugs provides better treatment by increasing the activity of the drugs due to mutual interactions. In this work, we used quantum theory of atoms in molecules and density functional theory to examine the structure and activity relationship between Favipiravir (FPV), a potential antiviral drug for ribonucleic acid viruses, and an inhibitor of influenza virus, and Ibuprofen (IBP), a nonsteroidal anti-inflammatory drug. The drug-likeness evaluations of IBP…FPV were also performed. It was observed that FPV and IBP interact in several ways via hydrogen bonding (HB) leading to changes in the activities of the combined drug and IBP…FPV is predicted to be orally bioavailable. It was also found that the highest occupied molecular orbital (HOMO) is located on IBP and the lowest unoccupied molecular orbital (LUMO) is distributed over FPV. The biological activity of the combined drug system is increased because of the small HOMO-LUMO energy gap value, partially covalent and non-covalent HB interactions.Communicated by Ramaswamy H. Sarma.

Structure and Computational Studies of New Sulfonamide Compound: {(4-nitrophenyl)sulfonyl}tryptophan.

Synthesis of sulfonamide through an indirect method that avoids contamination of the product with no need for purification has been carried out using the indirect process. Here, we report the synthesis of a novel sulfonamide compound, ({4-nitrophenyl}sulfonyl)tryptophan (DNSPA) from 4-nitrobenzenesulphonylchloride and L-tryptophan precursors. The slow evaporation method was used to form single crystals of the named compound from methanolic solution. The compound was characterized by X-ray crystallographic analysis and spectroscopic methods (NMR, IR, mass spectrometry, and UV-vis). The sulfonamide N-H NMR signal at 8.07-8.09 ppm and S-N stretching vibration at 931 cm-1 indicate the formation of the target compound. The compound crystallized in the monoclinic crystal system and P21 space group with four molecules of the compound in the asymmetric unit. Molecular aggregation in the crystal structure revealed a 12-molecule aggregate synthon sustained by O-H⋯O hydrogen bonds and stabilised by N-H⋯O intermolecular contacts. Experimental studies were complemented by DFT calculations at the B3LYP/6-311++G(d,p) level of theory. The computed structural and spectroscopic data are in good agreement with those obtained experimentally. The energies of interactions between the units making up the molecule were calculated. Molecular docking studies showed that DNSPA has a binding energy of -6.37 kcal/mol for E. coli DNA gyrase (5MMN) and -6.35 kcal/mol for COVID-19 main protease (6LU7).

Unveiling the non-polar [3+2] cycloaddition reactions of cyclic nitrones with strained alkylidene cyclopropanes within a molecular electron density theory study.

The role of cyclopropane substitution on the ethylene in zw-type [3+2] cycloaddition (32CA) reactions of cyclic nitrones has been studied within Molecular Electron Density Theory (MEDT) at the ωB97X-D/6-311G(d,p) computational level. Electron Localization Function (ELF) analysis of the ethylenes shows that the presence the cyclopropane only slightly increases the electron density in the C-C bonding region. Analysis of the Conceptual DFT reactivity indices indicates that the presence of the cyclopropane does not produce any remarkable change in the reactivity of these strained ethylenes. The marginal electrophilic character of ethylene makes the zw-type 32CA reactions of non-polar character. The presence of the cyclopropane in the ethylene decreases the activation enthalpy of the 32CA reactions by only 1.7 and 2.6 kcal mol-1, and also decreases the ortho regioselectivity. The loss of the strain present in the cyclopropane is responsible for the reduction of the activation enthalpy and the increase of the reaction enthalpy in these non-polar 32CA reactions. The presence of the cyclopropane does not cause any change, neither in the transition state structure (TS) geometries nor in their electronic structure. The very low global electron density transfer (GEDT) computed at the TSs confirms the non-polar character of these 32CA reactions. The ortho regioselectivity experimentally observed in these non-polar 32CA reactions is determined by the most favorable two-center interaction between the less electronegative C1 carbon of nitrone and the non-substituted methylene C5 carbon of the ethylenes.

Structural, vibrational and electronic properties of some tetrel-bonded complexes of the fluorinated methanes methyl fluoride, difluoromethane and fluoroform: an ab initio study.

A search has been conducted, by means of ab initio molecular orbital theory, for potential tetrel-bonded complexes formed between the fluorinated methanes methyl fluoride, difluoromethane and fluoroform, and the related hydrides ammonia, water, hydrogen fluoride, phosphine, hydrogen sulphide and hydrogen chloride. Eleven such complexes have been identified, six containing CH3F and five CH2F2. The complexes are typically less strongly bound than their hydrogen-bonded counterparts, and the interaction energies vary in a consistent way with the periodic trend of the electron donors. The intermolecular separations and changes of the relevant intramolecular bond lengths, the wavenumber shifts of the critical vibrational modes and the extents of charge transfer correlate, by and large, with the strengths of interaction.

Isolation, crystal structure and DFT study of 2,2,6,6-tetramethyl-4-oxo-piperidinium nitrate isolated from the stembark of Vitex doniana (Lamiaceae).

We present the results of the analysis of 2,2,6,6-tetramethyl-4-oxo-piperidinium nitrate isolated from the stem bark of Vitex doniana, a tree growing in Nigeria. The low-temperature molecular structure comprises one 2,2,6,6-tetramethyl-4-oxopiperidinium cation and one nitrate anion as one molecule in the asymmetric unit. The compound crystallizes in the monoclinic space group P21/n. A portion of the nitrate anion exhibits positional disorder with the main disorder component present 66.253(2) % of the time and the minor disorder component present 33.279(2) % of the time. In comparison with the previously reported room-temperature structure of C9H18N2O4 . The low-temperature structure shows similarity with the piperidinium ring adopting a slightly deformed chair conformation while the nitrate anion is disordered. DFT method was used to complement the experimental study.

Physicochemical analysis of wastewater generated from a coating industry in Mauritius.

The coating industry is one of the most important consumers of water and chemicals and consequently is a major water polluter in Mauritius. The focus of this study was to characterise wastewater generated by a coating industry in Mauritius. The objectives were to develop a wastewater sampling strategy and to analyse the pollutant parameters as per Mauritian regulations. The wastewater samples were analysed for physicochemical properties and metal abundances over a period of 6 months. The physicochemical parameters analysed were pH, electrical conductivity (EC), true colour, total suspended solids (TSS), biological oxygen demand (BOD5), chemical oxygen demand (COD), nitrate, phosphate, sulphate and free chlorine. The wastewater samples were also analysed for metal ions such as sodium, potassium, arsenic, cadmium, chromium, cobalt, copper, iron, lead, manganese, mercury, molybdenum, nickel and zinc. The results of the physicochemical parameters indicated the presence of biologically resistant organic matters in all the wastewater samples with elevated values of BOD5 and COD, and low biodegradability index, respectively. The coating industry wastewater samples were acidic and saline in nature. Moreover, they presented high concentrations of TSS, free chlorine and sodium ions compared to standard limits promulgated by the Mauritian Government. Spearman's rank correlation matrix with non-linear regression analysis showed significant associations among the measured parameters which were found to have a common origin in the coating industry wastewater. This research will be useful for regular monitoring and setting up an adequate coating industry wastewater treatment for the potential reuse in production processes in Mauritius.

Decoding the reaction mechanism of the cyclocondensation of ethyl acetate2-oxo-2-(4-oxo-4H-pyrido [1.2-a] pyrimidin-3-yl) polyazaheterocycle and ethylenediamine using bond evolution theory.

We investigated the flow of electron density along the cyclocondensation reaction between ethyl acetate 2-oxo-2-(4-oxo-4H-pyrido[1.2-a]pyrimidin-3-yl) polyazaheterocycle (1) and ethylenediamine (2) at the ωB97XD/6-311++G(d,p)computational method within of bond evolution theory (BET). The exploration of potential energy surface shows that this reaction has three channels (1-3) with the formation of product 3 via channel-2 (the most favorable one) as the main product and this is in good agreement with experimental observations. The BET analysis allows identifying unambiguously the main chemical events happening along channel-2. The mechanism along first step (TS2-a) is described by a series of four structural stability domains (SSDs), while five SSDs for the last two steps (TS2-b and TS2-c). The first and third steps can be summarized as follows, the formation of N1-C6 bond (SSD-II), then, the restoration of the nitrogen N1 lone pair (SSD-III), and finally, the formation of the last O1-H1 bond (SSD-IV). For the second step, the formation of hydroxide ion is noted, as a result of the disappearance of V(C6,O7) basin and the transformation of C6-N1 single bond into double one (SSD-IV). Finally, the appearance of V(O7,H2) basin lead to the elimination of water molecule within the last domain is observed.

Theoretical study of a derivative of chlorophosphine with aliphatic and aromatic Grignard reagents: SN2@P or the novel SN2@Cl followed by SN2@C?

The proposed SN2 reactions of a hindered organophosphorus reactant with aliphatic and aromatic nucleophiles [Ye et al., Org. Lett., 2017, 19, 5384-5387] were studied theoretically in order to explain the observed stereochemistry of the products. Our computations (using B3LYP as the functional) indicate that the reaction with the aliphatic nucleophile occurs through a backside SN2@P pathway while the reaction with the aromatic nucleophile proceeds through a novel SN2@Cl mechanism, followed by a frontside SN2@C mechanism.

Can the Antivirals Remdesivir and Favipiravir Work Better Jointly? In Silico Insights.

The proverb "Old is gold" is applicable in drug discovery and the proverb "All that Glitters is not Gold" is also appropriate. In the COVID-19 era, there has been a race for drugs to be effective against SARS-CoV-2. There are reports about the uses of Remdesivir and Favipiravir as existing antivirals against virus but none have been conclusive so far. In the attempts for innovations, the combination of drugs is also under trials. Therefore, we used the density functional theory method and quantum theory of atoms in molecules to investigate drug-drug interactions involving Remdesivir and Favipiravir. The computed parameters were related to the antiviral actions of both drugs together. The results indicate enhanced antiviral activity and it will be worthy to consider additional investigations with the combination of these two drugs.

Virtual screening, ADMET profiling, PASS prediction, and bioactivity studies of potential inhibitory roles of alkaloids, phytosterols, and flavonoids against COVID-19 main protease (Mpro).

The current research used a virtual screening method to study 57 isolated phytochemicals (alkaloids, phytosterols, and flavonoids) against the SARS-CoV-2 main protease (Mpro). The absorption, distribution, metabolism, excretion, and toxicity (ADMET) of the selected compounds were analysed using admetSAR tool while SwissADME and Molinspiration chemoinformatics tools were used to examine the oral bioavailability and drug-likeness properties. Parameters such as physicochemical properties, activity spectra for substances (PASS) prediction, bioactivity, binding mode, and molecular interactions were also analysed. Our results favoured Lupeol (-8.6 kcal/mol), Lupenone (-7.7 kcal/mol), Hesperetin (-7.4 kcal/mol), Apigenin (-7.3 kcal/mol) and Castasterone (-7.3 kcal/mol) as probable inhibitors of SARS-CoV-2. This is because of their good binding affinities, bioactivities, drug-likeness, ADMET properties, PASS properties, oral bioavailability, binding mode and their interactions with the active site of the target receptor compared to Remdesivir and Azithromycin. Therefore, these compounds could be explored towards the development of new therapeutic agents against SARS-CoV-2.

Interaction between favipiravir and hydroxychloroquine and their combined drug assessment: in silico investigations.

Hydroxychloroquine (HCQ) and favipiravir (FPV) are known to be effective antivirals, and there are reports about their use to fight the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) despite that these are not conclusive. The use of combined drugs is common in drug discovery, and thus, we investigated HCQ and FPV as a combined drug. The density functional theory method was used for the optimization of geometries, spectroscopic analysis and calculation of reactivity parameters. The quantum theory of atoms in molecules was applied to explain the nature of the hydrogen bonds and confirm the higher stability of the combined drug. We also evaluated the absorption, distribution, metabolism and excretion (ADME) parameters to assess their drug actions jointly using SwissADME. The preliminary findings of our theoretical study are promising for further investigations of more potent and selective antiviral drugs. Supplementary Information: The online version contains supplementary material available at 10.1007/s11696-021-01946-8.

Theoretical study of the interactions between peptide tyrosine tyrosine [PYY (1-36)], a newly identified modulator in type 2 diabetes pathophysiology, with receptors NPY1R and NPY4R.

PURPOSE: Diabetes mellitus is a common condition in the clinically obese. Bariatric surgery is one of the ways to put type 2 diabetes in remission. Recent findings propose the appetite-regulator peptide tyrosine tyrosine (PYY) as a therapeutic option for patients with type 2 diabetes. This novel gut hormone restores impaired insulin and glucagon secretion in pancreatic islets and is implicated in type 2 diabetes reversal after bariatric surgery. The current study elucidates the interactions between PYY and the NPY1R and NPY4R receptors using computational methods. METHODS: Protein structure prediction, molecular docking simulation, and molecular dynamics (MD) simulation were performed to elucidate the interactions of PYY with NPY1R and NPY4R. RESULTS: The predicted binding models of PYY-NPY receptors are in agreement with those described in the literature, although different interaction partners are presented for the C-terminal tail of PYY. Non-polar interactions are predicted to drive the formation of the protein complex. The calculated binding energies show that PYY has higher affinity for NPY4R (ΔGGBSA = -65.08 and ΔGPBSA = -87.62 kcal/mol) than for NPY1R (ΔGGBSA = -23.11 and ΔGPBSA = -50.56 kcal/mol). CONCLUSIONS: Based on the constructed models, the binding conformations obtained from docking and MD simulation for both the PYY-NPY1R and PYY-NPY4R complexes provide a detailed map of possible interactions. The calculated binding energies show a higher affinity of PYY for NPY4R. These findings may help to understand the mechanisms behind the improvement of diabetes following bariatric surgery.