Synthesis and structural studies of pentacycloundecane-based HIV-1 PR inhibitors: A hybrid 2D NMR and docking/QM/MM/MD approach

Molecular Acrobatics in Polycyclic Frames: Synthesis of "Kurmanediol" via Post-synthetic Modification of Cage Molecules by Olefinic Metathesis

We report late-stage ring-opening metathesis (ROM), ring-rearrangement metathesis (RRM), and ring-closing metathesis (RCM) approaches to generate expanded pentacycloundecane (PCUD) cage derivatives. These higher-order intricate polycyclic cage systems are aesthetically pleasing and structurally intriguing. Their assembly maintains molecular symmetry during the entire synthetic sequence. To this end, metathesis-based catalysts are used to execute the ROM, RRM, and RCM strategies. The synthetic approach to these cage polycycles involves the Diels-Alder reaction, [2 + 2] photocycloaddition, RRM, ROM, and RCM as key steps.

In-silico studies on wild orange (Citrus macroptera Mont.) compounds against COVID-19 pro-inflammation targets

One-fifth of COVID-19 patients suffer a severe course of COVID-19 (SARS-CoV-2) infection; however, the specific causes remain unclear. Despite numerous papers that have been flooded in different scientific journals clear clinical picture of COVID-19 aftermath persists to remain fuzzy. The survivors of severe COVID-19infection having defeated the virus are just the starting of an uncharted recovery path. Currently, there is no drug available that is safe to consume to combat this pandemic. However, researchers still struggling to find specific therapeutic solutions. The present study employed an in silico approach to assessing the inhibitory potential of the phytochemicals obtained from GC-MS analysis of Citrus macroptera against inflammatory proteins like COX-2, NMDAR and VCAM-1 which remains in a hyperactive state even after a patient is fully cured of this deadly mRNA virus. An extensive molecular docking investigation of the phyto-compounds at the active binding pockets of the inflammatory proteins revealed the promising inhibitory potential of the phytochemicals. Reasonable physicochemical attributes of the compounds following Lipinski's rule of five, VEBER and PAINS analysis further established them as potential therapeutic candidates against aforesaid inflammatory proteins. MM-GBSA binding free energy estimation revealed that Limonene was the most promising candidate displaying the highest binding efficacy with the concerned VCAM-1 protein included in the present analysis. An interesting finding is the phytochemicals exhibited better binding energy scores with the concerned COX-2, VCAM-1 and NMDA receptor proteins than the conventional drugs that are specifically targeted against them. Our in silico results suggest that all the natural phyto-compounds derived from C. macroptera could be employed in Post covid inflammation complexities after appropriate pre-clinical and clinical trials for further scientific validation.Communicated by Ramaswamy H. Sarma.

Understanding Drug Resistance of Wild-Type and L38HL Insertion Mutant of HIV-1 C Protease to Saquinavir

Acquired immunodeficiency syndrome (AIDS) is one of the most challenging infectious diseases to treat on a global scale. Understanding the mechanisms underlying the development of drug resistance is necessary for novel therapeutics. HIV subtype C is known to harbor mutations at critical positions of HIV aspartic protease compared to HIV subtype B, which affects the binding affinity. Recently, a novel double-insertion mutation at codon 38 (L38HL) was characterized in HIV subtype C protease, whose effects on the interaction with protease inhibitors are hitherto unknown. In this study, the potential of L38HL double-insertion in HIV subtype C protease to induce a drug resistance phenotype towards the protease inhibitor, Saquinavir (SQV), was probed using various computational techniques, such as molecular dynamics simulations, binding free energy calculations, local conformational changes and principal component analysis. The results indicate that the L38HL mutation exhibits an increase in flexibility at the hinge and flap regions with a decrease in the binding affinity of SQV in comparison with wild-type HIV protease C. Further, we observed a wide opening at the binding site in the L38HL variant due to an alteration in flap dynamics, leading to a decrease in interactions with the binding site of the mutant protease. It is supported by an altered direction of motion of flap residues in the L38HL variant compared with the wild-type. These results provide deep insights into understanding the potential drug resistance phenotype in infected individuals.

Improved Production of HIV-1 Subtype C Protease from Transgenic E. Coli

Background:

