Quantum Biochemistry Research Articles

Shedding light on cancer immunology at the molecular level: A quantum biochemistry study of representative PD-1/PD-L1 conformations

BackgroundProgrammed death 1 (PD-1) binding to PD-L1 is a potent mechanism used by immunogenic tumors to evade the immune system and the immune checkpoint PD-1PD-L1 has emerged as a promising target in the search for new drugs to improve cancer treatment. The crystallographic structure of humanPD-1humanPD-L1 shed light on the molecular characterization of this system and allowed computational studies to be carried out to characterize structural behaviors. MethodsThis study demonstrated the importance of analyzing the flexibility of protein systems through molecular dynamics simulations (MDS) and its impacts on the interaction energy obtained through quantum biochemistry. ResultsThe computational results obtained provide a description of the flexibility and energetic profile of the PD-1PD-L1 contact surface using representative conformations from MDS. Variations of up to 50 % in the total interaction energy values were detected depending on the scrutinized conformation, which can be mainly attributed to the flexibility of the CC' loop, FG loop and ASP85-GLN91 of PD-1 and the MET58-LYS62 segment of PD-L1. Quantum biochemistry revealed the three hot spots in PD-L1: ARG113L-ARG125L > ILE54L-VAL76L > ALA18L-ASP26L; and two energetic hot spots in PD-1: ALA125-ARG139 > VAL63-GLN88. Nonetheless, VAL63-GLN88 and GLY124-ARG139 exhibit significant variation in interaction energy between different conformations, while ARG113L-ARG125L is the only hot spot with high energetic fluctuation on the PD-L1 surface. ConclusionThis is the first application of MDS coupled to dimensionality reduction and density functional theory (DFT) demonstrating new structural and energetic features that might be useful in discovering/designing more potent PD-1PD-L1 inhibitors.

Alzheimer's Disease Immunotherapy and Mimetic Peptide Design for Drug Development: Mutation Screening, Molecular Dynamics, and a Quantum Biochemistry Approach Focusing on Aducanumab::Aβ2-7 Binding Affinity.

Seven treatments are approved for Alzheimer's disease, but five of them only relieve symptoms and do not alter the course of the disease. Aducanumab (Adu) and lecanemab are novel disease-modifying antiamyloid-β (Aβ) human monoclonal antibodies that specifically target the pathophysiology of Alzheimer's disease (AD) and were recently approved for its treatment. However, their administration is associated with serious side effects, and their use is limited to early stages of the disease. Therefore, drug discovery remains of great importance in AD research. To gain new insights into the development of novel drugs for Alzheimer's disease, a combination of techniques was employed, including mutation screening, molecular dynamics, and quantum biochemistry. These were used to outline the interfacial interactions of the Aducanumab::Aβ2-7 complex. Our analysis identified critical stabilizing contacts, revealing up to 40% variation in the affinity of the Adu chains for Aβ2-7 depending on the conformation outlined. Remarkably, two complementarity determining regions (CDRs) of the Adu heavy chain (HCDR3 and HCDR2) and one CDR of the Adu light chain (LCDR3) accounted for approximately 77% of the affinity of Adu for Aβ2-7, confirming their critical role in epitope recognition. A single mutation, originally reported to have the potential to increase the affinity of Adu for Aβ2-7, was shown to decrease its structural stability without increasing the overall binding affinity. Mimetic peptides that have the potential to inhibit Aβ aggregation were designed by using computational outcomes. Our results support the use of these peptides as promising drugs with great potential as inhibitors of Aβ aggregation.

Unveiling fructose and glucose binding to human serum albumin: fluorescence measurements and docking, molecular dynamics and quantum biochemistry computations

This research examines the interaction between human serum albumin (HSA) and various sugar forms (β-D-fructofuranose (FRC), α-D-glucopyranose (GLC), Keto-D-fructose (FRO), Aldehydo-D-glucose (GLO), and modified Aldehydo-D-glucose (GLOm)) using fluorescent spectroscopy, molecular docking simulations, molecular dynamics, protein conformational clusters (EnGens), molecular fractionation with conjugate caps (MFCC) and quantum biochemistry analysis. We analyze molecular and quantum aspects, uncovering interaction energies between sugar atoms and amino acids. Total interaction energy considers protein fragmentation, energetic decomposition, and interaction energy from a bottom-up perspective. Molecular dynamics reveal that unmodified Aldehydo-D-glucose (GLO) escapes HSA binding sites, explaining gradual glycation. We pioneer studying HSA’s binding mechanism with glucose and fructose in a 1:1 ratio using long molecular dynamics simulations. Results suggest the transitional GLOm form has a higher Sudlow I site propensity than unmodified glucose, crucial for K195 glycation. FRO and GLOm interaction tendencies move toward a deeper FA7 cavity, near its center. This approach effectively elucidates small molecule binding mechanisms, consistent with previous experimental results. Communicated by Ramaswamy H. Sarma

