Discovery of N-Containing (-)-Borneol Esters as Respiratory Syncytial Virus Fusion Inhibitors

Respiratory syncytial virus (RSV) causes acute respiratory infections, thus, posing a serious threat to the health of infants, children, and elderly people. In this study, we have discovered a series of potent RSV entry inhibitors with the (-)-borneol scaffold. The active compounds 3b, 5a, 5c, 7b, 9c, 10b, 10c, and 14b were found to exhibit activity against RSV A strain A2 in HEp-2 cells. The most active substances, 3b (IC50 = 8.9 μM, SI = 111) and 5a (IC50 = 5.0 μM, SI = 83), displayed more potency than the known antiviral agent Ribavirin (IC50 = 80.0 μM, SI = 50). Time-of-addition assay and temperature shift studies demonstrated that compounds 3b, 5a, and 6b inhibited RSV entry, probably by interacting with the viral F protein that mediated membrane fusion, while they neither bound to G protein nor inhibited RSV attachment to the target cells. Appling procedures of molecular modeling and molecular dynamics, the binding mode of compounds 3b and 5a was proposed. Taken together, the results of this study suggest (-)-borneol esters to be promising lead compounds for developing new anti-RSV agents.


Introduction
Increasing outbreaks of acute respiratory infections, mostly caused by viruses with epidemic or pandemic potential, pose a significant threat to public health. The most common viruses causing a large proportion of acute respiratory infections are influenza viruses, human rhinovirus, respiratory syncytial virus (RSV), adenovirus, and human coronavirus [1]. The COVID-19 pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has resulted in more than 518 million confirmed cases and more than six million deaths globally [2]. The influenza virus has raised alarming concerns worldwide. Millions of people have died from influenza pandemics (Spanish flu in 1918, Asian flu in 1957, and swine flu in 2009) [3]. Influenza viruses, coronaviruses, and other seasonal respiratory viruses tend to cause sporadic infections or seasonal epidemics in humans. However, a pandemic may occur when a new virus with sufficient transmissibility and pathogenicity emerges. Therefore, it is of great significance to find effective drugs against infections caused by respiratory viruses.
RSV is a negative-sense single-stranded RNA virus belonging to the family Paramyxoviridae. This virus is the most common cause of acute respiratory infections in infants and children [4]. Typically, RSV is mild and accompanied by upper respiratory symptoms, with the risk of severe RSV increased by factors such as young age (<6 months), congenital heart disease, chronic pulmonary diseases, and immunocompromised conditions [5]. Globally, infection with RSV is responsible for approximately 3.2 million hospital admissions and Table 1. Structure and antiviral activity of (-)-borneol-based esters.
fants and children [4]. Typically, RSV is mild and accompanied by upper respiratory symptoms, with the risk of severe RSV increased by factors such as young age (<6 months), congenital heart disease, chronic pulmonary diseases, and immunocompromised conditions [5]. Globally, infection with RSV is responsible for approximately 3.2 million hospital admissions and approximately 60,000 in-hospital deaths per year [6,7]. No specific etiotropic therapy is available, and the only compound (Synagis) designed to prevent this infection is practically inaccessible to the general population [8]. Several groups of researchers around the world are actively developing antiviral agents against respiratory syncytial infection, but none of them are yet suitable for wide use in clinical practice [9][10][11].
Recently, we found that (-)-borneol-based esters I-IV (Table 1) containing a saturated nitrogen ring or an acyclic tertiary nitrogen atom exhibited good anti-filovirus efficacy by affecting the glycoprotein (GP)-mediated membrane fusion process [12,13]. At the same time, in vitro results with an infectious virus and in a pseudovirus system, as well as theoretical calculations, allowed us to demonstrate that (-)-borneol esters IV with butyl substituents are likely to be fusion inhibitors targeting the S protein of SARS-CoV-2 [14]. Moreover, we found esters I-IV to effectively inhibit the replication of the influenza A virus [15] and esters III with a morpholine cycle to exhibit antiviral activity against orthopoxviruses [16]. These results indicate (-)-borneol-based esters to have broad-spectrum antiviral activity, acting on an early step of the virus life cycle, and that they should be studied as RSV inhibitors.  The selectivity index (SI, CC50/IC50 ratio) values of each compound were taken as the evaluation reference for antiviral potency. The antiviral activity against filoviruses and SARS-CoV-2 was studied using a "surrogate" system. This system included vesicular stomatitis virus (VSV) particles pseudotyped with Ebola virus GP (rVSV-ΔG-EboV-GP) or Marburg virus GP (rVSV-ΔG-MarV-GP), and lentiviral particles exhibiting the S protein of SARS-CoV-2 (pSARS-CoV-2-S). a SIMarV-GP = CC50/IC50MarV-GP, where IC50MarV-GP is the concentration of a compound required to inhibit rVSV-ΔG-MarV-GP infection of HEK293T cells by 50%. b SIEboV-GP = CC50/IC50EboV-GP, where IC50EboV-GP is the concentration of a compound required to inhibit rVSV-ΔG-EboV-GP infection of HEK293T cells by 50%. c SIH1N1 = CC50/IC50H1N1, where IC50H1N1: con-Virus n fants and children [4]. Typically, RSV is mild and accompanied by upper respiratory symptoms, with the risk of severe RSV increased by factors such as young age (<6 months), congenital heart disease, chronic pulmonary diseases, and immunocompromised conditions [5]. Globally, infection with RSV is responsible for approximately 3.2 million hospital admissions and approximately 60,000 in-hospital deaths per year [6,7]. No specific etiotropic therapy is available, and the only compound (Synagis) designed to prevent this infection is practically inaccessible to the general population [8]. Several groups of researchers around the world are actively developing antiviral agents against respiratory syncytial infection, but none of them are yet suitable for wide use in clinical practice [9][10][11].
Recently, we found that (-)-borneol-based esters I-IV (Table 1) containing a saturated nitrogen ring or an acyclic tertiary nitrogen atom exhibited good anti-filovirus efficacy by affecting the glycoprotein (GP)-mediated membrane fusion process [12,13]. At the same time, in vitro results with an infectious virus and in a pseudovirus system, as well as theoretical calculations, allowed us to demonstrate that (-)-borneol esters IV with butyl substituents are likely to be fusion inhibitors targeting the S protein of SARS-CoV-2 [14]. Moreover, we found esters I-IV to effectively inhibit the replication of the influenza A virus [15] and esters III with a morpholine cycle to exhibit antiviral activity against orthopoxviruses [16]. These results indicate (-)-borneol-based esters to have broad-spectrum antiviral activity, acting on an early step of the virus life cycle, and that they should be studied as RSV inhibitors.  The selectivity index (SI, CC50/IC50 ratio) values of each compound were taken as the evaluation reference for antiviral potency. The antiviral activity against filoviruses and SARS-CoV-2 was studied using a "surrogate" system. This system included vesicular stomatitis virus (VSV) particles pseudotyped with Ebola virus GP (rVSV-ΔG-EboV-GP) or Marburg virus GP (rVSV-ΔG-MarV-GP), and lentiviral particles exhibiting the S protein of SARS-CoV-2 (pSARS-CoV-2-S). a SIMarV-GP = CC50/IC50MarV-GP, where IC50MarV-GP is the concentration of a compound required to inhibit rVSV-ΔG-MarV-GP infection of HEK293T cells by 50%. b SIEboV-GP = CC50/IC50EboV-GP, where IC50EboV-GP is the concentration of a compound required to inhibit rVSV-ΔG-EboV-GP infection of HEK293T cells by 50%. c SIH1N1 = CC50/IC50H1N1, where IC50H1N1: con-fants and children [4]. Typically, RSV is mild and accompanied by upper respiratory symptoms, with the risk of severe RSV increased by factors such as young age (<6 months), congenital heart disease, chronic pulmonary diseases, and immunocompromised conditions [5]. Globally, infection with RSV is responsible for approximately 3.2 million hospital admissions and approximately 60,000 in-hospital deaths per year [6,7]. No specific etiotropic therapy is available, and the only compound (Synagis) designed to prevent this infection is practically inaccessible to the general population [8]. Several groups of researchers around the world are actively developing antiviral agents against respiratory syncytial infection, but none of them are yet suitable for wide use in clinical practice [9][10][11].
Recently, we found that (-)-borneol-based esters I-IV (Table 1) containing a saturated nitrogen ring or an acyclic tertiary nitrogen atom exhibited good anti-filovirus efficacy by affecting the glycoprotein (GP)-mediated membrane fusion process [12,13]. At the same time, in vitro results with an infectious virus and in a pseudovirus system, as well as theoretical calculations, allowed us to demonstrate that (-)-borneol esters IV with butyl substituents are likely to be fusion inhibitors targeting the S protein of SARS-CoV-2 [14]. Moreover, we found esters I-IV to effectively inhibit the replication of the influenza A virus [15] and esters III with a morpholine cycle to exhibit antiviral activity against orthopoxviruses [16]. These results indicate (-)-borneol-based esters to have broad-spectrum antiviral activity, acting on an early step of the virus life cycle, and that they should be studied as RSV inhibitors.   NA NA e SIpSARS-CoV-2-S = 18 e SIpSARS-CoV-2-S = 29 SISARS-CoV-2 = 6-102 * The selectivity index (SI, CC50/IC50 ratio) values of each compound were taken as the evaluation reference for antiviral potency. The antiviral activity against filoviruses and SARS-CoV-2 was studied using a "surrogate" system. This system included vesicular stomatitis virus (VSV) particles pseudotyped with Ebola virus GP (rVSV-ΔG-EboV-GP) or Marburg virus GP (rVSV-ΔG-MarV-GP), and lentiviral particles exhibiting the S protein of SARS-CoV-2 (pSARS-CoV-2-S). a SIMarV-GP = CC50/IC50MarV-GP, where IC50MarV-GP is the concentration of a compound required to inhibit rVSV-ΔG-MarV-GP infection of HEK293T cells by 50%. b SIEboV-GP = CC50/IC50EboV-GP, where IC50EboV-GP is the concentration of a compound required to inhibit rVSV-ΔG-EboV-GP infection of HEK293T cells by 50%. c SIH1N1 = CC50/IC50H1N1, where IC50H1N1: con-fants and children [4]. Typically, RSV is mild and accompanied by upper respiratory symptoms, with the risk of severe RSV increased by factors such as young age (<6 months), congenital heart disease, chronic pulmonary diseases, and immunocompromised conditions [5]. Globally, infection with RSV is responsible for approximately 3.2 million hospital admissions and approximately 60,000 in-hospital deaths per year [6,7]. No specific etiotropic therapy is available, and the only compound (Synagis) designed to prevent this infection is practically inaccessible to the general population [8]. Several groups of researchers around the world are actively developing antiviral agents against respiratory syncytial infection, but none of them are yet suitable for wide use in clinical practice [9][10][11].
Recently, we found that (-)-borneol-based esters I-IV (Table 1) containing a saturated nitrogen ring or an acyclic tertiary nitrogen atom exhibited good anti-filovirus efficacy by affecting the glycoprotein (GP)-mediated membrane fusion process [12,13]. At the same time, in vitro results with an infectious virus and in a pseudovirus system, as well as theoretical calculations, allowed us to demonstrate that (-)-borneol esters IV with butyl substituents are likely to be fusion inhibitors targeting the S protein of SARS-CoV-2 [14]. Moreover, we found esters I-IV to effectively inhibit the replication of the influenza A virus [15] and esters III with a morpholine cycle to exhibit antiviral activity against orthopoxviruses [16]. These results indicate (-)-borneol-based esters to have broad-spectrum antiviral activity, acting on an early step of the virus life cycle, and that they should be studied as RSV inhibitors.   NA NA e SIpSARS-CoV-2-S = 18 e SIpSARS-CoV-2-S = 29 SISARS-CoV-2 = 6-102 * The selectivity index (SI, CC50/IC50 ratio) values of each compound were taken as the evaluation reference for antiviral potency. The antiviral activity against filoviruses and SARS-CoV-2 was studied using a "surrogate" system. This system included vesicular stomatitis virus (VSV) particles pseudotyped with Ebola virus GP (rVSV-ΔG-EboV-GP) or Marburg virus GP (rVSV-ΔG-MarV-GP), and lentiviral particles exhibiting the S protein of SARS-CoV-2 (pSARS-CoV-2-S). a SIMarV-GP = CC50/IC50MarV-GP, where IC50MarV-GP is the concentration of a compound required to inhibit rVSV-ΔG-MarV-GP infection of HEK293T cells by 50%. b SIEboV-GP = CC50/IC50EboV-GP, where IC50EboV-GP is the concentration of a compound required to inhibit rVSV-ΔG-EboV-GP infection of HEK293T cells by 50%. c SIH1N1 = CC50/IC50H1N1, where IC50H1N1: con-fants and children [4]. Typically, RSV is mild and accompanied by upper respiratory symptoms, with the risk of severe RSV increased by factors such as young age (<6 months), congenital heart disease, chronic pulmonary diseases, and immunocompromised conditions [5]. Globally, infection with RSV is responsible for approximately 3.2 million hospital admissions and approximately 60,000 in-hospital deaths per year [6,7]. No specific etiotropic therapy is available, and the only compound (Synagis) designed to prevent this infection is practically inaccessible to the general population [8]. Several groups of researchers around the world are actively developing antiviral agents against respiratory syncytial infection, but none of them are yet suitable for wide use in clinical practice [9][10][11].
Recently, we found that (-)-borneol-based esters I-IV (Table 1) containing a saturated nitrogen ring or an acyclic tertiary nitrogen atom exhibited good anti-filovirus efficacy by affecting the glycoprotein (GP)-mediated membrane fusion process [12,13]. At the same time, in vitro results with an infectious virus and in a pseudovirus system, as well as theoretical calculations, allowed us to demonstrate that (-)-borneol esters IV with butyl substituents are likely to be fusion inhibitors targeting the S protein of SARS-CoV-2 [14]. Moreover, we found esters I-IV to effectively inhibit the replication of the influenza A virus [15] and esters III with a morpholine cycle to exhibit antiviral activity against orthopoxviruses [16]. These results indicate (-)-borneol-based esters to have broad-spectrum antiviral activity, acting on an early step of the virus life cycle, and that they should be studied as RSV inhibitors.   NA NA e SIpSARS-CoV-2-S = 18 e SIpSARS-CoV-2-S = 29 SISARS-CoV-2 = 6-102 * The selectivity index (SI, CC50/IC50 ratio) values of each compound were taken as the evaluation reference for antiviral potency. The antiviral activity against filoviruses and SARS-CoV-2 was studied using a "surrogate" system. This system included vesicular stomatitis virus (VSV) particles pseudotyped with Ebola virus GP (rVSV-ΔG-EboV-GP) or Marburg virus GP (rVSV-ΔG-MarV-GP), and lentiviral particles exhibiting the S protein of SARS-CoV-2 (pSARS-CoV-2-S). a SIMarV-GP = CC50/IC50MarV-GP, where IC50MarV-GP is the concentration of a compound required to inhibit rVSV-ΔG-MarV-GP infection of HEK293T cells by 50%. b SIEboV-GP = CC50/IC50EboV-GP, where IC50EboV-GP is the concentration of a compound required to inhibit rVSV-ΔG-EboV-GP infection of HEK293T cells by 50%. c SIH1N1 = CC50/IC50H1N1, where IC50H1N1: con- The selectivity index (SI, CC 50 /IC 50 ratio) values of each compound were taken as the evaluation reference for antiviral potency. The antiviral activity against filoviruses and SARS-CoV-2 was studied using a "surrogate" system. This system included vesicular stomatitis virus (VSV) particles pseudotyped with Ebola virus GP (rVSV-∆G-EboV-GP) or Marburg virus GP (rVSV-∆G-MarV-GP), and lentiviral particles exhibiting the S protein of SARS-CoV-2 (pSARS-CoV-2-S). a SI MarV-GP = CC 50 /IC 50MarV-GP , where IC 50MarV-GP is the concentration of a compound required to inhibit rVSV-∆G-MarV-GP infection of HEK293T cells by 50%. b SI EboV-GP = CC 50 /IC 50EboV-GP , where IC 50EboV-GP is the concentration of a compound required to inhibit rVSV-∆G- In this study, we expanded the library of (-)-borneol-based esters and evaluated their activity against RSV. The anti-RSV structure-activity relationship (SAR) and the antiviral mechanism of the key compounds were analyzed.

