Combined computational and experimental study about the incorporation of phosphorus into the structure of graphene oxide.

Phosphorus-containing graphene-based hybrids are materials with outstanding properties for diverse applications. In this work, an easy route to produce phosphorus-graphene oxide hybrid materials is described, involving the use of variable amounts of H3PO4 and H2SO4 during the reaction of oxidation of a graphitic precursor. The physical and chemical features of the hybrids change significantly with the variation in the acid amounts used in the syntheses. XPS and solid-state 13C and 31P NMR results show that the hybrids contain large amounts of oxygen functional groups, with the phosphorus incorporation proceeding mostly through the formation of phosphate-like linkages and other functions with C-O-P bonds. The experimental findings are supported by DFT calculations, which allow the assessment of the energetics and the geometry of the interaction between phosphate groups and graphene-based models; these calculations are also used to predict the chemical shifts in the 31P and 13C NMR spectra of the models, which show good agreement with the experimentally observed solid-state NMR spectra.

Modeling studies on the role of vitamins B1 (thiamin), B3 (nicotinamide), B6 (pyridoxamine), and caffeine as potential leads for the drug design against COVID-19.

In response to the COVID-19 pandemic, and the lack of effective and safe antivirals against it, we adopted a new approach in which food supplements with vital antiviral characteristics, low toxicity, and fast excretion have been targeted. The structures and chemical properties of the food supplements were compared to the promising antivirals against SARS-COV-2. Our goal was to exploit the food supplements to mimic the topical antivirals' functions but circumventing their severe side effects, which has limited the necessary dosage needed to exhibit the desired antiviral activity. On this line, after a comparative structural analysis of the chemicals mentioned above, and investigation of their potential mechanisms of action, we selected caffeine and some compounds of the vitamin B family and further applied molecular modeling techniques to evaluate their interactions with the RDB domain of the Spike protein of SARS-CoV-2 (SC2Spike) and its corresponding binding site on human ACE-2 (HssACE2). Our results pointed to vitamins B1 and B6 in the neutral form as potential binders to the HssACE2 RDB binding pocket that might be able to impair the SARS-CoV-2 mechanism of cell invasion, qualifying as potential leads for experimental investigation against COVID-19.

Improved Protocol for the Selection of Structures from Molecular Dynamics of Organic Systems in Solution: The Value of Investigating Different Wavelet Families.

Wavelets are mathematical tools used to decompose and represent another function described in the time domain, allowing the study of each component of the original function with a scale-compatible resolution. Thus, these transforms have been used to select conformations from molecular dynamics (MD) trajectories in systems of fundamental and technological interest. Recently, our research group has used wavelets to develop and validate a method, meant to select structures from MD trajectories, which we named OWSCA (optimal wavelet signal compression algorithm). Here, we moved forward on this project by demonstrating the efficacy of this method on the study of three different systems (non-flexible organic, flexible organic, and protein). For each system, 93 wavelets were investigated to verify which is the best one for a given organic system. The results show that the best wavelets were different for each system and, also, very close to the experimental values, with the wavelets db1, rbio 3.1, and bior1.1 being selected for the non-flexible, flexible organic, and protein systems, respectively. This reinforces our OWSCA as a very efficient and promising method for the selection of structures from MD trajectories of different classes of compounds. Our findings also point out that additional studies considering wavelet families are needed for defining the best wavelet for representing each system under study.

Are the current commercially available oximes capable of reactivating acetylcholinesterase inhibited by the nerve agents of the A-series?

The misuse of novichok agents in assassination attempts has been reported in the international media since 2018. These relatively new class of neurotoxic agents is claimed to be more toxic than the agents of the G and V series and so far, there is no report yet in literature about potential antidotes against them. To shed some light into this issue, we report here the design and synthesis of NTMGMP, a surrogate of A-242 and also the first surrogate of a novichok agent useful for experimental evaluation of antidotes. Furthermore, the efficiency of the current commercial oximes to reactivate NTMGMP-inhibited acetylcholinesterase (AChE) was evaluated. The Ellman test was used to confirm the complete inhibition of AChE, and to compare the subsequent rates of reactivation in vitro as well as to evaluate aging. In parallel, molecular docking, molecular dynamics and MM-PBSA studies were performed on a computational model of the human AChE (HssAChE)/NTMGMP complex to assess the reactivation performances of the commercial oximes in silico. Experimental and theoretical studies matched the exact hierarchy of efficiency and pointed to trimedoxime as the most promising commercial oxime for reactivation of AChE inhibited by A-242.