Tofacitinib and peficitinib inhibitors of Janus kinase for autoimmune disease treatment: a quantum biochemistry approach.

Autoimmune inflammatory diseases, such as rheumatoid arthritis (RA) and ulcerative colitis, are associated with an uncontrolled production of cytokines leading to the pronounced inflammatory response of these disorders. Their therapy is currently focused on the inhibition of cytokine receptors, such as the Janus kinase (JAK) protein family. Tofacitinib and peficitinib are JAK inhibitors that have been recently approved to treat rheumatoid arthritis. In this study, an in-depth analysis was carried out through quantum biochemistry to understand the interactions involved in the complexes formed by JAK1 and tofacitinib or peficitinib. Computational analyses provided new insights into the binding mechanisms between tofacitinib or peficitinib and JAK1. The essential amino acid residues that support the complex are also identified and reported. Additionally, we report new interactions, such as van der Waals; hydrogen bonds; and alkyl, pi-alkyl, and pi-sulfur forces, that stabilize the complexes. The computational results revealed that peficitinib presents a similar affinity to JAK1 compared to tofacitinib based on their interaction energies.

In silico approach revealed the membrane receptor PHO36 as a new target for synthetic anticandidal peptides.

Aim: Synthetic antimicrobial peptides (SAMPs) present the potential to fight systemic fungal infections. Here, the PHO36 receptor from Candida albicans was analyzed by in silico tools as a possible target for three anticandidal SAMPs: RcAlb-PepIII, PepGAT and PepKAA.Materials & methods: Molecular docking, dynamics and quantum biochemistry were employed to understand the individual contribution of amino acid residues in the interaction region.Results: The results revealed that SAMPs strongly interact with the PHO36 by multiple high-energy interactions. This is the first study to employ quantum biochemistry to describe the interactions between SAMPs and the PHO36 receptor.Conclusion: This work contributes to understanding and identifying new molecular targets with medical importance that could be used to discover new drugs against systemic fungal infections.

Quantum Biochemistry: Unveiling the Quantum Coherence in Enzyme Catalysis

Objective: This study investigates the influence of quantum coherence in enzyme catalysis, aiming to elucidate its role in biochemical processes. The primary objective is to unravel quantum effects within various enzymes (A-F) involved in crucial biochemical pathways. Methodology: Employing a multidisciplinary approach, advanced experimental techniques and computational methods were utilized. Quantum tunneling rates were measured through Reaction Progress Kinetic Analysis (RPKA) with High-Performance Liquid Chromatography (HPLC). Femtosecond-resolved spectroscopy captured quantum coherence times, while Two-Dimensional Infrared Spectroscopy (2D IR) probed vibrational coupling. Ultrafast Laser Spectroscopy provided insights into enzyme dynamics. Density Functional Theory (DFT) calculations and Ab Initio Simulations complemented experimental findings. Results: The results reveal distinct quantum signatures across all enzymes. Notably, Enzyme A demonstrates a quantum tunneling rate of 3.2 x 10^-2 s^-1. Quantum coherence times in Enzyme B showcase unprecedented femtosecond scales, while other enzymes exhibit diverse behaviors. DFT calculations for Enzyme E predict a 30% reduction in energy barriers. Ab Initio Simulations of Enzyme F unveil persistent entanglement states. Conclusion: The observed quantum phenomena suggest a profound interplay between quantum coherence and enzyme catalysis, emphasizing the enzyme-specific nature of quantum effects. The implications of energy barrier reduction and entanglement states provide insights into potential quantum-assisted catalytic mechanisms.

Noncovalent binding of carbofuran to acetylcholinesterase from Homo sapiens, Danio rerio, Apis mellifera and Caenorhabditis elegans: Homology modelling, molecular docking and dynamics, and quantum biochemistry description

Although various regulatory agencies have banned or severely restricted the use of carbofuran (CAR), recent reports indicate the presence of CAR residues in both cultivated and wild areas. This pesticide is a potent inhibitor of acetylcholinesterase (AChE), which acts by preventing the hydrolysis of acetylcholine (ACh). Given the critical role of AChE::ACh in the proper functioning of the nervous system, we thought it appropriate to investigate the binding of CAR to AChEs from Homo sapiens, Danio rerio, Apis mellifera, and Caenorhabditis elegans using homology modelling, molecular docking, molecular dynamics, and quantum biochemistry. Molecular docking and dynamics results indicated peculiar structural behavior in each AChE::CAR system. Quantum biochemistry results showed similar affinities for all complexes, confirming the description of carbofuran as a broad-spectrum pesticide, and have a limited correlation with IC50 values. We found the following decreasing affinity order of AChE species: H. sapiens > A. mellifera > C. elegans > D. rerio. The computational results suggest that CAR occupies different pockets in the AChEs studied. In addition, our results showed that CAR binds to hsAChE and ceAChE in a very similar manner: it has high affinities for the same subsites in both species and forms hydrogen bonds with residues (hsTYR124 and ceTRP107) occupying homologous positions in the peripheral site. This suggests that this nematode is a potential model to evaluate the toxicity of carbamates, even though the sequence identity between them is only 41 %. Interestingly, we also observed that the catalytic histidines of drAChE and amAChE exhibited favorable contacts with carbofuran, suggesting that the non-covalent binding of carbofuran to these proteins may promote faster carbamylation rates than the binding modes to human and worm acetylcholinesterases. Our computational results provide a better understanding of the binding mechanisms in these complexes, as well as new insights into the mechanism of carbamylation.

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