Chemistry
A series of (-)-borneol derivatives, 12 of them being novel, were obtained, as shown in Scheme 1. The synthesis of compounds 2-11a-c was reported previously [12]. The synthesis of the derivatives with thiomorpholine 12a-c, 4-benzylpiperidine 13a-c, and Nphenylpiperazine 14-15a-c cycles was carried out as described earlier with minor changes in the second stage. Acetonitrile was used as a solvent instead of methylene chloride. Potassium carbonate was used as a base and potassium iodine was added to activate alkyl halides 2a-c. The products were isolated for column chromatography, with hexane-ethyl acetate used as the eluent.
In this study, we expanded the library of (-)-borneol-based esters and evaluated their activity against RSV. The anti-RSV structure-activity relationship (SAR) and the antiviral mechanism of the key compounds were analyzed.

SAR for the Anti-RSV Activity of the Target Compounds
All target compounds were evaluated for their in vitro antiviral activity against the RSV A strain A2 in HEp-2 cells. The IC 50 and SI values of each compound were taken  Table 2. Table 2. In vitro anti-RSV activity, cytotoxicity, and logP values for compounds 3-15a-c. Compound 12 h, yields 35-53%; (iv) the corresponding amine (1.5 equiv), Et3N (1 equiv), CH2Cl2, room temperature or slight heating for 12-48 h, yields 32-91%. * Compound 1a was used as a commercial reagent.