Human Immunodeficiency Virus 1 (HIV-1) subtype C is responsible for the majority of infections of patients in Southern Africa. The HIV protease is a primary target for the development of highly efficient anti-retroviral pharmaceuticals because of its pivotal role in the maturation of the virus in the host cell. For target validation of novel HIV protease inhibitors, there is a need for the availability of an abundance of this protease.

Objective:

This study reports an optimized method to produce HIV-1 protease derived from HIV-1 subtype C.

Methods:

It involves the use of a transgenic E. coli strain that overexpresses the native form of the enzyme via inclusion bodies. A stringent method for the isolation, purification, and renaturation resulted in the production of highly pure active HIV-1 protease. In order to facilitate an increase in protease yields, an optimized growth strategy was developed. In this regard, a chemically defined medium with lower glucose content and devoid of essential amino acids of the TCA cycle was used as an alternative to the widely used nutrient-rich Luria Bertani (LB) medium.

Results:

Results indicated an increase in protease yield up to twice the amount, thereby making this medium an attractive alternative for increasing biomass and HIV protease production for future research.

Conclusion:

An optimized method for HIV-1 protease derived from HIV-1 subtype C production using chemically defined media was established. This was achieved using a known method to isolate and purify the enzyme with the use of a specialized feeding strategy.

Exploration of the effects of sequence variations between HIV-1 and HIV-2 proteases on their three-dimensional structures

HIV protease inhibitors (PIs) approved by the FDA (US Food and Drug Administration) are a major class of antiretroviral. HIV-2 protease (PR2) is naturally resistant to most of them as PIs were designed for HIV-1 protease (PR1). In this study, we explored the impact of amino-acid substitutions between PR1 and PR2 on the structure of protease (PR) by comparing the structural variability of 13 regions using 24 PR1 and PR2 structures complexed with diverse ligands. Our analyses confirmed structural rigidity of the catalytic region and highlighted the important role of three regions in the conservation of the catalytic region conformation. Surprisingly, we showed that the flap region, corresponding to a flexible region, exhibits similar conformations in PR1 and PR2. Furthermore, we identified regions exhibiting different conformations in PR1 and PR2, which could be explained by the intrinsic flexibility of these regions, by crystal packing, or by PR1 and PR2 substitutions. Some substitutions induce structural changes in the R2 and R4 regions that could have an impact on the properties of PI-binding site and could thus modify PI binding mode. Substitutions involved in structural changes in the elbow region could alter the flexibility of the PR2 flap regions relative to PR1, and thus play a role in the transition from the semi-open form to the closed form, and have an impact on ligand binding. These results improve the understanding of the impact of sequence variations between PR1 and PR2 on the natural resistance of HIV-2 to commercially available PIs.Communicated by Ramaswamy H. Sarma.

Exploring molecular mechanism of allosteric inhibitor to relieve drug resistance of multiple mutations in HIV-1 protease by enhanced conformational sampling

Recently, allosteric regulations of HIV-1 protease (PR) are suggested as a promising approach to relieve drug resistance of mutations toward inhibitors targeting the active site of PR. Replica-exchange molecular dynamics (REMD) simulations and normal mode analysis (NMA) are integrated to enhance conformational sampling of PR. Molecular mechanics generalized Born surface area (MM-GBSA) method was applied to calculate binding free energies of three inhibitors APV, DRV, and NIT to the wild-type (WT) and multidrug resistance (MDR) PRs. The results suggest that binding free energies of APV and DRV are decreased in the MDR PR relative to the WT PR, suggesting drug resistance of mutations on these two inhibitors. However, the binding ability of the allosteric inhibitor NIT is not impaired in the MDR PR. In addition, internal dynamics analysis based on REMD simulations proves that mutations hardly produce obvious effect on the conformation of the MDR PR in comparison to the WT PR. Scanning of hydrophobic contacts and hydrogen bond contacts of inhibitors with residues of PRs on the concatenated trajectories of REMD demonstrates that mutations change the symmetric interaction networks of APV and DRV with PR, but do not generate obvious influence on the asymmetric interaction network of NIT with PR. In summary, allosteric inhibitor NIT can adapt the MDR PR better than those inhibitors toward the active site of PR, thus allosteric inhibitors of PR may be a possible channel to overcome drug resistance of PR.

Ru(ii)-Catalyzed C6-selective C–H amidation of 2-pyridones

An efficient Ru-catalyzed C6 site-selective amidation of 2-pyridones has been accomplished with dioxazolone under mild conditions.