Plant Protease Inhibitors Against SARS-CoV-2 Main Protease: An In⟓ilico Approach

Aim: To evaluate using bioinformatics tools the interactions between plant protease inhibitors (PIs) and SARS-CoV-2 Mpro. Materials & methods: This was an in silico study based on molecular docking, molecular dynamics simulations and quantum biochemistry calculations. Results: The plant PIs pineapple cystatin and sesame cystatin interacted allosterically with Mpro, leading to significant structural alterations. These conformational changes lead to a reduction of the area and volume of the Mpro proteolytic site, likely affecting the protease activity and, consequently, reducing viral replication. Conclusion: This work highlights the therapeutic potential of plant PIs as candidates for future in vivo investigations into new therapeutics for COVID-19.

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.

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

Quantum Biochemistry Screening and In Vitro Evaluation of Leishmania Metalloproteinase Inhibitors.

Leishmanolysin, also known as major promastigote protease (PSP) or gp63, is the most abundant surface glycoprotein of Leishmania spp., and has been extensively studied and recognized as the main parasite virulence factor. Characterized as a metalloprotease, gp63 can be powerfully inactivated in the presence of a metal chelator. In this study, we first used the structural parameters of a 7-hydroxycoumarin derivative, L1 compound, to evaluate the theoretical–computational experiments against gp63, comparing it with an available metal chelator already described. The methodology followed was (i) analysis of the three-dimensional structure of gp63 as well as its active site, and searching the literature and molecular databases for possible inhibitors; (ii) molecular docking simulations and investigation of the interactions in the generated protein–ligand complexes; and (iii) the individual energy of the gp63 amino acids that interacted most with the ligands of interest was quantified by ab initio calculations using Molecular Fraction with Conjugated Caps (MFCC). MFCC still allowed the final quantum balance calculations of the protein interaction to be obtained with each inhibitor candidate binder. L1 obtained the best energy quantum balance result with −2 eV, followed by DETC (−1.4 eV), doxycycline (−1.3 eV), and 4-terpineol (−0.6 eV), and showed evidence of covalent binding in the enzyme active site. In vitro experiments confirmed L1 as highly effective against L. amazonensis parasites. The compound also exhibited a low cytotoxicity profile against mammalian RAW and 3T3 cells lines, presenting a selective index of 149.19 and 380.64 µM, respectively. L1 induced promastigote forms’ death by necrosis and the ultrastructural analysis revealed disruption in membrane integrity. Furthermore, leakage of the contents and destruction of the parasite were confirmed by Spectroscopy Dispersion analysis. These results together suggested L1 has a potential effect against L. amazonensis, the etiologic agent of diffuse leishmaniasis, and the only one that currently does not have a satisfactory treatment.

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.

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.

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 A2 (Lys49-PLA2) proteins can induce myonecrosis and are not effectively neutralized by current treatments. Thus, by taking advantage of crystallographic structures of Bothrops moojeni Lys49-PLA2 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-PLA2 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-PLA2 (AIN↔Lys49-PLA2) complex. Our simulations are a valuable tool to support the big challenge for neutralizing the damages in victims of snakebites.

Computational approach, scanning electron and fluorescence microscopies revealed insights into the action mechanisms of anticandidal peptide Mo-CBP3-PepIII

AimsThe Candida genus is composed of opportunistic pathogens that threaten public health. Given the increase in resistance to current drugs, it is necessary to develop new drugs to treat infections by these pathogens. Antimicrobial peptides are promising alternative molecules with low cost, broad action spectrum and low resistance induction. This study aimed to clarify the action mechanisms of synthetic peptides against Candida albicans. Main methodsThe mode of action of the anticandidal peptides Mo-CBP3-PepIII were analyzed through molecular dynamics and quantum biochemistry methods against Exo-β-1,3-glucanase (EXG), vital to cell wall metabolism. Furthermore, scanning electron (SEM) and fluorescence (FM) microscopies were employed to corroborate the in silico data. Key findingsMo-CBP3-PepIII strongly interacted with EXG (−122.2 kcal mol−1) at the active site, higher than the commercial inhibitor pepstatin. Also, molecular dynamics revealed the insertion of Mo-CBP3-PepIII into the yeast membrane. SEM analyses revealed that Mo-CBP3-PepIII induced cracks and scars of the cell wall and FM analyses confirmed the pore formation on the Candida membrane. SignificanceMo-CBP3-PepIII has strong potential as a new drug with a broad spectrum of action, given its different mode of action compared to conventional drugs.

