Quantum biochemical analysis of the TtgR regulator and effectors.

The recent expansion of multidrug-resistant (MDR) pathogens poses significant challenges in treating healthcare-associated infections. Although antibacterial resistance occurs by numerous mechanisms, active efflux of the drugs is a critical concern. A single species of efflux pump can produce a simultaneous resistance to several drugs. One of the best-studied efflux pumps is the TtgABC: a tripartite resistance-nodulation-division (RND) efflux pump implicated in the intrinsic antibiotic resistance in Pseudomonas putida DOT-T1E. The expression of the TtgABC gene is down-regulated by the HTH-type transcriptional repressor TtgR. In this context, by employing quantum chemistry methods based on the Density Functional Theory (DFT) within the Molecular Fragmentation with Conjugate Caps (MFCC) approach, we investigate the coupling profiles of the transcriptional regulator TtgR in complex with quercetin (QUE), a natural polyphenolic flavonoid, tetracycline (TAC), and chloramphenicol (CLM), two broad-spectrum antimicrobial agents. Our quantum biochemical computational results show the: [i] convergence radius, [ii] total binding energy, [iii] relevance (energetically) of the ligands regions, and [iv] most relevant amino acids residues of the TtgR-QUE/TAC/CLM complexes, pointing out distinctions and similarities among them. These findings improve the understanding of the binding mechanism of effectors and facilitate the development of new chemicals targeting TtgR, helping in the battle against the rise of resistance to antimicrobial drugs. These advances are crucial in the ongoing fight against rising antimicrobial drug resistance, providing hope for a future where healthcare-associated infections can be more beneficially treated.

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.

New ethionamide boosters and EthR2: structural and energetic analysis.

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

Blockade of the checkpoint PD-1 by its ligand PD-L1 and the immuno-oncological drugs pembrolizumab and nivolumab.

We investigate the interaction between the programmed cell death protein 1 (PD-1) and the programmed cell death ligand 1 (PD-L1), as well as the immuno-oncological drugs pembrolizumab (PEM), and nivolumab (NIV), through quantum chemistry methods based on the Density Functional Theory (DFT) and the molecular fractionation with conjugate caps (MFCC) scheme, in order to map their hot-spot regions. Our results showed that the total interaction energy order of the three complexes is in good agreement with the experimental binding affinity order: PD-1/PEM > PD-1/NIV > PD-1/PD-L1. Besides, a detailed investigation revealed the energetically most relevant residue-residue pairs-interaction for each complex. Our computational results give a better understanding of the interaction mechanism between the protein PD-1 and its ligands (natural and inhibitors), unleashing the immune surveillance to destroy the cancer cells by decreasing their immune evasion. They are also an efficient alternative towards the development of new small-molecules and antibody-based drugs, pointing out to new treatments for cancer therapy.

Exploring human porphobilinogen synthase metalloprotein by quantum biochemistry and evolutionary methods.

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

In silico approach of modified melanoma peptides and their immunotherapeutic potential.

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.

Critical properties of the SIS model on the clustered homophilic network.

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.

Molecular modelling and quantum biochemistry computations of a naturally occurring bioremediation enzyme: Alkane hydroxylase from Pseudomonas putida P1.

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.

Conformational, Optoelectronic and Vibrational Properties of the Entacapone Molecule: A Quantum Chemistry Study.

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.

Vibrational Spectroscopy and Phonon-Related Properties of the L-Aspartic Acid Anhydrous Monoclinic Crystal.

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

Sub-diffusive electronic transport in a DNA single-strand chain with electron-phonon coupling.

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.

Ising models on the regularized Apollonian network.

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.

Critical properties of a superdiffusive epidemic process.

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.

The two-dimensional site-diluted Ising model: a short-time-dynamics approach.

The site-diluted Ising ferromagnet is investigated on a square lattice, within short-time-dynamics numerical simulations, for different site concentrations. The dynamical exponents θ and z are obtained and it is shown that these exponents do depend strongly on the disorder, exhibiting a clear breakdown of universality, characterized by relative variations of nearly 100% in the range of site concentrations investigated. In what concerns the static exponents ß and ν, universality is preserved within the error bars.

Criticality of a contact process with coupled diffusive and non-diffusive fields.

We investigate the critical behavior of a model with two coupled critical densities, one of which is diffusive. The model simulates the propagation of an epidemic process in a population, which uses the underlying lattice to leave a track of the recent disease history. We determine the critical density of the population above which the system reaches an active stationary state with a finite density of active particles. We also perform a scaling analysis to determine the order parameter, the correlation length, and critical relaxation exponents. We show that the model does not belong to the usual directed percolation universality class and is compatible with the class of directed percolation with diffusive and conserved fields.

Majority-vote model on a random lattice.

The stationary critical properties of the isotropic majority vote model on random lattices with quenched connectivity disorder are calculated by using Monte Carlo simulations and finite size analysis. The critical exponents gamma and beta are found to be different from those of the Ising and majority vote on the square lattice model and the critical noise parameter is found to be q(c) =0.117+/-0.005 .

Optimizing the encounter rate in biological interactions: Lévy versus Brownian strategies.

An important application involving two-species reaction-diffusion systems relates to the problem of finding the best statistical strategy for optimizing the encounter rate between organisms. We investigate the general problem of how the encounter rate depends on whether organisms move in Lévy or Brownian random walks. By simulating a limiting generalized searcher-target model (e.g., predator-prey, mating partner, pollinator-flower), we find that Lévy walks confer a significant advantage for increasing encounter rates when the searcher is larger or moves rapidly relative to the target, and when the target density is low.

Critical behavior of a one-dimensional diffusive epidemic process.

We investigate the critical behavior of a one-dimensional diffusive epidemic propagation process by means of a Monte Carlo procedure. In the model, healthy (A) and sick (B) individuals diffuse on a lattice with diffusion constants D(A) and D(B), respectively. According to a Wilson renormalization calculation, the system presents a second-order phase transition between a steady reactive state and a vacuum state, with distinct universality classes for the cases D(A)=D(B) and D(A)D(B). In this work we perform a finite size scaling analysis of order parameter data at the vicinity of the critical point in dimension d=1. Our results show no signature of a first-order transition in the case of D(A)>D(B). A finite size scaling typical of second-order phase transitions fits well the data from all three regimes. We found that the correlation exponent nu=2 as predicted by field-theoretical arguments. Estimates for beta/nu are given for all relevant regimes.