Quantum biochemical analysis of the TtgR regulator and effectors

The recent expansion of multidrug-resistant (MDR) pathogens poses significant challenges in treating healthcare-associated infections. Although antibacterial resistance occurs by numerous mechanisms, active efflux of the drugs is a critical concern. A single species of efflux pump can produce a simultaneous resistance to several drugs. One of the best-studied efflux pumps is the TtgABC: a tripartite resistance-nodulation-division (RND) efflux pump implicated in the intrinsic antibiotic resistance in Pseudomonas putida DOT-T1E. The expression of the TtgABC gene is down-regulated by the HTH-type transcriptional repressor TtgR. In this context, by employing quantum chemistry methods based on the Density Functional Theory (DFT) within the Molecular Fragmentation with Conjugate Caps (MFCC) approach, we investigate the coupling profiles of the transcriptional regulator TtgR in complex with quercetin (QUE), a natural polyphenolic flavonoid, tetracycline (TAC), and chloramphenicol (CLM), two broad-spectrum antimicrobial agents. Our quantum biochemical computational results show the: [i] convergence radius, [ii] total binding energy, [iii] relevance (energetically) of the ligands regions, and [iv] most relevant amino acids residues of the TtgR-QUE/TAC/CLM complexes, pointing out distinctions and similarities among them. These findings improve the understanding of the binding mechanism of effectors and facilitate the development of new chemicals targeting TtgR, helping in the battle against the rise of resistance to antimicrobial drugs. These advances are crucial in the ongoing fight against rising antimicrobial drug resistance, providing hope for a future where healthcare-associated infections can be more beneficially treated.

Quantum biochemical analysis of the binding interactions between a potential inhibitory drug and the Ebola viral glycoprotein

Ebola virus disease (EVD) causes outbreaks and epidemics in West Africa that persist until today. The envelope glycoprotein of Ebola virus (GP) consists of two subunits, GP1 and GP2, and plays a key role in anchoring or fusing the virus to the host cell in its active form on the virion surface. Toremifene (TOR) is a ligand that mainly acts as an estrogen receptor antagonist; however, a recent study showed a strong and efficient interaction with GP. In this context, we aimed to evaluate the energetic affinity features involved in the interaction between GP and toremifene by computer simulation techniques using the Molecular Fractionation Method with Conjugate Caps (MFCC) scheme and quantum-mechanical (QM) calculations, as well as missense mutations to assess protein stability. We identified ASP522, GLU100, TYR517, THR519, LEU186, LEU515 as the most attractive residues in the EBOV glycoprotein structure that form the binding pocket. We divided toremifene into three regions and evaluated that region i was more important than region iii and region ii for the formation of the TOR-GP1/GP2 complex, which might control the molecular remodeling process of TOR. The mutations that caused more destabilization were ARG134, LEU515, TYR517 and ARG559, while those that caused stabilization were GLU523 and ASP522. TYR517 is a critical residue for the binding of TOR, and is highly conserved among EBOV species. Our results may help to elucidate the mechanism of drug action on the GP protein of the Ebola virus and subsequently develop new pharmacological approaches against EVD.Communicated by Ramaswamy H. Sarma.

Antiviral Activity and Molecular Dynamics Simulation of Hops Compounds against Oropouche Virus (Peribunyaviridae)

The Oropouche virus (OROV) is a member of the family Peribunyaviridae (order Bunyavirales) and the cause of a dengue-like febrile illness transmitted mainly by biting midges and mosquitoes. In this study, we aimed to explore acylphloroglucinols and xanthohumol from hops (Humulus lupulus L.) as a promising alternative for antiviral therapies. The evaluation of the inhibitory potential of hops compounds on the viral cycle of OROV was performed through two complementary approaches. The first approach applies cell-based assay post-inoculation experiments to explore the inhibitory potential on the latest steps of the viral cycle, such as genome translation, replication, virion assembly, and virion release from the cells. The second part covers in silico methods evaluating the ability of those compounds to inhibit the activity of the endonuclease domain, which is essential for transcription, binding, and cleaving RNA. In conclusion, the beta acids showed strongest inhibitory potential in post-treatment assay (EC50 = 26.7 µg/mL). Xanthohumol had the highest affinity for OROV endonuclease followed by colupulone and cohumulone. This result contrasts with that observed for docking and MM/PBSA analysis, where cohumulone was found to have a higher affinity. Finally, among the three tested ligands, Lys92 and Arg33 exhibited the highest affinity with the protein.

