Rescue of Rare CFTR Trafficking Mutants Highlights a Structural Location-Dependent Pattern for Correction.

Cystic Fibrosis (CF) is a genetic disease caused by mutations in the gene encoding the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) channel. Currently, more than 2100 variants have been identified in the gene, with a large number being very rare. The approval of modulators that act on mutant CFTR protein, correcting its molecular defect and thus alleviating the burden of the disease, revolutionized the field of CF. However, these drugs do not apply to all patients with CF, especially those with rare mutations-for which there is a lack of knowledge on the molecular mechanisms of the disease and the response to modulators. In this work, we evaluated the impact of several rare putative class II mutations on the expression, processing, and response of CFTR to modulators. Novel cell models consisting of bronchial epithelial cell lines expressing CFTR with 14 rare variants were created. The variants studied are localized at Transmembrane Domain 1 (TMD1) or very close to the signature motif of Nucleotide Binding Domain 1 (NBD1). Our data show that all mutations analyzed significantly decrease CFTR processing and while TMD1 mutations respond to modulators, those localized in NBD1 do not. Molecular modeling calculations confirm that the mutations in NBD1 induce greater destabilization of CFTR structure than those in TMD1. Furthermore, the structural proximity of TMD1 mutants to the reported binding site of CFTR modulators such as VX-809 and VX-661, make them more efficient in stabilizing the CFTR mutants analyzed. Overall, our data suggest a pattern for mutation location and impact in response to modulators that correlates with the global effect of the mutations on CFTR structure.

Unraveling the Aquaporin-3 Inhibitory Effect of Rottlerin by Experimental and Computational Approaches.

The natural polyphenolic compound Rottlerin (RoT) showed anticancer properties in a variety of human cancers through the inhibition of several target molecules implicated in tumorigenesis, revealing its potential as an anticancer agent. Aquaporins (AQPs) are found overexpressed in different types of cancers and have recently emerged as promising pharmacological targets. Increasing evidence suggests that the water/glycerol channel aquaporin-3 (AQP3) plays a key role in cancer and metastasis. Here, we report the ability of RoT to inhibit human AQP3 activity with an IC50 in the micromolar range (22.8 ± 5.82 µM for water and 6.7 ± 2.97 µM for glycerol permeability inhibition). Moreover, we have used molecular docking and molecular dynamics simulations to understand the structural determinants of RoT that explain its ability to inhibit AQP3. Our results show that RoT blocks AQP3-glycerol permeation by establishing strong and stable interactions at the extracellular region of AQP3 pores interacting with residues essential for glycerol permeation. Altogether, our multidisciplinary approach unveiled RoT as an anticancer drug against tumors where AQP3 is highly expressed providing new information to aquaporin research that may boost future drug design.

Optimization of an in Silico Protocol Using Probe Permeabilities to Identify Membrane Pan-Assay Interference Compounds.

Membrane pan-assay interference compounds (PAINS) are a class of molecules that interact nonspecifically with lipid bilayers and alter their physicochemical properties. An early identification of these compounds avoids chasing false leads and the needless waste of time and resources in drug discovery campaigns. In this work, we optimized an in silico protocol on the basis of umbrella sampling (US)/molecular dynamics (MD) simulations to discriminate between compounds with different membrane PAINS behavior. We showed that the method is quite sensitive to membrane thickness fluctuations, which was mitigated by changing the US reference position to the phosphate atoms of the closest interacting monolayer. The computational efficiency was improved further by decreasing the number of umbrellas and adjusting their strength and position in our US scheme. The inhomogeneous solubility-diffusion model (ISDM) used to calculate the membrane permeability coefficients confirmed that resveratrol and curcumin have distinct membrane PAINS characteristics and indicated a misclassification of nothofagin in a previous work. Overall, we have presented here a promising in silico protocol that can be adopted as a future reference method to identify membrane PAINS.

Alchemical Design of Pharmacological Chaperones with Higher Affinity for Phenylalanine Hydroxylase.

