Search

Inhibition of HIV-1 protease is a cornerstone of antiretroviral therapy. However, the notorious ability of HIV-1 to develop resistance to protease inhibitors (PIs), particularly darunavir (DRV), poses a major challenge. Using quantum chemistry and computer simulations, this study aims to investigate the interactions between two novel PIs, GRL-004 and GRL-063, as well as a wild-type (WT) HIV strain and a DRV-resistant mutant strain. To do this, we used molecular docking, molecular dynamics simulations, and quantum mechanical calculations to check how well GRL-004 and GRL-063 bound to both WT and DRV-resistant proteases. The results show that GRL-004 and GRL-063 bind very well to ASP29 in the WT strain. ASP29 is an important amino acid in the HIV protease dimer. Remarkably, amino acids such as ILE50 in the WT strains showed substantial binding energies to both drugs. Quantum energy calculations showed a slight reduction in the energy affinity of the interaction between the MUT strain and the ligand GRL-063, compared to the WT strain. GRL-004 showed similar interaction energy for both strains, suggesting that it has greater plasticity than GRL-063 despite its lower interaction affinity. Furthermore, GLY49B demonstrated strong binding energies regardless of mutations. Other relevant residues with strong binding energies include GLY49B, PHE82A, PRO81A, ASP29A, ASP25A and ALA28B. This study improves our understanding of receptor-ligand dynamics and the adaptability of new protease inhibitors (PIs), which has profound implications for the innovation of future antiretroviral drugs.

Exploring therapeutic options is crucial in the ongoing COVID-19 pandemic caused by SARS-CoV-2. Nirmatrelvir, which is a potent inhibitor that targets the SARS-CoV-2 M pro , shows promise as an antiviral treatment. Additionally, Ivermectin, which is a broad-spectrum antiparasitic drug, has demonstrated effectiveness against the virus in laboratory settings. However, its clinical implications are still debated. Using computational methods, such as molecular docking and 100 ns molecular dynamics simulations, we investigated how Nirmatrelvir and Ivermectin interacted with SARS-CoV-2 M pro(A) . Calculations using density functional theory were instrumental in elucidating the behavior of isolated molecules, primarily by analyzing the frontier molecular orbitals. Our analysis revealed distinct binding patterns: Nirmatrelvir formed strong interactions with amino acids, like MET49, MET165, HIS41, HIS163, HIS164, PHE140, CYS145, GLU166, and ASN142, showing stable binding, with a root-mean-square deviation (RMSD) of around 2.0 Å. On the other hand, Ivermectin interacted with THR237, THR239, LEU271, LEU272, and LEU287, displaying an RMSD of 1.87 Å, indicating enduring interactions. Both ligands stabilized M pro(A) , with Ivermectin showing stability and persistent interactions despite forming fewer hydrogen bonds. These findings offer detailed insights into how Nirmatrelvir and Ivermectin bind to the SARS-CoV-2 main protease, providing valuable information for potential therapeutic strategies against COVID-19.

Atomoxetine is a drug widely used for the treatment of the attention deficit hyperactivity disorder (ADHD) with reduced risk of adverse motor reactions and chemical dependence. However, the pharmacokinetics characteristics as well as the toxicological risk of atomoxetine deserves further investigation to comprehensively analyze the therapeutic and safety aspects of this drug. This study aimed to predict the physicochemical profile and medicinal chemistry characteristics of atomoxetine, alongside its pharmacokinetic properties-namely absorption, distribution, metabolism, and excretion-as well as its toxicology (ADMET) potential through the utilization of web-based in silico tools. This research emphasizes predicted physicochemical, medicinal chemistry, and absorption parameters of atomoxetine that could influence the efficacy and safety of this drug for ADHD treatment. Additionally, atomoxetine also presents noteworthy predicted risks of hepatotoxicity, cardiotoxicity, neurotoxicity, nephrotoxicity, respiratory system toxicity, skin toxicity, and carcinogenicity. These findings underscore the necessity for further assessments of atomoxetine's safety profile, particularly considering different patient populations and durations of drug treatment. The data reported here from in silico predictions suggest that closer monitoring is warranted when atomoxetine is administered to patients with ADHD. Moreover, controlled studies detailing reliable protocols for personalized dosing, considering the multifactorial variability in metabolism efficiency and toxicological potential, would enable a more comprehensive assessment of atomoxetine's safety profile., Competing Interests: No conflicts of interest declared concerning the publication of this article.