Quantum mechanics implementation in drug-design workflows: does it really help?

The pharmaceutical industry is progressively operating in an era where development costs are constantly under pressure, higher percentages of drugs are demanded, and the drug-discovery process is a trial-and-error run. The profit that flows in with the discovery of new drugs has always been the motivation for the industry to keep up the pace and keep abreast with the endless demand for medicines. The process of finding a molecule that binds to the target protein using in silico tools has made computational chemistry a valuable tool in drug discovery in both academic research and pharmaceutical industry. However, the complexity of many protein-ligand interactions challenges the accuracy and efficiency of the commonly used empirical methods. The usefulness of quantum mechanics (QM) in drug-protein interaction cannot be overemphasized; however, this approach has little significance in some empirical methods. In this review, we discuss recent developments in, and application of, QM to medically relevant biomolecules. We critically discuss the different types of QM-based methods and their proposed application to incorporating them into drug-design and -discovery workflows while trying to answer a critical question: are QM-based methods of real help in drug-design and -discovery research and industry?

Binding Free Energy Calculations of Nine FDA-approved Protease Inhibitors Against HIV-1 Subtype C I36T↑T Containing 100 Amino Acids Per Monomer

In this work, have investigated the binding affinities of nine FDA-approved protease inhibitor drugs against a new HIV-1 subtype C mutated protease, I36T↑T. Without an X-ray crystal structure, homology modelling was used to generate a three-dimensional model of the protease. This and the inhibitor models were employed to generate the inhibitor/I36T↑T complexes, with the relative positions of the inhibitors being superimposed and aligned using the X-ray crystal structures of the inhibitors/HIV-1 subtype B complexes as a reference. Molecular dynamics simulations were carried out on the complexes to calculate the average binding free energies for each inhibitor using the molecular mechanics generalized Born surface area (MM-GBSA) method. When compared to the binding free energies of the HIV-1 subtype B and subtype C proteases (calculated previously by our group using the same method), it was clear that the I36T↑T proteases mutations and insertion had a significant negative effect on the binding energies of the non-pepditic inhibitors nelfinavir, darunavir and tipranavir. On the other hand, ritonavir, amprenavir and indinavir show improved calculated binding energies in comparison with the corresponding data for wild-type C-SA protease. The computational model used in this study can be used to investigate new mutations of the HIV protease and help in establishing effective HIV drug regimes and may also aid in future protease drug design.

Non-B HIV-1 subtypes in sub-Saharan Africa: impact of subtype on protease inhibitor efficacy