SAR for the Anti-RSV Activity of the Target Compounds
All target compounds were evaluated for their in vitro antiviral activity against the RSV A strain A2 in НЕр-2 cells. The IC50 and SI values of each compound were taken as criteria for evaluating the antiviral activity. The structures and antiviral activities of compounds 3-15a-c are listed in Table 2. 12 h, yields 35-53%; (iv) the corresponding amine (1.5 equiv), Et3N (1 equiv), CH2Cl2, room temperature or slight heating for 12-48 h, yields 32-91%. * Compound 1a was used as a commercial reagent.

SAR for the Anti-RSV Activity of the Target Compounds
All target compounds were evaluated for their in vitro antiviral activity against the RSV A strain A2 in НЕр-2 cells. The IC50 and SI values of each compound were taken as criteria for evaluating the antiviral activity. The structures and antiviral activities of compounds 3-15a-c are listed in Table 2. perature or slight heating for 12-48 h, yields 32-91%. * Compound 1a was used as a commercial reagent.

SAR for the Anti-RSV Activity of the Target Compounds
All target compounds were evaluated for their in vitro antiviral activity against the RSV A strain A2 in НЕр-2 cells. The IC50 and SI values of each compound were taken as criteria for evaluating the antiviral activity. The structures and antiviral activities of compounds 3-15a-c are listed in Table 2.

SAR for the Anti-RSV Activity of the Target Compounds
All target compounds were evaluated for their in vitro antiviral activity against the RSV A strain A2 in НЕр-2 cells. The IC50 and SI values of each compound were taken as criteria for evaluating the antiviral activity. The structures and antiviral activities of compounds 3-15a-c are listed in Table 2. perature or slight heating for 12-48 h, yields 32-91%. * Compound 1a was used as a commercial reagent.

SAR for the Anti-RSV Activity of the Target Compounds
All target compounds were evaluated for their in vitro antiviral activity against the RSV A strain A2 in НЕр-2 cells. The IC50 and SI values of each compound were taken as criteria for evaluating the antiviral activity. The structures and antiviral activities of compounds 3-15a-c are listed in Table 2. perature or slight heating for 12-48 h, yields 32-91%. * Compound 1a was used as a commercial reagent.

SAR for the Anti-RSV Activity of the Target Compounds
All target compounds were evaluated for their in vitro antiviral activity against the RSV A strain A2 in НЕр-2 cells. The IC50 and SI values of each compound were taken as criteria for evaluating the antiviral activity. The structures and antiviral activities of compounds 3-15a-c are listed in Table 2.

SAR for the Anti-RSV Activity of the Target Compounds
All target compounds were evaluated for their in vitro antiviral activity against the RSV A strain A2 in НЕр-2 cells. The IC50 and SI values of each compound were taken as criteria for evaluating the antiviral activity. The structures and antiviral activities of compounds 3-15a-c are listed in Table 2.

SAR for the Anti-RSV Activity of the Target Compounds
All target compounds were evaluated for their in vitro antiviral activity against the RSV A strain A2 in НЕр-2 cells. The IC50 and SI values of each compound were taken as criteria for evaluating the antiviral activity. The structures and antiviral activities of compounds 3-15a-c are listed in Table 2.

SAR for the Anti-RSV Activity of the Target Compounds
All target compounds were evaluated for their in vitro antiviral activity against the RSV A strain A2 in НЕр-2 cells. The IC50 and SI values of each compound were taken as criteria for evaluating the antiviral activity. The structures and antiviral activities of compounds 3-15a-c are listed in Table 2.

SAR for the Anti-RSV Activity of the Target Compounds
All target compounds were evaluated for their in vitro antiviral activity against the RSV A strain A2 in НЕр-2 cells. The IC50 and SI values of each compound were taken as criteria for evaluating the antiviral activity. The structures and antiviral activities of compounds 3-15a-c are listed in Table 2.