Yellow fever outbreaks are prevalent, particularly in endemic regions. Given the lack of an established treatment for this disease, significant attention has been directed toward managing this arbovirus. In response, we developed a multiepitope vaccine designed to elicit an immune response, utilizing advanced immunoinformatic and molecular modeling techniques. To achieve this, we predicted B- and T-cell epitopes using the sequences from all structural (E, prM, and C) and nonstructural proteins of 196 YFV strains. Through comprehensive analysis, we identified 10 cytotoxic T-lymphocyte (CTL) and 5T-helper (Th) epitopes that exhibited overlap with B-lymphocyte epitopes. These epitopes were further evaluated for their affinity to a wide range of human leukocyte antigen system alleles and were rigorously tested for antigenicity, immunogenicity, allergenicity, toxicity, and conservation. These epitopes were linked to an adjuvant ( β -defensin) and to each other using ligands, resulting in a vaccine sequence with appropriate physicochemical properties. The 3D structure of this sequence was created, improved, and quality checked; then it was anchored to the Toll-like receptor. Molecular Dynamics and Quantum Mechanics/Molecular Mechanics simulations were employed to enhance the accuracy of docking calculations, with the QM portion of the simulations carried out utilizing the density functional theory formalism. Moreover, the inoculation model was able to provide an optimal codon sequence that was inserted into the pET-28a( +) vector for in silico cloning and could even stimulate highly relevant humoral and cellular immunological responses. Overall, these results suggest that the designed multi-epitope vaccine can serve as prophylaxis against the yellow fever virus., (© 2024. The Author(s).)

In this computational study, we advanced the understanding of the antigenic properties of the NADC-34-like isolate of the Porcine Reproductive and Respiratory Syndrome Virus (PRRSV), named YC-2020, relevant in veterinary pathology. We utilized sequence comparison analyses of the M and N proteins, comparing them with those of NADC34, identifying substantial amino acid homology that allowed us to highlight conserved epitopes and crucial variants. Through the application of Clustal Omega for multiple sequence alignment and platforms like Vaxijen and AllerTOP for predicting antigenic and allergenic potential, our analyses revealed important insights into the conservation and variation of epitopes essential for the development of effective diagnostic tools and vaccines. Our findings, aligned with initial experimental studies, underscore the importance of these epitopes in the development of targeted immunodiagnostic platforms and significantly contribute to the management and control of PRRSV. However, further studies are required to validate the computational predictions of antigenicity for this new viral isolate. This approach underscores the potential of computational models to enable ongoing monitoring and control of PRRSV evolution in swine. While this study provides valuable insights into the antigenic properties of the novel PRRSV isolate YC-2020 through computational analysis, it is important to acknowledge the limitations inherent to in silico predictions, specifically, the absence of laboratory validation., Competing Interests: Declaration of Competing Interest The authors declare no potential conflicts of interest with respect to the research, authorship, and/or publication of the project via internal resources., (Copyright © 2024 Elsevier B.V. All rights reserved.)

Oropouche virus (OROV) is a member of the Peribunyaviridae family and the causative agent of a dengue-like febrile illness transmitted by mosquitoes. Although mild symptoms generally occur, complications such as encephalitis and meningitis may develop. A lack of proper diagnosis, makes it a potential candidate for new epidemics and outbreaks like other known arboviruses such as Dengue, Yellow Fever and Zika virus. The study of natural molecules as potential antiviral compounds is a promising alternative for antiviral therapies. Wedelolactone (WDL) has been demonstrated to inhibit some viral proteins and virus replication, making it useful to target a wide range of viruses. In this study, we report the in silico effects of WDL on the OROV N-terminal polymerase and its potential inhibitory effects on several steps of viral infection in mammalian cells in vitro , which revealed that WDL indeed acts as a potential inhibitor molecule against OROV infection., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2024 The Authors. Published by Elsevier B.V.)

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., (© 2024. The Author(s).)

