Design and Optimization of Quinazoline Derivatives: New Non-nucleoside Inhibitors of Bovine Viral Diarrhea Virus

Bovine viral diarrhea virus (BVDV) belongs to the Pestivirus genus (Flaviviridae). In spite of the availability of vaccines, the virus is still causing substantial financial losses to the livestock industry. In this context, the use of antiviral agents could be an alternative strategy to control and reduce viral infections. The viral RNA-dependent RNA polymerase (RdRp) is essential for the replication of the viral genome and constitutes an attractive target for the identification of antiviral compounds. In a previous work, we have identified potential molecules that dock into an allosteric binding pocket of BVDV RdRp via a structure-based virtual screening approach. One of them, N-(2-morpholinoethyl)-2-phenylquinazolin-4-amine [1, 50% effective concentration (EC50) = 9.7 ± 0.5 μM], was selected to perform different chemical modifications. Among 24 derivatives synthesized, eight of them showed considerable antiviral activity. Molecular modeling of the most active compounds showed that they bind to a pocket located in the fingers and thumb domains in BVDV RdRp, which is different from that identified for other non-nucleoside inhibitors (NNIs) such as thiosemicarbazone (TSC). We selected compound 2-[4-(2-phenylquinazolin-4-yl)piperazin-1-yl]ethanol (1.9; EC50 = 1.7 ± 0.4 μM) for further analysis. Compound 1.9 was found to inhibit the in vitro replication of TSC-resistant BVDV variants, which carry the N264D mutation in the RdRp. In addition, 1.9 presented adequate solubility in different media and a high-stability profile in murine and bovine plasma.

The control of BVDV infections combines vaccination (which is not mandatory in some countries) with the removal of persistently infected animals, which are the main source of infective virus in the field (Hessman et al., 2009). The period taken to develop a protective immune response after vaccination provides a window for viral transmission. The use of antiviral agents could be an attractive strategy as a stopgap measure to prevent infection in naïve calves (Newcomer et al., 2014). Moreover, antivirals could also be used to treat valuable animals infected with pestiviruses in zoological collections and in in-breeding programs and in vitro embryo production, to cure established cell lines contaminated with pestiviruses, which are used, for example, to produce interferons and vaccines for medical use (Stringfellow et al., 2005;Newcomer and Givens, 2013;Newcomer et al., 2014;Musiu et al., 2020). However, no antivirals are currently available for controlling BVDV infections in either laboratories or farms. Over the last decades, several selective anti-BVDV compounds have been reported, including virus targeting and host targeting derivatives, e.g., polymerase inhibitors (Castro et al., 2019), protease inhibitors (Bukhtiyarova et al., 2001), human cellular enzyme inhibitors (Branza-Nichita et al., 2001), and entry inhibitors (Pascual et al., 2018). Furthermore, the design of effective antiviral therapies will likely require a combination of drugs with complementary mechanisms of action to decrease the probability of selecting resistant viral mutants (Newcomer and Givens, 2013). In the host cell, the viral replication is carried out by a large membrane-bound replication complex composed of viral RNA, viral polymerase, and other viral and cellular proteins. In this replication complex, a negative-strand RNA is synthesized from a positive-strand RNA genome, resulting in dsRNA intermediates from which numerous positive strands are produced for progeny viruses (Choi et al., 2006). The viral RNA-dependent RNA polymerase (NS5B RdRp) is responsible of viral RNA synthesis and has been used as a target of several nucleoside inhibitors (NIs) and non-nucleoside inhibitors (NNIs). NIs are widely used to treat hepatitis B virus (HBV), hepatitis C virus (HCV), HIV-1, and herpesvirus infections. They act by competing with nucleotide substrates for binding to the active site of the polymerase and induce termination of the nucleic acid chain synthesis (Finkielsztein et al., 2010). In turn, NNIs are highly specific and act allosterically by inhibiting the enzymatic activity of the polymerase. This group includes a series of very diverse chemical structures, such as N-propyl-N-[2-(2H-1,2,4triazino[5,6-b]indol-3-ylthio)ethyl]-1-propanamine (VP32947) (Baginski et al., 2000), bromophenyl-imidazo-pyridines (BPIP) (Paeshuyse et al., 2007), ethyl 2-methylimidazo[1,2a]pyrrolo[2,3-c]pyridin-8-carboxylate (AG110) (Paeshuyse et al., 2007), pyrazolo-triazolo-pyrimidinamine (LZ37) (Paeshuyse et al., 2009), 2-(2-benzimidazolyl)-5-[4-(2imidazolino)phenyl]furan (DB772) (Newcomer et al., 2012), 5,6-dimethoxy-1-indanone [thiosemicarbazone (TSC)] (Castro et al., 2011;Soraires Santacruz et al., 2017), 2phenylbenzimidazole , and arylazoenamine derivatives (Giliberti et al., 2010).
In a previous work, we identified a series of molecules that dock into the allosteric binding pocket of BVDV RdRp, employing a structure-based virtual screening approach. Five of them were found to be active against BVDV in vitro and displayed 50% effective concentration (EC 50 ) values in the sub-and low-micromolar range (Castro et al., 2019). From these series of structurally and functionally diverse compounds, Compound 1 (EC 50 = 9.7 ± 0.5 µM) was chosen for further optimization as it resulted an easily synthesizable and scalable compound (Figure 1).
The present work reports the synthesis and anti-BVDV activity of a new family of heterocyclic compounds derived from quinazolines. With the aim of characterizing their mode of action, we combined molecular dynamics (MD) simulation and the assessment of the antiviral activity against well-characterized BVDV mutant viruses that are resistant to NNIs of the viral RdRp. In addition, the compound cytotoxicity and in vitro pharmacokinetic properties were determined. Together, our results suggest that these molecules are suitable antivirals for BVDV and can be used as a good starting point to optimize antiviral potency.

