LDH-Indomethacin Nanoparticles Antitumoral Action: A Possible Coadjuvant Drug for Cancer Therapy.

Indomethacin (INDO) has a mechanism of action based on inhibiting fatty acids cyclooxygenase activity within the inflammation process. The action mechanism could be correlated with possible anticancer activity, but its high toxicity in normal tissues has made therapy difficult. By the coprecipitation method, the drug carried in a layered double hydroxides (LDH) hybrid matrix would reduce its undesired effects by promoting chemotherapeutic redirection. Therefore, different samples containing INDO intercalated in LDH were synthesized at temperatures of 50, 70, and 90 °C and synthesis times of 8, 16, 24, and 48 h, seeking the best structural organization. X-ray diffraction (XRD), vibrational Fourier transform infrared spectroscopy (FT-IR), scanning electron microscopy (SEM), spectrophotometric analysis in UV-VIS, and differential thermogravimetric analysis (TGA/DTA) were used for characterization. Our results indicate that higher temperatures and longer synthesis time through coprecipitation reduce the possibility of INDO intercalation. However, it was possible to establish a time of 16 h and a temperature of 50 °C as the best conditions for intercalation. In vitro results confirmed the cell viability potential and anticancer activity in the LDH-INDO sample (16 h and 50 °C) for gastric cancer (AGP01, ACP02, and ACP03), breast cancer (MDA-MB-231 and MCF-7), melanoma (SK-MEL-19), lung fibroblast (MRC-5), and non-neoplastic gastric tissue (MN01) by MTT assay. Cell proliferation was inhibited, demonstrating higher and lower toxicity against MDA-MB-231 and SK-MEL-19. Thus, a clinical redirection of INDO is suggested as an integral and adjunctive anticancer medication in chemotherapy treatment.

Role of UDP-N-acetylmuramic acid in the regulation of MurA activity revealed by molecular dynamics simulations.

The peptidoglycan biosynthesis pathway plays a vital role in bacterial cells, and facilitates peptidoglycan layer formation, a fundamental structural component of the bacterial cell wall. The enzymes in this pathway are candidates for antibiotic development, as most do not have mammalian homologues. The UDP-N-acetylglucosamine (UNAG) enolpyruvyl transferase enzyme (MurA) in the peptidoglycan pathway cytoplasmic step is responsible for the phosphoenolpyruvate (PEP)-UNAG catalytic reaction, forming UNAG enolpyruvate and inorganic phosphate. Reportedly, UDP-N-acetylmuramic acid (UNAM) binds tightly to MurA forming a dormant UNAM-PEP-MurA complex and acting as a MurA feedback inhibitor. MurA inhibitors are complex, owing to competitive binding interactions with PEP, UNAM, and UNAG at the MurA active site. We used computational methods to explore UNAM and UNAG binding. UNAM showed stronger hydrogen-bond interactions with the Arg120 and Arg91 residues, which help to stabilize the closed conformation of MurA, than UNAG. Binding free energy calculations using end-point computational methods showed that UNAM has a higher binding affinity than UNAG, when PEP is attached to Cys115. The unbinding process, simulated using τ-random acceleration molecular dynamics, showed that UNAM has a longer relative residence time than UNAG, which is related to several complex dissociation pathways, each with multiple intermediate metastable states. This prevents the loop from opening and exposing the Arg120 residue to accommodate UNAG and potential new ligands. Moreover, we demonstrate the importance of Cys115-linked PEP in closed-state loop stabilization. We provide a basis for evaluating novel UNAM analogues as potential MurA inhibitors. PUBLIC SIGNIFICANCE: MurA is a critical enzyme involved in bacterial cell wall biosynthesis and is involved in antibiotic resistance development. UNAM can remain in the target protein's active site for an extended time compared to its natural substrate, UNAG. The prolonged interaction of this highly stable complex known as the 'dormant complex' comprises UNAM-PEP-MurA and offers insights into antibiotic development, providing potential options against drug-resistant bacteria and advancing our understanding of microbial biology.

Computational Analysis of Triazole-Based Kojic Acid Analogs as Tyrosinase Inhibitors by Molecular Dynamics and Free Energy Calculations.