The rise of pyrazinamide (PZA)-resistant strains of Mycobacterium tuberculosis (MTB) poses a major challenge to conventional tuberculosis (TB) treatments. PZA, a cornerstone of TB therapy, must be activated by the mycobacterial enzyme pyrazinamidase (PZase) to convert its active form, pyrazinoic acid, which targets the ribosomal protein S1. Resistance, often associated with mutations in the RpsA protein, complicates treatment and highlights a critical gap in the understanding of structural dynamics and mechanisms of resistance, particularly in the context of the G97D mutation. This study utilizes a novel integration of computational techniques, including multiscale biomolecular and molecular dynamics simulations, physicochemical and medicinal chemistry predictions, quantum computations and virtual screening from the ZINC and Chembridge databases, to elucidate the resistance mechanism and identify lead compounds that have the potential to improve treatment outcomes for PZA-resistant MTB, namely ZINC15913786, ZINC20735155, Chem10269711, Chem10279789 and Chem10295790. These computational methods offer a cost-effective, rapid alternative to traditional drug trials by bypassing the need for organic subjects while providing highly accurate insight into the binding sites and efficacy of new drug candidates. The need for rapid and appropriate drug development emphasizes the need for robust computational analysis to justify further validation through in vitro and in vivo experiments., (© 2024 The Authors. Journal of Cellular and Molecular Medicine published by Foundation for Cellular and Molecular Medicine and John Wiley & Sons Ltd.)

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.

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Lima Neto, Bezerra, Barbosa, Araujo, Galvão, Lyra, Oliveira, Akash, Jardan, Nafidi, Bourhia and Fulco.)

Oropouche virus (OROV) is an emerging vector-borne arbovirus found in South America that causes Oropouche fever, a febrile infection similar to dengue fever. It has a high epidemic potential, causing illness in over 500,000 cases diagnosed since the virus was first discovered in 1955. Currently, the prevention of human viral infection depends on vaccination, but availability for many viruses is limited, and they are classified as neglected viruses. At present, there are no vaccines or antiviral treatments available. An alternative approach to limiting the spread of the virus is to selectively disrupt viral replication mechanisms. Here, we demonstrate the inhibitory effect of acridones, which efficiently inhibited viral replication by 99.9 % in vitro . To evaluate possible mechanisms of action, we conducted tests with dsRNA, an intermediate in virus replication, as well as MD simulations, docking, and binding free energy analysis. The results showed a strong interaction between FAC21 and the OROV endonuclease, which possibly limits the interaction of viral RNA with other proteins. Therefore, our results suggest a dual mechanism of antiviral action, possibly caused by ds-RNA intercalation. In summary, our findings demonstrate that a new generation of antiviral drugs could be developed based on the selective optimization of molecules., Competing Interests: The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Mauricio Lacerda Nogueira reports financial support was provided by The Coordinating Research on Emerging Arboviral Threats Encompassing the Neotropics., (© 2023 The Author(s).)

The Oropouche virus (OROV) is a member of the family Peribunyaviridae (order Bunyavirales ) and the cause of a dengue-like febrile illness transmitted mainly by biting midges and mosquitoes. In this study, we aimed to explore acylphloroglucinols and xanthohumol from hops ( Humulus lupulus L.) as a promising alternative for antiviral therapies. The evaluation of the inhibitory potential of hops compounds on the viral cycle of OROV was performed through two complementary approaches. The first approach applies cell-based assay post-inoculation experiments to explore the inhibitory potential on the latest steps of the viral cycle, such as genome translation, replication, virion assembly, and virion release from the cells. The second part covers in silico methods evaluating the ability of those compounds to inhibit the activity of the endonuclease domain, which is essential for transcription, binding, and cleaving RNA. In conclusion, the beta acids showed strongest inhibitory potential in post-treatment assay (EC 50 = 26.7 µg/mL). Xanthohumol had the highest affinity for OROV endonuclease followed by colupulone and cohumulone. This result contrasts with that observed for docking and MM/PBSA analysis, where cohumulone was found to have a higher affinity. Finally, among the three tested ligands, Lys92 and Arg33 exhibited the highest affinity with the protein.

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, Pf DHODH 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 Pf DHODH 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.

mRNA vaccines are a new class of vaccine that can induce potent and specific immune responses against various pathogens. However, the design of mRNA vaccines requires the identification and optimization of suitable antigens, which can be challenging and time consuming. Reverse vaccinology is a computational approach that can accelerate the discovery and development of mRNA vaccines by using genomic and proteomic data of the target pathogen. In this article, we review the advances of reverse vaccinology for mRNA vaccine design against SARS-CoV-2, the causative agent of COVID-19. We describe the steps of reverse vaccinology and compare the in silico tools used by different studies to design mRNA vaccines against SARS-CoV-2. We also discuss the challenges and limitations of reverse vaccinology and suggest future directions for its improvement. We conclude that reverse vaccinology is a promising and powerful approach to designing mRNA vaccines against SARS-CoV-2 and other emerging pathogens.