General Information
All reagents and solvents were obtained from commercial sources. Column chromatography was carried out by employing a silica gel (Kieselgel 60,Merck). Pre-coated F-254 silica gel plates were used for thin-layer analytical chromatography. 1 H and 13 C NMR spectra were recorded on Bruker BioSpin 600 MHz AVIII-600, Bruker Avance II 500 MHz, and Bruker 300 MHz spectrometers at room temperature. Chemical shifts (δ) are reported in ppm and coupling constants (J) in hertz. The mass spectrometer utilized was a Xevo G2-S QTOF (Waters Corporation, Manchester, UK) with an electrospray ionization (ESI) source. The mass spectrometer was operated in positive and negative ion modes with probe capillary voltages of 2.5 and 2.3 kV, respectively. The sampling cone voltage was 30 V. The source and desolvation gas temperatures were set to 120 and 350 • C, respectively. The nitrogen gas desolvation flow rate was 600 L h −1 , and the cone desolvation flow rate was 10 L h −1 . The mass spectrometer was calibrated across the range of m/z 50-1,200 using a 0.5 mM sodium formate solution prepared in 2-propanol/water (90:10 v/v). Data were drift-corrected during acquisition using a leucine enkephalin reference spray (LockSpray) infused at 2 µl min −1 . Data were acquired in the range of m/z 50-1,200, and the scan time was set from 22 to 1 s. Data acquisition and processing were carried out using MassLynx, version 4.1 (Waters Corporation, Milford, MA, USA). UV spectra were measured with a Shimadzu 3,600 UV/Vis/NIR spectrophotometer.
The synthesis of intermediate 4a-4g and 5a-5g are described in the Supplementary Material.

Molecular Modeling
For the molecular docking study, the RdRp protein structure was obtained from the RCSB Protein Data Bank (PDB: 1S48) (Choi et al., 2004). UCSF Chimera was used to mutate the MSEmodified amino acid to MET. Ligands and protein were then saved in a PDBQT file format for input into AutoDock Vina to carry out docking calculations (Trott and Olson, 2010). ADT was used to remove crystal waters, add polar hydrogens, and assign Gasteiger charges to the protein structure (Maier et al., 2015). For RdRp, the search space was a box with 30 × 30 × 30 Å XYZ dimensions, respectively, and was centered on −17, 177, and 19 XYZ coordinates, respectively.
MD simulations were performed using the Amber16 and AmberTools16 package (Case et al., 2016) with the Amber14SB force field (Maier et al., 2015). The docked poses of the compounds were employed as the initial configuration for the system. The TIP3P water model (Jorgensen et al., 1983) was used in a truncated octahedral box, extending 10 Å from the protein. A physiological salt concentration of 0.15 M was set by employing Na + and Cl − ions.
The simulations were performed using the PMEMD CUDA module of the Amber simulation suite and consisted of the following steps: an initial 1-ps run with a 0.01-fs timestep to eliminate bad contacts followed by an energy minimization; heating from 0 to 10 K over 10 ps with a 0.1-fs timestep with strong restraints (50 kcal mol −1 Å −2 ) on the protein residues and ligands; and then from 10 to 300 K over 90 ps with a 0.5-fs timestep and weaker restraints (10 kcal mol −1 Å −2 ). The system was then equilibrated for 400 ps at constant temperature and pressure with weak restraints on the ligands and on the CA atoms of the protein (1 kcal mol −1 Å −2 ). The Langevin thermostat was used with a collision frequency of 2.0 ps −1 (Davidchack et al., 2009); the SHAKE algorithm was used to constrain bonds, allowing a 2-fs timestep; and an 8.0-Å cutoff was used for nonbonded interactions. Finally, a 100-ns production simulation was performed under the NPT conditions described above.