Molecular docking, molecular dynamics (MD) simulations and the linear interaction energy (LIE) method were used here to predict binding modes and free energy for a set of 1,2,3-triazole-based KA analogs as potent inhibitors of Tyrosinase (TYR), a key metalloenzyme of the melanogenesis process. Initially, molecular docking calculations satisfactorily predicted the binding mode of evaluated KA analogs, where the KA part overlays the crystal conformation of the KA inhibitor into the catalytic site of TYR. The MD simulations were followed by the LIE method, which reproduced the experimental binding free energies for KA analogs with an r2 equal to 0.97, suggesting the robustness of our theoretical model. Moreover, the van der Waals contributions performed by some residues such as Phe197, Pro201, Arg209, Met215 and Val218 are responsible for the binding recognition of 1,2,3-triazole-based KA analogs in TYR catalytic site. Finally, our calculations provide suitable validation of the combination of molecular docking, MD, and LIE approaches as a powerful tool in the structure-based drug design of new and potent TYR inhibitors.

Computational Analysis of SAM Analogs as Methyltransferase Inhibitors of nsp16/nsp10 Complex from SARS-CoV-2.

Methyltransferases (MTases) enzymes, responsible for RNA capping into severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), are emerging important targets for the design of new anti-SARS-CoV-2 agents. Here, analogs of S-adenosylmethionine (SAM), obtained from the bioisosteric substitution of the sulfonium and amino acid groups, were evaluated by rigorous computational modeling techniques such as molecular dynamics (MD) simulations followed by relative binding free analysis against nsp16/nsp10 complex from SARS-CoV-2. The most potent inhibitor (2a) shows the lowest binding free energy (-58.75 Kcal/mol) and more potency than Sinefungin (SFG) (-39.8 Kcal/mol), a pan-MTase inhibitor, which agrees with experimental observations. Besides, our results suggest that the total binding free energy of each evaluated SAM analog is driven by van der Waals interactions which can explain their poor cell permeability, as observed in experimental essays. Overall, we provide a structural and energetic analysis for the inhibition of the nsp16/nsp10 complex involving the evaluated SAM analogs as potential inhibitors.

In Silico Study towards Repositioning of FDA-Approved Drug Candidates for Anticoronaviral Therapy: Molecular Docking, Molecular Dynamics and Binding Free Energy Calculations.

The SARS-CoV-2 targets were evaluated for a set of FDA-approved drugs using a combination of drug repositioning and rigorous computational modeling methodologies such as molecular docking and molecular dynamics (MD) simulations followed by binding free energy calculations. Six FDA-approved drugs including, Ouabain, Digitoxin, Digoxin, Proscillaridin, Salinomycin and Niclosamide with promising anti-SARS-CoV-2 activity were screened in silico against four SARS-CoV-2 proteins-papain-like protease (PLpro), RNA-dependent RNA polymerase (RdRp), SARS-CoV-2 main protease (Mpro), and adaptor-associated kinase 1 (AAK1)-in an attempt to define their promising targets. The applied computational techniques suggest that all the tested drugs exhibited excellent binding patterns with higher scores and stable complexes compared to the native protein cocrystallized inhibitors. Ouabain was suggested to act as a dual inhibitor for both PLpro and Mpro enzymes, while Digitoxin bonded perfectly to RdRp. In addition, Salinomycin targeted PLpro. Particularly, Niclosamide was found to target AAK1 with greater affinity compared to the reference drug. Our study provides comprehensive molecular-level insights for identifying or designing novel anti-COVID-19 drugs.

Metabolic Processing of Selenium-Based Bioisosteres of meso-Diaminopimelic Acid in Live Bacteria.

A primary component of all known bacterial cell walls is the peptidoglycan (PG) layer, which is composed of repeating units of sugars connected to short and unusual peptides. The various steps within PG biosynthesis are targets of potent antibiotics as proper assembly of the PG is essential for cellular growth and survival. Synthetic mimics of PG have proven to be indispensable tools to study the bacterial cell structure, growth, and remodeling. Yet, a common component of PG, meso-diaminopimelic acid (m-DAP) at the third position of the stem peptide, remains challenging to access synthetically and is not commercially available. Here, we describe the synthesis and metabolic processing of a selenium-based bioisostere of m-DAP (selenolanthionine) and show that it is installed within the PG of live bacteria by the native cell wall crosslinking machinery in mycobacterial species. This PG probe has an orthogonal release mechanism that could be important for downstream proteomics studies. Finally, we describe a bead-based assay that is compatible with high-throughput screening of cell wall enzymes. We envision that this probe will supplement the current methods available for investigating PG crosslinking in m-DAP-containing organisms.