Anabolic androgenic steroids (AAS) are substances with androgenic and anabolic characteristics. Among the many side effects of hormone therapy with AAS, the following stand out: heart problems, adrenal gland disorders, aggressive behavior, increased risk of prostate cancer, problems related to lack of libido and impotence. Such substances vary in the relationship between androgenic activity, and the activation of the androgen receptor (AR) is of fundamental importance for the singularity of the action of each AAS. In this sense, our study evaluates the aspects that comprise the interactions of testosterone agonists (TES), dihydrotestosterone (DHT) and tetrahydrogestrinone (THG) in complex with the AR. In addition, we also evaluated the impact of ligand-receptor affinity differences in a mutation model. We apply computational techniques based on density functional theory (DFT) and use, as methodology, Molecular Fractionation with Conjugate Caps (MFCC). The energetic specificities present in the interaction between the analyzed complexes attest that the highest affinity with the AR receptor is found for AR-THG, followed by AR-DHT, AR-TES and AR-T877A-DHT, respectively. Our results also show the differences and equivalences between the different agonists, in addition to evaluating the difference between the DHT ligand in complex with the wild-type and mutant receptor, presenting the main amino acid residues that involve the interaction with the ligands. The computational methodology used proves to be an operative and sophisticated choice to help in the search for pharmacological agents for various therapies that have androgen as a target.

Mayaro virus (MAYV) is an arbovirus found in the Americas that can cause debilitating arthritogenic disease. Although it is an emerging virus, the only current approach is vector control, as there are no approved vaccines to prevent MAYV infection nor therapeutics to treat it. In search of an effective vaccine candidate against MAYV, we used immunoinformatics and molecular modeling to attempt to identify promiscuous T-cell epitopes of the nonstructural polyproteins (nsP1, nsP2, nsP3, and nsP4) from 127 MAYV genomes sequenced in the Americas (08 Bolivia, 72 Brazil, 04 French Guiana, 05 Haiti, 20 Peru, 04 Trinidad and Tobago, and 14 Venezuela). For this purpose, consensus sequences of 360 proteins were used to identify short protein sequences that can bind to MHC I class (MHC II). Our analysis revealed 56 potential MHC-I/TCD8+ (29 MHC-II/TCD4+) epitopes, but only 6 (16) TCD8+ (TCD4+) epitopes showed high antigenicity and conservation, non-allergenicity, non-toxicity, and excellent population coverage. Finally, classical and quantum mechanical calculations (QM:MM) were used to improve the quality of the docking calculations, with the QM part of the simulations performed using the density functional theory formalism (DFT). These results provide insights for the advancement of diagnostic platforms, vaccine development, and immunotherapeutic interventions.Communicated by Ramaswamy H. Sarma.

The current COVID-19 pandemic, an infectious disease caused by the novel coronavirus (SARS-CoV-2), poses a threat to global health because of its high rate of spread and death. Currently, vaccination is the most effective method to prevent the spread of this disease. In the present study, we developed a novel multiepitope vaccine against SARS-CoV-2 containing Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1), Delta (B.1.617.2), and Omicron (BA.1) variants. To this end, we performed a robust immunoinformatics approach based on multiple epitopes of the four structural proteins of SARS-CoV-2 (S, M, N, and E) from 475 SARS-CoV-2 genomes sequenced from the regions with the highest number of registered cases, namely the United States, India, Brazil, France, Germany, and the United Kingdom. To investigate the best immunogenic epitopes for linear B cells, cytotoxic T lymphocytes (CTL), and helper T lymphocytes (HTL), we evaluated antigenicity, allergenicity, conservation, immunogenicity, toxicity, human population coverage, IFN-inducing, post-translational modifications, and physicochemical properties. The tertiary structure of a vaccine prototype was predicted, refined, and validated. Through docking experiments, we evaluated its molecular coupling to the key immune receptor Toll-Like Receptor 3 (TLR3). To improve the quality of docking calculations, quantum mechanics/molecular mechanics calculations (QM/MM) were used, with the QM part of the simulations performed using the density functional theory formalism (DFT). Cloning and codon optimization were performed for the successful expression of the vaccine in E. coli. Finally, we investigated the immunogenic properties and immune response of our SARS-CoV-2 multiepitope vaccine. The results of the simulations show that administering our prototype three times significantly increases the antibody response and decreases the amount of antigens. The proposed vaccine candidate should therefore be tested in clinical trials for its efficacy in neutralizing SARS-CoV-2., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Elsevier Ltd. All rights reserved.)

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., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

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

