In silico Prediction of Human Oral Absorption Based on QSAR Analyses of PAMPA Permeability

The parallel artificial membrane permeation assay (PAMPA) was developed as a model for the prediction of transcellular permeation in the process of drug absorption. Our research group has measured the PAMPA permeability of peptide‐related compounds, diverse drugs, and agrochemicals. This work led to a classical quantitative structure–activity relationship (QSAR) equation for PAMPA permeability coefficients of structurally diverse compounds based on simple physicochemical parameters such as lipophilicity at a particular pH (log P_oct and |pK_a?pH|), H‐bond acceptor ability (SA_HA), and H‐bond donor ability (SA_HD). Since the PAMPA permeability of lipophilic compounds decreased with their apparent lipophilicity due to the unstirred water layer (UWL) barrier on membrane surfaces and to membrane retention, a bilinear QSAR model was introduced to explain the permeability of a broader set of compounds using the same physicochemical parameters as those used for the linear model. We also compared PAMPA and Caco‐2 cell permeability coefficients of compounds transported by various absorption mechanisms. The compounds were classified according to their absorption pathway (passively transported compounds, actively transported compounds, and compounds excreted by efflux systems) in the plot of Caco‐2 vs. PAMPA permeability. Finally, based on the QSAR analyses of PAMPA permeability, an in silico prediction model of human oral absorption for possibly transported compounds was proposed, and the usefulness of the model was examined.     

Mechanism of Long-Chain Free Fatty Acid Protonation at the Membrane-Water Interface

Long-chain free fatty acids (FFAs) play an important role in several physiological and pathological processes such as lipid fusion, adjustments of membrane permeability and fluidity, and the regulation of enzyme and protein activities. FFA-facilitated membrane proton transport (flip-flop) and FFA-dependent proton transport by membrane proteins (e.g., mitochondrial uncoupling proteins) are governed by the difference between FFA’s intrinsic pK_a value and the pH in the immediate membrane vicinity. Thus far, a quantitative understanding of the process has been hampered, because the pK_a value shifts upon moving the FFA from the aqueous solution into the membrane. For the same FFA, pK_a values between 5 and 10.5 were reported. Here, we systematically evaluated the dependence of pK_a values on chain length and number of double bonds by measuring the ζ-potential of liposomes reconstituted with FFA at different pH values. The experimentally obtained intrinsic pK_a values (6.25, 6.93, and 7.28 for DOPC membranes) increased with FFA chain length (C16, C18, and C20), indicating that the hydrophobic energy of transfer into the bilayer is an important pK_a determinant. The observed pK_a decrease in DOPC with increasing number of FFA double bonds (7.28, 6.49, 6.16, and 6.13 for C20:0, C20:1, C20:2, and C20:4, respectively) is in line with a decrease in transfer energy. Molecular dynamic simulations revealed that the ionized carboxylic group of the FFAs occupied a fixed position in the bilayer independent of chain length, underlining the importance of Born energy. We conclude that pK_a is determined by the interplay between the energetic costs for 1) burying the charged moiety into the lipid bilayer and 2) transferring the hydrophobic protonated FFA into the bilayer.     

Permeability of lipid bilayer to anthracycline derivatives. Role of the bilayer composition and of the temperature

The uptake of three anthracycline derivatives: doxorubicin, daunorubicin and pirarubicin, into large unilamellar vesicles (LUV) in response to a driving force provided by DNA encapsulated inside the LUV has been investigated as a function of the temperature and of the bilayers lipid composition. The kinetics of the decay of the anthracycline fluorescence in the presence of DNA-containing liposome was used to follow the diffusion of the drug through the membrane. For the three drugs, the permeability coefficient of the neutral form of the drug (P⁰) decreases as the amount of negatively charged phospholipid in the bilayers increases. This can be explained by the fact that the kinetics of passive diffusion of the drugs depends on the amount of neutral form embedded in the polar head group region, which decreases as the quantity of negatively charged phospholipids increases. P⁰ also decreases as the amount of cholesterol, that makes the bilayer more rigid, increases. The activation energies, E_a, for the passage of the neutral form of these anthracyclines through the bilayers lie within 100±15 kJ·mol⁻¹, except for pirarubicin and doxorubicin through anionic phospholipid-rich membranes (E_a=57 kJ·mol⁻¹) and cholesterol-rich membranes (E_a=167 kJ·mol⁻¹).     

