Exploring the Binding Mechanism of GABAB Receptor Agonists and Antagonists through in Silico Simulations

Quantum mechanics insights into melatonin and analogs binding to melatonin MT1 and MT2 receptors

Melatonin receptors MT1 and MT2 are G protein-coupled receptors that mediate the effects of melatonin, a hormone involved in circadian rhythms and other physiological functions. Understanding the molecular interactions between these receptors and their ligands is crucial for developing novel therapeutic agents. In this study, we used molecular docking, molecular dynamics simulations, and quantum mechanics calculation to investigate the binding modes and affinities of three ligands: melatonin (MLT), ramelteon (RMT), and 2-phenylmelatonin (2-PMT) with both receptors. Based on the results, we identified key amino acids that contributed to the receptor-ligand interactions, such as Gln181/194, Phe179/192, and Asn162/175, which are conserved in both receptors. Additionally, we described new meaningful interactions with Gly108/Gly121, Val111/Val124, and Val191/Val204. Our results provide insights into receptor-ligand recognition's structural and energetic determinants and suggest potential strategies for designing more optimized molecules. This study enhances our understanding of receptor-ligand interactions and offers implications for future drug development.

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.

Investigation of protein-protein interactions and hotspot region on the NSP7-NSP8 binding site in NSP12 of SARS-CoV-2

Background: The RNA-dependent RNA polymerase (RdRp) complex, essential in viral transcription and replication, is a key target for antiviral therapeutics. The core unit of RdRp comprises the nonstructural protein NSP12, with NSP7 and two copies of NSP8 (NSP81 and NSP82) binding to NSP12 to enhance its affinity for viral RNA and polymerase activity. Notably, the interfaces between these subunits are highly conserved, simplifying the design of molecules that can disrupt their interaction. Methods: We conducted a detailed quantum biochemical analysis to characterize the interactions within the NSP12-NSP7, NSP12-NSP81, and NSP12-NSP82 dimers. Our objective was to ascertain the contribution of individual amino acids to these protein-protein interactions, pinpointing hotspot regions crucial for complex stability. Results: The analysis revealed that the NSP12-NSP81 complex possessed the highest total interaction energy (TIE), with 14 pairs of residues demonstrating significant energetic contributions. In contrast, the NSP12-NSP7 complex exhibited substantial interactions in 8 residue pairs, while the NSP12-NSP82 complex had only one pair showing notable interaction. The study highlighted the importance of hydrogen bonds and π-alkyl interactions in maintaining these complexes. Intriguingly, introducing the RNA sequence with Remdesivir into the complex resulted in negligible alterations in both interaction energy and geometric configuration. Conclusion: Our comprehensive analysis of the RdRp complex at the protein-protein interface provides invaluable insights into interaction dynamics and energetics. These findings can guide the design of small molecules or peptide/peptidomimetic ligands to disrupt these critical interactions, offering a strategic pathway for developing effective antiviral drugs.

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.

Quantum Biochemistry and MM-PBSA Description of the ZIKV NS2B-NS3 Protease: Insights into the Binding Interactions beyond the Catalytic Triad Pocket

The Zika virus protease NS2B-NS3 has a binding site formed with the participation of a H51-D75-S135 triad presenting two forms, active and inactive. Studies suggest that the inactive conformation is a good target for the design of inhibitors. In this paper, we evaluated the co-crystallized structures of the protease with the inhibitors benzoic acid (5YOD) and benzimidazole-1-ylmethanol (5H4I). We applied a protocol consisting of two steps: first, classical molecular mechanics energy minimization followed by classical molecular dynamics were performed, obtaining stabilized molecular geometries; second, the optimized/relaxed geometries were used in quantum biochemistry and molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) calculations to estimate the ligand interactions with each amino acid residue of the binding pocket. We show that the quantum-level results identified essential residues for the stabilization of the 5YOD and 5H4I complexes after classical energy minimization, matching previously published experimental data. The same success, however, was not observed for the MM-PBSA simulations. The application of quantum biochemistry methods seems to be more promising for the design of novel inhibitors acting on NS2B-NS3.

Exploring the Spike-hACE 2 Residue-Residue Interaction in Human Coronaviruses SARS-CoV-2, SARS-CoV, and HCoV-NL63

Coronaviruses (CoVs) have been responsible for three major outbreaks since the beginning of the 21st century, and the emergence of the recent COVID-19 pandemic has resulted in considerable efforts to design new therapies against coronaviruses. Thus, it is crucial to understand the structural features of their major proteins related to the virus-host interaction. Several studies have shown that from the seven known CoV human pathogens, three of them use the human Angiotensin-Converting Enzyme 2 (hACE-2) to mediate their host's cell entry: SARS-CoV-2, SARS-CoV, and HCoV-NL63. Therefore, we employed quantum biochemistry techniques within the density function theory (DFT) framework and the molecular fragmentation with conjugate caps (MFCC) approach to analyze the interactions between the hACE-2 and the spike protein-RBD of the three CoVs in order to map the hot-spot residues that form the recognition surface for these complexes and define the similarities and differences in the interaction scenario. The total interaction energy evaluated showed a good agreement with the experimental binding affinity order: SARS-2 > SARS > NL63. A detailed investigation revealed the energetically most relevant regions of hACE-2 and the spike protein for each complex, as well as the key residue-residue interactions. Our results provide valuable information to deeply understand the structural behavior and binding site characteristics that could help to develop antiviral therapeutics that inhibit protein-protein interactions between CoVs S protein and hACE-2.

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.