Assessment of host-guest molecular encapsulation of eugenol using ß-cyclodextrin.

Eugenol is a natural compound with well-known repellent activity. However, its pharmaceutical and cosmetic applications are limited, since this compound is highly volatile and thermolabile. Nanoencapsulation provides protection, stability, conservation, and controlled release for several compounds. Here, eugenol was included in ß-cyclodextrin, and the complex was characterized through X-ray diffraction analysis (XRD) and Fourier-transform infrared spectroscopy (FTIR). Additionally, we used molecular dynamics simulations to explore the eugenol-ß-cyclodextrin complex stability with temperature increases. Our computational result demonstrates details of the molecular interactions and conformational changes of the eugenol-ß-cyclodextrin complex and explains its stability between temperatures 27°C and 48°C, allowing its use in formulations that are subjected to varied temperatures.

QM/MM Study of the Fosfomycin Resistance Mechanism Involving FosB Enzyme.

Multidrug-resistant organisms contain antibiotic-modifying enzymes that facilitate resistance to a variety of antimicrobial compounds. Particularly, the fosfomycin (FOF) drug can be structurally modified by several FOF-modifying enzymes before it reaches the biological target. Among them, FosB is an enzyme that utilizes l-cysteine or bacillithiol in the presence of a divalent metal to open the epoxide ring of FOF and, consequently, inactivate the drug. Here, we have used hybrid quantum mechanics/molecular mechanics (QM/MM) and molecular dynamics (MD) simulations to explore the mechanism of the reaction involving FosB and FOF. The calculated free-energy profiles show that the cost to open the epoxide ring of FOF at the C2 atom is ∼3.0 kcal/mol higher than that at the C1 atom. Besides, our QM/MM MD results revealed the critical role of conformation change of Cys9 and Asn50 to release the drug from the active site. Overall, the present study provides insights into the mechanism of FOF-resistant proteins.

Drug repurposing and computational modeling for discovery of inhibitors of the main protease (Mpro) of SARS-CoV-2.

The main protease (Mpro or 3CLpro) is a conserved cysteine protease from the coronaviruses and started to be considered an important drug target for developing antivirals, as it produced a deadly outbreak of COVID-19. Herein, we used a combination of drug reposition and computational modeling approaches including molecular docking, molecular dynamics (MD) simulations, and the calculated binding free energy to evaluate a set of drugs in complex with the Mpro enzyme. Particularly, our results show that darunavir and triptorelin drugs have favorable binding free energy (-63.70 and -77.28 kcal mol-1, respectively) in complex with the Mpro enzyme. Based on the results, the structural and energetic features that explain why some drugs can be repositioned to inhibit Mpro from SARS-CoV-2 were exposed. These features should be considered for the design of novel Mpro inhibitors.

Exploring the Catalytic Mechanism of the RNA Cap Modification by nsp16-nsp10 Complex of SARS-CoV-2 through a QM/MM Approach.

The inhibition of key enzymes that may contain the viral replication of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have assumed central importance in drug discovery projects. Nonstructural proteins (nsps) are essential for RNA capping and coronavirus replication since it protects the virus from host innate immune restriction. In particular, nonstructural protein 16 (nsp16) in complex with nsp10 is a Cap-0 binding enzyme. The heterodimer formed by nsp16-nsp10 methylates the 5'-end of virally encoded mRNAs to mimic cellular mRNAs and thus it is one of the enzymes that is a potential target for antiviral therapy. In this study, we have evaluated the mechanism of the 2'-O methylation of the viral mRNA cap using hybrid quantum mechanics/molecular mechanics (QM/MM) approach. It was found that the calculated free energy barriers obtained at M062X/6-31+G(d,p) is in agreement with experimental observations. Overall, we provide a detailed molecular analysis of the catalytic mechanism involving the 2'-O methylation of the viral mRNA cap and, as expected, the results demonstrate that the TS stabilization is critical for the catalysis.

Facile Synthesis and Metabolic Incorporation of m-DAP Bioisosteres Into Cell Walls of Live Bacteria.