In silico description of the adsorption of cell signaling pathway proteins ovalbumin, glutathione, LC3, TLR4, ASC PYCARD, PI3K and NF-Kβ on 7.0 nm gold nanoparticles: obtaining their Lennard-Jones-like potentials through docking and molecular mechanics

The impact of vaccination on the world's population is difficult to calculate. For developing different types of vaccines, adjuvants are substances added to vaccines to increase the magnitude and durability of the immune response and the effectiveness of the vaccine. This work explores the potential use of spherical gold nanoparticles (AuNPs) as adjuvants. Thus, we employed docking techniques and molecular mechanics to describe how a AuNP 7.0 nm in diameter interacts with cell signaling pathway proteins. Initially, we used X-ray crystallization data of the proteins ovalbumin, glutathione, LC3, TLR4, ASC PYCARD, PI3K, and NF-Kβ to study the adsorption with an AuNP through molecular docking. Therefore, interaction energies were obtained for the AuNP complexes and individual proteins, as well as the AuNP and OVA complex (AuNP@OVA) with each cellular protein, respectively. Results showed that AuNPs had the highest affinity for OVA individually, followed by glutathione, ASC PYCARD domain, LC3, PI3K, NF-Kβ, and TLR4. Furthermore, when evaluating the AuNP@OVA complex, glutathione showed a greater affinity with more potent interaction energy when compared to the other studied systems.

Quantum biochemistry description of PI3Kα enzyme bound to selective inhibitors

The PI3K class I is composed of four PI3K isoforms that serve as regulatory enzymes governing cellular metabolism, proliferation, and survival. The hyperactivation of PI3Kα is observed in various types of cancer and is linked to poor prognosis. Unfortunately, the development inhibitors selectively targeting one of the isoforms remains challenging, with only few agents in clinical use. The main difficulty arises from the high conservation among residues at the ATP-binding pocket across isoforms, which also serves as target pocket for inhibitors. In this work, molecular dynamics and quantum calculations were performed to investigate the molecular features guiding the binding of selective inhibitors, alpelisib and GDC-0326, into the ATP-binding pocket of PI3Kα. While molecular dynamics allowed crystallographic coordinates to relax, the interaction eergy between each amino acid residues and inhibitors was obtained by combining the Molecular Fractionation with Conjugated Caps scheme with Density Functional Theory calculations. In addition, the atomic charge of ligands in the bound and unbound (free) was calculated. Results indicated that the most relevant residues for the binding of alpelisib are Ile932, Glu859, Val851, Val850, Tyr836, Met922, Ile800, and Ile848, while the most important residues for the binding of GDC-0326 are Ile848, Ile800, Ile932, Gln859, Glu849, and Met922. In addition, residues Trp780, Ile800, Tyr836, Ile848, Gln859 Val850, Val851, Ile932 and Met922 are common hotspots for both inhibitors. Overall, the results from this work contribute to improving the understanding of the molecular mechanisms controlling selectivity and highlight important interactions to be considered during the rational design of new agents.Communicated by Ramaswamy H. Sarma.

Losartan as an ACE inhibitor: a description of the mechanism of action through quantum biochemistry

Losartan (LST) is a potent and selective angiotensin II (Ang II) type 1 (AT1) receptor antagonist widely used in the treatment of hypertension. The formation of Ang II is catalyzed by the angiotensin I-converting enzyme (ACE) through proteolytic cleavage of angiotensin I (Ang I), which is involved in the control of blood pressure. Despite the vast literature on the relationship of losartan with the renin-angiotensin system (RAS), the actions of losartan on the sACE enzyme are so far poorly understood. In view of this, we investigated how losartan can interact with the sACE enzyme to block its activity and intracellular signaling. After performing docking assays following quantum biochemistry calculations using losartan and sACE crystallographic data, we report that their interaction results reveal a new mechanism of action with important implications for understanding its effects on hypertension.

Computational methods for calculation of protein-ligand binding affinities in structure-based drug design