A Permeability Study of O2 and the Trace Amine p-Tyramine through Model Phosphatidylcholine Bilayers

We study here the permeability of the hydrophobic O₂ molecule through a model DPPC bilayer at 323K and 350K, and of the trace amine p-tyramine through PC bilayers at 310K. The tyramine results are compared to previous experimental work at 298K. Nonequilibrium work methods were used in conjunction to simultaneously obtain both the potential of mean force (PMF) and the position dependent transmembrane diffusion coefficient, D(z), from the simulations. These in turn were used to calculate the permeability coefficient, P, through the inhomogeneous solubility-diffusion model. The results for O₂ are consistent with previous simulations, and agree with experimentally measured P values for PC bilayers. A temperature dependence in the permeability of O₂ through DPPC was obtained, with P decreasing at higher temperatures. Two relevant species of p-tyramine were simulated, from which the PMF and D(z) were calculated. The charged species had a large energetic barrier to crossing the bilayer of ~ 21 kcal/mol, while the uncharged, deprotonated species had a much lower barrier of ~ 7 kcal/mol. The effective in silico permeability for p-tyramine was calculated by applying three approximations, all of which gave nearly identical results (presented here as a function of the pK_a). As the permeability value calculated from simulation was highly dependent on the pK_a of the amine group, a further pK_a study was performed that also varied the fraction of the uncharged and zwitterionic p-tyramine species. Using the experimental P value together with the simulated results, we were able to label the phenolic group as responsible for the pK_a1 and the amine for the pK_a2, that together represent all of the experimentally measured pK_a values for p-tyramine. This agrees with older experimental results, in contrast to more recent work that has suggested there is a strong ambiguity in the pK_a values.     

Prediction of Secondary Ionization of the Phosphate Group in Phosphotyrosine Peptides

A computational approach, based on a continuum molecular electrostatics model, for the calculation of the pK_a values of secondary ionization of the phosphate group in phenyl phosphate derivatives is described. The method uses the ESP atomic charges of the mono-anionic and di-anionic forms of the ionizable phosphate group, computed with the use of the density functional method, and applies a new concept of the model group, being the reference state for the pK_a calculations. Both conformational flexibility and tautomeric degrees of freedom are taken into account in the calculations. The method was parameterized using experimentally available pK_a values of four derivatives of phenyl phosphates, and phosphotyrosine. Subsequently this parameterization was used to predict pK_a of the phosphate group in a short peptide Gly-Gly-Tyr(P)-Ala, and in a longer peptide consisting of 12 residues, the latter in water, and in a complex with a protein—phospholipase. The agreement between the computed and the experimental pK_a values is better than ±0.3 pH units for the optimized solute dielectric constant of 11–13. This approach is promising and its extension to other phospho-amino acids is in progress.     

A PAMPA Study of the Permeability‐Enhancing Effect of New Ceramide Analogues

There is a major need in drug discovery for quick, precise, and cost‐effective high‐throughput screening (HTS) systems in the early stages of drug research. The Parallel Artificial Membrane Permeability Assay (PAMPA) aims at predicting the passive membrane properties of drugs. Since 1998, model membranes have been developed to predict gastro‐intestinal absorption or transport through the blood–brain barrier. This paper presents recent results in a project aiming to improve the prediction of transdermal penetration. Using the PAMPA system, we investigated the effect of four newly synthetized ceramide analogues (certramides) on the permeability of three model compounds (ciprofloxacin, nifedipine, and verapamil). The certramides differ in the length of one alkyl chain, while the other alkyl chain and the head group remained the same. A relationship between the membrane concentration of certramides (from 0 to 100%) and the permeability of compounds was found, and the results of different certramides were compared. The strongest effect on permeability was caused by the ceramide analogue CTR(C₁₂–C₁₆). The reproducibility of the experiments and the impact of presence or absence of organic solvents (dodecane and CHCl₃) in the membrane were also investigated.     

Molecular interactions involved in proton-dependent gating in KcsA potassium channels