Bacterial cell walls contain peptidoglycan (PG), a scaffold that provides proper rigidity to resist lysis from internal osmotic pressure and a barrier to protect cells against external stressors. It consists of repeating sugar units with a linkage to a stem peptide that becomes cross-linked by cell wall transpeptidases (TP). While synthetic PG fragments containing l-lysine in the third position on the stem peptide are easier to access, those with meso-diaminopimelic acid (m-DAP) pose a severe synthetic challenge. Herein, we describe a solid phase synthetic scheme based on widely available building blocks to assemble meso-cystine (m-CYT), which mimics key structural features of m-DAP. To demonstrate proper mimicry of m-DAP, cell wall probes were synthesized with m-CYT in place of m-DAP and evaluated for their metabolic processing in live bacterial cells. We found that m-CYT-based cell wall probes were properly processed by TPs in various bacterial species that endogenously contain m-DAP in their PG. Additionally, we have used hybrid quantum mechanical/molecular mechanical (QM/MM) and molecular dynamics (MD) simulations to explore the influence of m-DAP analogs on the PG cross-linking. The results showed that the cross-linking mechanism of transpeptidases occurred through a concerted process. We anticipate that this strategy, which is based on the use of inexpensive and commercially available building blocks, can be widely adopted to provide greater accessibility of PG mimics for m-DAP containing organisms.

Pellucidin A promotes antinociceptive activity by peripheral mechanisms inhibiting COX-2 and NOS: In vivo and in silico study.

Peperomia pellucida (PP) belongs to the Peperomia genus, which has a pantropic distribution. PP is used to treat a wide range of symptoms and diseases, such as pain, inflammation, and hypertension. Intriguingly, PP extract is used by different tropical countries for its anti-inflammatory and antinociceptive effects. In fact, these outcomes have been shown in animal models, though the exact bioactive products of PP that exert such results are yet to be discovered. To determine and elucidate the mechanism of action of one of these compounds, we evaluated the antinociceptive effect of the novel dimeric ArC2 compound, Pellucidin A by using in vivo and in silico models. Animals were then subjected to chemical, biphasic and thermal models of pain. Pellucidin A induced an antinociceptive effect against chemical-induced pain in mice, demonstrated by the decrease of the number of writhes, reaching a reduction of 43% and 65% in animals treated with 1 and 5 mg/kg of Pellucidin A, respectively. In the biphasic response (central and peripheral), animals treated with Pellucidin A showed a significant reduction of the licking time exclusively during the second phase (inflammatory phase). In the hot-plate test, Pellucidin A did not have any impact on the latency time of the treated animals. Moreover, in vivo and in silico results show that Pellucidin A's mechanism of action in the inflammatory pain occurs most likely through interaction with the nitric oxide (NO) pathway. Our results demonstrate that the antinociceptive activities of Pellucidin A operate under mechanism(s) of peripheral action, involving inflammatory mediators. This work provides insightful novel evidence of the biological properties of Pellucidin A, and leads to a better understanding of its mechanism of action, pointing to potential pharmacological use.

Evaluating the Performance of a Non-Bonded Cu2+ Model Including Jahn-Teller Effect into the Binding of Tyrosinase Inhibitors.

Tyrosinase (TYR) is a metalloenzyme classified as a type-3 copper protein, which is involved in the synthesis of melanin through a catalytic process beginning with the conversion of the amino acid l-Tyrosine (l-Tyr) to l-3,4-dihydroxyphenylalanine (l-DOPA). It plays an important role in the mechanism of melanogenesis in various organisms including mammals, plants, and fungi. Herein, we used a combination of computational molecular modeling techniques including molecular dynamic (MD) simulations and the linear interaction energy (LIE) model to evaluate the binding free energy of a set of analogs of kojic acid (KA) in complex with TYR. For the MD simulations, we used a dummy model including the description of the Jahn-Teller effect for Cu2+ ions in the active site of this enzyme. Our results show that the LIE model predicts the TYR binding affinities of the inhibitor in close agreement to experimental results. Overall, we demonstrate that the classical model provides a suitable description of the main interactions between analogs of KA and Cu2+ ions in the active site of TYR.

Assessment of the Cruzain Cysteine Protease Reversible and Irreversible Covalent Inhibition Mechanism.