In 2012, 25 million people [71% of global human immunodeficiency virus (HIV) infection] were estimated to be living with HIV in sub-Saharan Africa. Of these, approximately 1.6 million were new infections and 1.2 million deaths occurred. South Africa alone accounted for 31% of HIV/acquired immunodeficiency syndrome (AIDS) deaths in sub-Saharan Africa. This disturbing statistic indicates that South Africa remains the epicenter of the HIV/AIDS pandemic, compounded by the fact that only 36% of HIV-positive patients in South Africa have access to antiretroviral (ARV) treatment. Drug resistance mutations have emerged, and current ARVs show reduced efficacy against non-B subtypes. In addition, several recent studies have shown an increased prevalence of non-B African HIV strains in the Americas and Europe. Therefore, the use of ARVs in a non-B HIV-1 subtype context requires further investigation. HIV-1 subtype C protease, found largely in sub-Saharan Africa, has been under-investigated when compared with the subtype B protease, which predominates in North America and Europe. This review, therefore, focuses on HIV-1 proteases from B and C subtypes.

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 IC₅₀ 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 (IC₅₀) 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 (IC₅₀ = 0.076 μM) and the PCU-lactone-CO-EAIS peptide inhibitors (IC₅₀ = 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 (∆G_solv) 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 IC₅₀ data.

Electronic supplementary material

The online version of this article (doi:10.1186/s12929-015-0115-5) contains supplementary material, which is available to authorized users.

Recent patents and emerging therapeutics for HIV infections: a focus on protease inhibitors

The inclusion of protease inhibitors (PIs) in highly active antiretroviral therapy has significantly improved clinical outcomes in HIV-1-infected patients. To date, PIs are considered to be the most important therapeutic agents for the treatment of HIV infections. Despite high anti-HIV-1 potency, poor oral bioavailability of PIs has been a major concern. For achieving therapeutic concentrations, large doses of PIs are administered, which results in unacceptable systemic toxicities. Such severe and long-term toxicities necessitate the development of safer and potentially promising PIs. Recently, considerable attention has been paid to the development of newer compounds capable of inhibiting wild-type and resistant HIV-1 protease. Some of these PIs have displayed potent HIV-1 protease inhibitory activity. In this review, we have made an attempt to provide an overview on clinically approved and newly developing PIs, and related recent patents in the development of novel PIs.

Antiproliferative and antiviral activity of three libraries of adamantane derivatives

Three libraries of adamantane derivatives were synthesized and evaluated for antiviral and antiproliferative activities against a broad variety of DNA and RNA viruses. Whereas none of the compounds exhibit antiviral activity at subtoxic concentrations, antiproliferative activity was found against murine leukemia cells (L1210), human T-lymphocyte cells (CEM), and cervix carcinoma cells (HeLa) for 4, 8, and 10.

Combined quantum mechanics/molecular mechanics (QM/MM) methods in computational enzymology

Computational enzymology is a rapidly maturing field that is increasingly integral to understanding mechanisms of enzyme-catalyzed reactions and their practical applications. Combined quantum mechanics/molecular mechanics (QM/MM) methods are important in this field. By treating the reacting species with a quantum mechanical method (i.e., a method that calculates the electronic structure of the active site) and including the enzyme environment with simpler molecular mechanical methods, enzyme reactions can be modeled. Here, we review QM/MM methods and their application to enzyme-catalyzed reactions to investigate fundamental and practical problems in enzymology. A range of QM/MM methods is available, from cheaper and more approximate methods, which can be used for molecular dynamics simulations, to highly accurate electronic structure methods. We discuss how modeling of reactions using such methods can provide detailed insight into enzyme mechanisms and illustrate this by reviewing some recent applications. We outline some practical considerations for such simulations. Further, we highlight applications that show how QM/MM methods can contribute to the practical development and application of enzymology, e.g., in the interpretation and prediction of the effects of mutagenesis and in drug and catalyst design.

Synthesis, screening and computational investigation of pentacycloundecane-peptoids as potent CSA-HIV PR inhibitors

Herein, we present the first pentacycloundecane (PCU) diol peptoid derived HIV protease inhibitors with IC(50) values ranging from 6.5 to 0.075 μM. Five derivatives were synthesized in an attempt to understand the structure activity relationship of this class of compounds for HIV protease inhibition. NMR spectroscopy (new Efficient Adiabatic Symmetrized Rotating Overhauser Effect Spectroscopy, EASY-ROESY) was employed to determine the predominant conformation of the active compound. In this study docking studies and MD simulations provided insight into the binding theme of this class of peptoid inhibitors to the CSA-HIV PR active site. Conserved and stable hydrogen bonding between the hydroxyl groups of the inhibitors and the active site Asp25/Asp25' residues were observed from the docking and along the MD trajectories.

Synthesis, 2D-NMR and molecular modelling studies of pentacycloundecane lactam-peptides and peptoids as potential HIV-1 wild type C-SA protease inhibitors

In this study, eight non-natural peptides and peptoids incorporating the pentacycloundecane (PCU) lactam were designed and synthesized as potential inhibitors of the wild type C-SA HIV-protease. Five of these inhibitors gave IC(50) values ranging from 0.5 up to 0.75 µM against the resistance-prone wild type C-South African HIV-protease. NMR EASY-ROESY studies enabled us to describe the secondary structure of three of these compounds in solution. The 3D structures of the selected cage peptides were also modelled in solution using QM/MM/MD simulations. Satisfactory agreement between the NMR observations and the low energy calculated structures exists. Only one of these inhibitors (11 peptoid), which showed the best IC(50)(0.5 µM), exhibited a definable 3-D structure in solution. Autodock4 and AutodockVina were used to model the potential interaction between these inhibitors and the HIV-PR. It appears that the docking results are too crude to be correlated with the relative narrow range of experimental IC(50) values (0.5-10 µM). The PCU-peptides and peptoides were several orders less toxic (145 μM for 11 and 102 μM for 11 peptoid) to human MT-4 cells than lopinavir (0.025 μM). This is the first example of a polycyclic cage framework to be employed as an HIV-PR transition state analogue inhibitor and can potentially be utilized for other diseases related proteases. [Figure: see text].