Corrigendum: Assessing potency and binding kinetics of soluble adenylyl cyclase (sAC) inhibitors to maximize therapeutic potential
Erratum: Assessing potency and binding kinetics of soluble adenylyl cyclase (sAC) inhibitors to maximize therapeutic potential
- 1Department of Pharmacology, Weill Cornell Medicine, New York, NY, United States
- 2Department of Biochemistry, University of Bayreuth, Bayreuth, Germany
- 3Tri-Institutional Therapeutics Discovery Institute, New York, NY, United States
In mammalian cells, 10 different adenylyl cyclases produce the ubiquitous second messenger, cyclic adenosine monophosphate (cAMP). Amongst these cAMP-generating enzymes, bicarbonate (HCO3−)-regulated soluble adenylyl cyclase (sAC; ADCY10) is uniquely essential in sperm for reproduction. For this reason, sAC has been proposed as a potential therapeutic target for non-hormonal contraceptives for men. Here, we describe key sAC-focused in vitro assays to identify and characterize sAC inhibitors for therapeutic use. The affinity and binding kinetics of an inhibitor can greatly influence in vivo efficacy, therefore, we developed improved assays for assessing these efficacy defining features.
Introduction
Adenylyl cyclases (ACs) are the enzymes responsible for cyclizing adenosine triphosphate (ATP) into cyclic adenosine monophosphate (cAMP), a ubiquitous second messenger that plays an essential role in a wide variety of cellular signaling pathways. For many years following its initial discovery (Sutherland and Rall, 1958), it was widely thought that cAMP was exclusively produced at the plasma membrane of mammalian cells by nine different G-protein regulated, transmembrane-domain containing ACs (tmACs; ADCY1-9). However, in 1975, seemingly in opposition to this paradigm, a G-protein insensitive cytosolic AC activity was identified in a mammalian testis lysate (Braun and Dods, 1975). This cytosolic AC activity could utilize Mn2+ ions as a cofactor for the cyclization reaction, biochemically differentiating it from other class III adenylyl cyclases (like tmACs), which primarily relied on Mg2+ ions for AC activity (Braun and Dods, 1975; Braun et al., 1977). Together, these unique features indicated that this cytosolic AC activity was distinct from classical tmACs. Over 20 years later, a 10th mammalian AC gene (ADCY10) was cloned and its protein product, soluble AC (sAC), was identified as the enzyme responsible for the cytosolic AC activity (Buck et al., 1999). sAC is now recognized as one of 10 mammalian AC isoforms.
While sAC is the only mammalian AC isoform that lacks transmembrane spanning segments (Kamenetsky et al., 2006), like the other mammalian ACs, its enzymatic activity is attributed to two heterologous catalytic domains (C1 and C2) (Buck et al., 1999; Chen et al., 2000). RNA expression profiling studies predict the existence of sAC isoforms containing only the second (C2) catalytic domain (Farrell et al., 2008; Chen et al., 2013), thus far, only C1 and C2 containing isoforms of sAC have been biochemically characterized and demonstrated to be abundantly expressed in testes and sperm (Buck et al., 1999; Jaiswal and Conti, 2001; Hess et al., 2005). Further studies are needed to confirm the expression of, and understand the contribution to sAC biology of, C2-only isoforms. The experiments in this study focus exclusively on the C1-C2 containing isoforms.
Biochemically, both Ca2+ and HCO3− directly and synergistically stimulate the activity of sAC; Ca2+ lowers the substrate Km of sAC, thus increasing its apparent affinity for ATP, while HCO3− increases the Vmax of the reaction, thereby enhancing the rate at which sAC generates product cAMP (Litvin et al., 2003; Steegborn et al., 2005; Kleinboelting et al., 2014). Through these physiological activators, sAC plays a role in various biological processes (Wiggins et al., 2018; Rossetti et al., 2021) including reproduction (Buffone et al., 2014). Prior to ejaculation, sperm stored in the epididymis are morphologically mature, yet dormant, and lack the ability to fertilize an egg (Yanagimachi, 1994). Upon ejaculation, sperm begin to swim and gain the “capacity” to fertilize as they transit through the female reproductive tract in a molecular process called capacitation (Visconti et al., 1998). Pharmacological (Hess et al., 2005; Ramos-Espiritu et al., 2016; Balbach et al., 2021) and genetic (Esposito et al., 2004; Hess et al., 2005; Xie et al., 2006; Akbari et al., 2019) evidence, in both mice and humans, indicate that sAC generated cAMP is an essential component of the signal transduction pathway that mediates sperm motility and capacitation. Due to the essential role of sAC in these processes, there is great interest in developing sAC inhibitors for use as novel, on-demand, non-hormonal male contraceptives (Balbach et al., 2020; Ferreira et al., 2022). Despite its widespread expression (Wiggins et al., 2018; Rossetti et al., 2021), there are numerous strategies to safely and effectively target sAC for contraception (recently reviewed in Ferreira et al., 2022).
To develop sAC inhibitors for therapeutic use, we employed structure-based drug design to improve the potency and properties of a well-characterized, small molecule sAC inhibitor, LRE1 (Fushimi et al., 2021). LRE1 was identified in a high-throughput screen (Ramos-Espiritu et al., 2016). It is an allosteric inhibitor that binds to sAC with micromolar potency (IC50 = 3.2 μM, Table 1) at the dimer interface of C1 and C2 and occupies the same site where the activator HCO3− binds (Kleinboelting et al., 2014). To assist in our drug discovery efforts, we developed multiple sAC-focused assays to thoroughly characterize sAC inhibitors. These assays are described in detail below. Assays for assessing inhibitor selectivity for sAC over tmACs (Ramos-Espiritu et al., 2016; Balbach et al., 2021) and inhibitor potency in physiologically relevant systems (Ramos-Espiritu et al., 2016; Balbach et al., 2021) were included in our drug discovery program, but they are not discussed as this paper focuses specifically on the in vitro characterization of inhibitor potency and binding kinetics.
Materials and methods
Materials
Human sAC protein used in all assays was purified from Sf9 cells via baculovirus expression. Biochemical potency and jump dilution assays used N-terminal tagged GST-sACt (Litvin et al., 2003); SPR used highly pure C-terminal tagged sACt-His6 that was prepared via sequential chromatography that included, in this order, affinity (Ni-Sepharose HP from GE Healthcare), ion exchange (Mono Q from GE Healthcare), and size exclusion (Superdex 200 from GE Healthcare) chromatographies. The synthesis of TDI-10229 is described in Fushimi et al., 2021 and the synthesis of other compounds is described in Miller et al., 2022.
In vitro adenylyl cyclase activity assay with purified sAC protein (biochemical potency assay)
All in vitro adenylyl cyclase activity assays performed with purified human sAC protein utilized the classical “two-column” method, developed by Salomon (Salomon, 1979). In brief, conversion of α-32P labeled ATP into 32P labeled cAMP is quantitated following sequential Dowex and Alumina chromatography to purify the generated 32P labeled cAMP. Two variations of this assay were employed: the “standard” assay and the “subnanomolar” assay. In the standard assay, ∼5 nM of sACt protein was incubated in assay buffer containing 4 mM MgCl2, 2 mM CaCl2, 1 mM ATP, and 40 mM NaHCO3. The reaction was initiated by addition of ∼1,000,000 counts per minute of α-32P labeled ATP. In the subnanomolar assay, ∼0.25 nM of sACt protein was incubated in buffer containing 10 mM MnCl2 and 2 mM ATP, and the reaction was initiated by ∼3,000,000 counts per minute of α-32P labeled ATP Thus, this assay included three-fold higher specific activity of the substrate ATP. For both assay types, the buffer also contained 50 mM Tris-HCl pH 7.5, 3 mM DTT, and 0.03% BSA. Reactions were performed using a final volume of 100 μL at 30°C, and cAMP was quantitated after 30 min. For concentration-response curves, purified sACt protein was preincubated with the indicated inhibitors or vehicle (1% v/v DMSO) for 15 min. GraphPad Prism (www.graphpad.com, version 9.3.1) was used for curve fitting of the concentration-response data and determination of IC50 values.
Cellular adenylyl cyclase activity assay in sAC overexpressing (4–4) cells (cellular potency assay)
To measure sAC-dependent cAMP accumulation in cells, we used human “4-4 cells”, which are HEK293 cells that stably overexpress sACt (Zippin et al., 2013). For the assay, 1 × 105 4-4 cells were seeded in each well of a 24-well plate and incubated for 24 h in DMEM+10% FBS media at 37°C/5% CO2. One hour before the assay the media was aspirated and replaced with 300 μL fresh media. Cells were preincubated with sAC inhibitor at the indicated concentrations or vehicle (0.7% DMSO) for 10 min at 37°C. Intracellular cAMP accumulation was initiated by addition of 500 μM IBMX (Sigma-Aldrich), and cAMP was allowed to accumulate for 5-min. Intracellular cAMP was quantitated using the Direct cAMP Elisa kit (Enzo Life Sciences) following manufacturer’s instructions. GraphPad Prism (www.graphpad.com, version 9.3.1) was used for curve fitting of the concentration-response data and determination of IC50 values.
Surface plasmon resonance with purified sAC protein
On-rate (kon), off-rate (koff), and absolute affinity (Ki) values for inhibitor binding to sAC were obtained at 25 or 37°C with a Biacore 8 K instrument (Cytiva) using a single cycle kinetics protocol. In PBS-P+ buffer (1 mM KH2PO4, 150 mM NaCl, 6 mM Na2HPO4, 0.05% (w/v) P20 Surfactant), 50 μg/ml of recombinant purified His-tagged sACt protein was covalently immobilized on a Series S Sensor NTA chip (Cytiva) using Ni2+-His tag chelation followed by amine coupling with a 1:1 mixture of 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide and N-hydroxysuccinimide. After coupling, 1 M ethanolamine (Cytiva) followed by 350 mM EDTA (Cytiva) were, respectively, used to block any remaining reactive groups on the surface of the chip and to strip the Ni2+. Following chip preparation, TBS-P+ running buffer (50 mM Tris-HCl pH 7.5, 150 mM NaCl, and 0.05% (w/v) P20 Surfactant supplemented with 1% DMSO) was flowed over the surface of the chip. After a stable baseline was obtained, five different concentrations of the indicated inhibitors dissolved in TBS-P+ running buffer with a final DMSO concentration of 1%, were sequentially injected into a single channel for 120 s at a flowrate of 50 μl/min, followed by 60 s of TBS-P+ running buffer +1% DMSO alone. Concentrations used for TDI-10229 were 5,000, 1,250, 312.5, 78.1, and 19.5 nM. For all other inhibitors, the concentrations used were: 1,000, 250, 62.5, 15.6 and 4 nM. Subsequent to the highest concentration, compounds were allowed to dissociate for 600 s in the presence of TBS-P+ running buffer containing 1% DMSO. All experiments were performed in parallel in otherwise identically prepared reference channels lacking immobilized protein. To process the collected data, responses from the reference channels were subtracted from the responses from the active channels. From the reference-subtracted data, fitted curves, kon, koff, and Ki values were determined with the Biacore 8 K Insight Evaluation Software Version 2.0 (Cytiva) using a 1:1 binding kinetics model.
In vitro jump dilution assay with purified sAC protein
Jump dilution adenylyl cyclase assays were performed using a modified version of the subnanomolar in vitro adenylyl cyclase activity assay described above. All assays were performed at 30°C, in the presence of 0.03% BSA, and each had a final DMSO concentration of 0.01%. At the start of each assay, ∼25 nM recombinant purified human sACt protein was preincubated with the indicated inhibitors for 15 min. Each inhibitor was used at an initial concentration 10-fold above their IC50 values. Following the preincubation period, 1 μl of the enzyme-inhibitor solution was diluted by addition of 99 μl of a reaction solution containing: 2 mM ATP, 10 mM Mn2+, 50 mM Tris pH 7.5, 3 mM DTT, 0.03% BSA and ∼4,000,000 counts per minute of α-32P labeled ATP. After this 100-fold dilution step, each inhibitor was present at a concentration 10-fold below their IC50 values to minimize inhibitor rebinding during the reaction. Parallel reactions were stopped every 6 min over the course of 1 h and generated cAMP was measured using the “two-column” method referenced above. Data were fit to the equation: % Total cAMP Formed = vst+((v0-vs) (1-e(−kobst)))/kobs, with kobs being an estimate of the dissociation rate constant (koff) (Kumar and Lowery, 2017) using GraphPad Prism (www.graphpad.com, version 9.3.1). To calculate vs (uninhibited enzyme velocity) and v0 (inhibited enzyme velocity), the jump dilution assay was performed in the presence of only DMSO or an excess concentration of inhibitor, respectively.
Results
Assessing the biochemical and cellular potency of sAC inhibitors
The lead optimization phase of drug development is initially focused on efforts to improve the biochemical potency (IC50 for an inhibitor) of a lead compound via successive rounds of medicinal chemistry. Thus, a successful drug discovery program requires a robust potency assay that can efficiently and accurately determine the IC50 values of newly synthesized compounds. For our program, we utilized a highly reproducible, in vitro assay to determine the biochemical potencies of our sAC inhibitors on purified recombinant sAC protein. This in vitro biochemical potency assay utilized a radioactivity-based method to measure AC activity, which was pioneered by Yarom Salomon (Salomon, 1979). This method requires the use of α-32P labelled ATP that, in the presence of an active AC, is converted to 32P labelled cAMP. The radioactive cAMP reaction product is purified and separated from the unconverted radioactive ATP substrate using sequential Dowex and Alumina chromatography (Salomon et al., 1974). We routinely use this method to reliably measure the in vitro activity of several different ACs, including sAC (Levin et al., 1992; Chen et al., 2000; Litvin et al., 2003; Ramos-Espiritu et al., 2016).
For our “standard” biochemical potency assay, the AC activity of ∼5 nM purified recombinant sAC protein is measured in a reaction that includes the physiologically relevant activators Ca2+ and HCO3− (Litvin et al., 2003). For this assay, substrate ATP is kept at its Km of 1 mM (Litvin et al., 2003), which closely mimics intracellular concentrations of ATP, to identify inhibitors that are either competitive, uncompetitive, or noncompetitive with substrate. Using our standard assay, we generated concentration-response curves for sAC inhibitors to calculate an accurate IC50 value via non-linear curve fitting (Figure 1A). Driven primarily by this standard assay to assess potency, we described the development and characterization of two sAC inhibitors, TDI-10229 (IC50 = 158.6 nM) (Balbach et al., 2021; Fushimi et al., 2021) and TDI-11861 (IC50 ≤ 2.5 nM) (Table 1) (Miller et al., 2022), that are useful for in vivo interrogation of sAC’s therapeutic potential.
FIGURE 1. Measuring Biochemical IC50 Values with In Vitro sAC Activity Assays. Concentration-response curves of indicated inhibitors on purified recombinant human soluble adenylyl cyclase (sAC) protein. Data was collected at 30°C in the presence of either (A) 1 mM ATP, 2 mM Ca2+, 4 mM Mg2+, 40 mM HCO3− and ∼5 nM sAC protein or (B) 2 mM ATP, 10 mM Mn2+, and ∼0.25 nM sAC protein. Data is normalized to respective DMSO-treated controls and shown as mean ± SEM (n ≥ 3) (C) An overlay of the fitted IC50 curves from the standard (dashed lines) and subnanomolar (solid lines) assay. Colored arrow indicates the shift in IC50 for TDI-11891.
In addition to TDI-11861, our medicinal chemistry efforts yielded numerous sAC inhibitors with biochemical potencies in the low nanomolar range (Miller et al., 2022) (Figure 1A; Table 1). When determining in vitro IC50 values, IC50 values of less than half the active target protein concentration cannot be measured due to the tight binding phenomenon (Williams and Morrison, 1979; Cook and Cleland, 2007). Since the concentration of sAC protein in our standard assay is ∼5 nM, the assay reached its theoretical limit and is not suitable for assessing inhibitor potencies of TDI-11861, TDI-11891, and TDI-11893, whose calculated IC50 values fell at or below 2.5 nM. To increase the potency range of our biochemical assay and gain the ability to measure low pM IC50 values, we tested lower concentrations of sAC protein in our cyclase assay (Supplementary Figure S1). We found 0.25 nM sAC produced sufficient cAMP over 30 min to assess inhibitor potencies. Thus, our new conditions decreased the amount of sAC protein 20-fold, from ∼5 to ∼0.25 nM. Because Mn2+-dependent in vitro sAC activity is ∼20-fold higher than Mg2+-dependent in vitro sAC activity (Litvin et al., 2003), we measured sAC activity in the presence of Mn2+ as the sole divalent cation. sAC activity was first detected in the presence of Mn2+ (Braun and Dods, 1975), and Mn2+ remains the most potent in vitro stimulator of sAC activity. TDI-10229 (Fushimi et al., 2021) and subsequent inhibitors (Miller et al., 2022) extend into the ATP binding site, so they inhibit both the physiologically-stimulated and Mn2+ stimulated activities. In parallel, we tripled the specific activity of the substrate ATP (i.e., increasing the radioactive ATP added to the assay from ∼1,000,000 to ∼3,000,000 counts per minute per assay). We refer to these assay conditions, with 20-fold reduced sAC protein levels, as a “subnanomolar” sAC activity assay.
We re-assessed potencies of several sAC inhibitors using this newly developed subnanomolar assay (Figure 1B). As expected, the inhibitors with IC50 values above 2.5 nM, TDI-10229 and TDI-11155, whose potencies are within the measurable range of the standard assay, had similar IC50 values in both assays (Figure 1C; Table 1). For the three sAC inhibitors with IC50 values at or near the theoretical lower limit of the standard assay (i.e., TDI-11861, TDI-11891, and TDI-11893), the subnanomolar assay revealed enhanced potencies. While TDI-11861 and TDI-11893 each exhibited IC50 values of 1.7 nM, TDI-11891, which is the enantiomer of TDI-11893 (Table 1), was ∼7.5-fold more potent with an IC50 of 0.33 nM in the subnanomolar assay (Figure 1C; Table 1). Thus, while the standard assay remains relevant for assessing potencies in the presence of the physiological activators Ca2+ and HCO3−, the newly developed assay conditions are suitable for distinguishing sAC inhibitors with subnanomolar potencies.
While the standard and subnanomolar assays are useful tools for determining the potency of our inhibitors against purified sAC protein (i.e., biochemical potency), we are ultimately interested in inhibitor efficacy in a cellular context. In addition to how tightly an inhibitor binds its target, efficacy against an intracellular enzyme target is also determined by factors that are present only in the context of cells, such as membrane permeability, non-specific protein binding and compound stability in media and inside cells. To determine the cellular potency of our inhibitors, we developed a cell-based assay using engineered HEK293 cells that stably overexpress sAC, referred to as 4-4 cells (Zippin et al., 2013). sAC is the predominant source of cAMP inside 4-4 cells, and in the presence of a pan-selective phosphodiesterase (PDE) inhibitor (i.e., IBMX) that prevents cAMP degradation, cAMP accumulates within 4-4 cells in a sAC-dependent manner (Bitterman et al., 2013). Measuring cAMP accumulation in the presence of increasing doses of sAC inhibitors reveals their cellular potency. sAC inhibitors such as LRE1 (IC50 = 5.3 μM), TDI-10229 (IC50 = 113.5 nM), TDI-11155 (IC50 = 15.7 nM), and TDI-11861 (IC50 = 5.1 nM) had similar biochemical and cellular IC50 values while other sAC inhibitors, such as TDI-11893 (IC50 = 19.4 nM) and TDI-11891 (IC50 = 2.3 nM), had shifted cellular IC50 values compared to their biochemical IC50 values (Figure 2; Table 1), indicating cell-context specific influences on their potencies.
FIGURE 2. Measuring Cellular IC50 Values with a Cellular sAC Activity Assay. Concentration-response curves of indicated inhibitors in sAC-overexpressing 4-4 cells. Cells were preincubated with sAC inhibitors then treated with 500 μM 3-isobutyl-1-methylxanthine (IBMX) for 5 min. The amount of accumulated cAMP at 37°C was then measured. Data is normalized to respective DMSO-treated controls and shown as mean ± SEM (n ≥ 3).
Using surface plasmon resonance to measure affinities and binding kinetics of sAC inhibitors
A major goal of lead optimization is to obtain highly potent compounds. As detailed above, we achieved this goal by using sAC-focused in vitro potency assays to identify sAC inhibitors that had picomolar biochemical potencies and favorable low nanomolar cellular potencies. We were also extremely interested in inhibitor binding kinetics, especially the inhibitor dissociation rates. As an example, for sAC inhibitors to exhibit maximal efficacy as male contraceptives, compounds which bind sAC in sperm in the man will have to remain engaged with their target post-ejaculation, after the sperm enter the inhibitor-free female reproductive tract.
The binding kinetics of an inhibitor involves the various rates that define the interactions between the inhibitor and its target protein; those being the bimolecular rate at which the inhibitor binds to its target (i.e., kon, association constant, or on-rate) and the rate at which the target-inhibitor complex dissociates (i.e., koff, dissociation constant, or off-rate). The on-rate and off-rate of an inhibitor define its absolute biochemical affinity, as Ki = koff/kon. Consequently, a high intrinsic affinity is an indication of either a slow off-rate, a fast on-rate, or both. Another useful metric that can be evaluated and enhanced during drug development is the residence time (τ), which is defined as the reciprocal of the off-rate (1/koff) and is a estimation of the total lifetime of the inhibitor-target complex (Copeland, 2016; Bernetti et al., 2019). It is well precedented that a significant percentage of drug candidates have strong preclinical data but ultimately still fail in clinical trials due to a lack of efficacy (Kola and Landis, 2004). Recent evidence from a variety of drug discovery efforts indicates that residence time can be a major determinant of in vivo efficacy, defining it as an important inhibitor property to be optimized during lead optimization (Copeland et al., 2006; Swinney, 2008; Lu and Tonge, 2010; Schuetz et al., 2017).
To incorporate residence time assessment into our drug discovery program and define the absolute affinity of the inhibitors, we employed surface plasmon resonance (SPR), a biophysical technique that is widely utilized to measure both the on-rates and off-rates of biomolecular interactions (Zhang et al., 2016). For SPR, we sequentially flowed increasing concentrations of each sAC inhibitor over a sensor chip containing covalently immobilized sAC protein. Inhibitors were applied for 120 s to determine the on-rate (kon), followed by 60 s of buffer. After the final, highest concentration was bound, the sAC inhibitor was allowed to dissociate for 600 s to measure the off-rate (koff). The on-rate, off-rate, Ki, and residence time values were calculated by evaluation software using a 1:1 binding kinetics model. Residence times were found to increase with the biochemical potency of the inhibitors (i.e., the more potent sAC inhibitors had longer residence times) (Figures 3A–E; Table 2). For example, at room temperature (25°C) a sAC inhibitor with mid-nanomolar biochemical potency (TDI-10229, IC50 = 158.6 nM) had the shortest residence time (25 s) of the compounds tested (Figure 3A), while a sAC inhibitor with a mid-picomolar biochemical potency (TDI-11891, biochemical IC50 = 0.33 nM) had a significantly longer residence time of ∼3,000 s (Figure 3E). Among the sAC inhibitors tested, the on-rates, which are lower than diffusion controlled, did not change with respect to biochemical potency; the assessed sAC inhibitors had on-rates of ∼200,000 (Ms)−1 regardless of their biochemical potency (Table 2). Thus, the observed correlation between measured absolute affinity (Ki) (Table 2) and biochemical potencies (Table 1) reflects differences in inhibitor residence times.
FIGURE 3. Assessing sAC Inhibitor Binding Kinetics with Surface Plasmon Resonance. Representative single-cycle SPR sensorgrams of (A,F) TDI-10229, (B,G) TDI-11155, (C,H) TDI-11861 (D,I) TDI-11893, and (E,J) TDI-11891 associating to and dissociating from immobilized sAC protein. Traces (dotted lines) are shown along with best fits generated via a 1:1 binding model (colored lines). Experiments (n ≥ 4) were conducted at either (A–E) 25°C or (F–J) 37°C.
These SPR experiments were performed at room temperature (25°C) which facilitates measuring off-rates of compounds with shorter residence times. For compounds with long residence times, we repeated the SPR experiments at the physiologically relevant temperature of 37°C. As expected, the elevated temperature speeds up the molecular interactions occurring within the system, leading to accelerated on-rates and shortened residence times. Increasing the temperature significantly reduced sAC inhibitor residence times by 3- to 6.5-fold, increased on-rates by 1.6- to 3-fold, and increased Ki values by approximately 2-fold (Figures 3F–J; Table 2). However, at 37°C, the measured residence times and absolute affinities still positively correlated with biochemical potency (Table 1) (i.e., the most potent sAC inhibitors exhibited the longest residence times and the strongest absolute affinities).
An in vitro jump dilution assay to measure sAC inhibitor residence times
SPR is a valuable technique to directly compare the binding kinetics and absolute affinities of sAC inhibitors, but it is a purely biophysical technique that studies the inhibitor/target interaction with the enzyme immobilized. Therefore, we also developed an in vitro “jump dilution” adenylyl cyclase assay. Jump dilution assays are an experimental format frequently used as an alternative/complement to SPR to measure the residence time of small molecules (Walkup et al., 2015; Zhang et al., 2016). Our in vitro jump dilution cyclase assay was a modified version of the subnanomolar biochemical potency assay. Each sAC inhibitor was preincubated with sAC protein at an inhibitor concentration 10-fold above its biochemical IC50 and a sAC concentration 100-fold above its level in the subnanomolar sAC activity assay. After a set preincubation time, the protein-inhibitor mixture underwent a 100-fold “jump dilution,” such that the final inhibitor concentration is 10-fold below its biochemical IC50 (Figure 4A). At this minimal inhibitory concentration, sAC activity recovered as the inhibitor dissociated from the protein over the course of the 60-min experiment (Figure 4B). The off-rate of each inhibitor is the rate at which sAC activity recovered to uninhibited levels; the reciprocal of this rate is the inhibitor residence time (Kumar and Lowery, 2017).
FIGURE 4. In Vitro Jump Dilution Assay for Determining sAC Inhibitor Residence Times. (A) Schematic diagram of the jump dilution assay (figure adapted from BellBrook Labs “A Guide to Measuring Drug Target Residence Times with Biochemical Assays”) (B) In vitro jump dilution curves of indicated inhibitors. All assays were done at 30°C in the presence of 2 mM ATP, 10 mM Mn2+ and ∼0.25 nM of purified recombinant human sAC protein. Following a 100-fold dilution, sAC cyclase activity was measured every 6 min for 60 min. Data is normalized to respective DMSO-treated controls and is shown as mean ± SEM (n ≥ 4).
As expected, the rank order of residence times determined with the newly developed in vitro jump dilution assay correlated well with the order of compound residence times measured via SPR (Table 2). However, the exact length of the residence time differed in the two assays, and this difference, either an increase or decrease in residence time, varied from inhibitor to inhibitor. For example, TDI-11861 had a 2-fold shorter residence time in the jump dilution assay (∼2,200 s vs. ∼1,000 s) while TDI-11891 had a 2-fold longer residence time in the jump dilution assay (∼3,100 s vs. ∼5,700 s) (Table 2). The maximum length of our in vitro jump dilution assay is 60 min (3,600 s), and thus the overall length of the assay was shorter than the measured residence time of TDI-11891. Since the length of a jump dilution assay should exceed the longest measured residence time (Zhang et al., 2016), we increased the jump dilution assay runtime to 180 min for TDI-11891, and measured sAC activity every 18 min (Supplementary Figure S2A). In this extended jump dilution assay, the measured residence time of TDI-11891 increased only slightly to ∼6,500 s (Supplementary Figure S2B).
Discussion
sAC generated cAMP-signaling has been linked to various physiological processes, including lysosomal acidification (Rahman et al., 2016), regulation of intraocular eye pressure (Lee et al., 2011; Gandhi et al., 2017), and the motility and capacitation of sperm (Hess et al., 2005; Xie et al., 2006; Buffone et al., 2014). Due to these varied biological functions, sAC inhibition is a novel approach with several potential therapeutic applications (Wiggins et al., 2018; Ferreira et al., 2022). To develop sAC inhibitors for therapeutic use, we initiated drug discovery efforts that sought to bring a compound identified in a high-throughput screen (LRE1) through the lead optimization phase to ultimately produce sAC inhibitors with “drug-like” properties. In this paper, we detail several sAC-focused in vitro assays that we used to thoroughly characterize efficacy defining properties of sAC inhibitors and, as a result, establish vital tools in sAC-directed drug discovery research.
For in vitro potency assessment, we use two assays; our previously described standard assay, which measures potency in a reaction containing physiologically relevant sAC activators (Balbach et al., 2021; Fushimi et al., 2021; Miller et al., 2022), and the newly developed subnanomolar assay, which is used to determine the biochemical potency of our highest affinity sAC inhibitors. The subnanomolar assay overcomes the protein concentration-mediated limitations of the standard assay, affording us the capability to reliably determine IC50 values for sAC inhibitors whose biochemical potencies are in the high-to mid-picomolar range. The new dynamic range achieved with the subnanomolar assay allows us to continue our structure-guided medicinal chemistry efforts to further improve the biochemical potencies of sAC inhibitors. However, the subnanomolar assay does have limitations because it includes Mn2+ as the sole divalent cation in the assay and does not include allosteric modulators, such as HCO3−. The compounds used in this study, in addition to allosterically inhibiting sAC via binding in the bicarbonate binding site (Ramos-Espiritu et al., 2016), interfere with ATP binding (Fushimi et al., 2021; Miller et al., 2022). To complement and expand the standard and subnanomolar sAC activity assays, we have also incorporated a cellular potency assay into our drug discovery program. We use this assay to identify sAC inhibitors that function in a cellular context, eliminating compounds with poor cellular target engagement from further development.
Our lead optimization program also utilizes two separate assays, SPR and in vitro jump dilution, to assess the binding kinetics (i.e., on-rate and/or off-rate) of our sAC inhibitors. We conduct our SPR experiments at two different temperatures, 25°C and the more physiologically relevant 37°C. As expected, sAC inhibitor residence times are much shorter at 37°C, providing a more realistic evaluation of functional residence times during in vivo experiments. For this chemical series of sAC inhibitors, the prime determinant for absolute affinity and biochemical potency appears to be residence time, which generally increased concurrently with both properties; on-rates of our sAC inhibitors did not change in relation to potency or affinity. To complement our SPR experiments, which were done in absence of any AC reaction cofactors and/or substrate on immobilized protein, we developed an in vitro jump dilution assay to determine if sAC inhibitor residence times were altered in the presence of an active, cAMP-producing enzyme. The residence times measured via in vitro jump dilution assays are consistent with the results from SPR experiments; the rank order of sAC inhibitor residence times is similar between the two types of assays. For some sAC inhibitors, we observe up to 2-fold differences between in vitro jump dilution residence times and SPR residence times, which could be due to different assay temperatures, inter-assay variability or altered inhibitor binding due to the presence of ATP/Mn2+ and the allosteric changes that occur during the enzymatic reaction. Finally, we increased the length of the in vitro jump dilution assay to reliably determine the residence time for the tightest binding sAC inhibitor, TDI-11891, whose residence time is longer than 3,600 s. Although lengthening the assay did not reveal an appreciable difference for TDI-11891, the extended jump dilution assay allows for a reliable assessment of inhibitors with significantly longer residence times. This may be required if, in the future, our drug discovery program generates sAC inhibitors with residence times even greater than TDI-11891.
The assays described in this paper are important research tools that can be utilized to characterize sAC inhibitors. They comprise the initial stages of a workflow designed to produce sAC inhibitors with optimized therapeutic properties (i.e., high potency and long residence times). The standard assay is used for initial assessment of the biochemical potency of a sAC inhibitor; if the measured IC50 value is less than or equal to the theoretical limit of the standard assay (∼2.5 nM) then it is followed by the subnanomolar assay to determine whether the true IC50 is lower than originally measured. Once sAC inhibitors with the desired in vitro potencies are identified, the cellular potency (4-4 cell) assay is used to exclude sAC inhibitors with cell-specific limitations from further assessment and development. In parallel, residence times and absolute affinity are determined via SPR at room temperature, and if sufficiently tight binding, at 37°C, to best mimic binding kinetics that will occur in vivo. The jump dilution assay is then used as an independent check to confirm the measured residence times and to determine whether the residence time may be altered due to allosteric changes during the enzymatic reaction. For inhibitors whose measured residence time exceeds 3,600 s, the extended jump dilution will reliably determine a significantly long residence time. At this point in the process, promising sAC inhibitors should be assessed for sAC-selectivity. The mammalian enzymes most evolutionarily and biochemically similar to sAC are the family of 9 transmembrane adenylyl cyclases (tmACs); thus, we examine inhibitory activity against several tmACs (Bitterman et al., 2013; Ramos-Espiritu et al., 2016; Fushimi et al., 2021). Potent and selective inhibitor candidates are then subjected to a pharmacodynamic efficacy assay. As mentioned previously, sAC is essential for sperm capacitation (Esposito et al., 2004; Hess et al., 2014; Akbari et al., 2019; Balbach et al., 2021), so for this application we measure concentration dependence in several sperm capacitation-specific assays (Ramos-Espiritu et al., 2016; Balbach et al., 2021). These assays are directly applicable if sAC inhibitors are being developed as contraceptives. If sAC inhibitors are being developed for another therapeutic use, then other function-specific assays need to be conducted to determine inhibitor potency in physiologically relevant systems. Ultimately, when the ideal candidate sAC inhibitor is identified, the in vivo efficacy of the inhibitor must be determined in an appropriate in vivo experiment to provide proof-of-concept evidence for the therapeutic benefit of sAC inhibition.
In summary, what we report here is a workflow designed for vetting sAC inhibitors for drug discovery programs where binding kinetics, specifically residence times, and affinity are efficacy defining features.
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Author contributions
JB, LL, RM, CS, and PM conceived and designed research; TR, JF, LG, and HB performed experiments; TR prepared figures; JB, LL, and TR wrote the paper. All authors approved final version of manuscript.
Funding
This work was supported by the NIH; NIA F31 AG069501 (TR), NICHD F31 HD105363 (JF), NICHD P50 HD100549 (JB and LL) and by generous funding from the Male Contraceptive Initiative (MCI) (to JB and LL).
Acknowledgments
We thank Michael B. Murphy at Cytiva as well as J. Fraser Glickman, Lavoisier Ramos-Espiritu, and other staff at Rockefeller University’s Fisher Drug Discovery Resource Center for technical assistance with the Biacore 8 K instrument and SPR experiments. We gratefully acknowledge the support provided by the Tri-Institutional Therapeutics Discovery Institute (TDI), a 501(c) (3) organization. TDI receives financial support from Takeda Pharmaceutical Company, TDI’s parent institutes (Memorial Sloan Kettering Cancer Center, The Rockefeller University and Weill Cornell Medicine) and from a generous contribution from Mr. Lewis Sanders and other philanthropic sources.
Conflict of interest
JB and LL have licensed commercialization of a panel of monoclonal antibodies directed against sAC. CS, JB, LL, and employees of the Tri-Institutional Therapeutics Discovery Institute (TDI) may benefit from the further development and licensure of molecules described herein.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2022.1013845/full#supplementary-material
References
Akbari, A., Pipitone, G. B., Anvar, Z., Jaafarinia, M., Ferrari, M., Carrera, P., et al. (2019). ADCY10 frameshift variant leading to severe recessive asthenozoospermia and segregating with absorptive hypercalciuria. Hum. Reprod. 34, 1155–1164. doi:10.1093/humrep/dez048
Balbach, M., Fushimi, M., Huggins, D. J., Steegborn, C., Meinke, P. T., Levin, L. R., et al. (2020). Optimization of lead compounds into on-demand, nonhormonal contraceptives: Leveraging a public-private drug discovery institute collaboration. Biol. Reprod. 103, 176–182. doi:10.1093/biolre/ioaa052
Balbach, M., Ghanem, L., Rossetti, T., Kaur, N., Ritagliati, C., Ferreira, J., et al. (2021). Soluble adenylyl cyclase inhibition prevents human sperm functions essential for fertilization. Mol. Hum. Reprod. 27, gaab054. doi:10.1093/molehr/gaab054
Bernetti, M., Masetti, M., Rocchia, W., and Cavalli, A. (2019). Kinetics of drug binding and residence time. Annu. Rev. Phys. Chem. 70, 143–171. doi:10.1146/annurev-physchem-042018-052340
Bitterman, J. L., Ramos-Espiritu, L., Diaz, A., Levin, L. R., and Buck, J. (2013). Pharmacological distinction between soluble and transmembrane adenylyl cyclases. J. Pharmacol. Exp. Ther. 347, 589–598. doi:10.1124/jpet.113.208496
Braun, T., and Dods, R. F. (1975). Development of a Mn-2+-sensitive, “soluble” adenylate cyclase in rat testis. Proc. Natl. Acad. Sci. U. S. A. 72, 1097–1101. doi:10.1073/pnas.72.3.1097
Braun, T., Frank, H., Dods, R., and Sepsenwol, S. (1977). Mn2+-sensitive, soluble adenylate cyclase in rat testis Differentiation from other testicular nucleotide cyclases. Biochim. Biophys. Acta 481, 227–235. doi:10.1016/0005-2744(77)90155-3
Buck, J., Sinclair, M. L., Schapal, L., Cann, M. J., and Levin, L. R. (1999). Cytosolic adenylyl cyclase defines a unique signaling molecule in mammals. Proc. Natl. Acad. Sci. U. S. A. 96, 79–84. doi:10.1073/pnas.96.1.79
Buffone, M. G., Wertheimer, E. V., Visconti, P. E., and Krapf, D. (2014). Central role of soluble adenylyl cyclase and cAMP in sperm physiology. Biochim. Biophys. Acta 1842, 2610–2620. doi:10.1016/j.bbadis.2014.07.013
Chen, J., Martinez, J., Milner, T. A., Buck, J., and Levin, L. R. (2013). Neuronal expression of soluble adenylyl cyclase in the mammalian brain. Brain Res. 1518, 1–8. doi:10.1016/j.brainres.2013.04.027
Chen, Y., Cann, M. J., Litvin, T. N., Iourgenko, V., Sinclair, M. L., Levin, L. R., et al. (2000). Soluble adenylyl cyclase as an evolutionarily conserved bicarbonate sensor. Science 289, 625–628. doi:10.1126/science.289.5479.625
Cook, P. F., and Cleland, W. W. (2007). Enzyme kinetics and mechanism. Garl. Sci. 2007, 416. doi:10.4324/9780203833575
Copeland, R. A., Pompliano, D. L., and Meek, T. D. (2006). Drug-target residence time and its implications for lead optimization. Nat. Rev. Drug Discov. 5, 730–739. doi:10.1038/nrd2082
Copeland, R. A. (2016). The drug-target residence time model: A 10-year retrospective. Nat. Rev. Drug Discov. 15, 87–95. doi:10.1038/nrd.2015.18
Esposito, G., Jaiswal, B. S., Xie, F., Krajnc-Franken, M. A. M., Robben, T. J. A. A., Strik, A. M., et al. (2004). Mice deficient for soluble adenylyl cyclase are infertile because of a severe sperm-motility defect. Proc. Natl. Acad. Sci. U. S. A. 101, 2993–2998. doi:10.1073/pnas.0400050101
Farrell, J., Ramos, L., Tresguerres, M., Kamenetsky, M., Levin, L. R., and Buck, J. (2008). Somatic “soluble” adenylyl cyclase isoforms are unaffected in Sacy tm1Lex/Sacy tm1Lex “knockout” mice. PLoS One 3, e3251. doi:10.1371/journal.pone.0003251
Ferreira, J., Levin, L. R., and Buck, J. (2022). Strategies to safely target widely expressed soluble adenylyl cyclase for contraception. Front. Pharmacol. 13, 953903. doi:10.3389/fphar.2022.953903
Fushimi, M., Buck, H., Balbach, M., Gorovyy, A., Ferreira, J., Rossetti, T., et al. (2021). Discovery of TDI-10229: A potent and orally bioavailable inhibitor of soluble adenylyl cyclase (sAC, ADCY10). ACS Med. Chem. Lett. 12, 1283–1287. doi:10.1021/acsmedchemlett.1c00273
Gandhi, J. K., Roy Chowdhury, U., Manzar, Z., Buck, J., Levin, L. R., Fautsch, M. P., et al. (2017). Differential intraocular pressure measurements by tonometry and Direct cannulation after treatment with soluble adenylyl cyclase inhibitors. J. Ocul. Pharmacol. Ther. 33, 574–581. doi:10.1089/jop.2017.0027
Hess, K. C., Jones, B. H., Marquez, B., Chen, Y., Ord, T. S., Kamenetsky, M., et al. (2005). The “soluble” adenylyl cyclase in sperm mediates multiple signaling events required for fertilization. Dev. Cell 9, 249–259. doi:10.1016/j.devcel.2005.06.007
Hess, K. C., Liu, J., Manfredi, G., Mühlschlegel, F. A., Buck, J., Levin, L. R., et al. (2014). A mitochondrial CO2-adenylyl cyclase-cAMP signalosome controls yeast normoxic cytochrome c oxidase activity. FASEB J. 28, 4369–4380. doi:10.1096/fj.14-252890
Jaiswal, B. S., and Conti, M. (2001). Identification and functional analysis of splice variants of the germ cell soluble adenylyl cyclase. J. Biol. Chem. 276, 31698–31708. doi:10.1074/jbc.M011698200
Kamenetsky, M., Middelhaufe, S., Bank, E. M., Levin, L. R., Buck, J., and Steegborn, C. (2006). Molecular details of cAMP generation in mammalian cells: A tale of two systems. J. Mol. Biol. 362, 623–639. doi:10.1016/j.jmb.2006.07.045
Kleinboelting, S., Diaz, A., Moniot, S., van den Heuvel, J., Weyand, M., Levin, L. R., et al. (2014). Crystal structures of human soluble adenylyl cyclase reveal mechanisms of catalysis and of its activation through bicarbonate. Proc. Natl. Acad. Sci. U. S. A. 111, 3727–3732. doi:10.1073/pnas.1322778111
Kola, I., and Landis, J. (2004). Can the pharmaceutical industry reduce attrition rates? Nat. Rev. Drug Discov. 3, 711–715. doi:10.1038/nrd1470
Kumar, M., and Lowery, R. G. (2017). A high-throughput method for measuring drug residence time using the transcreener ADP assay. SLAS Discov. 22, 915–922. doi:10.1177/2472555217695080
Lee, Y. S., Tresguerres, M., Hess, K., Marmorstein, L. Y., Levin, L. R., Buck, J., et al. (2011). Regulation of anterior chamber drainage by bicarbonate-sensitive soluble adenylyl cyclase in the ciliary body. J. Biol. Chem. 286, 41353–41358. doi:10.1074/jbc.M111.284679
Levin, L. R., Han, P. L., Hwang, P. M., Feinstein, P. G., Davis, R. L., and Reed, R. R. (1992). The Drosophila learning and memory gene rutabaga encodes a Ca2+/Calmodulin-responsive adenylyl cyclase. Cell 68, 479–489. doi:10.1016/0092-8674(92)90185-f
Litvin, T. N., Kamenetsky, M., Zarifyan, A., Buck, J., and Levin, L. R. (2003). Kinetic properties of “soluble” adenylyl cyclase. Synergism between calcium and bicarbonate. J. Biol. Chem. 278, 15922–15926. doi:10.1074/jbc.M212475200
Lu, H., and Tonge, P. J. (2010). Drug-target residence time: Critical information for lead optimization. Curr. Opin. Chem. Biol. 14, 467–474. doi:10.1016/j.cbpa.2010.06.176
Miller, M., Rossetti, T., Ferreira, J., Ghanem, L., Balbach, M., Kaur, N., et al. (2022). Design, synthesis, and pharmacological evaluation of second-generation soluble Adenylyl Cyclase (sAC, ADCY10) inhibitors with slow dissociation rates. Manuscript Submitted for Publication.
Rahman, N., Ramos-Espiritu, L., Milner, T. A., Buck, J., and Levin, L. R. (2016). Soluble adenylyl cyclase is essential for proper lysosomal acidification. J. Gen. Physiol. 148, 325–339. doi:10.1085/jgp.201611606
Ramos-Espiritu, L., Kleinboelting, S., Navarrete, F. A., Alvau, A., Visconti, P. E., Valsecchi, F., et al. (2016). Discovery of LRE1 as a specific and allosteric inhibitor of soluble adenylyl cyclase. Nat. Chem. Biol. 12, 838–844. doi:10.1038/nchembio.2151
Rossetti, T., Jackvony, S., Buck, J., and Levin, L. R. (2021). Bicarbonate, carbon dioxide and pH sensing via mammalian bicarbonate-regulated soluble adenylyl cyclase. Interface Focus 11, 20200034. doi:10.1098/rsfs.2020.0034
Salomon, Y., Londos, C., and Rodbell, M. (1974). A highly sensitive adenylate cyclase assay. Anal. Biochem. 58, 541–548. doi:10.1016/0003-2697(74)90222-X
Schuetz, D. A., de Witte, W. E. A., Wong, Y. C., Knasmueller, B., Richter, L., Kokh, D. B., et al. (2017). Kinetics for drug discovery: An industry-driven effort to target drug residence time. Drug Discov. Today 22, 896–911. doi:10.1016/j.drudis.2017.02.002
Steegborn, C., Litvin, T. N., Levin, L. R., Buck, J., and Wu, H. (2005). Bicarbonate activation of adenylyl cyclase via promotion of catalytic active site closure and metal recruitment. Nat. Struct. Mol. Biol. 12, 32–37. doi:10.1038/nsmb880
Sutherland, E. W., and Rall, T. W. (1958). Fractionation and characterization of a cyclic adenine ribonucleotide formed by tissue particles. J. Biol. Chem. 232, 1077–1091. doi:10.1016/s0021-9258(19)77423-7
Swinney, D. C. (2008). Applications of binding kinetics to drug discovery. Pharm. Med. 22, 23–34. doi:10.1007/BF03256679
Visconti, P. E., Galantino-Homer, H., Moore, G. D., Bailey, J. L., Ning, X., Fornes, M., et al. (1998). The molecular basis of sperm capacitation. J. Androl. 19, 242–248.
Walkup, G. K., You, Z., Ross, P. L., Allen, E. K. H., Daryaee, F., Hale, M. R., et al. (2015). Translating slow-binding inhibition kinetics into cellular and in vivo effects. Nat. Chem. Biol. 11, 416–423. doi:10.1038/nchembio.1796
Wiggins, S. V., Steegborn, C., Levin, L. R., and Buck, J. (2018). Pharmacological modulation of the CO2/HCO3-/pH-calcium-and ATP-sensing soluble adenylyl cyclase. Pharmacol. Ther. 190, 173–186. doi:10.1016/j.pharmthera.2018.05.008
Williams, J. W., and Morrison, J. F. (1979). “[17] the kinetics of reversible tight-binding inhibition,” in Enzyme Kinetics and mechanism Part A initial Rate and inhibitor methods methods in enzymology (Amsterdam, Netherlands: Elsevier), 437–467. doi:10.1016/0076-6879(79)63019-7
Xie, F., Garcia, M. A., Carlson, A. E., Schuh, S. M., Babcock, D. F., Jaiswal, B. S., et al. (2006). Soluble adenylyl cyclase (sAC) is indispensable for sperm function and fertilization. Dev. Biol. 296, 353–362. doi:10.1016/j.ydbio.2006.05.038
Yanagimachi, R. (1994). “Mammalian fertilization,” in The Physiology of reproduction. Editors E. Knobil, and J. D. Neill (New York: Raven Press), 189–317.
Zhang, R., Barbieri, C. M., Garcia-Calvo, M., Myers, R. W., McLaren, D., and Kavana, M. (2016). Moderate to high throughput in vitro binding kinetics for drug discovery. Front. Biosci. Sch. Ed. 8, 278–297. doi:10.2741/s462
Keywords: soluble adenylyl cyclase, male contraceptive, residence time, drug development, picomolar potency, binding kinetics, lead optimization, SPR
Citation: Rossetti T, Ferreira J, Ghanem L, Buck H, Steegborn C, Myers RW, Meinke PT, Levin LR and Buck J (2022) Assessing potency and binding kinetics of soluble adenylyl cyclase (sAC) inhibitors to maximize therapeutic potential. Front. Physiol. 13:1013845. doi: 10.3389/fphys.2022.1013845
Received: 07 August 2022; Accepted: 01 September 2022;
Published: 28 September 2022.
Edited by:
Muhammad Aslam, University of Giessen, GermanyReviewed by:
Domenico De Rasmo, Institute of Biomembranes, ItalyKonstantinos Lefkimmiatis, University of Pavia, Italy
Sarah L Sayner, University of South Alabama, United States
Copyright © 2022 Rossetti, Ferreira, Ghanem, Buck, Steegborn, Myers, Meinke, Levin and Buck. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lonny R. Levin, bGxldmluQG1lZC5jb3JuZWxsLmVkdQ==