Phenylketonuria (PKU) is a rare metabolic disease caused by variations in a human gene, PAH, encoding phenylalanine hydroxylase (PAH), and the enzyme converting the essential amino acid phenylalanine into tyrosine. Many PKU-causing variations compromise the conformational stability of the encoded enzyme, decreasing or abolishing its catalytic activity, and leading to an elevated concentration of phenylalanine in the blood, which is neurotoxic. Several therapeutic approaches have been developed to treat the more severe manifestations of the disorder, but they are either not entirely effective or difficult to adhere to throughout life. In a search for novel pharmacological chaperones to treat PKU, a lead compound was discovered (compound IV) that exhibited promising in vitro and in vivo chaperoning activity on PAH. The structure of the PAH-IV complex has been reported. Here, using alchemical free energy calculations (AFEC) on the structure of the PAH-IV complex, we design a new generation of compound IV-analogues with a higher affinity for the enzyme. Seventeen novel analogues were synthesized, and thermal shift and isothermal titration calorimetry (ITC) assays were performed to experimentally evaluate their stabilizing effect and their affinity for the enzyme. Most of the new derivatives bind to PAH tighter than lead compound IV and induce a greater thermostabilization of the enzyme upon binding. Importantly, the correspondence between the calculated alchemical binding free energies and the experimentally determined ΔΔGb values is excellent, which supports the use of AFEC to design pharmacological chaperones to treat PKU using the X-ray structure of their complexes with the target PAH enzyme.

Evaluation of EGCG Loading Capacity in DMPC Membranes.

Catechins are molecules with potential use in different pathologies such as diabetes and cancer, but their pharmaceutical applications are often hindered by their instability in the bloodstream. This issue can be circumvented using liposomes as their nanocarriers for in vivo delivery. In this work, we studied the molecular details of (-)-epigallocatechin-3-gallate (EGCG) interacting with 1,2-dimyristoyl- sn-glycero-3-phosphocholine (DMPC) monolayer/bilayer systems to understand the catechin loading ability and liposome stability, using experimental and computational techniques. The molecular dynamics simulations show the EGCG molecules deep inside the lipid bilayer, positioned below the lipid ester groups, generating a concentration-dependent lipid condensation. This effect was also inferred from the surface pressure isotherms of DMPC monolayers. In the polarization-modulated infrared reflection absorption spectra assays, the predominant effect at higher concentrations of EGCG (e.g., 20 mol %) was an increase in lipid tail disorder. The steady-state fluorescence data confirmed this disordered state, indicating that the catechin-induced liposome aggregation outweighs the condensation effects. Therefore, by adding more than 10 mol % EGCG to the liposomes, a destabilization of the vesicles occurs with the ensuing release of entrapped catechins. The loading capacity for DMPC seems to be limited by its disordered lipid arrangements, typical of a fluid phase. To further increase the clinical usefulness of liposomes, lipid bilayers with more stable and organized assemblies should be employed to avoid aggregation at large concentrations of catechin.

Insights on the Mechanism of Action of INH-C10 as an Antitubercular Prodrug.

Tuberculosis remains one of the top causes of death worldwide, and combating its spread has been severely complicated by the emergence of drug-resistance mutations, highlighting the need for more effective drugs. Despite the resistance to isoniazid (INH) arising from mutations in the katG gene encoding the catalase-peroxidase KatG, most notably the S315T mutation, this compound is still one of the most powerful first-line antitubercular drugs, suggesting further pursuit of the development of tailored INH derivatives. The N'-acylated INH derivative with a long alkyl chain (INH-C10) has been shown to be more effective than INH against the S315T variant of Mycobacterium tuberculosis, but the molecular details of this activity enhancement are still unknown. In this work, we show that INH N'-acylation significantly reduces the rate of production of both isonicotinoyl radical and isonicotinyl-NAD by wild type KatG, but not by the S315T variant of KatG mirroring the in vivo effectiveness of the compound. Restrained and unrestrained MD simulations of INH and its derivatives at the water/membrane interface were performed and showed a higher preference of INH-C10 for the lipidic phase combined with a significantly higher membrane permeability rate (27.9 cm s-1), compared with INH-C2 or INH (3.8 and 1.3 cm s-1, respectively). Thus, we propose that INH-C10 is able to exhibit better minimum inhibitory concentration (MIC) values against certain variants because of its better ability to permeate through the lipid membrane, enhancing its availability inside the cell. MIC values of INH and INH-C10 against two additional KatG mutations (S315N and D735A) revealed that some KatG variants are able to process INH faster than INH-C10 into an effective antitubercular form (wt and S315N), while others show similar reaction rates (S315T and D735A). Altogether, our results highlight the potential of increased INH lipophilicity for treating INH-resistant strains.