Test Compounds
All compounds were solubilized in DMSO at a concentration of 5 and 1 mg/ml and subsequently serially diluted in a cell culture medium. The highest concentration evaluated of each compound contained DMSO at a final concentration <1% (non-cytotoxic concentration). The thiosemicarbazone of 5,6dimethoxy-1-indanone (TSC) was employed as a reference inhibitor of BVDV RdRp (Finkielsztein et al., 2008;Castro et al., 2011) (EC 50 = 1.8 ± 0.6 µM). Finally, the stability of compounds in DMSO and water:DMSO (98:2) was evaluated as described in the Supplementary Material.

Cytotoxicity Assay
MDBK cells were seeded in 96-well plates at a density of 1 × 10 4 cells per well in GM. After 24 h, the culture medium was removed, monolayers were washed twice with phosphatebuffered saline (PBS), and serial dilutions of the test compounds were made in E-MEM supplemented with donor horse serum (DHS, Gibco, IM) and added. Cultures were incubated for 3 days at 37 • C, and then, cell viability was determined by the MTS/PMS method according to the manufacturer's instructions (Promega). After incubation for 3 h at 37 • C and 5% CO 2 , the absorbance was determined at 490 nm. The 50% cytotoxic concentration (CC 50 ) was defined as the concentration of compound that inhibited the proliferation of exponentially growing cells by 50% and was calculated by interpolation in dose-response curves.

Anti-BVDV Assay
MDBK cells were seeded in 96-well microplates at a density of 1 × 10 4 cells per well in GM. After incubation for 24 h at 37 • C, the culture medium was removed, and the monolayers were washed twice with PBS and infected with BVDV with an inoculum that resulted in a >80% cytopathic effect (CPE) after 3 days of incubation at 37 • C (MOI = ∼0.01 PFU/cell). Then, serial dilutions of the test compounds in IM were added. Mockinfected cells and mock-treated infected cells were included as controls in each plate. After incubation for 3 days at 37 • C and 5% CO 2 , the culture medium was removed, cell monolayers were washed twice with PBS, and IM supplemented with MTS/PMS was added to each well. The absorbance of each well was read at 490 nm. The EC 50 was defined as the concentration of compound that prevented 50% of the virus-induced CPE and was calculated by interpolation in dose-response curves.

Cross-Resistance Assay
In order to determine the mechanism of antiviral action of this class of compounds, anti-BVDV assays were performed as stated above with different BVDV NADLs that are resistant to TSC: BVDV R1 (NS5B A392E) and BVDV R2 (NS5B N264D). Compound 1.9 was used at its maximum non-cytotoxic concentration (12.5 µM). As positive control, TSC was used at 34.0 µM (20 × EC 50 ). Mock-infected cells and mock-treated infected cells were included as controls in each assay plate. After 3 days, the culture medium was removed, cell monolayers were washed twice with PBS, and IM supplemented with MTS/PMS solution was added to each well. The percentage of CPE inhibition was calculated for each experimental condition.

Shake-Flask Method for Drug Solubility Testing
The shake-flask method was used as a reference assay (Alsenz and Kansy, 2007;Baka et al., 2008;Bollini et al., 2013). Compound 1.9 (1.0 mg) was measured in simulated gastric fluid (SGF, pH 1.2), simulated intestinal fluid (SIF, pH 6.8), and PBS solution (pH 7.4) by UV-Vis-NIR spectrophotometry. Verapamil was used as internal control. More detailed information is available in the Supplementary Material.