SAR for the Anti-RSV Activity of the Target Compounds
All target compounds were evaluated for their in vitro antiviral activity against the RSV A strain A2 in НЕр-2 cells. The IC50 and SI values of each compound were taken as criteria for evaluating the antiviral activity. The structures and antiviral activities of compounds 3-15a-c are listed in Table 2. a CC50 is the median cytotoxic concentration, i.e., the concentration causing 50% cell death. b IC50 is the 50% inhibitory concentration, i.e., the concentration causing a 50% decrease in virus replication. c SI is the selectivity index, i.e., the CC50/IC50 ratio. CC50 and IC50 values are presented as the mean value ± standard deviation. All the values were calculated based on three independent experiments. logP, the octanol-water partition coefficient, was calculated based on the molecular descriptors of compounds using the plugin QikProp. Green means the activity of the compound is more high.
In the series of (-)-borneol esters 3-4a-c containing alicyclic substituents, derivative 3b, carrying methyl groups at the nitrogen atom was found to be the most potent, with an IC50 value of 8.9 µ M and an SI value of 111. The analogs of derivative 3b with different linker lengths (compounds 3a and 3c) and the analogs with ethyl groups (compounds 4ac) turned out to be more toxic and exhibited lower activity than derivative 3b. Within the subfamily of (-)-borneol esters 5-15a-c with N-containing heterocycles, derivatives 5a-c containing morpholine moieties showed high activity with IC50s of 5.0, 8.5, and 5.8 µ M depending on linker length, with n = 1, 2, and 3, respectively. The morpholine derivatives 5a and 5c, with n = 1 and n = 3, demonstrated low toxicity, resulting in high SI values of 83 and 53, respectively. The introduction of heterocycle that is less polar than morpholine, such as thiomorpholine (derivatives 12a-c) with the value of logPs in the range of 3-4, resulted in potency loss. Meanwhile, among the piperidine derivatives 6a-c, with logP values also in the range of 3-4, the agent 6b with linker length n = 2 showed moderate anti-RSV activity, with an IC50 value of 14.6 µ M and an SI value of 22. The introduction of a methyl group at the 4-position of the piperidine cycle yielded compounds with higher logPs. Compound 7b with linker length n = 2 showed moderate antiviral potency, with an IC50 value of 6.2 µ M and an SI value of 26. Notably, compounds 8a-c, with the azepane cycle differing from the piperidine derivatives by an additional methylene group in the cycle, turned out to be not active. The 4-methylpiperazine derivative 9c with linker length n = 2 had high antiviral activity, with an IC50 value of 5.1 µ M. Additionally, the 4-ethylpiperazine derivatives 10b and 10c (n = 2, n = 3) showed micromolar IC50 values  a CC50 is the median cytotoxic concentration, i.e., the concentration causing 50% cell death. b IC50 is the 50% inhibitory concentration, i.e., the concentration causing a 50% decrease in virus replication. c SI is the selectivity index, i.e., the CC50/IC50 ratio. CC50 and IC50 values are presented as the mean value ± standard deviation. All the values were calculated based on three independent experiments. logP, the octanol-water partition coefficient, was calculated based on the molecular descriptors of compounds using the plugin QikProp. Green means the activity of the compound is more high.
In the series of (-)-borneol esters 3-4a-c containing alicyclic substituents, derivative 3b, carrying methyl groups at the nitrogen atom was found to be the most potent, with an IC50 value of 8.9 µ M and an SI value of 111. The analogs of derivative 3b with different linker lengths (compounds 3a and 3c) and the analogs with ethyl groups (compounds 4ac) turned out to be more toxic and exhibited lower activity than derivative 3b. Within the subfamily of (-)-borneol esters 5-15a-c with N-containing heterocycles, derivatives 5a-c containing morpholine moieties showed high activity with IC50s of 5.0, 8.5, and 5.8 µ M depending on linker length, with n = 1, 2, and 3, respectively. The morpholine derivatives 5a and 5c, with n = 1 and n = 3, demonstrated low toxicity, resulting in high SI values of 83 and 53, respectively. The introduction of heterocycle that is less polar than morpholine, such as thiomorpholine (derivatives 12a-c) with the value of logPs in the range of 3-4, resulted in potency loss. Meanwhile, among the piperidine derivatives 6a-c, with logP values also in the range of 3-4, the agent 6b with linker length n = 2 showed moderate anti-RSV activity, with an IC50 value of 14.6 µ M and an SI value of 22. The introduction of a methyl group at the 4-position of the piperidine cycle yielded compounds with higher logPs. Compound 7b with linker length n = 2 showed moderate antiviral potency, with an IC50 value of 6.2 µ M and an SI value of 26. Notably, compounds 8a-c, with the azepane cycle differing from the piperidine derivatives by an additional methylene group in the cycle, turned out to be not active. The 4-methylpiperazine derivative 9c with linker length n = 2 had high antiviral activity, with an IC50 value of 5.1 µ M. Additionally, the 4-ethylpiperazine derivatives 10b and 10c (n = 2, n = 3) showed micromolar IC50 values  a CC50 is the median cytotoxic concentration, i.e., the concentration causing 50% cell death. b IC50 is the 50% inhibitory concentration, i.e., the concentration causing a 50% decrease in virus replication. c SI is the selectivity index, i.e., the CC50/IC50 ratio. CC50 and IC50 values are presented as the mean value ± standard deviation. All the values were calculated based on three independent experiments. logP, the octanol-water partition coefficient, was calculated based on the molecular descriptors of compounds using the plugin QikProp. Green means the activity of the compound is more high.
In the series of (-)-borneol esters 3-4a-c containing alicyclic substituents, derivative 3b, carrying methyl groups at the nitrogen atom was found to be the most potent, with an IC50 value of 8.9 µ M and an SI value of 111. The analogs of derivative 3b with different linker lengths (compounds 3a and 3c) and the analogs with ethyl groups (compounds 4ac) turned out to be more toxic and exhibited lower activity than derivative 3b. Within the subfamily of (-)-borneol esters 5-15a-c with N-containing heterocycles, derivatives 5a-c containing morpholine moieties showed high activity with IC50s of 5.0, 8.5, and 5.8 µ M depending on linker length, with n = 1, 2, and 3, respectively. The morpholine derivatives 5a and 5c, with n = 1 and n = 3, demonstrated low toxicity, resulting in high SI values of 83 and 53, respectively. The introduction of heterocycle that is less polar than morpholine, such as thiomorpholine (derivatives 12a-c) with the value of logPs in the range of 3-4, resulted in potency loss. Meanwhile, among the piperidine derivatives 6a-c, with logP values also in the range of 3-4, the agent 6b with linker length n = 2 showed moderate anti-RSV activity, with an IC50 value of 14.6 µ M and an SI value of 22. The introduction of a methyl group at the 4-position of the piperidine cycle yielded compounds with higher logPs. Compound 7b with linker length n = 2 showed moderate antiviral potency, with an IC50 value of 6.2 µ M and an SI value of 26. Notably, compounds 8a-c, with the azepane cycle differing from the piperidine derivatives by an additional methylene group in the cycle, turned out to be not active. The 4-methylpiperazine derivative 9c with linker length n = 2 had high antiviral activity, with an IC50 value of 5.1 µ M. Additionally, the 4-ethylpiperazine derivatives 10b and 10c (n = 2, n = 3) showed micromolar IC50 values  a CC50 is the median cytotoxic concentration, i.e., the concentration causing 50% cell death. b IC50 is the 50% inhibitory concentration, i.e., the concentration causing a 50% decrease in virus replication. c SI is the selectivity index, i.e., the CC50/IC50 ratio. CC50 and IC50 values are presented as the mean value ± standard deviation. All the values were calculated based on three independent experiments. logP, the octanol-water partition coefficient, was calculated based on the molecular descriptors of compounds using the plugin QikProp. Green means the activity of the compound is more high.
In the series of (-)-borneol esters 3-4a-c containing alicyclic substituents, derivative 3b, carrying methyl groups at the nitrogen atom was found to be the most potent, with an IC50 value of 8.9 µ M and an SI value of 111. The analogs of derivative 3b with different linker lengths (compounds 3a and 3c) and the analogs with ethyl groups (compounds 4ac) turned out to be more toxic and exhibited lower activity than derivative 3b. Within the subfamily of (-)-borneol esters 5-15a-c with N-containing heterocycles, derivatives 5a-c containing morpholine moieties showed high activity with IC50s of 5.0, 8.5, and 5.8 µ M depending on linker length, with n = 1, 2, and 3, respectively. The morpholine derivatives 5a and 5c, with n = 1 and n = 3, demonstrated low toxicity, resulting in high SI values of 83 and 53, respectively. The introduction of heterocycle that is less polar than morpholine, such as thiomorpholine (derivatives 12a-c) with the value of logPs in the range of 3-4, resulted in potency loss. Meanwhile, among the piperidine derivatives 6a-c, with logP values also in the range of 3-4, the agent 6b with linker length n = 2 showed moderate anti-RSV activity, with an IC50 value of 14.6 µ M and an SI value of 22. The introduction of a methyl group at the 4-position of the piperidine cycle yielded compounds with higher logPs. Compound 7b with linker length n = 2 showed moderate antiviral potency, with an IC50 value of 6.2 µ M and an SI value of 26. Notably, compounds 8a-c, with the azepane cycle differing from the piperidine derivatives by an additional methylene group in the cycle, turned out to be not active. The 4-methylpiperazine derivative 9c with linker length n = 2 had high antiviral activity, with an IC50 value of 5.1 µ M. Additionally, the 4-ethylpiperazine derivatives 10b and 10c (n = 2, n = 3) showed micromolar IC50 values  a CC50 is the median cytotoxic concentration, i.e., the concentration causing 50% cell death. b IC50 is the 50% inhibitory concentration, i.e., the concentration causing a 50% decrease in virus replication. c SI is the selectivity index, i.e., the CC50/IC50 ratio. CC50 and IC50 values are presented as the mean value ± standard deviation. All the values were calculated based on three independent experiments. logP, the octanol-water partition coefficient, was calculated based on the molecular descriptors of compounds using the plugin QikProp. Green means the activity of the compound is more high.
In the series of (-)-borneol esters 3-4a-c containing alicyclic substituents, derivative 3b, carrying methyl groups at the nitrogen atom was found to be the most potent, with an IC50 value of 8.9 µ M and an SI value of 111. The analogs of derivative 3b with different linker lengths (compounds 3a and 3c) and the analogs with ethyl groups (compounds 4ac) turned out to be more toxic and exhibited lower activity than derivative 3b. Within the subfamily of (-)-borneol esters 5-15a-c with N-containing heterocycles, derivatives 5a-c containing morpholine moieties showed high activity with IC50s of 5.0, 8.5, and 5.8 µ M depending on linker length, with n = 1, 2, and 3, respectively. The morpholine derivatives 5a and 5c, with n = 1 and n = 3, demonstrated low toxicity, resulting in high SI values of 83 and 53, respectively. The introduction of heterocycle that is less polar than morpholine, such as thiomorpholine (derivatives 12a-c) with the value of logPs in the range of 3-4, resulted in potency loss. Meanwhile, among the piperidine derivatives 6a-c, with logP values also in the range of 3-4, the agent 6b with linker length n = 2 showed moderate anti-RSV activity, with an IC50 value of 14.6 µ M and an SI value of 22. The introduction of a methyl group at the 4-position of the piperidine cycle yielded compounds with higher logPs. Compound 7b with linker length n = 2 showed moderate antiviral potency, with an IC50 value of 6.2 µ M and an SI value of 26. Notably, compounds 8a-c, with the azepane cycle differing from the piperidine derivatives by an additional methylene group in the cycle, turned out to be not active. The 4-methylpiperazine derivative 9c with linker length n = 2 had high antiviral activity, with an IC50 value of 5.1 µ M. Additionally, the 4-ethylpiperazine derivatives 10b and 10c (n = 2, n = 3) showed micromolar IC50 values and SI values of 20 and 21. It is worth mentioning that the logP values of piperazine derivatives 9-10a-c are similar to those of active morpholine derivatives 5a-c.
As mentioned above, (-)-borneol esters act in the early stage of viral life, with surface a CC50 is the median cytotoxic concentration, i.e., the concentration causing 50% cell death. b IC50 is the 50% inhibitory concentration, i.e., the concentration causing a 50% decrease in virus replication. c SI is the selectivity index, i.e., the CC50/IC50 ratio. CC50 and IC50 values are presented as the mean value ± standard deviation. All the values were calculated based on three independent experiments. logP, the octanol-water partition coefficient, was calculated based on the molecular descriptors of compounds using the plugin QikProp. Green means the activity of the compound is more high.
In the series of (-)-borneol esters 3-4a-c containing alicyclic substituents, derivative 3b, carrying methyl groups at the nitrogen atom was found to be the most potent, with an IC50 value of 8.9 µ M and an SI value of 111. The analogs of derivative 3b with different linker lengths (compounds 3a and 3c) and the analogs with ethyl groups (compounds 4ac) turned out to be more toxic and exhibited lower activity than derivative 3b. Within the subfamily of (-)-borneol esters 5-15a-c with N-containing heterocycles, derivatives 5a-c containing morpholine moieties showed high activity with IC50s of 5.0, 8.5, and 5.8 µ M depending on linker length, with n = 1, 2, and 3, respectively. The morpholine derivatives 5a and 5c, with n = 1 and n = 3, demonstrated low toxicity, resulting in high SI values of 83 and 53, respectively. The introduction of heterocycle that is less polar than morpholine, such as thiomorpholine (derivatives 12a-c) with the value of logPs in the range of 3-4, resulted in potency loss. Meanwhile, among the piperidine derivatives 6a-c, with logP values also in the range of 3-4, the agent 6b with linker length n = 2 showed moderate anti-RSV activity, with an IC50 value of 14.6 µ M and an SI value of 22. The introduction of a methyl group at the 4-position of the piperidine cycle yielded compounds with higher logPs. Compound 7b with linker length n = 2 showed moderate antiviral potency, with an IC50 value of 6.2 µ M and an SI value of 26. Notably, compounds 8a-c, with the azepane cycle differing from the piperidine derivatives by an additional methylene group in the cycle, turned out to be not active. The 4-methylpiperazine derivative 9c with linker length n = 2 had high antiviral activity, with an IC50 value of 5.1 µ M. Additionally, the a CC50 is the median cytotoxic concentration, i.e., the concentration causing 50% cell death. b IC50 is the 50% inhibitory concentration, i.e., the concentration causing a 50% decrease in virus replication. c SI is the selectivity index, i.e., the CC50/IC50 ratio. CC50 and IC50 values are presented as the mean value ± standard deviation. All the values were calculated based on three independent experiments. logP, the octanol-water partition coefficient, was calculated based on the molecular descriptors of compounds using the plugin QikProp. Green means the activity of the compound is more high.
In the series of (-)-borneol esters 3-4a-c containing alicyclic substituents, derivative 3b, carrying methyl groups at the nitrogen atom was found to be the most potent, with an IC50 value of 8.9 µ M and an SI value of 111. The analogs of derivative 3b with different linker lengths (compounds 3a and 3c) and the analogs with ethyl groups (compounds 4ac) turned out to be more toxic and exhibited lower activity than derivative 3b. Within the subfamily of (-)-borneol esters 5-15a-c with N-containing heterocycles, derivatives 5a-c containing morpholine moieties showed high activity with IC50s of 5.0, 8.5, and 5.8 µ M depending on linker length, with n = 1, 2, and 3, respectively. The morpholine derivatives 5a and 5c, with n = 1 and n = 3, demonstrated low toxicity, resulting in high SI values of 83 and 53, respectively. The introduction of heterocycle that is less polar than morpholine, such as thiomorpholine (derivatives 12a-c) with the value of logPs in the range of 3-4, resulted in potency loss. Meanwhile, among the piperidine derivatives 6a-c, with logP values also in the range of 3-4, the agent 6b with linker length n = 2 showed moderate anti-RSV activity, with an IC50 value of 14.6 µ M and an SI value of 22. The introduction of a methyl group at the 4-position of the piperidine cycle yielded compounds with higher logPs. Compound 7b with linker length n = 2 showed moderate antiviral potency, with an IC50 value of 6.2 µ M and an SI value of 26. Notably, compounds 8a-c, with the azepane cycle differing from the piperidine derivatives by an additional methylene group in the cycle, turned out to be not active. The 4-methylpiperazine derivative 9c with linker length n = 2 had high antiviral activity, with an IC50 value of 5.1 µ M. Additionally, the  a CC50 is the median cytotoxic concentration, i.e., the concentration causing 50% cell death. b IC50 is the 50% inhibitory concentration, i.e., the concentration causing a 50% decrease in virus replication. c SI is the selectivity index, i.e., the CC50/IC50 ratio. CC50 and IC50 values are presented as the mean value ± standard deviation. All the values were calculated based on three independent experiments. logP, the octanol-water partition coefficient, was calculated based on the molecular descriptors of compounds using the plugin QikProp. Green means the activity of the compound is more high.
In the series of (-)-borneol esters 3-4a-c containing alicyclic substituents, derivative 3b, carrying methyl groups at the nitrogen atom was found to be the most potent, with an IC50 value of 8.9 µ M and an SI value of 111. The analogs of derivative 3b with different linker lengths (compounds 3a and 3c) and the analogs with ethyl groups (compounds 4ac) turned out to be more toxic and exhibited lower activity than derivative 3b. Within the subfamily of (-)-borneol esters 5-15a-c with N-containing heterocycles, derivatives 5a-c containing morpholine moieties showed high activity with IC50s of 5.0, 8.5, and 5.8 µ M depending on linker length, with n = 1, 2, and 3, respectively. The morpholine derivatives 5a and 5c, with n = 1 and n = 3, demonstrated low toxicity, resulting in high SI values of 83 and 53, respectively. The introduction of heterocycle that is less polar than morpholine, such as thiomorpholine (derivatives 12a-c) with the value of logPs in the range of 3-4, resulted in potency loss. Meanwhile, among the piperidine derivatives 6a-c, with logP values also in the range of 3-4, the agent 6b with linker length n = 2 showed moderate anti-RSV activity, with an IC50 value of 14.6 µ M and an SI value of 22. The introduction of a methyl group at the 4-position of the piperidine cycle yielded compounds with higher logPs. Compound 7b with linker length n = 2 showed moderate antiviral potency, with an IC50 value of 6.2 µ M and an SI value of 26. Notably, compounds 8a-c, with the azepane cycle differing from the piperidine derivatives by an additional methylene group in the cycle, turned out to be not active. The 4-methylpiperazine derivative 9c with linker length n = 2 had high antiviral activity, with an IC50 value of 5.1 µ M. Additionally, the 3 5.213 27.8 ± 7.4 14.7 ± 2.1 2 Ribavirin <4000 80.1 ± 13.5 50 a CC 50 is the median cytotoxic concentration, i.e., the concentration causing 50% cell death. b IC 50 is the 50% inhibitory concentration, i.e., the concentration causing a 50% decrease in virus replication. c SI is the selectivity index, i.e., the CC 50 /IC 50 ratio. CC 50 and IC 50 values are presented as the mean value ± standard deviation. All the values were calculated based on three independent experiments. logP, the octanol-water partition coefficient, was calculated based on the molecular descriptors of compounds using the plugin QikProp. Green means the activity of the compound is more high.
In the series of (-)-borneol esters 3-4a-c containing alicyclic substituents, derivative 3b, carrying methyl groups at the nitrogen atom was found to be the most potent, with an IC 50 value of 8.9 µM and an SI value of 111. The analogs of derivative 3b with different linker lengths (compounds 3a and 3c) and the analogs with ethyl groups (compounds 4a-c) turned out to be more toxic and exhibited lower activity than derivative 3b. Within the subfamily of (-)-borneol esters 5-15a-c with N-containing heterocycles, derivatives 5a-c containing morpholine moieties showed high activity with IC 50 s of 5.0, 8.5, and 5.8 µM depending on linker length, with n = 1, 2, and 3, respectively. The morpholine derivatives 5a and 5c, with n = 1 and n = 3, demonstrated low toxicity, resulting in high SI values of 83 and 53, respectively. The introduction of heterocycle that is less polar than morpholine, such as thiomorpholine (derivatives 12a-c) with the value of logPs in the range of 3-4, resulted in potency loss. Meanwhile, among the piperidine derivatives 6a-c, with logP values also in the range of 3-4, the agent 6b with linker length n = 2 showed moderate anti-RSV activity, with an IC 50 value of 14.6 µM and an SI value of 22. The introduction of a methyl group at the 4-position of the piperidine cycle yielded compounds with higher logPs. Compound 7b with linker length n = 2 showed moderate antiviral potency, with an IC 50 value of 6.2 µM and an SI value of 26. Notably, compounds 8a-c, with the azepane cycle differing from the piperidine derivatives by an additional methylene group in the cycle, turned out to be not active. The 4-methylpiperazine derivative 9c with linker length n = 2 had high antiviral activity, with an IC 50 value of 5.1 µM. Additionally, the 4-ethylpiperazine derivatives 10b and 10c (n = 2, n = 3) showed micromolar IC 50 values and SI values of 20 and 21. It is worth mentioning that the logP values of piperazine derivatives 9-10a-c are similar to those of active morpholine derivatives 5a-c.
As mentioned above, (-)-borneol esters act in the early stage of viral life, with surface glycoproteins being the likely targets of these compounds. Therefore, the fusion (F) glycoprotein on the RSV virion surface appears to be a likely target for synthesized (-)borneol derivatives. The active binding site of the F protein of known inhibitors is enriched with amino acid residues containing aromatic fragments [17]. In this regard, we focused our attention on N-substituted piperazines 11a-c, 14-15a-c, and 4-substituted piperidines 13a-c with aromatic rings. However, as shown in Table 2, aromatic substituents in the piperazine or piperidine cycle led to increased toxicity, with the CC 50 values between 20.1 and 116.4 µM, excluding compound 14b with a CC 50 value of 365.7 µM. Additionally, compound 14b exhibited the best antiviral activity, with an IC 50 value of 8.1 µM and an SI value of 45. Other derivatives in this subseries of (-)-borneol derivatives showed lower SI values. It is worth noting that the high toxicity of (-)-borneol esters 11a-c was observed earlier in HEK293T cells during a study of anti-filovirus activity, with CC 50 values ranging from 30.1 to 59.4 µM [12]. The high toxicity may be due to higher lipophilicity. Highly lipophilic (logP > 5) compounds tend to bind to hydrophobic targets other than the desired one; therefore, there is an increased risk of toxicity [18].
Overall, the morpholine cycle and methyl substitutions at the nitrogen atom were the most favored structural fragments for the best antiviral activity. The introduction of a more hydrophobic moiety directly next to the aliphatic linker, such as ethylene groups, thiomorpholine, and the piperidine cycle, resulted in analogs with low potencies. Installing a piperazine cycle allowed introducing various substituents at the nitrogen atom, and several derivatives with good antiviral activity were found in this structure. Notably, the logP values for highly active compounds were in the range of 2.1-3.9, which is similar to the logP values for most F protein inhibitors in phase II clinical trials [19].