Hereditary nonspherocytic hemolytic anemia caused by red cell glucose-6-phosphate isomerase (GPI) deficiency in two Portuguese patients: Clinical features and molecular study.

Glucose-6-phosphate isomerase (GPI) deficiency cause hereditary nonspherocytic hemolytic anemia (HNSHA) of variable severity in individuals homozygous or compound heterozygous for mutations in GPI gene. This work presents clinical features and genotypic results of two patients of Portuguese origin with GPI deficiency. The patients suffer from a mild hemolytic anemia (Hb levels ranging from 10 to 12.7g/mL) associated with macrocytosis, reticulocytosis, hyperbilirubinemia, hyperferritinemia and slight splenomegaly. Genomic DNA sequencing revealed in one patient homozygosity for a new missense mutation in exon 3, c.260G>C (p.Gly87Ala), and in the second patient compound heterozygosity for the same missense mutation (p.Gly87Ala), along with a frameshift mutation resulting from a single nucleotide deletion in exon 14, c.1238delA (p.Gln413Arg fs*24). Mutation p.Gln413Arg fs*24 is the first frameshift null mutation to be described in GPI deficiency. Molecular modeling suggests that the structural change induced by the p.Gly87Ala pathogenic variant has direct impact in the structural arrangement of the region close to the active site of the enzyme.

Self-assembly molecular dynamics simulations shed light into the interaction of the influenza fusion Peptide with a membrane bilayer.

Influenza virus is one of the most devastating human pathogens. In order to infect host cells, this virus fuses its membrane with the host membrane in a process mediated by the glycoprotein hemagglutinin. During fusion, the N-terminal region of hemagglutinin, which is known as the fusion peptide (FP), inserts into the host membrane, promoting lipid mixing between the viral and host membranes. Therefore, this peptide plays a key role in the fusion process, but the exact mechanism by which it promotes lipid mixing is still unclear. To shed light into this matter, we performed molecular dynamics (MD) simulations of the influenza FP in different environments (water, dodecylphosphocholine (DPC) micelles, and a dimyristoylphosphatidylcholine (DMPC) membrane). While in pure water the peptide lost its initial secondary structure, in simulations performed in the presence of DPC micelles it remained stable, in agreement with previous experimental observations. In simulations performed in the presence of a preassembled DMPC bilayer, the peptide became unstructured and was unable to insert into the membrane as a result of technical limitations of the method used. To overcome this problem, we used a self-assembly strategy, assembling the membrane together with the peptide. These simulations revealed that the peptide can adopt a membrane-spanning conformation, which had not been predicted by previous MD simulation studies. The peptide insertion had a strong effect on the membrane, lowering the bilayer thickness, disordering nearby lipids, and promoting lipid tail protrusion. These results contribute to a better understanding of the role of the FP in the fusion process.

Decoding the secrets: how conformational and structural regulators inhibit the human 20S proteasome.

Acquired resistance to drugs that modulate specific protein functions, such as the human proteasome, presents a significant challenge in targeted therapies. This underscores the importance of devising new methodologies to predict drug binding and potential resistance due to specific protein mutations. In this work, we conducted an extensive computational analysis to ascertain the effects of selected mutations (Ala49Thr, Ala50Val, and Cys52Phe) within the active site of the human proteasome. Specifically, we sought to understand how these mutations might disrupt protein function either by altering protein stability or by impeding interactions with a clinical administered drug. Leveraging molecular dynamics simulations and molecular docking calculations, we assessed the effect of these mutations on protein stability and ligand affinity. Notably, our results indicate that the Cys52Phe mutation critically impacts protein-ligand binding, providing valuable insights into potential proteasome inhibitor resistance.

Probing the Interactions of Thiazole Abietane Inhibitors with the Human Serine Hydrolases ABHD16A and ABHD12.

12-Thiazole abietanes are highly selective reversible inhibitors of hABHD16A that could potentially alleviate neuroinflammation. In this study, we used synthetic chemistry, competitive activity-based protein profiling, and computational methodologies to try to establish relevant structural determinants of activity and selectivity of this class of compounds for inhibiting ABHD16A over ABHD12. Five compounds significantly inhibited hABHD16A but also very efficiently discriminated between inhibition of hABHD16A and hABHD12, with compound 35 being the most effective, at 100 µM (55.1 ± 8.7%; p < 0.0001). However, an outstanding switch in the selectivity toward ABHD12 was observed in the presence of a ring A ester, if the C2' position of the thiazole ring possessed a 1-hydroxyethyl group, as in compound 28. Although our data were inconclusive as to whether the observed enzyme inhibition is allosteric or not, we anticipate that the structure-activity relationships presented herein will inspire future drug discovery efforts in this field.