SGF, SIF, PBS, Mouse, and Bovine Plasma Stability
Stability assays were performed according to the procedure described in a previous work (Leal et al., 2019). Samples were taken at different incubation times and analyzed by high-performance liquid chromatography-mass spectrometry (HPLC-MS) using a Waters Alliance e2695 system (Waters Corporation, Milford, MA, USA) fitted with a Phenomenex Kinetex R XB-C18 column (150 × 4.6 mm, 5-µm particle size, Phenomenex Inc., Torrance, CA, USA) coupled to a Waters SQD2 single quadrupole mass spectrometer with an ESI source. A gradient elution was applied using 40% water and 0.1% acetic acid (mobile phase A) and 40% methanol (mobile phase B), with the following program: 0-3 min 40% B; 3-7 min 40-90% B; and 7-20 min 90% B. The flow rate was kept constant at 0.35 ml min −1 . After each sample injection, the gradient was returned to its initial conditions in 16 min. The injection volume was 5 µl. The column temperature was 35 • C. The mass spectrometer was operated in positive ion mode with a probe capillary voltage of 2.5 kV. The sampling cone voltage was set to 35.0 V. The source and desolvation gas temperatures were set at 150 and 350 • C, respectively. The nitrogen gas desolvation flow rate was 600 L h −1 , and the cone gas flow rate was 10 L h −1 . The mass spectrometer was calibrated across the range of m/z 20-2,023 with a sodium and cesium iodide solution. Data were acquired in scan mode with a scan duration of 0.2 s and in SIR mode with unit resolution. Data acquisition and processing were carried out using MassLynx software, version 4.1. Enalapril and salicylaldehyde isonicotinoyl hydrazine (SIH) (Kovaríková et al., 2008) were used as positive controls for murine and bovine plasma, respectively. SIH was synthesized according to the procedures described by Gaonkar et al. (2006) and Trzesowska-Kruszynska (2013).

Synthesis of Quinazoline Derivatives
Compound 1 emerged as one of the most promising scaffolds from our structure-based virtual screening approach toward the development of RdRp inhibitors (Castro et al., 2019). To improve the antiviral activity, we prepared a range of quinazoline analogs introducing modification at positions 4, 2, and 7 of the quinazoline scaffold (1.1-1.24).

Antiviral Activity Against wtBVDV
The synthesized compounds (1.1-1.24) were evaluated for their antiviral potency against BVDV in cultured MDBK cells employing TSC as antiviral activity positive control. The values of EC 50 , CC 50 , and selectivity index (SI, CC 50 /EC 50 ratio) of the target compounds are summarized in Figure 1 and Table 1.
Results demonstrated that the antiviral activity of compounds 1.7-1.10, 1.12, 1.18, 1.21, and 1.23 was higher than that of the lead compound 1, with EC 50 values ranging from 1.3 to 4.6 µM (Figure 1). Among them, 1.8 (EC 50 = 1.4 ± 0.1 µM) and 1.9 (EC 50 = 1.8 ± 0.4 µM) proved to be the most potent inhibitors with good SIs. A high in vitro SI suggests a specific inhibition of a viral target, which would result in an in vivo higher therapeutic index.
Preliminary SAR studies have revealed that different substituents located at position 4 on the quinazoline core have significant effects on the antiviral potency (Figure 2). For instance, aromatic groups like benzyl, methyl pyridine, and methyl pyridine N-oxide sharply reduced the antiviral potency (1.1-1.5 and 1.13). However, the presence of charges at physiological pH on the R 4 substituent conferred the compound higher bioactivity (1.7-1.10 and 1.12). Interestingly, the replacement of the N-hexyl substituent (1.6) for the polar N,N-dimethyl propane diamine (1.7) contributed greatly to the activity. Mock-treated cells were added as control. Cells were incubated and allowed to proliferate for 3 days at 37 • C and 5% CO 2 . Cell viability was then determined by the MTS/PMS method. (B) Antiviral activity. Sub-confluent MDBK cell monolayers grown in 96-well plates were infected with BVDV NADL (MOI = 0.01) and treated with different concentrations of the compounds. Mock-infected cells and mock-treated infected cells were included as controls in each plate. After 3 days at 37 • C and 5% CO 2 , cell viability was determined by the MTS/PMS method, and the percentage of the viral CPE inhibition was calculated. SCHEME 1 | General procedure for the synthesis of quinazoline derivatives.
The impact of the substitution in the aromatic ring at R 2 of compounds 1.8 and 1.9 on the antiviral activity was also assessed. The introduction of electron-donating substituents in the phenyl ring, such as N,N-dimethyl (1.18, EC 50 = 3.5 ± 0.6 µM), methoxy (1.19, EC 50 = 12.8 ± 0.2 µM), or methyl group (1.20, EC 50 = 8.5 ± 3.4 µM) generated analogs with either similar or reduced antiviral activity, as compared to the lead compounds 1.8 and 1.9 (EC 50 = 1.4 ± 0.1 and 1.8 ± 0.4 µM, respectively). However, the incorporation of a nitro group (1.14-1.16) as a strong electron-withdrawing substituent or two halogen atoms (1.17) in the phenyl group gave rise to either less active or inactive compounds.
Furthermore, four derivatives (1.21-1.24) lacking the phenyl substitution at position 2 of the quinazoline core were synthesized. Except compound 1.10, all the derivatives showed lower potency than their corresponding analogs (1 and 1.8-1.10). In addition, this series of compounds showed considerably low cytotoxicity.

