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.

Deciphering Interactions between Potential Inhibitors and the Plasmodium falciparum DHODH Enzyme: A Computational Perspective

Malaria is a parasitic disease that, in its most severe form, can even lead to death. Insect-resistant vectors, insufficiently effective vaccines, and drugs that cannot stop parasitic infestations are making the fight against the disease increasingly difficult. It is known that the enzyme dihydroorotate dehydrogenase (DHODH) is of paramount importance for the synthesis of pyrimidine from the Plasmodium precursor, that is, for its growth and reproduction. Therefore, its blockade can lead to disruption of the parasite's life cycle in the vertebrate host. In this scenario, PfDHODH inhibitors have been considered candidates for a new therapy to stop the parasitic energy source. Given what is known, in this work, we applied molecular fractionation with conjugated caps (MFCC) in the framework of the quantum formalism of density functional theory (DFT) to evaluate the energies of the interactions between the enzyme and the different triazolopyrimidines (DSM483, DMS557, and DSM1), including a complex carrying the mutation C276F. From these results, it was possible to identify the main features of each system, focusing on the wild-type and mutant PfDHODH and examining the major amino acid residues that are part of the four complexes. Our analysis provides new information that can be used to develop new drugs that could prove to be more effective alternatives to present antimalarial drugs.

Identification of 3, 4-dihydroxy complexes as potential antiviral via DFT, molecular docking, molecular dynamics and MM/PBSA against rabies and dengue receptors

The quest to identify antiviral drug candidates for dengue and rabies viral diseases is a great challenge for the researchers. While different research is being conducted on the repurposed drugs against these two viruses, no drug compound has gained success in treating them. Therefore, in this study, 3, 4-dihydroxy complexes have been virtually designed to investigate their antiviral properties and analyze their efficiency in interaction with the concerned viral diseases. DFT calculations are carried out to study the electronic and thermodynamic properties to understand the stability and reactivity of the reported compounds. These compounds were subjected to molecular docking studies to understand the binding interactions with NS5 Dengue virus mRNA 2'-O-methyltransferase and phosphoprotein C-terminal domain of Rabies virus. MD simulation, hydrogen bond analysis, and MM/PBSA were performed at 100 ns to support the obtained docking results.Communicated by Ramaswamy H. Sarma.

Assessment of Covalently Binding Warhead Compounds in the Validation of the Cytomegalovirus Nuclear Egress Complex as an Antiviral Target

Herpesviral nuclear egress is a regulated process of viral capsid nucleocytoplasmic release. Due to the large capsid size, a regular transport via the nuclear pores is unfeasible, so that a multistage-regulated export pathway through the nuclear lamina and both leaflets of the nuclear membrane has evolved. This process involves regulatory proteins, which support the local distortion of the nuclear envelope. For human cytomegalovirus (HCMV), the nuclear egress complex (NEC) is determined by the pUL50-pUL53 core that initiates multicomponent assembly with NEC-associated proteins and capsids. The transmembrane NEC protein pUL50 serves as a multi-interacting determinant that recruits regulatory proteins by direct and indirect contacts. The nucleoplasmic core NEC component pUL53 is strictly associated with pUL50 in a structurally defined hook-into-groove complex and is considered as the potential capsid-binding factor. Recently, we validated the concept of blocking the pUL50-pUL53 interaction by small molecules as well as cell-penetrating peptides or an overexpression of hook-like constructs, which can lead to a pronounced degree of antiviral activity. In this study, we extended this strategy by utilizing covalently binding warhead compounds, originally designed as binders of distinct cysteine residues in target proteins, such as regulatory kinases. Here, we addressed the possibility that warheads may likewise target viral NEC proteins, building on our previous crystallization-based structural analyses that revealed distinct cysteine residues in positions exposed from the hook-into-groove binding surface. To this end, the antiviral and NEC-binding properties of a selection of 21 warhead compounds were investigated. The combined findings are as follows: (i) warhead compounds exhibited a pronounced anti-HCMV potential in cell-culture-based infection models; (ii) computational analysis of NEC primary sequences and 3D structures revealed cysteine residues exposed to the hook-into-groove interaction surface; (iii) several of the active hit compounds exhibited NEC-blocking activity, as shown at the single-cell level by confocal imaging; (iv) the clinically approved warhead drug ibrutinib exerted a strong inhibitory impact on the pUL50-pUL53 core NEC interaction, as demonstrated by the NanoBiT assay system; and (v) the generation of recombinant HCMV ∆UL50-ΣUL53, allowing the assessment of viral replication under conditional expression of the viral core NEC proteins, was used for characterizing viral replication and a mechanistic evaluation of ibrutinib antiviral efficacy. Combined, the results point to a rate-limiting importance of the HCMV core NEC for viral replication and to the option of exploiting this determinant by the targeting of covalently NEC-binding warhead compounds.

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).