Time-of-Addition Assay
After a detailed SAR study, compounds 3b, 5a, and 6b with high SI values were selected as key compounds to further investigate the mechanism of RSV replication inhibition. To identify the likely stage of the RSV life cycle that is blocked by the key compounds, we first evaluated their antiviral activity using a time-of-addition assay in HEp-2 cells. The viral life cycle of RSV follows the usual basic stages of the enveloped virus. Within the first hour after binding to the receptor, the virion penetrates the cell cytoplasm, with the subsequent release of the genome. Furthermore, within 2-20 h after penetration into the cell, the intensive replication and assembly of new viral particles occurs. After 20 h, the release of new viral particles on the cell surface begins. The time-of-addition (TOA) approach determines how long the addition of a compound can be postponed before it loses its antiviral activity. The dynamics of virus growth depending on the time of drug addition is shown in Figure 1. When added at 1 h post infection, compounds 3b, 5a, and 6b were found to exhibit a significant reduction in viral titer. The antiviral activity was maintained when they were added to the cells at up to 2 h post infection. When the compounds tested were added at later time points, they lost their inhibitory activity. These data suggest that the anti-RSV activity of compounds 3b, 5a, and 6b is manifested at the viral entry step, probably due to a fusion protein (F) or G-surface glycoprotein being inhibited.