Theoretical assessment of the performances of commercial oximes on the reactivation of acetylcholinesterase inhibited by the nerve agent A-242 (novichok).

The nerve agents of the A-series are relatively recent chemical weapons with no antidote available yet. Once inside the human body, those chemicals act similarly to the classic nerve agents, by binding to the catalytic residue Serine 203 (Ser203) of human acetylcholinesterase (HssAChE) and thus preventing the proper function of this enzyme. However, there is no experimental evidence yet if the current antidotes for intoxication by nerve agents are also capable of restoring AChE inhibited by the nerve agents of the A-series. In order to launch some light on this issue, we used computational techniques (molecular docking, molecular dynamics and MM-PBSA interaction energy calculations) to assess the performances of the four currently available commercial oximes (2-PAM, HI-6, obidoxime and trimedoxime) when in contact with HssAChE inhibited by the agent A-242. Based on the near-attack conformation (NAC) criterion, our results suggest that the commercial oximes would have limited efficacy to reactivate the enzyme since they are not able to properly approach the adduct Ser203-A-242. Among those oximes, trimedoxime seems to be the most promising, since it showed lower values of energy in the MM-PBSA calculations, a higher stability inside the catalytic anionic center (CAS) of HssAChE, and was able to adopt a position closer to the NAC that could enable the reactivation mechanism.

Docking and molecular dynamics studies of potential new leads against DBL3x derived from chondroitin sulfate A (CSA): a new approach for the treatment of malaria.

In this work the DBL3x domain of the erythrocyte membrane protein from Plasmodium Falciparum (PfEMP1), was revisited as a potential molecular target for the development of new drugs against malaria. This protein interacts with chondroitin sulfate A (CSA), a glycosaminoglycan present in the substance fundamental for connective tissues of vertebrates and is implicated in malaria complications in pregnant women. We performed molecular docking and molecular dynamic studies of DBL3x complexed with CSA and five analogues, where the sulfate group was replaced by phosphate, in order to evaluate if the better electrostatic interactions provided by phosphate groups could afford better binders capable of preventing the binding of CSA to DBL3x. Results suggest that all proposed compounds have high affinity towards DBL3x and could bind better to the DBL3x domain of PfEMP1 than CSA, qualifying as potential inhibitors of this protein and, therefore, new potential leads for the drug design against malaria.Communicated by Ramaswamy H. Sarma.

Identification of novel potential ricin inhibitors by virtual screening, molecular docking, molecular dynamics and MM-PBSA calculations: a drug repurposing approach.

Ricin is a potent cytotoxin with no available antidote. Its catalytic subunit, RTA, damages the ribosomal RNA (rRNA) of eukaryotic cells, preventing protein synthesis and eventually leading to cell death. The combination between easiness of obtention and high toxicity turns ricin into a potential weapon for terrorist attacks, urging the need of discovering effective antidotes. On this context, we used computational techniques, in order to identify potential ricin inhibitors among approved drugs. Two libraries were screened by two different docking algorithms, followed by molecular dynamics simulations and MM-PBSA calculations in order to corroborate the docking results. Three drugs were identified as potential ricin inhibitors: deferoxamine, leucovorin and plazomicin. Our calculations showed that these compounds were able to, simultaneously, form hydrogen bonds with residues of the catalytic site and the secondary binding site of RTA, qualifying as potential antidotes against intoxication by ricin.Communicated by Ramaswamy H. Sarma.

Biodegradation of Organophosphorus Compounds Predicted by Enzymatic Process Using Molecular Modelling and Observed in Soil Samples Through Analytical Techniques and Microbiological Analysis: A Comparison.

