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.

Identification of promiscuous T cell epitopes on Mayaro virus structural proteins using immunoinformatics, molecular modeling, and QM:MM approaches.

The Mayaro virus (MAYV) belongs to genus Alphavirus (family Togaviridae) and has been reported in several countries, especially in tropical regions of America. Due to its outbreaks and potential lack of medication, an effective vaccine formulation is strongly required. This study aimed to predict promiscuous T cell epitopes from structural polyproteins of MAYV using an immunoinformatics approach. For this purpose, consensus sequences were used to identify short protein sequences capable of binding to MHC class I and class II alleles. Our analysis pointed out 4 MHC-I/TCD8+ and 21 MHC-II/TCD4+ epitopes on capside (1;3), E1 (2;5), E2 (1;10), E3 (0;2), and 6 K (0;1) proteins. These predicted epitopes were characterized by high antigenicity, immunogenicity, conservancy, non-allergenic, non-toxic, and good population coverage rate values for North and South American geographical areas. Afterwards, we used the crystal structure of human toll-like receptor 3 (TLR3) ectodomain as a template to predict, through docking essays, the placement of a vaccine prototype at the TLR3 receptor binding site. Finally, classical and quantum mechanics/molecular mechanics (QM:MM) computations were employed to improve the quality of docking calculations, with the QM part of the simulations being accomplished by using the density functional theory (DFT) formalism. These results provide important insights into the advancement of diagnostic platforms, the development of vaccines, and immunotherapeutic interventions.

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

The emergent Zika virus (ZIKV) infection has become a threat to global health due to its association with severe neurological abnormalities, namely Guillain-Barré Syndrome (GBS) in adults and Congenital Zika virus Syndrome (CZS) in neonates. Many studies are nowadays being conducted to find an effective antiviral drug against ZIKV. In particular, NS2B-NS3 protease is an attractive drug target due to its essential function in viral replication, although a drug is not yet commercially available. In this context, we present here a comparative structural study, based on quantum chemistry calculations, to analyze the intermolecular binding energies between the crystallographic structure of NS2B-NS3 protease and dipeptide boronic acid (cn-716) and aldehyde (acyl-KR-aldehyde) peptidomimetic inhibitors, by using the molecular fractionation with conjugate caps (MFCC) scheme within the density functional theory (DFT) formalism. Most intermolecular interactions in cn-716/NS2B-NS3 (acyl-KR-aldehyde/NS2B-NS3) are due to the amino acid residues Asp83*, His51, Asp129, Ser81*, Gly133, Ala132, Tyr161, Asn152 and Asp75 (Asp83*, Asp129, His51, Asn152, Tyr161, Tyr130, Gly153, Gly151, Asp75, Pro131, and Gly82). Additionally, we have considered missense mutation analysis of these residues to evaluate the destabilization and the increase of the flexibility of the protease, showing that mutation of the residues Tyr161 and Tyr130 causes more impact. Our simulations are a valuable tool for a better understanding of the binding mechanism of recognized inhibitors of NS2B-NS3 protease, and can lead to the rational design and development of novel anti-Zika drugs with improved efficiency.

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

GABAB is a G protein-coupled receptor that functions as a constitutive heterodimer composed of the GABAB1a/b and GABAB2 subunits. It mediates slow and prolonged inhibitory neurotransmission in the nervous system, representing an attractive target for the treatment of various disorders. However, the molecular mechanism of the GABAB receptor is not thoroughly understood. Therefore, a better description of the binding of existing agonists and antagonists to this receptor is crucial to improve our knowledge about G protein-coupled receptor structure as well as for helping the development of new potent and more selective therapeutic agents. In this work, we used the recent X-ray cocrystallization data of agonists (GABA and baclofen) and antagonists (2-hydroxysaclofen, SCH50911, and CGP54626) bound to the GABAB orthosteric site together with quantum biochemistry and the molecular fractionation with conjugate caps (MFCC) scheme to describe the individual contribution of each amino acid residue involved in the GABAB-ligand interaction, pointing out differences and similarities among the compounds. Our quantum biochemical computational results show that the total binding energy of the ligands to the GABAB ligand pocket, with radius varying from 2.0 to 9.0 Å, is well-correlated with the experimental binding affinity. In addition, we found that the binding site is very similar for agonists or antagonists, showing small differences in the importance of the most significant amino acids. Finally, we predict the energetic relevance of the regions of the five ligands as well as the influence of each protein lobe on GABAB-ligand binding. These results provide important new information on the binding mechanism of the GABAB receptor and should facilitate the development of new chemicals targeting this receptor.

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.