Temperature Shift Assay
The F protein plays an essential role in virus entry into host cells together with the surface glycoprotein G. The antiviral assay was performed with temperature modification to exclude the inhibition of surface glycoprotein G as a possible target of action for compounds 3b, 5a, and 6b. RSV entry into host cells involves two separate events: viruscell fusion, with glycoprotein G being the primary functional component, and cell-cell fusion, with F protein promoting the fusion of virus and host cell membranes. While RSV attachment to cells can occur at 4 °C, its fusion to the target cell membrane occurs at temperatures above 18 °C. Therefore, a compound to block virus attachment should have antiviral activity when added during virus incubation at 4 °C, followed by washing until the incubation temperature is shifted to 37 °C. Thus, to investigate the effect of compounds 3b, 5a, and 6b on virus attachment, we tested the antiviral activity under experimental conditions (A) and (B) (Figure 2a). The conditions (A) were as follows: НЕр-2 cells were incubated with RSV at 4 °C for 1 h to ensure the virus binding while preventing virus internalization; unbound virus and the tested compounds were washed away, and the cells were then shifted to 37 °C for the rest of the experiment. The conditions (B) involved the standard protocol for anti-RSV activity measures. As illustrated in Figure  2b, compounds 3b, 5a, and 6b showed no anti-RSV activity under the experimental conditions (A). The antiviral activity was maintained in experiment (B) when the unbound virus was washed away, and the compounds tested were added to the cells. These findings confirm the possibility of post-binding events at the viral entry stage being inhibited by the compounds tested.

Temperature Shift Assay
The F protein plays an essential role in virus entry into host cells together with the surface glycoprotein G. The antiviral assay was performed with temperature modification to exclude the inhibition of surface glycoprotein G as a possible target of action for compounds 3b, 5a, and 6b. RSV entry into host cells involves two separate events: virus-cell fusion, with glycoprotein G being the primary functional component, and cell-cell fusion, with F protein promoting the fusion of virus and host cell membranes. While RSV attachment to cells can occur at 4 • C, its fusion to the target cell membrane occurs at temperatures above 18 • C. Therefore, a compound to block virus attachment should have antiviral activity when added during virus incubation at 4 • C, followed by washing until the incubation temperature is shifted to 37 • C. Thus, to investigate the effect of compounds 3b, 5a, and 6b on virus attachment, we tested the antiviral activity under experimental conditions (A) and (B) (Figure 2a). The conditions (A) were as follows: HEp-2 cells were incubated with RSV at 4 • C for 1 h to ensure the virus binding while preventing virus internalization; unbound virus and the tested compounds were washed away, and the cells were then shifted to 37 • C for the rest of the experiment. The conditions (B) involved the standard protocol for anti-RSV activity measures. As illustrated in Figure 2b, compounds 3b, 5a, and 6b showed no anti-RSV activity under the experimental conditions (A). The antiviral activity was maintained in experiment (B) when the unbound virus was washed away, and the compounds tested were added to the cells. These findings confirm the possibility of post-binding events at the viral entry stage being inhibited by the compounds tested.

Binding Site Analysis
Given the results of the time-of-addition assay, temperature shift study, and previously reported data, the F protein was suggested to be the likely target of the compounds under test. The RSV F protein is a main surface glycoprotein that mediates the fusion of the viral membrane with a host cell membrane. It is structurally similar to F proteins from other Paramyxoviridae (hydrophobic domains, heptad repeats, cysteine residues, etc.) and has the same manner of proteolytic activation, resulting in the exposition of a hydrophobic fusion peptide. The inactive precursor F0 is cleaved by the furin-like protease into an N-terminal F2 subunit and a C-terminal membrane-anchored F1 subunit carrying the fusion peptide. The fusion peptide is unique in that it is capable of causing membrane fusion even without G glycoprotein.
The RSV F protein is used as a major target for antivirals and vaccine develop-ment because of its importance in the viral replication cycle, its conserved sequence and structure, its exposed position in the virion, and its strong immunogenicity.
Therefore, we performed molecular docking studies of synthesized compounds in the active binding site of known inhibitors of the F protein. The known inhibitors of the F protein, Rilematovir [20], Sisunatovir (RV521) [19] and Presatovir [21], were found to target the same RSV F hydrophobic cavity DS-Cav1. The cavity DS-Cav1 [22] was located in the region of a fusion peptide at the N-terminus of the F1-subunit and the α-helix of the heptad repeat (HRB) at the C-terminus of the F1-subunit. To localize the binding site, the complex inhibitor RV521 bound to the F protein (PDB code 7KQD) was used [19]. The symmetrical cavity of the binding site of F protein inhibitors was located inside the trimer ( Figure 3A), at the N-terminal of the F1 or fusion peptide (amino acid residues: Phe137, Phe140, and Leu141) and the beginning of HBR (amino acid residues: Met396, Thr400, Asp486, Glu487, and Phe488). The RV521 reference inhibitor was located symmetrically at the binding site. Thus, the hydrophobic portion of the molecule was located close to the fusion peptide, in the hydrophobic cavity of the site ( Figure 3B). The benzimidazole