The bacterial potassium channel KcsA is gated open by the binding of protons to amino acids on the intracellular side of the channel. We have identified, via channel mutagenesis and x-ray crystallography, two pH-sensing amino acids and a set of nearby residues involved in molecular interactions that influence gating. We found that the minimal mutation of one histidine (H25) and one glutamate (E118) near the cytoplasmic gate completely abolished pH-dependent gating. Mutation of nearby residues either alone or in pairs altered the channel’s response to pH. In addition, mutations of certain pairs of residues dramatically increased the energy barriers between the closed and open states. We proposed a Monod–Wyman–Changeux model for proton binding and pH-dependent gating in KcsA, where H25 is a “strong” sensor displaying a large shift in pK_a between closed and open states, and E118 is a “weak” pH sensor. Modifying model parameters that are involved in either the intrinsic gating equilibrium or the pK_a values of the pH-sensing residues was sufficient to capture the effects of all mutations.     

Discovery of isatin and 1H-indazol-3-ol derivatives as d-amino acid oxidase (DAAO) inhibitors

d-Amino acid oxidase (DAAO) is a potential target in the treatment of schizophrenia as its inhibition increases brain d-serine level and thus contributes to NMDA receptor activation. Inhibitors of DAAO were sought testing [6+5] type heterocycles and identified isatin derivatives as micromolar DAAO inhibitors. A pharmacophore and structure-activity relationship analysis of isatins and reported DAAO inhibitors led us to investigate 1H-indazol-3-ol derivatives and nanomolar inhibitors were identified. The series was further characterized by pK_a and isothermal titration calorimetry measurements. Representative compounds exhibited beneficial properties in in vitro metabolic stability and PAMPA assays. 6-fluoro-1H-indazol-3-ol (37) significantly increased plasma d-serine level in an in vivo study on mice. These results show that the 1H-indazol-3-ol series represents a novel class of DAAO inhibitors with the potential to develop drug candidates.     

Reaction of Local Anesthetics with Phospholipids : A possible chemical basis for anesthesia

Local anesthetics (LA) have been found to interact with phospholipids and lipids extracted from nerve and muscle. This reaction is demonstrated by: (a) Inhibition by LA of phospholipid (and tissue lipid) facilitated transport of calcium from a methanol: water phase into chloroform. This action is dependent upon the cationic form of the LA. (b) LA increase the electrical resistance of "membranes" prepared by impregnating Millipore filters with cephalin:cholesterol or tissue lipid extracts and bathed with NaCl or KCl solutions. (c) LA coagulate aqueous dispersions of cephalin, phosphatidyl serine, phosphatidyl ethanolamine, and inositide, an action shared by calcium. The order of potency in coagulating cephalin sols is tetracaine > calcium > butacaine > procaine. Na⁺ and K⁺ do not coagulate phospholipid dispersions at 0.1 M concentration and antagonize the effect of Ca²⁺. (d) LA produce a marked fall in the pH of cephalin sols equivalent to that produced by calcium, (e) Ca²⁺ and LA form 1:2 molar complexes with phospholipids probably by ion-ion and ion-induced polar type of binding at the phosphate groups of the lipid. It is suggested that such reactions with cell membrane phospholipids may underlie inhibitory effects of LA on cellular ion fluxes and provide a chemical basis for anesthetic action.     

CELL CYCLE EVENTS IN THE HYDROCORTISONE REGULATION OF ALKALINE PHOSPHATASE IN HELA S3 CELLS

The increase in alkaline phosphatase in asynchronous cultures of HeLa S₃ cells grown in medium supplemented with hydrocortisone is characterized by a lag period of 10–12 hr. Present studies utilizing synchronous cell populations indicate: (a) a minimum of 8–10 hr of incubation with hydrocortisone is necessary for maximum induction of alkaline phosphatase; (b) the increase in enzyme activity produced by hydrocortisone is initiated exclusively in the synthetic phase of the cell cycle; (c) alkaline phosphatase activity does not vary appreciably over a normal control cell cycle. Radioactive hydrocortisone is rapidly distributed into HeLa cells irrespective of their position in the cell cycle, indicating that inductive effects are not governed by selective permeability during the cell cycle. Hydrocortisone-1,2-[³H] diffuses back from the cell into the medium when the cells are incubated in fresh medium containing no hydrocortisone, and the alkaline phosphatase induction, under these conditions, is completely reversible.     

The role of lipid-protein interactions in NADH-cytochrome c reductase (rotenone-insensitive) of rat liver mitochondria