Organophosphorus compounds (OP) are chemicals widely used as pesticides in different applications such as agriculture and public health (vector control), and some of the highly toxic forms have been used as chemical weapons. After application of OPs in an environment, they persist for a period, suffering a degradation process where the biotic factors are considered the most relevant forms. However, to date, the biodegradation of OP compounds is not well understood. There are a plenty of structure-based biodegradation estimation methods, but none of them consider enzymatic interaction in predicting and better comprehending the differences in the fate of OPs in the environment. It is well known that enzymatic processes are the most relevant processes in biodegradation, and that hydrolysis is the main pathway in the natural elimination of OPs in soil samples. Due to this, we carried out theoretical studies in order to investigate the interactions of these OPs with a chosen enzyme-the phosphotriesterase. This one is characteristic of some soils' microorganisms, and has been identified as a key player in many biodegradation processes, thanks to its capability for fast hydrolyzing of different OPs. In parallel, we conducted an experiment using native soil in two conditions, sterilized and not sterilized, spiked with specific amounts of two OPs with similar structure-paraoxon-ethyl (PXN) and O-(4-nitrophenyl) O-ethyl methylphosphonate (NEMP). The amount of OP present in the samples and the appearance of characteristic hydrolysis products were periodically monitored for 40 days using analytical techniques. Moreover, the number of microorganisms present was obtained with plate cell count. Our theoretical results were similar to what was achieved in experimental analysis. Parameters calculated by enzymatic hydrolysis were better for PXN than for NEMP. In soil, PXN suffered a faster hydrolysis than NEMP, and the cell count for PXN was higher than for NEMP, highlighting the higher microbiological toxicity of the latter. All these results pointed out that theoretical study can offer a better comprehension of the possible mechanisms involved in real biodegradation processes, showing potential in exploring how biodegradation of OPs relates with enzymatic interactions.

Molecular Modeling and In Vitro Studies of a Neutral Oxime as a Potential Reactivator for Acetylcholinesterase Inhibited by Paraoxon.

The present work aimed to compare the small, neutral and monoaromatic oxime, isatin-3-oxime (isatin-O), to the commercial ones, pralidoxime (2-PAM) and obidoxime, in a search for a new potential reactivator for acetylcholinesterase (AChE) inhibited by the pesticide paraoxon (AChE/POX) as well as a novel potential scaffold for further synthetic modifications. The multicriteria decision methods (MCDM) allowed the identification of the best docking poses of those molecules inside AChE/POX for further molecular dynamic (MD) studies, while Ellman's modified method enabled in vitro inhibition and reactivation assays. In corroboration with the theoretical studies, our experimental results showed that isatin-O have a reactivation potential capable of overcoming 2-PAM at the initial moments of the assay. Despite not achieving better results than obidoxime, this molecule is promising for being an active neutral oxime with capacity of crossing the blood⁻brain barrier (BBB), to reactivate AChE/POX inside the central and peripheral nervous systems. Moreover, the fact that isatin-O can also act as anticonvulsant makes this molecule a possible multipotent reactivator. Besides, the MCDM method showed to be an accurate method for the selection of the best docking poses generated in the docking studies.

Docking and molecular dynamics studies of peripheral site ligand-oximes as reactivators of sarin-inhibited human acetylcholinesterase.

In the present work, we performed docking and molecular dynamics simulations studies on two groups of long-tailored oximes designed as peripheral site binders of acetylcholinesterase (AChE) and potential penetrators on the blood brain barrier. Our studies permitted to determine how the tails anchor in the peripheral site of sarin-inhibited human AChE, and which aminoacids are important to their stabilization. Also the energy values obtained in the docking studies corroborated quite well with the experimental results obtained before for these oximes.

Molecular modeling toward selective inhibitors of dihydrofolate reductase from the biological warfare agent Bacillus anthracis.

In the present work, we applied docking and molecular dynamics techniques to study 11 compounds inside the enzymes dihydrofolate reductase (DHFR) from the biological warfare agent Bacillus anthracis (BaDHFR) and Homo sapiens sapiens (HssDHFR). Six of these compounds were selected for a study with the mutant BaF96IDHFR. Our results corroborated with experimental data and allowed the proposition of a new molecule with potential activity and better selectivity for BaDHFR.

Design of inhibitors of thymidylate kinase from Variola virus as new selective drugs against smallpox.

Recently we constructed a homology model of the enzyme thymidylate kinase from Variola virus (VarTMPK) and proposed it as a new target to the drug design against smallpox. In the present work, we used the antivirals cidofovir and acyclovir as reference compounds to choose eleven compounds as leads to the drug design of inhibitors for VarTMPK. Docking and molecular dynamics (MD) studies of the interactions of these compounds inside VarTMPK and human TMPK (HssTMPK) suggest that they compete for the binding region of the substrate and were used to propose the structures of ten new inhibitors for VarTMPK. Further docking and MD simulations of these compounds, inside VarTMPK and HssTMPK, suggest that nine among ten are potential selective inhibitors of VarTMPK.