Molecular Modeling 2.4.1. Binding Site Analysis
Given the results of the time-of-addition assay, temperature shift study, and previously reported data, the F protein was suggested to be the likely target of the compounds under test. The RSV F protein is a main surface glycoprotein that mediates the fusion of the viral membrane with a host cell membrane. It is structurally similar to F proteins from other Paramyxoviridae (hydrophobic domains, heptad repeats, cysteine residues, etc.) and has the same manner of proteolytic activation, resulting in the exposition of a hydrophobic fusion peptide. The inactive precursor F0 is cleaved by the furin-like protease into an Nterminal F2 subunit and a C-terminal membrane-anchored F1 subunit carrying the fusion peptide. The fusion peptide is unique in that it is capable of causing membrane fusion even without G glycoprotein.
The RSV F protein is used as a major target for antivirals and vaccine develop-ment because of its importance in the viral replication cycle, its conserved sequence and structure, its exposed position in the virion, and its strong immunogenicity.
Therefore, we performed molecular docking studies of synthesized compounds in the active binding site of known inhibitors of the F protein. The known inhibitors of the F protein, Rilematovir [20], Sisunatovir (RV521) [19] and Presatovir [21], were found to target the same RSV F hydrophobic cavity DS-Cav1. The cavity DS-Cav1 [22] was located in the region of a fusion peptide at the N-terminus of the F 1 -subunit and the α-helix of the heptad repeat (HRB) at the C-terminus of the F 1 -subunit. To localize the binding site, the complex inhibitor RV521 bound to the F protein (PDB code 7KQD) was used [19]. The symmetrical cavity of the binding site of F protein inhibitors was located inside the trimer ( Figure 3A), at the N-terminal of the F1 or fusion peptide (amino acid residues: Phe137, Phe140, and Leu141) and the beginning of HBR (amino acid residues: Met396, Thr400, Asp486, Glu487, and Phe488). The RV521 reference inhibitor was located symmetrically at the binding site. Thus, the hydrophobic portion of the molecule was located close to the fusion peptide, in the hydrophobic cavity of the site ( Figure 3B). The benzimidazole fragment was surrounded by three phenylalanines 488, with which it formed pi-pi stacking interactions ( Figure 4A). The protonated -NH3+ group was located in the donor-acceptor region of the site (Figure 3B), forming a salt bridge with Asp486 ( Figure 4A). The RMSD value for the redocking procedure was 0.04 Å2, indicating the correctness of the docking protocol.
Pharmaceuticals 2022, 15, x FOR PEER REVIEW 10 of 21 fragment was surrounded by three phenylalanines 488, with which it formed pi-pi stacking interactions ( Figure 4A). The protonated -NH3+ group was located in the donor-acceptor region of the site (Figure 3B), forming a salt bridge with Asp486 ( Figure 4A). The RMSD value for the redocking procedure was 0.04 Å 2, indicating the correctness of the docking protocol.

Docking Study
The DS-Cav1 binding site was considered a probable site for the binding of (-)-borneol derivatives. Ten docking positions were set as the maximum possible within the framework of the IFD protocol. For some ligands, fewer than 10 positions were implemented (Table S1). Most of the compounds could form hydrogen bonds or salt bridges with the amino acid HRB (Table S1). Salt bridges were formed between the protonated nitrogen atom and the negatively charged amino acids, most often with Asp and Glu. The

Docking Study
The DS-Cav1 binding site was considered a probable site for the binding of (-)-borneol derivatives. Ten docking positions were set as the maximum possible within the framework of the IFD protocol. For some ligands, fewer than 10 positions were implemented (Table S1). Most of the compounds could form hydrogen bonds or salt bridges with the amino acid HRB (Table S1). Salt bridges were formed between the protonated nitrogen atom and the negatively charged amino acids, most often with Asp and Glu. The presence of an aromatic fragment in the compounds (11a-c, 13-15a-c) was characterized by the formation of pi-pi stacking interactions with Phe of the fusion peptide and/or heptad repeat. In general, a correlation was observed between the molecular docking results and in vitro antiviral activity ( Figure S3). Compounds that effectively inhibited RSV were characterized by increased affinity for the binding site. Conversely, inactive compounds showed worse binding affinity to the binding site.
The analysis of the binding mode of lead compounds 3b and 5a in the DS-Cav1 active site revealed that the arrangement was similar to that of the reference inhibitor RV521. The hydrophobic part of compound 3b was located in the hydrophobic cavity near the N-terminus of the F 1 subunit ( Figure 3C). The positively charged nitrogen of compound 3b was located in the donor-acceptor region, next to the negatively charged amino acid residues Glu487 and Asp486. The hydrophobic 1,7,7-trimethylbicyclo [2.2.1] heptane fragment was surrounded by the amino acid residues of Phe488. In contrast to the RV521 reference inhibitor, molecule 3b was located closer to the HBR. The key interactions were the hydrogen bond with Asp486 and the salt bridge with Glu487. Given that molecule 3b was smaller, it did not occupy the entire binding site ( Figure 4B). An interesting fact is that molecule 3a with a shorter linker (n = 1) and molecule 3c with a longer linker (n = 3) were characterized by similar key interactions, but worse energy parameters (Table S1). We attribute this to the possibility of protonation. For example, deprotonated molecule 3a was characterized only by hydrophobic interactions with hydrophobic phenylalanine at the binding site ( Figure S1), which most likely affects the energy parameters of binding (Table S1). At the same time, protonated compound 3a formed additional symmetrical salt bridges with Asp486 A and B chains.
A similar situation was observed for compounds 5a-c. Compound 5b protonated more easily (pKa = 7.91) than compounds 5a (pKa = 7.05) and 5c (pKa = 6.21) (Table S1). Compounds 5a-c were characterized by low inhibitory IC 50 values. However, the low toxicity of 5a allowed this compound to be considered as a lead compound. The location in the binding site of deprotonated and protonated forms of compound 5a was characterized primarily by the formation of hydrophobic interactions. In the case of the deprotonated form, the formation of hydrogen bridges between the oxygen atom of the morpholine fragment and the amino acid residues Phe137 and Arg339 of one of the protomers ( Figure S2A) was observed. The binding site arrangement of the protonated form 5a was similar to that of 3b ( Figure 4C): the hydrophobic 1,7,7-trimethylbicyclo [2.2.1] heptane fragment was surrounded by Phe488 residues of each protomer, with the protonated nitrogen atom of the morpholine fragment forming a hydrogen bond with Asp486 and the Asp486 salt bridge of the adjacent protomer ( Figure S2B). Carbonyl oxygen formed hydrogen bridges with Phe488 and Asp489 ( Figure 4C).

Molecular Dynamics (MD) Simulation of Compounds 3b and 5a
In order to evaluate the location of ligands 3b and 5a within the active binding site, MD simulation was conducted. The model system was created in such a way that the transmembrane domain (amino acids 525-550) of the F protein was immersed in the viral membrane ( Figure S4). We believe that in MD simulations of such proteins, the presence of a membrane is necessary for the correct interpretation of the secondary structure of the protein during the entire simulation time. The root mean square deviation (RMSD) was taken as an important basis to measure the stability of the system. According to the RMSD analysis, both systems equalized to 60 ns of simulation ( Figures S5 and S6). Ligands 3b and 5a were located inside a symmetrical binding site throughout the simulation time. The result of the population analysis indicates minor movements of ligands within the binding site ( Figure 5A).
MD simulation was conducted. The model system was created in such a way that the transmembrane domain (amino acids 525-550) of the F protein was immersed in the viral membrane ( Figure S4). We believe that in MD simulations of such proteins, the presence of a membrane is necessary for the correct interpretation of the secondary structure of the protein during the entire simulation time. The root mean square deviation (RMSD) was taken as an important basis to measure the stability of the system. According to the RMSD analysis, both systems equalized to 60 ns of simulation ( Figures S5 and S6). Ligands 3b and 5a were located inside a symmetrical binding site throughout the simulation time. The result of the population analysis indicates minor movements of ligands within the binding site ( Figure 5A).  Figure 5B shows an analysis of the location of ligands 3b and 5a in the binding site starting from 60 ns of simulation. In both cases, the hydrophobic part of the molecules was located close to the fusion peptide, next to Phe137 in the case of 3b and next to Phe140 in the case of 5a. The positively charged parts of the molecule formed hydrogen bridges with negatively charged amino acids. Similar hydrogen bonds with Asp486 and Glu487 were observed for the binding of JNJ-240868 [23] and TMC-353121 inhibitors [17] containing a protonated amino group and a protonated morpholine fragment, respectively. Moreover, it was reported that the mechanism of inhibitory action of known inhibitors of the F protein might be due to the effect of ligands on the conformation of amino acids Phe488, Asp489, and Thr400 [17]. Compounds 3b and 5a formed intermolecular interactions with at least two of these amino acids ( Figure 5B).
Thus, the results of molecular docking and MD simulations indicate the following: (1) the probable binding site of (-)-borneol derivatives is a symmetrical hydrophobic cavity located inside the F protein between the fusion peptide and the heptad repeat B; (2) lead compounds 3b and 5a form strong intermolecular interactions with functional amino acids that persist for a long period of simulation. At the same time, despite the significant structural differences between the compounds being studied and known inhibitors, similar intermolecular interactions (hydrogen and salt bridges, pi-cation stacking bridges) with key amino acids have been recorded.

Chemistry
Reagents and solvents were purchased from commercial suppliers and used as received. Reactions were monitored by gas chromatography-mass spectrometry (GC-MS), using a 7820A gas chromatograph (Agilent Technologies, Santa Clara, CA, USA); flame- and coupling constants (J) were measured in Hertz. Splitting patterns were designated as follows: s, singlet; d, doublet; t, triplet; m, multiplet. High-resolution mass spectra (HR-MS) were recorded on a DFS Thermo Scientific spectrometer (Agilent Tech., Santa Clara, CA, USA) in full scan mode (0-500 m/z, 70 eV electron impact ionization, and direct sample administration). The purity of the target compounds was determined by a gas chromatography method. All of the target compounds reported in this paper had a purity of ≥95%. Compounds 2-11a-c were obtained according to previously described methods [12].