Drug design is an expensive and time-consuming process. Any method that allows reducing the time the costs of the drug development project can have great practical value for the pharmaceutical industry. In structure-based drug design, affinity prediction methods are of great importance. The majority of methods used to predict binding free energy in protein-ligand complexes use molecular mechanics methods. However, many limitations of these methods in describing interactions exist. An attempt to go beyond these limits is the application of quantum-mechanical description for all or only part of the analyzed system. However, the extensive use of quantum mechanical (QM) approaches in drug discovery is still a demanding challenge. This chapter briefly reviews selected methods used to calculate protein-ligand binding affinity applied in virtual screening (VS), rescoring of docked poses, and lead optimization stage, including QM methods based on molecular simulations.

Characterization of the binding interaction between atrazine and human serum albumin: Fluorescence spectroscopy, molecular dynamics and quantum biochemistry

Atrazine (ATR), one of the most used herbicides worldwide, causes persistent contamination of water and soil due to its high resistance to degradation. ATR is associated with low fertility and increased risk of prostate cancer in humans, as well as birth defects, low birth weight and premature delivery. Describing ATR binding to human serum albumin (HSA) is clinically relevant to future studies about pharmacokinetics, pharmacodynamics and toxicity of ATR, as albumin is the most abundant carrier protein in plasma and binds important small biological molecules. In this work we characterize, for the first time, the binding of ATR to HSA by using fluorescence spectroscopy and performing simulations using molecular docking, classical molecular dynamics and quantum biochemistry based on density functional theory (DFT). We determine the most likely binding sites of ATR to HSA, highlighting the fatty acid binding site FA8 (located between subdomains IA-IB-IIA and IIB-IIIA-IIIB) as the most important one, and evaluate each nearby amino acid residue contribution to the binding interactions explaining the fluorescence quenching due to ATR complexation with HSA. The stabilization of the ATR/FA8 complex was also aided by the interaction between the atrazine ring and SER454 (hydrogen bond) and LEU481(alkyl interaction).

Preclinical Investigations to Explore the Difference between the Diastereomers [177Lu]Lu-SibuDAB and [177Lu]Lu-RibuDAB toward Prostate Cancer Therapy

[¹⁷⁷Lu]Lu-Ibu-DAB-PSMA, a radioligand modified with ibuprofen as the albumin binder, showed higher accumulation in PSMA-positive tumors of mice than the clinically used [¹⁷⁷Lu]Lu-PSMA-617 but lower retention in non-targeted tissues than previously developed albumin-binding PSMA radioligands. The aim of this study was to investigate whether the stereochemistry of the incorporated ibuprofen affects the radioligand's in vitro and in vivo properties and to select the more favorable radioligand for further development. For this purpose, SibuDAB and RibuDAB containing (S)- and (R)-ibuprofen, respectively, were synthesized and labeled with lutetium-177. In vitro, the two isomers had similar properties; however, [¹⁷⁷Lu]Lu-SibuDAB showed increased binding to mouse and human plasma proteins (91 ± 1 and 88 ± 2%, respectively) compared to [¹⁷⁷Lu]Lu-RibuDAB (75 ± 2 and 79 ± 2%, respectively). In vivo, [¹⁷⁷Lu]Lu-SibuDAB was metabolically more stable than [¹⁷⁷Lu]Lu-RibuDAB with ∼90 vs ∼67% intact radioligand detected in the blood at 4 h post injection (p.i.). In line with the lower albumin-binding affinity, the blood clearance of [¹⁷⁷Lu]Lu-RibuDAB in mice was considerably faster [27% of injected activity (% IA), 1 h p.i.] than for [¹⁷⁷Lu]Lu-SibuDAB (50% IA, 1 h p.i.). Time-dependent biodistribution studies performed in tumor-bearing athymic nude mice showed high PSMA-specific tumor uptake for both isomers. A twofold increased area under the curve (AUC_0→8d) of the blood retention was determined for [¹⁷⁷Lu]Lu-SibuDAB as compared to [¹⁷⁷Lu]Lu-RibuDAB, whereas the kidney AUC_0→8d value of [¹⁷⁷Lu]Lu-SibuDAB was only half as high as for [¹⁷⁷Lu]Lu-RibuDAB. As a result, a more favorable tumor-to-kidney AUC_0→8d ratio was obtained for [¹⁷⁷Lu]Lu-SibuDAB, which was also visualized on SPECT/CT images. Based on its improved kidney clearance and higher metabolic stability, [¹⁷⁷Lu]Lu-SibuDAB was selected as the more favorable radioligand. Therapy studies performed with [¹⁷⁷Lu]Lu-SibuDAB (5 MBq/mouse) demonstrated the anticipated therapeutic superiority over the current gold-standard [¹⁷⁷Lu]Lu-PSMA-617 (5 MBq/mouse). The significantly increased survival time of mice treated with [¹⁷⁷Lu]Lu-SibuDAB as compared to those injected with [¹⁷⁷Lu]Lu-PSMA-617 justifies further development of this novel radioligand toward clinical application.