Molecular Modeling
To characterize the interactions between the compounds and the allosteric site of the RdRp, we performed docking and MD simulation analyses. The docked poses were used as initial system configurations, and all protein-ligand complexes remained stable along each 100-ns MD simulations (Supplementary Figure 1).
In particular, the phenyl quinazoline 1.9 presented two favorable hydrogen bond interactions between the HN atoms   of the ligand and the O atom in R295 and E675. In both cases, interatomic distances and angles were very favorable, with average values of 2 Å and 150 • , respectively (Figure 3 and Supplementary Figure 2). This ligand also presented a paralleldisplaced π-π interaction with Y289. The 7-chloroquinazoline series (1.21-1.24) presented different binding modes with respect to phenyl quinazoline analogs, as these compounds move slightly inside the binding pocket during the interaction. The most active compound of this series (1.21) established a hydrogen bond of moderate strength between the oxygen of the morpholine ring and the NH of K319.
In a previous work (Castro et al., 2019), we reported a quinoline compound [2-(4-ethylphenyl)-N-isopentylquinolin-4amine] with the potential to inhibit essential processes during viral RNA synthesis. A computational model showed the importance of the interaction of this compound with Y674, R295, and N217 located in the fingers and thumb domains of BVDV RdRp. Interestingly, these new compounds (1.9 and 1.21) accommodate their quinazoline moiety into the deep hydrophobic pocket, similarly to quinoline derivatives. The BVDV RdRp structure has the shape of a right hand, comprising the fingers, palm, and thumb domains. In the RdRp structure, R295, Y289, and V306 are located in the fingers domain and are part of the motif II of the polymerase, which is conserved between different RNA polymerases of positive-strand RNA viruses. This motif was shown to establish contact with the phosphate backbone and bases; and it was reported that the substitution at the residue R295 completely abolished RNA synthesis. On the other hand, Y674 and E675 are located within the C-terminal loop in the thumb domain. The β-thumb region interacts with the fingers and palm domains through this long C-terminal loop (residues 670-679). Together with a loop between α-20 and α-21 in the thumb domain, the β-thumb reduces the volume of the template channel, and the C-terminal loop, in combination with the fingertip region, facilitates the translocation of the template and product RNA. Interestingly, similar to compounds 1.9 and 1.21, docking analysis indicates that active compounds interact with R295 in the fingertip domain of the polymerase, providing a possible explanation to their mode of action (Supplementary Figure 3).
In the present work, a new series of potential BVDV NNIs were identified. According to our observations, these new quinazoline derivatives would interact with different hydrophobic pockets. Our results are in line with previously reported studies demonstrating that the quinoxaline derivatives FIGURE 4 | Antiviral activity of 1.9 against TSC BVDV-resistant mutants. The anti-BVDV assay was performed as stated in Materials and Methods. For infection, wtBVDV NADL and two different TSC-resistant BVDV NADLs, BVDV R1 (NS5B A392E) and BVDV R2 (NS5B N264D), were used. Compound 1.9 was used at its maximum non-cytotoxic concentration (12.5 µM). As positive control, TSC was used at 34 µM. Mock-infected cells and non-treated infected cells (NTC) were included as controls in each plate. After 3 days, cell viability was determined by the MTS/PMS method, and the percentage of the viral CPE relative to the corresponding NTC was calculated. *p < 0.05 relative to NTC.
that are active against BVDV would act as NNIs of viral RdRp and showing their main interactions with residues R295 and Y674 (Carta et al., 2016). In that work, the authors identified two hydrophobic cavities HC1 (A221, I261, I287, and Y289) and HC2 (V216, Y303, V306, K307, P408, and A412) as pharmacophoric requirements. Interestingly, our compounds showed similar requirements, and while compound 1.9 interacts with residues from HC1 and HC2, the 7-chloroquinazoline derivative 1.21 moved slightly inside the binding pocket during the interaction with the RdRp structure and showed interaction with residues from HC2.   Antiviral Activity Against BVDV Resistant to NNI TSC As stated above, compound 1.9 would interact with a hydrophobic pocket other than the one identified for other BVDV RdRp NNIs such as TSC (Castro et al., 2011), suggesting that 1.9 could be active even against TSC-resistant BVDV. To test this hypothesis, we evaluated the antiviral activity of compound 1.9 against TSC-resistant BVDV (BVDV R1 and BVDV R2) ( Figure 4). Interestingly, when TSC R -BVDV was treated with compound 1.9, different results were observed. BVDV R1, which carries the NS5B A392E mutation, was inhibited 50% by compound 1.9, whereas BVDV R2, which carries the N264D mutation, was completely sensitive to compound 1.9 at 12.5 µM (100% inhibition). The partial resistance observed for BVDV R1 supports RdRp as a target for compound 1.9 since the A392E mutation partially blocks its antiviral activity. On the other hand, the sensitivity observed for BVDV R2 to compound 1.9 is in agreement with our previous hypothesis and also reinforces the binding of this compound to a different hydrophobic pocket as shown by our molecular modeling studies.
The analysis of MD simulations demonstrated that A392 is located in the vicinity of the compound 1.9 binding site, thus explaining the partially resistant phenotype of BVDV R1 to the compound (Figure 5). In addition, this mutation was reported to disturb the local conformation of the protein. We speculate that this change would prevent the access of compounds to the predicted binding site Castro et al., 2011).