Reversible and irreversible covalent ligands are advanced cysteine protease inhibitors in the drug development pipeline. K777 is an irreversible inhibitor of cruzain, a necessary enzyme for the survival of the Trypanosoma cruzi (T. cruzi) parasite, the causative agent of Chagas disease. Despite their importance, irreversible covalent inhibitors are still often avoided due to the risk of adverse effects. Herein, we replaced the K777 vinyl sulfone group with a nitrile moiety to obtain a reversible covalent inhibitor (Neq0682) of cysteine protease. Then, we used advanced experimental and computational techniques to explore details of the inhibition mechanism of cruzain by reversible and irreversible inhibitors. The isothermal titration calorimetry (ITC) analysis shows that inhibition of cruzain by an irreversible inhibitor is thermodynamically more favorable than by a reversible one. The hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) and Molecular Dynamics (MD) simulations were used to explore the mechanism of the reaction inhibition of cruzain by K777 and Neq0682. The calculated free energy profiles show that the Cys25 nucleophilic attack and His162 proton transfer occur in a single step for a reversible inhibitor and two steps for an irreversible covalent inhibitor. The hybrid QM/MM calculated free energies for the inhibition reaction correspond to -26.7 and -5.9 kcal mol-1 for K777 and Neq0682 at the MP2/MM level, respectively. These results indicate that the ΔG of the reaction is very negative for the process involving K777, consequently, the covalent adduct cannot revert to a noncovalent protein-ligand complex, and its binding tends to be irreversible. Overall, the present study provides insights into a covalent inhibition mechanism of cysteine proteases.

Exploring Chloride Selectivity and Halogenase Regioselectivity of the SalL Enzyme through Quantum Mechanical/Molecular Mechanical Modeling.

The catalytic mechanism of SalL chlorinase has been investigated by combining quantum mechanical/molecular mechanical (QM/MM) techniques and umbrella sampling simulations to compute free energy profiles. Our results shed light on the interesting fact that the substitution of chloride with fluorine in SalL chlorinase leads to a loss of halogenase activity. The potential of mean force based on DFTB3/MM analysis shows that fluorination corresponds to a barrier 13.5 kcal·mol-1 higher than chlorination. Additionally, our results present a molecular description of SalL acting as a chlorinase instead of a methyl-halide transferase.

Mycobacterium abscessus l,d-Transpeptidases Are Susceptible to Inactivation by Carbapenems and Cephalosporins but Not Penicillins.

As a growing number of clinical isolates of Mycobacterium abscessus are resistant to most antibiotics, new treatment options that are effective against these drug-resistant strains are desperately needed. The majority of the linkages in the cell wall peptidoglycan of M. abscessus are synthesized by nonclassical transpeptidases, namely, the l,d-transpeptidases. Emerging evidence suggests that these enzymes represent a new molecular vulnerability in this pathogen. Recent studies have demonstrated that inhibition of these enzymes by the carbapenem class of ß-lactams determines their activity against Mycobacterium tuberculosis Here, we studied the interactions of ß-lactams with two l,d-transpeptidases in M. abscessus, namely, LdtMab1 and LdtMab2, and found that both the carbapenem and cephalosporin, but not penicillin, subclasses of ß-lactams inhibit these enzymes. Contrary to the commonly held belief that combination therapy with ß-lactams is redundant, doripenem and cefdinir exhibit synergy against both pansusceptible M. abscessus and clinical isolates that are resistant to most antibiotics, which suggests that dual-ß-lactam therapy has potential for the treatment of M. abscessus Finally, we solved the first crystal structure of an M. abscessus l,d-transpeptidase, LdtMab2, and using substitutions of critical amino acids in the catalytic site and computational simulations, we describe the key molecular interactions between this enzyme and ß-lactams, which provide an insight into the molecular basis for the relative efficacy of different ß-lactams against M. abscessus.

Targeting the cell wall of Mycobacterium tuberculosis: a molecular modeling investigation of the interaction of imipenem and meropenem with L,D-transpeptidase 2.

The single crystal X-ray structure of the extracellular portion of the L,D-transpeptidase (ex-LdtMt2 - residues 120-408) enzyme was recently reported. It was observed that imipenem and meropenem inhibit activity of this enzyme, responsible for generating L,D-transpeptide linkages in the peptidoglycan layer of Mycobacterium tuberculosis. Imipenem is more active and isothermal titration calorimetry experiments revealed that meropenem is subjected to an entropy penalty upon binding to the enzyme. Herein, we report a molecular modeling approach to obtain a molecular view of the inhibitor/enzyme interactions. The average binding free energies for nine commercially available inhibitors were calculated using MM/GBSA and Solvation Interaction Energy (SIE) approaches and the calculated energies corresponded well with the available experimentally observed results. The method reproduces the same order of binding energies as experimentally observed for imipenem and meropenem. We have also demonstrated that SIE is a reasonably accurate and cost-effective free energy method, which can be used to predict carbapenem affinities for this enzyme. A theoretical explanation was offered for the experimental entropy penalty observed for meropenem, creating optimism that this computational model can serve as a potential computational model for other researchers in the field.