New ethionamide boosters and EthR2: structural and energetic analysis

Ethionamide (ETH) is a high-profile drug for the treatment of patients with multidrug-resistant Mycobacterium tuberculosis and, in order to produce its inhibitory effects, it needs to be bioactivated by monooxygenase EthA. This process is under the control of the transcriptional repressors EthR and EthR2, so that their inhibition results in the boosting of ethionamide activation. Herein, through crystallographic data and computer simulations, we calculated the interaction binding energies of four inhibitors with improved in vitro potency, namely BDM76060 (PDB ID: 6HS1), BDM72201 (PDB ID: 6HRX), BDM76150 (PDB ID: 6HS2) and BDM72719 (PDB ID: 6HRY), in complexes with the transcriptional repressor EthR2, using density functional theory (DFT) within the molecular fractionation with conjugated caps (MFCC) approach. It was observed that these ligands share the same binding site within a 10.0 Å radius of the EthR2 protein; however, their structural particularities have a significant impact on the global energies of systems. The BDM72201 and BDM72719 components are weakly attached to the binding site, while BDM76060 and BDM76150 components produce stronger bonds, corroborating with experimental studies demonstrating that BDM76060 and BDM76150 are more successful in producing inhibitory effects. BDM76060 and BDM76150 have many functional groups that increase the contact surface with the protein and attract a more significant number of amino acid residues, being able to produce polarities that generate stronger interactions. In the current scenario of a growing number of cases of bacterial resistance, the obtained data can be used to guide clinical trials of these inhibitors and other inhibitors that act on the alternative EthR2 pathway, focusing on improving the activity of ethionamide, its effectiveness, a reduction in the treatment time and exposure to cytotoxic effects.

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.

Exploring human porphobilinogen synthase metalloprotein by quantum biochemistry and evolutionary methods

Previous studies have shown the porphobilinogen synthase (PBGS) zinc-binding mechanism and its conservation among the living cells. However, the precise molecular interaction of zinc with the active center of the enzyme is unknown. In particular, quantum chemistry techniques within the density functional theory (DFT) framework have been the key methodology to describe metalloproteins, when one is looking for a compromise between accuracy and computational feasibility. Considering this, we used DFT-based models within the molecular fractionation with conjugate caps scheme to evaluate the binding energy features of zinc interacting with the human PBGS. Besides, phylogenetic and clustering analyses were successfully employed in extracting useful information from protein sequences to identify groups of conserved residues that build the ions-binding site. Our results also report a conservative assessment of the relevant amino acids, as well as the benchmark analysis of the calculation models used. The most relevant intermolecular interactions in Zn2+-PBGS are due to the amino acids CYS0122, CYS0124, CYS0132, ASP0169, SER0168, ARG0221, HIS0131, ASP0120, GLY0133, VAL0121, ARG0209, and ARG0174. Among these residues, we highlighted ASP0120, GLY0133, HIS0131, SER0168, and ARG0209 by co-occurring in all clusters generated by unsupervised clustering analysis. On the other hand, the triple cysteines at 2.5 Å from zinc (CYS0122, CYS0124, and CYS0132) have the highest energy attraction and are absent in the taxa Viridiplantae, Sar, Rhodophyta, and some Bacteria. Additionally, the performance of the DFT-based models shows that the processing time-dependence is more associated with the choice of the basis set than the exchange-correlation functional.

Quantum binding energy features of the drug olmesartan bound to angiotensin type-1 receptors in the therapeutics of stroke

We investigated the binding energies of 105 residues within a 10 Å pocket radius, predicted the energetic relevance of olmesartan regions, and the influence of individual protein segments on OLM -AT1 binding.

The GABAB Receptor-Structure, Ligand Binding and Drug Development

The γ-aminobutyric acid (GABA) type B receptor (GABA_B-R) belongs to class C of the G-protein coupled receptors (GPCRs). Together with the GABA_A receptor, the receptor mediates the neurotransmission of GABA, the main inhibitory neurotransmitter in the central nervous system (CNS). In recent decades, the receptor has been extensively studied with the intention being to understand pathophysiological roles, structural mechanisms and develop drugs. The dysfunction of the receptor is linked to a broad variety of disorders, including anxiety, depression, alcohol addiction, memory and cancer. Despite extensive efforts, few compounds are known to target the receptor, and only the agonist baclofen is approved for clinical use. The receptor is a mandatory heterodimer of the GABA_B1 and GABA_B2 subunits, and each subunit is composed of an extracellular Venus Flytrap domain (VFT) and a transmembrane domain of seven α-helices (7TM domain). In this review, we briefly present the existing knowledge about the receptor structure, activation and compounds targeting the receptor, emphasizing the role of the receptor in previous and future drug design and discovery efforts.

mTOR–mLST8 interaction: hot spot identification through quantum biochemistry calculations

Quantum calculation of mTOR–mLST8 interaction.

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.