Physicochemical and in vitro Pharmacokinetic
Properties of 1.9 Compound 1.9 was selected to perform in vitro physicochemical studies (solubility and stability). The solubility (S) was evaluated at SGF (pH 1.2), SIF (pH 6.8), and PBS solution (pH 7.4), employing the shake-flask method with UV spectroscopy (Bollini et al., 2013). Compound 1.9 presented low solubility values ( Table 2), which were within the range normally observed for oral drugs (http://pharmeuropa.edqm.eu/home/). The stability was evaluated by calculating the percentage of remaining drug (HPLC-MS) after contact with different incubation media (SIF, SGF, PBS, mouse, and bovine plasma) for 120 min. Figure 6 shows the high-stability profile of 1.9 under all the conditions employed. Neither compound modification nor degradation was detected after the different incubation periods. The stability profiles of enalapril and SIH, used as positive controls for mouse and bovine plasma, respectively, are also shown in Figure 5B. Additionally, no significant differences were found in the concentration of compound 1.9 between mouse and bovine plasma samples.

CONCLUSIONS
From a series of 24 derivatives of N-(2-morpholinoethyl)-2-phenylquinazolin-4-amine (1), eight compounds showed improved anti-BVDV activity. A molecular modeling analysis of the most active compound 1.9 showed that this compound would bind to a hydrophobic pocket of the BVDV RdRp that is different from the one previously described for most BVDV RdRp NNIs. Our results indicate that compound 1.9 inhibits the replication of TSC-resistant BVDV variants in vitro, which carry a mutation in the viral polymerase that is selected by almost all the NNIs described for BVDV (N264D). Taking into account the main residues with which compound 1.9 would interact, its binding to BVDV RdRp could affect the template channel, through the interaction with both the fingers domain (motif II) and the β-thumb region, and prevent the synthesis of viral RNA. In addition, compound 1.9 presented adequate solubility in different media and a high-stability profile in murine and bovine plasma samples. These parameters are important for the development of oral drugs. Taken together, these results suggest that compound 1.9 is an NNI of BVDV that can be considered a promising candidate to treat BVDV infection and encourages further studies assessing the mechanism of action and potential combination with other previously reported antivirals.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding authors.

AUTHOR CONTRIBUTIONS
MB and LC designed and supervised the study. NA, LB, and MB performed the docking and computational simulations. GF, AB, and MB conducted chemical synthesis. EC, RR, ME, and MF performed biological experiments. DF performed the experimental stability study of compound 1.9, and GF performed solubility assays. All authors were involved in the preparation of the manuscript and approved the final version.