In vitro antiplasmodial activity, pharmacokinetic profiles and interference in isoprenoid pathway of 2-aniline-3-hydroxy-1.4-naphthoquinone derivatives.

BACKGROUND: Plasmodium falciparum has shown multidrug resistance, leading to the necessity for the development of new drugs with novel targets, such as the synthesis of isoprenic precursors, which are excellent targets because the pathway is different in several steps when compared with the human host. Naphthoquinone derivatives have been described as potentially promising for the development of anti-malarial leader molecules. In view of that, the focus in this work is twofold: first, evaluate the in vitro naphthoquinone antiplasmodial activity and cytotoxicity; secondly, investigate one possible action mechanism of two derivatives of hydroxy-naphthoquinones. RESULTS: The two hydroxy-naphthoquinones derivatives have been tested against P. falciparum in vitro, using strains of parasites chloroquine-sensitive (3D7) and chloroquine-resistant (Dd2), causing 50% inhibition of parasite growth with concentrations that varied from 7 to 44.5 µM. The cell viability in vitro against RAW Cell Line displayed IC50 = 483.5 and 714.9 µM, whereas, in primary culture tests using murine macrophages, IC50 were 315.8 and 532.6 µM for the two selected compounds, causing no haemolysis at the doses tested. The in vivo acute toxicity assays exhibited a significant safety margin indicated by a lack of systemic and behavioural toxicity up to 300 mg/kg. It is suggested that this drug seems to inhibit the biosynthesis of isoprenic compounds, particularly the menaquinone and tocopherol. CONCLUSIONS: These derivatives have a high potential for the development of new anti-malarial drugs since they showed low toxicity associated to a satisfactory antiplasmodial activity and possible inhibition of a metabolic pathway distinct from the pathways found in the mammalian host.

Supramolecular aggregates of oligosaccharides with co-solvents in ternary systems for the solubilizing approach of triamcinolone.

A second compound is generally associated with oligosaccharides as a strategy to maximize the solubilizing effect for nonpolar compounds. This study elucidated the role and the mechanism whereby liquid compounds interact in these supramolecular aggregates in the solubilization of triamcinolone. Three different oligosaccharides (beta-cyclodextrin, 2-hydroxipropil-beta-cyclodextrin, and randomly methylated beta-cyclodextrin) and two potent co-solvents (triethanolamine and N-methyl pyrrolidone) were carefully evaluated by using three distinct experimental approaches. Incredibly stable complexes were formed with cyclodextrins (CDs). The structure of the complexes was elucidated by magnetic resonance spectra 2D-ROESY. The interactions of the protons of ring "A" of the drug with H(3) and H(5) protons of the CD cavity observed in the binary complexes remained in both ternary complexes. Unlike the observed ternary associations with triethanolamine, N-methyl pyrrolidone competed with the triamcinolone CD cavity and considerably decreased the stability of the complex and the solubility of the drug. The molecular dynamics (MD) and quantum mechanics:molecular mechanics (QM:MM) calculations supported that triethanolamine stabilized the drug-CD interactions for the conformer identified in the 2D-ROESY experiments, improving the quality and uniformity of the formed complex. The role played by the co-solvent in the ternary complexes depends on its specific ability to interact with the CD cavity in the presence of the drug, which can be predicted in theoretical studies to select the best candidate.