General Procedure for the Synthesis of Compounds 12-15a-c
The general procedure for the synthesis of compounds 12a-15a was as follows. (1S,2R,4S)-1,7,7-Trimethylbicyclo [2.2.1]heptan-2-yl 2-chloroacetate 2a (1 equiv) was dissolved in CH 3 CN and the appropriate amine (1.2 equiv), K 2 CO 3 (4 equiv) and KI (1 equiv) were added. The reaction mixture was stirred at room temperature for 12 h. After the reaction was completed, the organic solvent was evaporated. The residue was diluted with brine, extracted twice with CH 2 Cl 2 , dried over anhydrous sodium sulfate, and evaporated to dryness. The obtained product was purified via silica gel column chromatography (hexane-ethyl acetate eluent).
The general procedure for the synthesis of compounds 12-15b-c was as follows.
to a concentration of 1000 µg/mL, and a series of 2-fold dilutions was prepared. The final concentration of DMSO was more than twice lower than its 50% cytotoxic con-centration, and thus, could not affect cell viability. A one-day culture of HEp2 cells, grown in 96-well plates, with cell concentration 3 × 10 5 /well of the plate, was checked visually in an inverted microscope for the integrity of the monolayer. Plates were selected for work, where the cell closure was 60-80%. Dilutions of the compounds at the appropriate concentration were added to the plate in a volume of 100 µL in each well in 2 replicates for each tested concentration. The plates were incubated for 24 h at 37 • C in the presence of 5% CO 2 . Cell viability was assessed using the MTT test. The MTT solution was prepared on maintenance medium at a concentration of 0.5 mg/mL. Then, 0.1 mL of MTT solution was added to each well. After 1.5 h of MTT contact at 37 • C at a CO 2 concentration of 5%, MTT was discarded with the cells of the well, and 0.1 mL of 96% ethyl alcohol was poured, after which the optical density in the wells was measured at 535 nm. Based on the data obtained, the CC 50 was calculated.

Antiviral Assay
The respiratory syncytial virus (strain A2) and cell culture HEp-2 were obtained from the working collection of the Laboratory of Chemotherapy for Viral Infections of Smorodintsev Research Institute of Influenza. The antiviral activity against the respiratory syncytial virus (RSV A strain A2) was assessed in a series of 3-fold dilutions of test compounds, starting from 1 2 CC 50 , which were added to HEp-2 cell culture at a double concentration, 100 µL per well, followed by an addition of 100 µL of the virus in a series of 10-fold dilutions. The cells were incubated at 37 • C and 5% CO 2 for 1 h. Then, the virus was washed out, and the compounds were again added at a single concentration and incubated at 37 • C and 5% CO 2 for 6 days. For the enzyme-linked immunosorbent assay (ELISA), the cell culture was fixed with cold 80% acetone at −20 • C for 15 min and then washed with phosphate-buffered saline containing 0.05% Tween 20. Next, a solution of primary mouse anti-RSV F protein antibodies was added to the culture and incubated at room temperature under continuous stirring for 2 h. Then, the cells were again washed with buffer, secondary anti-mouse antibodies were added, and the cells were incubated under continuous stirring for 2 h. Then, the antibodies were washed off, and a substrate-chromogenic mixture with tetramethylbenzidine was added. After 5 min, the reaction was stopped with 0.1 M of sulfuric acid, and the optical density of the solution was measured at a wavelength of 450 nm. Wells with absorbance values of two-fold or greater than the cell control were considered contaminated. The virus titer was calculated using the Reed and Muench method. All experiments were performed in triplicate.

Time-of-Addition Assay
Compounds 3b, 5a, and 6b were added at different time points before, after, or simultaneously with the introduction of the virus. The time of addition of the compound was counted from point 0-the time of the entry of the virus into the cell. During the period (-1)-0, the cells together with the virus were incubated at 40 • C. All other experiments were carried out at 37 • C. RSV virus A (2 mL) was added to the cells at a time that was conventionally designated as point -1, after which the cells were kept for an hour at a temperature of 40 • C. Then, at point 0, the virus was unbound. The cells were transferred to a thermostat at 37 • C, where they were incubated for 25 h. After this period, the medium was taken from each well, and a series of ten-fold dilutions were made on fresh cell culture and incubated for 6 days. For each compound, 2 repetitions were made by different operators. The virus titer was estimated by ELISA. The compounds were added at the following times relative to the addition of the virus: point -2-the compound was introduced one hour before cell infection (prophylactic regimen); point 0-at the moment of temperature change; points 1, 2, 4, 6, 24, at 1, 2, 4, 6, and 24 h after the temperature change, respectively. In the wells marked (-2)-(25), the compound was kept throughout the experiment, starting from point -2 and until the end of the experiment-25 h. No compound was added to the control wells. Instead, a similar volume of medium was added.

Ligands and Protein Preparations
The geometric parameters of the ligands were optimized by the OPLS4 force field method [24], considering possible conformers. To estimate the ADME parameters of the ligands, we used the QikProp plugin. The ADME estimation method is based on the molecular description of structures. Acidity constants for some compounds were calculated using quantum-chemical software Jaguar [25]. The most popular DFT method, B3LYP [26] with Danning's basis set cc-pVTZ(+) [27], was used for the optimization procedure. The contribution of protonated and deprotonated molecules in the molecular docking procedure was considered on the basis of the Henderson-Hasselbalch equation. The geometrical parameters of the full-size F protein (PDB code 7LVW [28]) of human RVS A2 were downloaded from the noncommercial base data Protein Data Bank [29]. To localize the binding site complex, the protomer-inhibitor RV521 was used (PDB code 7KQD [19]). Model structures were prepared using the plugin Protein Preparation Wizard: hydrogen atoms were added and minimized, side chains of amino acids were edited, multiple chemical bonds were restored, and water molecules were removed. Geometric parameters were optimized in the OPLS4 force field [24].

Molecular Docking
All theoretical calculations were carried out using Schrodinger Small Molecule Drug Discovery Release 2022-1 software. The molecular docking of compounds was carried out using the induced fit docking (IFD) protocol with the following conditions: flexible protein and ligand, grid matrix size of 15 Å, and amino acids within a radius of 5 Å from the ligand were optimized considering the influence of the ligand. The ranking of docking solutions was carried out by evaluating the following calculation parameters: the docking score (based on GlideScore minus fines), the ligand efficiency (LE) and the parameter of the model energy value (Emodel), including the GlideScore value, the energy of unrelated interactions, and the parameters of the energy spent for the formation of the binding site stacking.

Molecular Dynamics
Trimer-ligand complexes for compounds 3b and 5a were used for subsequent molecular dynamics simulations. Given the structure of the full-length F protein, the transmembrane domain (525-550 amino acids) was placed in the POPC membrane (phosphatidylcholine). Phosphatidylcholine is a part of most cell membranes of viruses [30,31]. The complexes were placed in an orthorhombic box with a size of 15 × 35 × 55 Å, filled with 0.15 M of an aqueous solution of NaCl. The solvent model was TIP3P. Counterions were added to the system to maintain neutrality. The thermodynamic ensemble was NPγT. The period of recorded dynamics simulation was 100 nanoseconds at a temperature of 310 K (37 • C). The protocol for preparing the system for simulation included a preliminary minimization and the balancing of system components.
The population analysis was performed using the VolMap plug-in implemented in the VMD program [32]. A map based on the weighted atomic density and weighted atomic population was considered. To quantify protein structure changes before and after simulation, we used the MultiSeq plug-in implemented in the VMD program. The RMSD parameter was chosen as a change descriptor. The changes were visually represented as a color gradient (Blue-White-Red).

Conclusions
A series of thirty-nine esters based on (-)-borneol were synthesized by a simple threestep procedure and evaluated for antiviral activity against RSV. The cytotoxicity of the compounds tested was studied in HEp2 cells. Two of the (-)-borneol derivatives synthesized, 3b (with methyl substitutions at the nitrogen atom) and 5a (with morpholine cycle), demonstrated the best antiviral activities. Moreover, the derivatives with piperidine (6b), 4-methylpiperidine (7b), 4-metylpiperazine (9c), and 4-ethylpiperazine (10b,c) were found to exhibit significant antiviral activity. The time-of-addition experiments demonstrated that compounds 3b, 5a, and 6b were likely to exhibit inhibitory activities against RSV infection by binding to the viral F protein and blocking the fusion of the viral and host cell membranes. Docking studies and molecular dynamics simulation on a representative target (F protein) revealed that the tested compounds fit well into RSV F hydrophobic cavity DS-Cav1. To conclude, (-)-borneol ester with an N-containing heterocycle has been identified as a promising therapeutic candidate for developing a broad-spectrum antiviral agent.