Simulating the inhibition reaction of Mycobacterium tuberculosis L,D-transpeptidase 2 by carbapenems.

A theoretical free energy study describes the inactivation of a new tuberculosis target, the l,d-transpeptidase 2 enzyme. A new reaction mechanism of two carbapenem inhibitors is proposed and their molecular features are determined using QM/MM and PMF approaches. The theoretical findings with the new proposed mechanism agree in principle with the experimental data.

Insights into the mechanism of oxidation of dihydroorotate to orotate catalysed by human class 2 dihydroorotate dehydrogenase: a QM/MM free energy study.

The dihydroorotate dehydrogenase (DHOD) enzyme catalyzes the unique redox reaction in the de novo pyrimidine biosynthesis pathway. In this reaction, the oxidation of dihydroorotate (DHO) to orotate (OA) and reduction of the flavin mononucleotide (FMN) cofactor is catalysed by DHOD. The class 2 DHOD, to which the human enzyme belongs, was experimentally shown to follow a stepwise mechanism but the data did not allow the determination of the order of bond-breaking in a stepwise oxidation of DHO. The goal of this study is to understand the reaction mechanism at the molecular level of class 2 DHOD, which may aid in the design of inhibitors that selectively impact the activity of only certain members of the enzyme family. In this paper, the catalytic mechanism of oxidation of DHO to OA in human DHOD was studied using a hybrid Quantum Mechanical/Molecular Mechanical (QM/MM) approach and Molecular Dynamics (MD) simulations. The free energy barriers calculated reveal that the mechanism in human DHOD occurs via a stepwise reaction pathway. In the first step, a proton is abstracted from the C5 of DHO to the deprotonated Ser215 side chain. Whereas, in the second step, the transfer of the hydride or hydride equivalent from the C6 of DHO to the N5 of FMN, where free energy barrier calculated by the DFT/MM level is 10.84 kcal mol(-1). Finally, a residual decomposition analysis was carried out in order to elucidate the influence of the catalytic region residues during DHO oxidation.

Pentacycloundecane lactam vs lactone norstatine type protease HIV inhibitors: binding energy calculations and DFT study.

BACKGROUND: Novel pentacycloundecane (PCU)-lactone-CO-EAIS peptide inhibitors were designed, synthesized, and evaluated against wild-type C-South African (C-SA) HIV-1 protease. Three compounds are reported herein, two of which displayed IC50 values of less than 1.00 µM. A comparative MM-PB(GB)SA binding free energy of solvation values of PCU-lactam and lactone models and their enantiomers as well as the PCU-lactam-NH-EAIS and lactone-CO-EAIS peptide inhibitors and their corresponding diastereomers complexed with South African HIV protease (C-SA) was performed. This will enable us to rationalize the considerable difference between inhibitory concentration (IC50) of PCU-lactam-NH-EAIS and PCU-lactone-CO-EAIS peptides. RESULTS: The PCU-lactam model exhibited more negative calculated binding free energies of solvation than the PCU-lactone model. The same trend was observed for the PCU-peptide inhibitors, which correspond to the experimental activities for the PCU-lactam-NH-EAIS peptide (IC50 = 0.076 µM) and the PCU-lactone-CO-EAIS peptide inhibitors (IC50 = 0.850 µM). Furthermore, a density functional theory (DFT) study on the natural atomic charges of the nitrogen and oxygen atoms of the three PCU-lactam, PCU-lactim and PCU-lactone models were performed using natural bond orbital (NBO) analysis. Electrostatic potential maps were also used to visualize the electron density around electron-rich regions. The asymmetry parameter (η) and quadrupole coupling constant (χ) values of the nitrogen and oxygen nuclei of the model compounds were calculated at the same level of theory. Electronic molecular properties including polarizability and electric dipole moments were also calculated and compared. The Gibbs theoretical free solvation energies of solvation (∆Gsolv) were also considered. CONCLUSIONS: A general trend is observed that the lactam species appears to have a larger negative charge distribution around the heteroatoms, larger quadrupole constant, dipole moment and better solvation energy, in comparison to the PCU-lactone model. It can be argued that these characteristics will ensure better eletronic interaction between the lactam and the receptor, corresponding to the observed HIV protease activities in terms of experimental IC50 data.