Quantum biochemistry, molecular docking, and dynamics simulation revealed synthetic peptides induced conformational changes affecting the topology of the catalytic site of SARS-CoV-2 main protease

The recent outbreak caused by SARS-CoV-2 continues to threat and take many lives all over the world. The lack of an efficient pharmacological treatments are serious problems to be faced by scientists and medical staffs worldwide. In this work, an in silico approach based on the combination of molecular docking, dynamics simulations, and quantum biochemistry revealed that the synthetic peptides RcAlb-PepI, PepGAT, and PepKAA, strongly interact with the main protease (Mpro) a pivotal protein for SARS-CoV-2 replication. Although not binding to the proteolytic site of SARS-CoV-2 Mpro, RcAlb-PepI, PepGAT, and PepKAA interact with other protein domain and allosterically altered the protease topology. Indeed, such peptide-SARS-CoV-2 Mpro complexes provoked dramatic alterations in the three-dimensional structure of Mpro leading to area and volume shrinkage of the proteolytic site, which could affect the protease activity and thus the virus replication. Based on these findings, it is suggested that RcAlb-PepI, PepGAT, and PepKAA could interfere with SARS-CoV-2 Mpro role in vivo. Also, unlike other antiviral drugs, these peptides have no toxicity to human cells. This pioneering in silico investigation opens up opportunity for further in vivo research on these peptides, towards discovering new drugs and entirely new perspectives to treat COVID-19. Communicated by Ramaswamy H. Sarma

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.

ACE2-derived peptides interact with the RBD domain of SARS-CoV-2 spike glycoprotein, disrupting the interaction with the human ACE2 receptor

Vaccines could be the solution to the current SARS-CoV-2 outbreak. However, some studies have shown that the immunological memory only lasts three months. Thus, it is imperative to develop pharmacological treatments to cope with COVID-19. Here, the in silico approach by molecular docking, dynamic simulations and quantum biochemistry revealed that ACE2-derived peptides strongly interact with the SARS-CoV-2 RBD domain of spike glycoprotein (S-RBD). ACE2-Dev-PepI, ACE2-Dev-PepII, ACE2-Dev-PepIII and ACE2-Dev-PepIV complexed with S-RBD provoked alterations in the 3D structure of S-RBD, leading to disruption of the correct interaction with the ACE2 receptor, a pivotal step for SARS-CoV-2 infection. This wrong interaction between S-RBD and ACE2 could inhibit the entry of SARS-CoV-2 in cells, and thus virus replication and the establishment of COVID-19 disease. Therefore, we suggest that ACE2-derived peptides can interfere with recognition of ACE2 in human cells by SARS-CoV-2 in vivo. Bioinformatic prediction showed that these peptides have no toxicity or allergenic potential. By using ACE2-derived peptides against SARS-CoV-2, this study points to opportunities for further in vivo research on these peptides, seeking to discover new drugs and entirely new perspectives to treat COVID-19. Communicated by Ramaswamy H. Sarma

Quantum biochemistry in cancer immunotherapy: New insights about CTLA-4/ipilimumab and design of ipilimumab-derived peptides with high potential in cancer treatment

Cancer is a group of diseases involving disordered growth of abnormal cells with the potential to invade and spread to other parts of the body. Today, immunotherapy is the most efficient treatment, with fewer side effects. Notably, the employment of monoclonal antibodies to inhibit checkpoint proteins, such as CTLA-4, has caused much excitement among cancer immunotherapy researchers. Thus, in-depth analysis through quantum biochemistry and molecular dynamics simulations was performed to understand the complex formed by ipilimumab and its target CTLA-4. Our computational results provide a better understanding of the binding mechanisms and new insights about the CTLA-4: ipilimumab interaction, identifying essential amino acid residues to support the complex. Additionally, we report new interactions such as aromatic-aromatic, aromatic-sulfur, and cation-pi interactions to stabilize the CTLA-4:ipilimumab complex. Finally, quantum biochemistry analyses reveal the most important amino acid residues involved in the CTLA-4:ipilimumab interface, which were used to design synthetic peptides to inhibit CTLA-4. The computational results presented here provide a better understanding of the CTLA-4:ipilimumab binding mechanisms, and can support the development of alternative antibody-based drugs with high relevance in cancer immunotherapy.