Blockade of the checkpoint PD-1 by its ligand PD-L1 and the immuno-oncological drugs pembrolizumab and nivolumab

We investigate the interaction between the programmed cell death protein 1 (PD-1) and the programmed cell death ligand 1 (PD-L1), as well as the immuno-oncological drugs pembrolizumab (PEM), and nivolumab (NIV), through quantum chemistry methods based on the Density Functional Theory (DFT) and the molecular fractionation with conjugate caps (MFCC) scheme, in order to map their hot-spot regions. Our results showed that the total interaction energy order of the three complexes is in good agreement with the experimental binding affinity order: PD-1/PEM > PD-1/NIV > PD-1/PD-L1. Besides, a detailed investigation revealed the energetically most relevant residue-residue pairs-interaction for each complex. Our computational results give a better understanding of the interaction mechanism between the protein PD-1 and its ligands (natural and inhibitors), unleashing the immune surveillance to destroy the cancer cells by decreasing their immune evasion. They are also an efficient alternative towards the development of new small-molecules and antibody-based drugs, pointing out to new treatments for cancer therapy.

Probe into the Molecular Mechanism of Ibuprofen Interaction with Warfarin Bound to Human Serum Albumin in Comparison to Ascorbic and Salicylic Acids: Allosteric Inhibition of Anticoagulant Release

The release of anticoagulant drugs such as warfarin from human serum albumin (HSA) has been important not only mechanistically but also clinically for patients who take multiple drugs simultaneously. In this study, the role of some commonly used drugs, including s-ibuprofen, ascorbic acid, and salicylic acid, was investigated in the release of warfarin bound to HSA in silico. The effects of the aforementioned drugs on the HSA-warfarin complex were investigated with molecular dynamics (MD) simulations using two approaches; in the first perspective, molecular docking was used to model the interaction of each drug with the HSA-warfarin complex, and in the second approach, drugs were positioned randomly and distant from the binary complex (HSA-warfarin) in a physiologically relevant concentration. The results obtained from both approaches indicated that s-ibuprofen and ascorbic acid both displayed allosteric effects on the release of warfarin from HSA. Although ascorbic acid aided in warfarin release, leading to destabilization of HSA, ibuprofen demonstrated a stabilizing effect on releasing the anticoagulant drug through several noncovalent interactions, including hydrophobic, electrostatic, and hydrogen-bonding interactions with the protein. The calculated binding free energy and energy contribution of involved residues using the molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) method, along with root mean square deviation (RMSD) values, protein gyration, and free energy surface (FES) mapping of the protein, provided valuable details on the nature of the interactions of each drug on the release of warfarin from HSA. These results can provide important information on the mechanisms of anticoagulant release that has not been revealed in molecular details previously.

Applicability of capillary electrophoresis-frontal analysis for displacement studies: Effect of several drugs on l-tryptophan and lidocaine binding to human serum albumin

