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.

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.

A Bayesian Analysis of Plant DNA Length Distribution via κ-Statistics.

We report an analysis of the distribution of lengths of plant DNA (exons). Three species of Cucurbitaceae were investigated. In our study, we used two distinct κ distribution functions, namely, κ-Maxwellian and double-κ, to fit the length distributions. To determine which distribution has the best fitting, we made a Bayesian analysis of the models. Furthermore, we filtered the data, removing outliers, through a box plot analysis. Our findings show that the sum of κ-exponentials is the most appropriate to adjust the distribution curves and that the values of the κ parameter do not undergo considerable changes after filtering. Furthermore, for the analyzed species, there is a tendency for the κ parameter to lay within the interval (0.27;0.43).

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.

Binding energies of the drugs capreomycin and streptomycin in complex with tuberculosis bacterial ribosome subunits.

Despite advances, tuberculosis remains a significant infectious disease, whose mortality presents alarming numbers. Although it can be cured, the number of cases of antimicrobial resistant strains is increasing, requiring the use of less efficient second-line drugs. Capreomycin and streptomycin are part of this group, being antibiotics whose mechanism of action is the inhibition of protein synthesis when interacting with the tuberculosis bacterial ribosome. Their binding mechanisms are distinct: capreomycin is able to bind to both ribosomal (30S and 50S) subunits, whereas streptomycin binds only to the smaller one (30S). In this context, the biochemical characterization of these binding sites for a proper understanding of their complex interactions is of crucial importance to increase their efficacy. Through crystallographic data and computer simulations, in this work we calculated the interaction binding energies of capreomycin and streptomycin in complex with the tuberculosis bacterial ribosome subunits, by using density functional theory (DFT) within the molecular fractionation with conjugated caps (MFCC) approach. For capreomycin in the 30S (50S) subunit, we investigated the binding energies of 44 (30) residues presented within a pocket radius of 14 Å (30 Å). Regarding streptomycin, 60 nucleotide (25 amino acid) residues distributed up to 12.5 Å (15 Å) away from the drug in the 30S subunit (S12 protein) were taken into account. We also identify the contributions of hydrogen bonds and hydrophobic interactions in the drug-receptor complex, and the regions of the drugs that most contributed to the anchorages of them in their binding sites, as well as identify residues that are most associated with mutations.

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.

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.

DFT Calculations with van der Waals Interactions of Hydrated Calcium Carbonate Crystals CaCO3·(H2O, 6H2O): Structural, Electronic, Optical, and Vibrational Properties.

The role of hydration on the structural, electronic, optical, and vibrational properties of monohydrated (CaCO3·H2O, hexagonal, P31, Z = 9) and hexahydrated (CaCO3·6H2O, monoclinic, C2/c, Z = 4) calcite crystals is assessed with the help of published experimental and theoretical data applying density functional theory within the generalized gradient approximation and a dispersion correction scheme. We show that the presence of water increases the main band gap of monohydrocalcite by 0.4 eV relative to the anhydrous structure, although practically not changing the hexahydrocalcite band gap. The gap type, however, is modified from indirect to direct as one switches from the monohydrated to the hexahydrated crystal. A good agreement was obtained between the simulated vibrational infrared and Raman spectra and the experimental data, with an infrared signature of hexahydrocalcite relative to monohydrocalcite being observed at 837 cm(-1). Other important vibrational signatures of the lattice, water molecules, and CO3(2-) were identified as well. Analysis of the phonon dispersion curves shows that, as the hydration level of calcite increases, the longitudinal optical-transverse optical phonon splitting becomes smaller. The thermodynamics properties of hexahydrocalcite as a function of temperature resemble closely those of calcite, while monohydrocalcite exhibits a very distinct behavior.

A quantum biochemistry investigation of willardiine partial agonism in AMPA receptors.

We employ quantum biochemistry methods based on the Density Functional Theory (DFT) approach to unveil the detailed binding energy features of willardiines co-crystallized with the AMPA receptor. Our computational results demonstrate that the total binding energies of fluorine-willardiine (FW), hydrogen-willardiine (HW), bromine-willardiine (BrW) and iodine-willardiine (IW) to the iGluR2 ligand-pocket correlate with the agonist binding energies, whose experimental sequential data match our computational counterpart, excluding the HW case. We find that the main contributions to the total willardiine-iGluR2 binding energy are due to the amino acid residues in decreasing order Glu705 > Arg485 > Ser654 > Tyr450 > T655. Furthermore, Met708, which is positioned close to the 5-substituent, attracts HW and FW, but repels BrW and IW. Our results contribute significantly to an improved understanding of the willardiine-iGluR2 binding mechanisms.