In silico analysis of energy interactions between nociceptin/orfanin FQ receptor and two antagonists with potential antidepressive action

This study addresses the binding energies of NOPR-ligand complexes and presents the main amino acid residues involved in the interaction between these complexes.

Quantum Biochemical Investigation of Lys49-PLA2 from Bothrops moojeni

Envenomation via snakebites occurs largely in areas where it is harder to access the hospital. Its mortality rate and sequelae acquired by the survivors symbolize a big challenge for antivenom therapy. In particular, the homologous phospholipase A₂ (Lys49-PLA₂) proteins can induce myonecrosis and are not effectively neutralized by current treatments. Thus, by taking advantage of crystallographic structures of Bothrops moojeni Lys49-PLA₂ complexed with VRD (varespladib) and AIN (aspirin), a quantum biochemistry study based on the molecular fractionation with conjugate cap scheme within the density functional theory formalism is performed to unveil these complexes' detailed interaction energies. The calculations revealed that important interactions between ligands and the Lys49-PLA₂ pocket could occur up to a pocket radius of r = 6.5 (5.0 Å) for VRD (AIN), with the total interaction energy of the VRD ligand being higher than that of the AIN ligand, which is well-correlated with the experimental binding affinity. Furthermore, we have identified the role played by the amino acids LYS0069, LYS0049, LEU0005, ILE0009, CYS0029, GLY0030, HIS0048, PRO0018, ALA0019, CYS0045, TYR0052, TYR0022, PRO0125*, and PHE0126* (LYS0069, LYS0049, GLY0032, LEU0002, and LEU0005) in the VRD↔Lys49-PLA₂ (AIN↔Lys49-PLA₂) complex. Our simulations are a valuable tool to support the big challenge for neutralizing the damages in victims of snakebites.

In silico studies of the inhibition mechanism of dengue with papain

Dengue virus is becoming a major global disease; the envelope protein is the major target for vaccine development against Dengue. Nowadays, the attention has focused on developing inhibitors based on Papain is a promising target for treating Dengue. In the present work, the theoretical studies of E-protein(Cys74-Glu79;Lys110)…Papain(Cys25, Asn175 and His159) complexes are analysed by Density Functional Theory (M06-2X/cc-pVDZ) method. Among the E-protein(Cys74-Glu79;Lys110)…Papain(Cys25, Asn175 and Hys159) complexes, E-protein(Glu76)…Papain(Cys25) complex has the highest interaction value of -352.22 kcal/mol. Moreover, the natural bond orbital analysis also supports the above results. The 100 ns Molecular Dynamics simulation reveals that, E-protein(Ala54-Ile129)…Papain(Cys25) complex had the lowest root mean square deviation value of 1 Å compared to the E-protein(Ala54-Ile129)… Papain(Asn175 & His159) complexes. The salt bridge formation between the Asp103 and Lys110 residues are the important stabilizing factor in E-protein(Ala54-Ile129)…Papain(Cys25) complex. This result can extend our knowledge of the functional behaviour of Papain and provides structural insight to target Envelope protein as forthcoming drug targets in Dengue.

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.

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.

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 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.

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.