Structural determinants for the membrane insertion of the transmembrane peptide of hemagglutinin from influenza virus.

Membrane fusion is a process involved in a high range of biological functions, going from viral infections to neurotransmitter release. Fusogenic proteins increase the slow rate of fusion by coupling energetically downhill conformational changes of the protein to the kinetically unfavorable fusion of the membrane lipid bilayers. Hemagglutinin is an example of a fusogenic protein, which promotes the fusion of the membrane of the influenza virus with the membrane of the target cell. The N-terminus of the HA2 subunit of this protein contains a fusion domain described to act as a destabilizer of the target membrane bilayers, leading eventually to a full fusion of the two membranes. On the other hand, the C-terminus of the same subunit contains a helical transmembrane domain which was initially described to act as the anchor of the protein to the membrane of the virus. However, in recent years the study of this peptide segment has been gaining more attention since it has also been described to be involved in the membrane fusion process. Yet, the structural characterization of the interaction of such a protein domain with membrane lipids is still very limited. Therefore, in this work, we present a study of this transmembrane peptide domain in the presence of DMPC membrane bilayers, and we evaluate the effect of several mutations, and the effect of peptide oligomerization in this interaction process. Our results allowed us to identify and confirm amino acid residue motifs that seem to regulate the interaction between the segment peptide and membrane bilayers. Besides these sequence requirements, we have also identified length and tilt requirements that ultimately contribute to the hydrophobic matching between the peptide and the membrane. Additionally, we looked at the association of several transmembrane peptide segments and evaluated their direct interaction and stability inside a membrane bilayer. From our results we could conclude that three independent TM peptide segments arrange themselves in a parallel arrangement, very similarly to what is observed for the C-terminal regions of the hemagglutinin crystallographic structure of the protein, to where the segments are attached.

Major Achievements in the Design of Quadruplex-Interactive Small Molecules.

Organic small molecules that can recognize and bind to G-quadruplex and i-Motif nucleic acids have great potential as selective drugs or as tools in drug target discovery programs, or even in the development of nanodevices for medical diagnosis. Hundreds of quadruplex-interactive small molecules have been reported, and the challenges in their design vary with the intended application. Herein, we survey the major achievements on the therapeutic potential of such quadruplex ligands, their mode of binding, effects upon interaction with quadruplexes, and consider the opportunities and challenges for their exploitation in drug discovery.

In Silico Characterization of African Swine Fever Virus Nucleoprotein p10 Interaction with DNA.

African swine fever virus (ASFV) is the etiological agent of a highly contagious, hemorrhagic infectious swine disease, with a tremendous sanitary and economic impact on a global scale. Currently, there are no globally available vaccines or treatments. The p10 protein, a structural nucleoprotein encoded by ASFV, has been previously described as capable of binding double-stranded DNA (dsDNA), which may have implications for viral replication. However, the molecular mechanism that governs this interaction is still unknown, mostly due to the lack of a structural model for this protein. In this work, we have generated an ab initio model of the p10 protein and performed extensive structural characterization, using molecular dynamics simulations to identify the motifs and residues regulating DNA recognition. The helix-turn-helix motif identified at the C-terminal region of the protein was shown to be crucial to the dsDNA-binding efficiency. As with other DNA-binding proteins, two distinct serine and lysine-rich regions found in the two helices were identified as key players in the binding to DNA, whose importance was later validated using experimental binding assays. Altogether, these findings may contribute to a better understanding of the p10 function in ASFV replication.

Structural analysis of Thermus thermophilus HB27 mannosyl-3-phosphoglycerate synthase provides evidence for a second catalytic metal ion and new insight into the retaining mechanism of glycosyltransferases.