Evidence of ternary complex formation in Trypanosoma cruzi trans-sialidase catalysis.

Trypanosoma cruzi trans-sialidase (TcTS) is a key target protein for Chagas disease chemotherapy. In this study, we investigated the implications of active site flexibility on the biochemical mechanism of TcTS. Molecular dynamics studies revealed remarkable plasticity in the TcTS catalytic site, demonstrating, for the first time, how donor substrate engagement with the enzyme induces an acceptor binding site in the catalytic pocket that was not previously captured in crystal structures. Furthermore, NMR data showed cooperative binding between donor and acceptor substrates, supporting theoretical results. In summary, our data put forward a coherent dynamic framework to understand how a glycosidase evolved its highly efficient trans-glycosidase activity.

Molecular modeling studies of the structural, electronic, and UV absorption properties of benzophenone derivatives.

Benzophenone derivatives (BZP), an important class of organic UV filters, are widely used in sunscreen products due to their ability to absorb in the UVA and UVB ranges. The structural, electronic, and spectral properties of BZP derivatives have been studied by density functional theory (DFT) and time-dependent DFT (TD-DFT) methods. DFT/B3LYP with the 6-31G(d) basis set is an accurate method for optimizing the geometry of BZPs. The absorption maxima obtained from the TD-DFT calculations in a vacuum were in agreement with the experimental absorption bands and showed that the main electronic transitions in the UVA/UVB range present π → π* character, the major transition being HOMO → LUMO. The oscillator strength seems to increase in the presence of disubstitution at the para position. For protic substituents, the position appears to be related to the absorption band. Absorption in the UVB range occurs in the presence of para substitution, whereas ortho substitution leads to absorption in the UVA spectral region. The obtained results provide some features for BZP derivatives that can be useful for customizing absorption properties (wavelengths and intensities) and designing new BZP derivatives as sunscreens.

Structure-based design and biological profile of (E)-N-(4-Nitrobenzylidene)-2-naphthohydrazide, a novel small molecule inhibitor of IκB kinase-ß.

In this study, we describe the rational design, molecular modeling and pharmacological profile of a novel IKK-ß inhibitor (E)-N-(4-nitrobenzylidene)-2-naphthohydrazide (LASSBio-1524). The design based on the IKK-ß active site, and a privileged structure template yielded a novel IKK-ß inhibitor scaffold with significant selectivity over IKK-α and CHK2, as assessed by an in vitro kinase assay. For a better understanding of the structural requirements of IKK-ß inhibition, molecular dynamics simulations of LASSBio-1524 (3) were performed. The NAH derivative LASSBio-1524 (3), was able to suppress arachidonic acid-induced edema formation in a dose-dependent manner, demonstrating an in vivo anti-inflammatory effect. The molecular architecture of this novel, low-molecular weight IKK-ß inhibitor is encouraging for further lead optimization toward the development of innovative anti-inflammatory drug candidates.

A new class of mechanism-based inhibitors for Trypanosoma cruzi trans-sialidase and their influence on parasite virulence.

One of the most interesting aspects of Trypanosoma cruzi is its adaptation to obtain sialic acid from its host, fulfilling this need exclusively through the reaction catalyzed by enzymatically active trans-sialidase (aTS), thought to play an important role in the pathogenesis of Chagas' disease. Herein, we report that 2-difluoromethyl-4-nitrophenyl-3,5-dideoxy-d-glycero-alpha-d-galacto-2-nonulopyranosid acid (NeuNAcFNP) inactivates aTS time- and dose-dependently, and this inhibition was not relieved removing the inhibitor. Also, NeuNAcFNP causes a decrease in infection of mammalian cells. Characterization of labeled aTS by matrix-assisted laser desorption/ionization time-of-flight/time-of-flight mass spectrometry revealed that inactivation of the enzyme occurs through formation of a covalent bond between Arg245 and Asp247 and the inhibitor aglycone. Participation of Asp247 in the catalytic mechanism was proved by constructing a TSD247A mutant, which presents only residual activity. Molecular dynamic simulations indicate that the D247A mutation results in a more open catalytic cleft. In summary, NeuNAcFNP is the first reported mechanism-based inhibitor of aTS, representing a new template for drug design and opening new possibilities for chemotherapy of Chagas' disease, as well as for the elucidation of aTS function in T. cruzi pathogenesis and biology.