The phospholipid depletion of rat liver mitochondria, induced by acetone-extraction or by digestion with phospholipase A₂ or phospholipase C, greatly inhibited the activity of NADH-cytochrome c reductase (rotenone-insensitive). A great decrease of the reductase activity also occurred in isolated outer mitochondrial membranes after incubation with phospholipase A₂. The enzyme activity was almost completely restored by the addition of a mixture of mitochondrial phospholipids to either lipid-deficient mitochondria, or lipid-deficient outer membranes. The individual phospholipids present in the outer mitochondrial membrane induced little or no stimulation of the reductase activity. Egg phosphatidylcholine was the most active phospholipid, but dipalmitoyl phosphatidylcholine was almost ineffective. The lipid depletion of mitochondria resulted in the disappearance of the non-linear Arrhenius plot which characterized the native reductase activity. A non-linear plot almost identical to that of the native enzyme was shown by the enzyme reconstituted with mitochondrial phospholipids. Triton X-100, Tween 80 or sodium deoxycholate induced only a small activation of NADH-cytochrome c reductase (rotenone-insensitive) in lipiddeficient mitochondria. The addition of cholesterol to extracted mitochondrial phospholipids at a 1 : 1 molar ratio inhibited the reactivation of NADH-cytochrome c reductase (rotenone-insensitive) but not the binding of phospholipids to lipid-deficient mitochondria or lipid-deficient outer membranes.These results show that NADH-cytochrome c reductase (rotenone-insensitive) of the outer mitochondrial membrane requires phospholipids for its activity. A mixture of phospholipids accomplishes this requirement better than individual phospholipids or detergents. It also seems that the membrane fluidity may influence the reductase activity.     

A role for phospholipids in the binding and metabolism of drugs by hepatic microsomes. Use of the fluorescent hydrophobic probe 1-anilinonaphthalene-8-sulphonate

1. The pretreatment of rat liver microsomes with phospholipase C or D decreased the N-demethylation of (+)-benzphetamine. The hydroxylation of aniline was essentially unchanged by pretreatment of microsomes with phospholipase C. 2. Some components of the microsomal mixed-function oxidase system were impaired by phospholipases. 3. The fluorescence of 1-anilinonaphthalene-8-sulphonate (ANS) was greatly enhanced by microsomes. Phospholipase C or D markedly decreased ANS–microsome fluorescence. Quantum yield of ANS–microsome fluorescence appeared to be related directly to phospholipid content of microsomes. 4. Most of the drugs studied enhanced ANS–microsome fluorescence. Warfarin, however, displaced ANS fluorescence competitively from microsomes. The latter effect was postulated as being due to warfarin competing with ANS for the cationic site on microsomal phosphatidylcholine. 5. ANS fluorescence was also increased by the presence of phospholipid micelles. The fluorescence of ANS–phosphatidylcholine micelles was modified by warfarin and (+)-benzphetamine in a manner similar to that observed with microsomes. Warfarin decrease of fluorescence was absent when ANS was bound to phosphatidic acid, which lacks a cationic site. 6. Trypsin pretreatment of microsomes did not modify ANS–microsome fluorescence, including drug-induced changes. 7. It was postulated that phospholipids have a permissive role in the metabolism of most drugs by hepatic microsomes and that the ANS probe might reflect interactions of compounds with microsomal membrane phospholipids.     

Stable Water Use Efficiency of Tibetan Alpine Meadows in Past Half Century: Evidence from Wool δ13C Values

Understanding the influences of climatic changes on water use efficiency (WUE) of Tibetan alpine meadows is important for predicting their long-term net primary productivity (NPP) because they are considered very sensitive to climate change. Here, we collected wool materials produced from 1962 to 2010 and investigated the long-term WUE of an alpine meadow in Tibet on basis of the carbon isotope values of vegetation (δ ¹³C_veg). The values of δ ¹³C_veg decreased by 1.34‰ during 1962–2010, similar to changes in δ ¹³C values of atmospheric CO₂. Carbon isotope discrimination was highly variable and no trend was apparent in the past half century. Intrinsic water use efficiency (W _i) increased by 18 μmol·mol^–1 (approximately 23.5%) during 1962–2010 because the increase in the intercellular CO₂ concentration (46 μmol·mol^–1) was less than that in the atmospheric CO₂ concentration (C _a, 73 μmol·mol^–1). In addition, W _i increased significantly with increasing growing season temperature and C _a. However, effective water use efficiency (W _e) remained relatively stable, because of increasing vapor pressure deficit. C _a, precipitation, and growing season temperature collectively explained 45% of the variation of W _e. Our findings indicate that the W _e of alpine meadows in the Tibetan Plateau remained relatively stable by physiological adjustment to elevated C _a and growing season temperature. These findings improve our understanding and the capacity to predict NPP of these ecosystems under global change scenarios.     