Mannosyl-3-phosphoglycerate synthase is a glycosyltransferase involved in the two-step synthetic pathway of mannosylglycerate, a compatible solute that accumulates in response to salt and/or heat stresses in many microorganisms thriving in hot environments. The three-dimensional structure of mannosyl-3-phosphoglycerate synthase from Thermus thermophilus HB27 in its binary complex form, with GDP-alpha-D-mannose and Mg(2+), shows a second metal binding site, about 6 A away from the mannose moiety. Kinetic and mutagenesis studies have shown that this metal site plays a role in catalysis. Additionally, Asp(167) in the DXD motif is found within van der Waals contact distance of the C1' atom in the mannopyranose ring, suggesting its action as a catalytic nucleophile, either in the formation of a glycosyl-enzyme intermediate according to the double-displacement S(N)2 reaction mechanism or in the stabilization of the oxocarbenium ion-like intermediate according to the D(N)*A(Nss) (S(N)i-like) reaction mechanism. We propose that either mechanism may occur in retaining glycosyltransferases with a GT-A fold, and, based on the gathered structural information, we identified an extended structural signature toward a common scaffold between the inverting and retaining glycosyltransferases.

Identification of Pan-Assay INterference compoundS (PAINS) Using an MD-Based Protocol.

Pan-assay interference compounds (PAINS) are promiscuous molecules with multiple behaviors that interfere with assay readouts. Membrane PAINS are a subset of these compounds that influence the function of membrane proteins by nonspecifically perturbing the lipid membranes that surround them. Here, we describe a computational protocol to identify potential membrane PAINS molecules by calculating the effect that a given compound has on the bilayer deformation propensity.

New (Iso)quinolinyl-pyridine-2,6-dicarboxamide G-Quadruplex Stabilizers. A Structure-Activity Relationship Study.

G-quadruplex (G4)-interactive small molecules have a wide range of potential applications, not only as drugs, but also as sensors of quadruplex structures. The purpose of this work is the synthesis of analogues of the bis-methylquinolinium-pyridine-2,6-dicarboxamide G4 ligand 360A, to identify relevant structure-activity relationships to apply to the design of other G4-interactive small molecules bearing bis-quinoline or bis-isoquinoline moieties. Thermal denaturation experiments revealed that non-methylated derivatives with a relative 1,4 position between the amide linker and the nitrogen of the quinoline ring are moderate G4 stabilizers, with a preference for the hybrid h-Telo G4, a 21-nt sequence present in human telomeres. Insertion of a positive charge upon methylation of quinoline/isoquinoline nitrogen increases compounds' ability to selectively stabilize G4s compared to duplex DNA, with a preference for parallel structures. Among these, compounds having a relative 1,3-position between the charged methylquinolinium/isoquinolinium nitrogen and the amide linker are the best G4 stabilizers. More interestingly, these ligands showed different capacities to selectively block DNA polymerization in a PCR-stop assay and to induce G4 conformation switches of hybrid h-Telo G4. Molecular dynamic simulations with the parallel G4 formed by a 21-nt sequence present in k-RAS gene promoter, showed that the relative spatial orientation of the two methylated quinoline/isoquinoline rings determines the ligands mode and strength of binding to G4s.

In Silico End-to-End Protein-Ligand Interaction Characterization Pipeline: The Case of SARS-CoV-2.

SARS-CoV-2 triggered a worldwide pandemic disease, COVID-19, for which an effective treatment has not yet been settled. Among the most promising targets to fight this disease is SARS-CoV-2 main protease (Mpro), which has been extensively studied in the last few months. There is an urgency for developing effective computational protocols that can help us tackle these key viral proteins. Hence, we have put together a robust and thorough pipeline of in silico protein-ligand characterization methods to address one of the biggest biological problems currently plaguing our world. These methodologies were used to characterize the interaction of SARS-CoV-2 Mpro with an α-ketoamide inhibitor and include details on how to upload, visualize, and manage the three-dimensional structure of the complex and acquire high-quality figures for scientific publications using PyMOL (Protocol 1); perform homology modeling with MODELLER (Protocol 2); perform protein-ligand docking calculations using HADDOCK (Protocol 3); run a virtual screening protocol of a small compound database of SARS-CoV-2 candidate inhibitors with AutoDock 4 and AutoDock Vina (Protocol 4); and, finally, sample the conformational space at the atomic level between SARS-CoV-2 Mpro and the α-ketoamide inhibitor with Molecular Dynamics simulations using GROMACS (Protocol 5). Guidelines for careful data analysis and interpretation are also provided for each Protocol.