Competing Interests: Declaration of interests The authors declare no competing interests.

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.

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

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.

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.

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., (Copyright © 2021. Published by Elsevier B.V.)

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., (© The Author(s) 2021. Published by Oxford University Press.)

Melanoma is a type of skin cancer with increasing incidence worldwide and high lethality. Conventional forms of treatment are not effective in advanced cancer stages. Hence, immunotherapeutic approaches have been tested to modulate immune response against tumor cells. Some vaccine models using tumor-associated antigens (TAAs) such as glycoprotein 100 (gp100) have been studied, but their expected effectiveness has not been shown until now. Antigen immunogenicity is a crucial point to improve the immune response, and therefore mutations are inserted in peptide sequences. It is possible to understand the interactions which occur between peptides and immune system molecules through computer simulation, and this is essential in order to guide efficient vaccine models. In this work, we have calculated the interaction binding energies of crystallographic data based on modified gp100 peptides and HLA-A*0201 using density functional theory (DFT) and the molecular fractionation with conjugated caps (MFCC) approach. Our results show the most relevant residue-residue interactions, the impact of three mutations in their binding sites, and the main HLA-A*0201 amino acids for peptide-HLA binding.

The spreading of epidemics in complex networks has been a subject of renewed interest of several scientific branches. In this regard, we have focused our attention on the study of the susceptible-infected-susceptible (SIS) model, within a Monte Carlo numerical simulation approach, representing the spreading of epidemics in a clustered homophilic network. The competition between infection and recovery that drives the system either to an absorbing or to an active phase is analyzed. We estimate the static critical exponents β ∕ ν , 1 ∕ ν and γ ∕ ν , through finite-size scaling (FSS) analysis of the order parameter ρ and its fluctuations, showing that they differ from those associated with the contact process on a scale-free network, as well as those predicted by the heterogeneous mean-field theory., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2020 Elsevier B.V. All rights reserved.)

The worldwide dependence of population on fossil fuels continues to have several harmful implications for the environment. Bioethanol is an excellent option for renewable fuel to replace the current greenhouse gas emitters. In addition, its production by enzymatic route has gained space among the industrial processes because it replaces the traditional acid treatment. Due to its high versatility, the xylanase family is used in this process as an accessory enzyme for degrading the lignocellulosic substrate of biomass. A chimera built by a xylanolytic domain (Xyl) and a xylose-binding protein (XBP) showed an experimentally improved catalytic efficiency and interdomain allosteric modulation after xylose binding. In this context, we performed a quantum biochemistry characterization of the interactions between these domains and dynamic cross-correlation (DCC) analysis after performing molecular dynamics (DM) simulations of the systems in the presence and absence of xylose in the XBP active site. We used the density functional theory (DFT) within the molecular fractionation with the conjugated caps (MFCC) approach to describe the pair energies, and the corresponding energy difference between the chimera domains responsible for the allosteric effect and amino acid DCC to evaluate the interdomain coupling differences between the energy states. The detailed energetic investigation together with the related structural and dynamics counterparts revealed the molecular mechanisms of chimeric improvement of the xylanase activity observed experimentally. This mechanism was correlated with greater stability and high connectivity at the interdomain interface in the xylose bound relative to the free chimera. We identify the contributions of hydrogen bonds, hydrophobic interactions and water-mediated interactions in the interdomain region responsible for stability together with the structural and dynamical elements related to the allosteric effect. Taken together, these observations led to a comprehensive understanding of the chimera's modulatory action that occurs through the formation of a highly connected interface that makes the essential movements related to xylanolytic activity in xylanase correlated to those of the xylose-binding protein.

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.

GABA B is a G protein-coupled receptor that functions as a constitutive heterodimer composed of the GABA B1a/b and GABA B2 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 GABA B 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 GABA B 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 GABA B -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 GABA B 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 GABA B -ligand binding. These results provide important new information on the binding mechanism of the GABA B receptor and should facilitate the development of new chemicals targeting this receptor.

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.

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.

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 IC 50  = 483.5 and 714.9 μM, whereas, in primary culture tests using murine macrophages, IC 50 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.

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.