Comparison of phospholipid effects on insulin-sensitive low Km cyclic AMP phosphodiesterase in adipocyte plasma membranes and microsomes

Both adipocyte plasma membranes and microsomes possess insulin-sensitive low K_m cyclic AMP phosphodiesterase activity. The activity of the enzyme from both sources was susceptible to activation by several anionic phospholipids. Activators of the plasma membrane enzyme were lysophosphatidylglycerol > lysophosphatidylcholine > lysophosphatidylserine > phosphatidylserine > phosphatidylglycerol. These same phospholipids activated the microsomal enzyme but the extent of activation by each phospholipid was reversed. Neutral phospholipids and other anionic phospholipids were without effect. The phospholipids had no effect on high K_m cAMP phosphodiesterase in either membrane. The results suggest that the phospholipid headgroup was an important determinant for enzyme activation by phospholipid. The increased susceptibility of the plasma membrane enzyme to lysophospholipid may be attributed to a difference in the plasma membrane enzyme compared to the microsomal membrane enzyme or to differences in plasma membrane and microsomal membrane phospholipid composition and their ability to regulate low K_m cAMP phosphodiesterase activity.     

Excitation Energy Transfer Dynamics and Excited-State Structure in Chlorosomes of Chlorobium phaeobacteroides

The excited-state relaxation within bacteriochlorophyll (BChl) e and a in chlorosomes of Chlorobium phaeobacteroides has been studied by femtosecond transient absorption spectroscopy at room temperature. Singlet-singlet annihilation was observed to strongly influence both the isotropic and anisotropic decays. Pump intensities in the order of 10¹¹ photons × pulse⁻¹ × cm⁻² were required to obtain annihilation-free conditions. The most important consequence of applied very low excitation doses is an observation of a subpicosecond process within the BChl e manifold (~200–500 fs), manifesting itself as a rise in the red part of the Q_y absorption band of the BChl e aggregates. The subsequent decay of the kinetics measured in the BChl e region and the corresponding rise in the baseplate BChl a is not single-exponential, and at least two components are necessary to fit the data, corresponding to several BChl e→BChl a transfer steps. Under annihilation-free conditions, the anisotropic kinetics show a generally slow decay within the BChl e band (10–20 ps) whereas it decays more rapidly in the BChl a region (~1 ps). Analysis of the experimental data gives a detailed picture of the overall time evolution of the energy relaxation and energy transfer processes within the chlorosome. The results are interpreted within an exciton model based on the proposed structure.     

Phospholipid Membrane-Mediated Hemozoin Formation: The Effects of Physical Properties and Evidence of Membrane Surrounding Hemozoin

Phospholipid membranes are thought to be one of the main inducers of hemozoin formation in Plasmodia and other blood-feeding parasites. The “membrane surrounding hemozoin” has been observed in infected cells but has not been observed in in vitro experiments. This study focused on observing the association of phospholipid membranes and synthetic β-hematin, which is chemically identical to hemozoin, and on a further exploration into the mechanism of phospholipid membrane-induced β-hematin formation. Our results showed that β-hematin formation was induced by phospholipids in the fluid phase but not in the gel phase. The ability of phospholipids to induce β-hematin formation was inversely correlated with gel-to-liquid phase transition temperatures, suggesting an essential insertion of heme into the hydrocarbon chains of the phospholipid membrane to form β-hematin. For this study, a cryogenic transmission electron microscope was used to achieve the first direct observation of the formation of a monolayer of phospholipid membrane surrounding β-hematin.     

The N-terminal Domain Tethers the Voltage-gated Calcium Channel β2e-subunit to the Plasma Membrane via Electrostatic and Hydrophobic Interactions