The effective concentration of a drug in the blood, i.e. the concentration of a free drug in the blood, is influenced by the strength of drug binding onto plasma proteins. Besides its efficacy, these interactions subsequently influence the liberation, absorption, distribution, metabolism, excretion, and toxicological properties of the drug. It is important to not only determine the binding strength and stoichiometry, but also the binding site of a drug on the plasma protein molecule, because the co-administration of drugs with the same binding site can affect the above-mentioned concentration and as a result the pharmacological behavior of the drugs and lead to side effects caused by the change in free drug concentration, its toxicity. In this study, the binding characteristics of six drugs with human serum albumin, the most abundant protein in human plasma, were determined by capillary electrophoresis-frontal analysis, and the obtained values of binding parameters were compared with the literature data. The effect of several drugs and site markers on the binding of l-tryptophan and lidocaine to human serum albumin was investigated in subsequent displacement studies which thus demonstrated the usability of capillary electrophoresis as an automated high-throughput screening method for drug-protein binding studies.

Periodic DFT Calculations-Review of Applications in the Pharmaceutical Sciences

In the introduction to this review the complex chemistry of solid-state pharmaceutical compounds is summarized. It is also explained why the density functional theory (DFT) periodic calculations became recently so popular in studying the solid APIs (active pharmaceutical ingredients). Further, the most popular programs enabling DFT periodic calculations are presented and compared. Subsequently, on the large number of examples, the applications of such calculations in pharmaceutical sciences are discussed. The mentioned topics include, among others, validation of the experimentally obtained crystal structures and crystal structure prediction, insight into crystallization and solvation processes, development of new polymorph synthesis ways, and formulation techniques as well as application of the periodic DFT calculations in the drug analysis.

Interaction energies between two antiandrogenic and one androgenic agonist receptor in the presence of a T877A mutation in prostate cancer: a quantum chemistry analysis

We provide a detailed picture of the interactions between the androgen receptor T877A mutation and three different compounds used in the treatment of prostate cancer.

The urokinase plasminogen activator binding to its receptor: a quantum biochemistry description within an in/homogeneous dielectric function framework with application to uPA–uPAR peptide inhibitors

DFT calculations using the MFCC fragment-based model considering a spatial-dependent dielectric function based on the Poisson–Boltzmann approximation were performed to describe the uPA–uPAR interactions.

Intermolecular interactions of cn-716 and acyl-KR-aldehyde dipeptide inhibitors against Zika virus

Structural representation and graphic panel showing the most relevant residues that contribute to the ZIKV NS2B–NS3–ligand complexes.

Role of Herborn (K240E) and Milano Slow (D375H) human serum albumin variants towards binding of phenylbutazone and ibuprofen

Human serum albumin (HSA) is the binding cargo in blood plasma. The binding of drugs to HSA determines the pharmacokinetics and pharmacodynamics of the drugs. There are 67 natural genetic variants of HSA were reported in literature. Studying the effect of albumin modifications on drug binding helps to treat the patients with proper medication. In the present study, we have aimed to understand the effect of two natural variants of HSA, such as Herborn (K240E) and Milano Slow (D375H) on the binding of phenylbutazone and ibuprofen. For this, we have generated K240E and D375H mutants and also double mutant (K240E/D375H) of HSA using site directed mutagenesis. The recombinant HSA and its variants were expressed in Pichia pastoris. The interaction of HSA and its variants to phenylbutazone and ibuprofen was studied using fluorescence spectroscopy. Our results showed that there is no significant effect of K240E and D375H mutations on phenylbutazone and ibuprofen binding. But the effect is significant when both the mutations were there in a single protein (K240E/D375H). Further, the CD spectroscopy data showed that there is no effect of phenylbutazone and ibuprofen binding on the conformation of protein, except in case of D375H, where there is a conformational change in the binding pocket with the ibuprofen binding.

Ribosomal RNA-Aminoglycoside Hygromycin B Interaction Energy Calculation within a Density Functional Theory Framework

We intend to investigate the drug-binding energy of each nucleotide inside the aminoglycoside hygromycin B (hygB) binding site of 30S ribosomal RNA (rRNA) subunit by using the molecular fractionation with conjugate caps (MFCC) strategy based on the density functional theory (DFT), considering the functional LDA/PWC, OBS, and the dielectric constant parametrization. Aminoglycosides are bactericidal antibiotics that have high affinity to the prokaryotic rRNA, inhibiting the synthesis of proteins by acting on the main stages of the translation mechanism, whereas binding to rRNA 16S, a component of the 30S ribosomal subunit in prokaryotes. The identification of the nucleotides presenting the most negative binding energies allows us to stabilize hygB in a suitable binding pocket of the 30S ribosomal subunit. In addition, it should be highlighted that mutations in these residues may probably lead to resistance to ribosome-targeting antibiotics. Quantum calculations of aminoglycoside hygromycin B-ribosome complex might contribute to further quantum studies with antibiotics like macrolides and other aminoglycosides.