Many species of bacteria involved in degradation of n-alkanes have an important constitutional metabolic enzyme, the alkane hydroxylase called AlkB, specialized in the conversion of hydrocarbons molecules that can be used as carbon and/or energy source. This enzyme plays an important role in the microbial degradation of oil, chlorinated hydrocarbons, fuel additives, and many other compounds. A number of these enzymes has been biochemically characterized in detail because the potential of alkane hydroxylases to catalyse high added-value reactions is widely recognized. Nevertheless, the industrial and process bioremediation application of them is restricted, owing to their complex biochemistry, challenging process requirements, and the limited number of their three-dimensional structures. Furthermore, AlkB has great potential as biocatalysts for selective transformation of a wide range of chemically inert unreactive alkanes into reactive chemical precursors that can be used as tools for bioremediation and bioprocesses. Aiming to understand the possible ways the AlkB enzyme Pseudomonas putida P1 interacts with octane, octanol and 1-octyne, we consider its suitable biochemical structure taking into account a 3-D homology modelling. Besides, by using a quantum chemistry computational model based on the density functional theory (DFT), we determine possible protein-substrate interaction regions measured by means of its binding energy simulated throughout the Molecular Fractionation with Conjugated Caps (MFCC) approach., (Copyright © 2017 Elsevier Inc. All rights reserved.)

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 chemistry study were carried out looking for the conformational, optoelectronic and vibrational properties of the entacapone molecule, an efficient drug used in the Parkinson's disease treatment. Classical annealing was performed to explore the entacapone's molecular configurations, searching for optimal geometries. The quantum optimization calculations were made using three different functional combination levels of the density functional theory (DFT). The structural data (bond length, bond and torsion angles), charge population analysis (absorption spectra) and molecular orbital study (HOMO and LUMO) were obtained considering the lower energy optimized conformation of the entacapone molecule. Furthermore, a complete assignment of the harmonic vibrational frequencies were achieved through their infrared (IR) and Raman spectra.

The infrared absorption and Raman scattering spectra of the monoclinic P21 l-aspartic acid anhydrous crystal were recorded and interpreted with the help of density functional theory (DFT) calculations. The effect of dispersive forces was taken into account, and the optimized unit cells allowed us to obtain the vibrational normal modes. The computed data exhibits good agreement with the measurements for low wavenumbers, allowing for a very good assignment of the infrared and Raman spectral features. The vibrational spectra of the two lowest energy conformers of the l-aspartic molecule were also evaluated using the hybrid B3LYP functional for the sake of comparison, showing that the molecular calculations give a limited description of the measured IR and Raman spectra of the l-aspartic acid crystal for wavenumbers below 1000 cm(-1). The results obtained reinforce the need to use solid-state calculations to describe the vibrational properties of molecular crystals instead of calculations for a single isolated molecule picture even for wavenumbers beyond the range usually associated with lattice modes (200 cm(-1) < ω < 1000 cm(-1)).

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.

We investigate the electronic wavepacket dynamics in a finite segment of a DNA single-strand chain considering the electron-phonon coupling. Our theoretical approach makes use of an effective tight-binding Hamiltonian to describe the electron dynamics, together with a classical harmonic Hamiltonian to treat the intrinsic DNA vibrations. An effective time-dependent Schrödinger equation is then settled up and solved numerically for an initially localized wave-packet using the standard Dormand-Prince eighth-order Runge-Kutta method. Our numerical results indicate the presence of a sub-diffusive electronic wavepacket spread mediated by the electron-phonon interaction.

We investigate the critical properties of Ising models on a regularized Apollonian network (RAN), here defined as a kind of Apollonian network in which the connectivity asymmetry associated with its corners is removed. Different choices for the coupling constants between nearest neighbors are considered and two different order parameters are used to detect the critical behavior. While ordinary ferromagnetic and antiferromagnetic models on a RAN do not undergo a phase transition, some antiferrimagnetic models show an interesting infinite-order transition. All results are obtained by an exact analytical approach based on iterative partial tracing of the Boltzmann factor as an intermediate step for the calculation of the partition function and the order parameters.

We introduce a superdiffusive one-dimensional epidemic process model on which infection spreads through a contact process. Healthy (A) and infected (B) individuals can jump with distinct probabilities D(A) and D(B) over a distance ℓ distributed according to a power-law probability P(ℓ)[proportionality]1/ℓ(μ). For μ≥3 the propagation is equivalent to diffusion, while μ<3 corresponds to Lévy flights. In the D(A)>D(B) diffusion regime, field-theoretical results have suggested a first-order transition, a prediction not supported by several numerical studies. An extensive numerical study of the critical behavior in both the diffusive (μ≥3) and superdiffusive (μ<3) D(A)>D(B) regimes is also reported. We employed a finite-size scaling analysis to obtain the critical point as well as the static and dynamic critical exponents for several values of μ. All data support a second-order phase transition with continuously varying critical exponents which do not belong to the directed percolation universality class.