The β-subunit associates with the α₁ pore-forming subunit of high voltage-activated calcium channels and modulates several aspects of ion conduction. Four β-subunits are encoded by four different genes with multiple splice variants. Only two members of this family, β_2a and β_2e, associate with the plasma membrane in the absence of the α₁-subunit. Palmitoylation on a di-cysteine motif located at the N terminus of β_2a promotes membrane targeting and correlates with the unique ability of this protein to slow down inactivation. In contrast, the mechanism by which β_2e anchors to the plasma membrane remains elusive. Here, we identified an N-terminal segment in β_2e encompassing a cluster of positively charged residues, which is strictly required for membrane anchoring, and when transferred to the cytoplasmic β_1b isoform it confers membrane localization to the latter. In the presence of negatively charged phospholipid vesicles, this segment binds to acidic liposomes dependently on the ionic strength, and the intrinsic fluorescence emission maxima of its single tryptophan blue shifts considerably. Simultaneous substitution of more than two basic residues impairs membrane targeting. Coexpression of the fast inactivating R-type calcium channels with wild-type β_2e, but not with a β_2e membrane association-deficient mutant, slows down inactivation. We propose that a predicted α-helix within this domain orienting parallel to the membrane tethers the β_2e-subunit to the lipid bilayer via electrostatic interactions. Penetration of the tryptophan side chain into the lipidic core stabilizes the membrane-bound conformation. This constitutes a new mechanism for membrane anchoring among the β-subunit family that also sustains slowed inactivation.     

QSAR study on permeability of hydrophobic compounds with artificial membranes

We previously reported a classical quantitative structure-activity relationship (QSAR) equation for permeability coefficients (P(app-pampa)) by parallel artificial membrane permeation assay (PAMPA) of structurally diverse compounds with simple physicochemical parameters, hydrophobicity at a particular pH (logP(oct) and |pK(a)-pH|), hydrogen-accepting ability (SA(HA)), and hydrogen-donating ability (SA(HD)); however, desipramine, imipramine, and testosterone, which have high logP(oct) values, were excluded from the derived QSAR equation because their measured P(app-pampa) values were lower than calculated. In this study, for further investigation of PAMPA permeability of hydrophobic compounds, we experimentally measured the P(app-pampa) of more compounds with high hydrophobicity, including several pesticides, and compared the measured P(app-pampa) values with those calculated from the QSAR equation. As a result, compounds having a calculated logP(app-pampa)>-4.5 showed lower measured logP(app-pampa) than calculated because of the barrier of the unstirred water layer and the membrane retention of hydrophobic compounds. The bilinear QSAR model explained the PAMPA permeability of the whole dataset of compounds, whether hydrophilic or hydrophobic, with the same parameters as the equation in the previous study. In addition, PAMPA permeability coefficients correlated well with Caco-2 cell permeability coefficients. Since Caco-2 cell permeability is effective for the evaluation of human oral absorption of compounds, the proposed bilinear model for PAMPA permeability could be useful for not only effective screening for several drug candidates but also the risk assessment of chemicals and agrochemicals absorbed by humans.     

In Depth Analysis on the Binding Sites of Adamantane Derivatives in HCV (Hepatitis C Virus) p7 Channel Based on the NMR Structure

Background

The recently solved solution structure of HCV (hepatitis C virus) p7 ion channel provides a solid structure basis for drug design against HCV infection. In the p7 channel the ligand amantadine (or rimantadine) was determined in a hydrophobic pocket. However the pharmocophore (−NH₂) of the ligand was not assigned a specific binding site.

Results

The possible binding sites for amino group of adamantane derivatives is studied based on the NMR structure of p7 channel using QM calculation and molecular modeling. In the hydrophobic cavity and nearby three possible binding sites are proposed: His17, Phe20, and Trp21. The ligand binding energies at the three binding sites are studied using high level QM method CCSD(T)/6–311+G(d,p) and AutoDock calculations, and the interaction details are analyzed. The potential application of the binding sites for rational inhibitor design are discussed.

Conclusions

Some useful viewpoints are concluded as follows. (1) The amino group (−NH₂) of adamantane derivatives is protonated (−NH₃ ⁺), and the positively charged cation may form cation-π interactions with aromatic amino acids. (2) The aromatic amino acids (His17, Phe20, and Trp21) are the possible binding sites for the protonated amino group (−NH₃ ⁺) of adamantane derivatives, and the cation-π bond energies are 3 to 5 times stronger than the energies of common hydrogen bonds. (3) The higher inhibition potent of rimantadine than amantadine probably because of its higher pK_a value (pK_a = 10.40) and the higher positive charge in the amino group. The potential application of p7 channel structure for inhibitor design is discussed.