A quantum biochemistry approach to investigate checkpoint inhibitor drugs for cancer

We investigate the coupling profiles of the receptor PD-1 in complex with its natural ligand PDL1 and two inhibitor drugs.

Bezerra K, I. N. Oliveira J, Albuquerque EL, Fulco UL. Binding energies of the drugs capreomycin and streptomycin in complex with tuberculosis bacterial ribosome subunits

We calculated the interaction binding energies of the drugs capreomycin and streptomycin in complex with the tuberculosis bacterial ribosome subunits.

Inhibition of the checkpoint protein PD-1 by the therapeutic antibody pembrolizumab outlined by quantum chemistry

Much of the recent excitement in the cancer immunotherapy approach has been generated by the recognition that immune checkpoint proteins, like the receptor PD-1, can be blocked by antibody-based drugs with profound effects. Promising clinical data have already been released pointing to the efficiency of the drug pembrolizumab to block the PD-1 pathway, triggering the T-lymphocytes to destroy the cancer cells. Thus, a deep understanding of this drug/receptor complex is essential for the improvement of new drugs targeting the protein PD-1. In this context, by employing quantum chemistry methods based on the Density Functional Theory (DFT), we investigate in silico the binding energy features of the receptor PD-1 in complex with its drug inhibitor. Our computational results give a better understanding of the binding mechanisms, being also an efficient alternative towards the development of antibody-based drugs, pointing to new treatments for cancer therapy.

Computational biochemical investigation of the binding energy interactions between an estrogen receptor and its agonists

We present the energy profiles of estrogen receptor–agonist ligand interactions in atomic detail using a quantum biochemical approach.

Computational investigation of the α2β1 integrin–collagen triple helix complex interaction

In this paper, quantum biochemistry methods have been used to describe important protein–protein interactions for the complex integrin–collagen.

Explaining urokinase type plasminogen activator inhibition by amino-5-hydroxybenzimidazole and two naphthamidine-based compounds through quantum biochemistry

Urokinase plasminogen activator (uPA) is a biomarker and therapeutic target for several cancer types whose inhibition has been shown to slow tumor growth and metastasis.

Outlining migrainous through dihydroergotamine–serotonin receptor interactions using quantum biochemistry

We present the electronic structure of the complex dihydroergotamine–serotonin receptor to unveil new medications to treat migraine and related diseases.

Energetic description of cilengitide bound to integrin

We employ quantum chemistry methods to investigate the binding energy features of the cyclic RGD pentapeptide cilengitide interacting with the integrin receptor α_Vβ₃.

An improved quantum biochemistry description of the glutamate–GluA2 receptor binding within an inhomogeneous dielectric function framework

A new methodology to define the inhomogeneous dielectric constant of protein residues, to apply to the calculation of protein–ligand properties such as the electrostatic interaction.

Quantum binding energy features of the T3-785 collagen-like triple-helical peptide

Structural representation of the T3-785 collagen-like triple-helical peptide depicting the 15 most and fewest energetically significant amino acids.

Explaining RANKL inhibition by OPG through quantum biochemistry computations and insights into peptide-design for the treatment of osteoporosis

Quantum biochemistry computations are applied to precisely describe important protein–protein interactions, providing a basis for the design of inhibitory peptides against osteoporosis.

A quantum chemistry investigation of a potential inhibitory drug against the dengue virus

The total interaction energy of the inhibitor Bz-nKRR-H bound to a serine protease of the dengue virus is mainly due to the action of Asn152, Met49, Tyr161, Asp129 and Gly151 (Met84, Met75, Asp81, Asp79 and Asp80) residues at the NS3 (NS2B) subunit.

Studies of intermolecular interactions in solid dispersions using advanced surface chemical analysis

The aim of this study is to utilise an advanced surface chemical analysis based on X-ray photoelectron spectroscopy (XPS) to determine and characterise drug/polymer interactions in solid dispersions manufactured via hot melt extrusion (HME).