Dioxygen and nitric oxide pathways and affinity to the catalytic site of rubredoxin:oxygen oxidoreductase from Desulfovibrio gigas.

Rubredoxin:oxygen oxidoreductase (ROO) is the terminal oxidase of a soluble electron transfer chain found in Desulfovibrio gigas. This protein belongs to the flavodiiron family and was initially described as an oxygen reductase, converting this substrate to water and acting as an oxygen-detoxifying system. However, more recent studies evidenced also the ability for this protein to act as a nitric oxide reductase, suggesting an alternative physiological role. To clarify the apparent bifunctional nature of this protein, we performed molecular dynamics simulations of the protein, in different redox states, together with O(2) and NO molecules in aqueous solution. The two small molecules were parameterized using free-energy calculations of the hydration process. With these simulations we were able to identify specific protein paths that allow the diffusion of both these molecules through the protein towards the catalytic centers. Also, we have tried to characterize the preference of ROO towards the presence of O(2) and/or NO at the active site. By using free-energy simulations, we did not find any significant preference for ROO to accommodate both O(2) and NO. Also, from our molecular dynamics simulations we were able to identify similar diffusion profiles for both O(2) and NO molecules. These two conclusions are in good agreement with previous experimental works stating that ROO is able to catalyze both O(2) and NO.

Protein thermal stabilization by charged compatible solutes: Computational studies in rubredoxin from Desulfovibrio gigas.

Molecular dynamics simulation studies of rubredoxin from Desulfovibrio gigas (RDG) were used to characterize the molecular mechanism of thermal stabilization by the compatible solute (CS) diglycerol-phospate (DGP). DGP is a negatively charged CS that accumulates under salt stress in Archaeoglobus fulgidus. Experimental results show that a 100 mM DGP solution exerts a strong protection effect in the half-life of RDG at 363 K (Lamosa et al., Appl Environ Microbiol 2000;66:1974-1979). RDG was simulated in four aqueous solutions at 300 and 363 K: pure aqueous media, 100 mM DGP, 100 mM NaCl, and 500 mM DGP. Our work shows that the 100 mM DGP solution is able to maintain the average protein structure when the temperature is increased, preventing the occurrence of large-scale deviation of a mobile loop involved in the first steps of RDG unfolding. The molecular mechanism of thermal denaturation protection by DGP seems to involve the direct interaction between the protein and the CS by hydrogen bond interactions near the mobile loop. Several clusters of DGP molecules are formed and preferentially localized at neutral electrostatic regions of the surface. The increase of DGP concentration to 500 mM did not yield better stabilization of the protein suggesting that the thermal protective role of this charged CS is achieved at low concentrations, as shown experimentally.

Proactive response to tackle the threat of emerging drugs: Synthesis and toxicity evaluation of new cathinones.

The emergence of potentially dangerous new psychoactive substances (NPS) imposes enormous challenges on forensic laboratories regarding their rapid and unambiguous identification. Access to comprehensive databases is essential for a quick characterization of these substances, allowing them to be categorized according to national and international legislations. In this work, it is reported the synthesis and structural characterization by NMR and MS of a library encompassing 21 cathinones, 4 of which are not yet reported in the literature, but with structural characteristics that make them a target for clandestine laboratories. This in-house library will be an important tool accessible to forensic laboratories, for the quick identification of seized NPS. The in vitro cytotoxicity of all cathinones was assessed in HepG2 cells, to have a preliminary but effective indication of their human hepatotoxicity potential. The two new cathinones DMB (8) and DMP (9) were the more cytotoxic, followed by the already seized mephedrone (2), 3,4-DMMC (3), 4-MDMC (7), NEB (12) with EC50 values ranging from 0.81mM for (3) to 1.28mM for (2). Results suggest an increase of cytotoxicity with the increase of the chain length of the acyl moiety and with the substitution (with one or two methyl groups) in the aromatic ring. The nature of the amine moiety seems to play only a minor role in the cytotoxic effect. Molecular dynamics simulations were performed to evaluate the molecular details related with the observed cytotoxicities. Although these studies indicated that cathinones are able to cross lipid bilayers with relative ease, when in their neutral forms, it was observed only a partial correlation between lipophilicity and cytotoxicity, indicating that membrane trafficking may not be the only key factor influencing the bioactivity of these compounds. This work is a valuable contribution to the forensic science field since a quick identification of novel cathinones is urgent to match their rapid increase in the market.