Modelling of the Usefulness of Carbon Nanotubes as Antiviral Compounds for Treating Alzheimer Disease ()
1. Introduction
Research in nanobiotechnology has rapidly increased since the development of molecular dynamic simulations (MDSs) in 1980, X-ray crystallography and scanning tunneling microscope in 1982. Nanobiotechnology is an interdisciplinary area combining the molecular biological approach with the micronanotechnology. This combination aims to design and develop new devices by manipulating at the nanoscale their size, shape, distinct properties. It also presents revolutionary opportunities by implying the creation of nano-materials, against the infection, cancer cells and cardiac disorder, designed to interact with targeted sites in the body at sub-cellular scales with a high degree of certainty. This has led the scientific researchers to find the lightest, strongest and most conductive carbon nano-materials that capable of transporting different biomolecules through their surfaces. Carbon nanodevices are a family of very small tubes which are wholly composed of carbon atoms, having diameter measured by nanometre level (one-billionth scale) which is about to ten-thousand times smaller than the human hair, such as peptides, fullerenes, nano-rods, nano-buds, graphenes and cylindrical carbon nanotubes. Carbon nanotubes (CNTs) are selective nanoparticle because of their huge potential, low solubility and toxicity, outstanding properties, maximum loading capability and extraordinary thermal conductivity [1] [2] . They are classified into two main sub-groups; multi-walled and single-walled CNTs (MWCNTS and SWCNTs). These nanotubes can be widely functionalized with proteins, bio-active peptides and drugs, and used to deliver their loads to the targeted cells, such as for treating the infected sites, inhibiting the growth of pathogens or attacking the cancer cells [3] [4] .
CNTs have attracted a lot of interesting researches since their discovery [5] [6] . Due to their unusual properties that can be practically modified with different bio-molecules through several techniques [4] [7] . There have been several studies addressed the ability of CNTs to conjugate with a wide variety of drugs for treatment purposes. CNTs have been explored for disease therapy applications, especially for cancer treatment [8] . In addition, they are also used to build up smaller, lighter and more efficient nano-sensors as scaffolds for tissues repair and cell growth [9] [10] [11] . To drive the drug delivery systems into the lymph-node cancer cells, synthesized MWCNTs with magnetite nano-particles were functionalized folic-acid and successfully loaded with Cisplatin (an anti-cancer drug) [12] [13] . Furthermore, CNTs have also been employed as drug delivery agents for treatment of hypertension (carvedilol) [14] , asthma (theophylline) [15] , human immunodeficiency virus (HIV) [16] and inflammation (dapsone, dexamethasone and ibuprofen) [17] [18] [19] [20] as well as some of brain diseases [21] [22] , e.g., Alzheimer disease (AD) [23] . SWCNTs have recently been used as safe carriers to increase the effectiveness of AD treatment (against Mitchidoria in the brain) [23] .
AD is the most common type of dementia, progressive and irreversible brain disorder that causes problems with behavior and slowly destroys thinking skills and memory. Its symptoms often develops gradually and gets worse overtime and classified into three stages; mild, moderate then becomes severe, which can be noticeably seen in the X-ray image shown in Figure 1, interfering with daily life tasks. AD has no current and effective cure so far to get rid of the progressive form of Alzheimer. Nowadays, treatment of AD symptoms are available and researchers still paying more attentions to find new techniques to slow worsening
Figure 1. Schematic geametry (X-ray image) expose the difference between the health brain, and the mild and severe stages of Alzheimer disease.
of losing memory, improve quality of daily tasks and prevent Alzheimer’s symptoms from developing. Yang et al. [23] who have shown that SWCNTs with diameter 0.8 to 1.6 nm and variant lengths 5 - 300 nm were successfully carried out with Acetylcholine (ACh) then directly delivered into brain, lysosomes are the targeted organelles of SWCNTs not mitochondria (MTCHD). This technique used the SWCNT as a catalyst due to its lack in the brain areas to increase the effectiveness of ACh drug. ACh can not be delivered directly into the brain with low doses because of its poor lipophilicity, this has recently been administrated that lack of ACh can be overcome by using SWCNT as carrier with high bio-safety [13] . Inside the brain, SWCNTs enter directly into the lysosomes, not to the MTCHD, which are the targeted organelles, but the high doses rise up the opportunity of SWCNTs to enter to the MTCHD. SWCNTs can be successfully exploited to deliver ACh into the targeted lysosomes to achieve therapeutic effects without undesired toxic effects [23] .
The proposed model is designed to investigate the mechanism of encapsulation of two different antiviral compounds, Acetylcholine (used for the early stages) and Ravastigmine (RAV) (used for the moderate and severe stages) with chemical formulas
and
as shown in Figure 2(a) and Figure 2(b), respectively, inside the SWCNTs of variant radii rc . Firstly, we obtain the mathematical geometry for each antiviral compound and also evaluate the total energy arising form each antiviral compound encapsulated inside a SWCNT of radius rc . ACh structure comprised as two-connected sphere with different radii rs , while RAV modelled as two possible structres; an ellipsoid and cylinder, each interacting with the interior wall of SWCNTs with variant radii rc . This paper is structured as, in the first section, we briefly outline the significance of usefulness of SWCNTs in many modern applications at nano-scales. Next, we apply the van der Waals force, Lennard-Jones potential and a discrete-continuum approximation to investigate the medical application which describes the possibility of conjugation between a SWCNT and antiviral compounds specialized to treat AD by inhibiting the growth of MTCDH (infected site in the brain). We also determine the magnitude of the potential energy arising from the interaction between antiviral compounds and nanotubes with various radii rc . Followed by a discussion and analysis our results in section 3. Finally, conclusions and remarks are given in the last section.
Figure 2. Schematic geometry for: (a) ACh antiviral compound as an interior atom interacting with a SWCNT of radius rc (sb:single bond and db:double bond); (b) RAV antiviral compound as an interior atom interacting with a SWCNT of radius rc (sb:single bond and db:double bond); (c) Antiviral compound as an interior atom at point P off-setting from the central-axis by a distance δ.
2. Mathematical Model
Here, we obtain the bio-physical model which describes the adsorption of two different drugs into SWCNTs with variant radii rc as a mathematical model by using van der Waals forces. We use the Lennard-Jones potential and the continuum approach together to model the encapsulation of these drugs inside a SWCNT. Next, we use the Cartesian coordinate
as a reference system to model the two interacting molecules, the certain biomolecule and the cylindrical nanotube. We assume that the point at atom P has coordinates
and the SWCNT defined as a cylindrical tube parameterized by
, where
,
,
and
as shown in Figure 2(c). Therefore, the distance
from point P to typical surface element of the cylindrical tube given by
(1)
The Lennard-Jones potential given as
(2)
where
is the potential function,
denotes the distance between two molecular structures, and A and B are the attractive and repulsive constants. The physical parameters,
and
, are calculated by using the empirical combining laws given by given by
,
and
, where
is the well depth,
is the van der Waals diameter and
is the non-bond energy [24] [25] . Here, we apply continuum approximation, atoms are assumed to be uniformly distributed over the surfaces of the two interacting molecules, to evaluate the interaction energy between two well-defined molecules by performing double integral over the surface of each molecule. From the work of Thamwattana et al. [26] , the interaction energy is given by
(3)
where
is the atomic surface densities of atoms on the nanotube and
is a typical surface element located on the interacting molecule. The integral
(
) is defined by
(4)
we may re-write the equation 4 by using the hypergeometric function as
(5)
Next, we assume that the atom at point P is within the volume element of each biomolecule. Thus, we can determine the molecular interaction arising from the certain drug by performing the volume integral of
over the volume of the certain drug, namely
(6)
where
is the distance from the nanotube axis to a typical point of the certain biomolecule and
is the mean volume density of the biomolecule, which depends on the assumed configuration of the interacting biomolecule and
can be given as
(7)
2.1. Insertion of ACh as Two-Connected Spheres into SWCNT
Here, we assume ACh structure modelled as two-connected spheres, the larger sphere centred at the origin point with radius
and the samller sphere with radius
located on the left side of the origin point. Each sphere assumed to be as a spherical shell parameterized
, where
,
,
and
is the radius of the spherical shell as shown in Figure 3(a). Further, the distance is given by
and the spherical volume element is
. From the work of Thamwattana et al. [26] , the interaction energy between a spherical molecule and a cylindrical nanotube is given as
(8)
where
and
are the atomic volume densities of the cylindrical nanotube and spheroidal molecule, respectively. So, the integral
(
) is given by
(9)
by using the relation between the beta and hypergeometric functions,
can be expressed in terms of
(10)
2.2. Insertion of RAV into SWCNT
To evaluate the total energy arising from the RAV drug interaction with SWCNT of radius rc , we consider two possible structures as models for RAV molecule which are an ellipsoid and cylinder as shown in Figure 3(b) and Figure 3(c), respectively.
2.2.1. An Ellipsoid Model
The RAV molecule assumed to be as a spheroidal structure, parameterized by
, where
,
,
, and a and b are the equatorial semi-axis length and polar semi-axis length (along the z-axis) of spheroidal structure, respectively, as shown in Figure 3(b). Further, the distance is given by
and the spheroidal volume element is
. From the work of Thamwattana et al. [26] , the interaction energy between a spheroidal molecule and a cylindrical nanotube is given as
(11)
where
is the mean volume density of the spheroidal molecule, respectively. So, the integral
(
) is given by
(12)
The Integral
can be expressed in terms of
Figure 3. The possible configuration for each antiviral compound: (a) ACh molecule splitted as two-connected sphere located on the left and right sides of the origin point; (b) RAV molecule as an ellipsoid structure; (c) RAV molecule as a perfect cylinder, each configuration interacting with a SWCNT of radius rc .
(13)
To obtain and evaluate the interaction energy for each configuration as shown in Figure 3, we need to determine the potential energy arising from the specific atom at point P inside the cylindrical nanotube as shown in Figure 2(c) (this atom is withing the volume of RAV).
2.2.2. Cylindrical Model
Here, we model the RAV molecule modelled as a perfect cylinder located at the origin (centered) with radius a and length
as shown in Figure 3(c). A typical point in the cylinder can be parameterized by
, where
,
and
. Therefore, the distance
is given by
and the volume element of cylinder is
. From Thamwattana’s work et al. [26] , the interaction energy between a cylindrical molecule and a cylindrical nanotube given as
(14)
where
is the mean volume density of the cylindrical molecule. So, the integral
(
) is given by
(15)
So,
is given as
(16)
3. Results and Discussion
In this section, we apply Lennard-Jones potential and the discrete-continuum approach to evaluate the interaction energy of each drug interacting inside SWCNTs with variant radii rc . The non-bond energy, well-depth
and van der Waals diameter
are shown in Table 1. The physical parameters and illustrated radii rc of CNTs are given in Table 2. The attractive and repulsive constants are calculated by using the combining laws (
and
) and are given in Table 3. The volume density for each configuration calculated as the total number of atoms that are containing the specific molecule are divided by the volume of the molecule structure, spherical shape (
), spheroidal structure (
) and cylindrical shell (
), which are
,
and
, respectively. Next, we evaluate and plot the minimum energies (for all configurations) which are arising from the ACh-SWCNT and RAV-SWCNT interactions. We also deduce the critical radius of SWCNT that will accept both antiviral compounds (ACh-SWCNT and RAV-SWCNT). The minimum energies for all configurations are obtained based on the equilibrium position of each molecule being away from the interior wall of CNT and its radius rc along the range of z-axis. In this model, we observe the encapsulation of ACh and RAV inside the nanotubes with radius in the range 3.204 Å < rc < 7.551 Å and the minimum energies obtained for both configurations when rc greater than 3.391 Å. The lowest interaction energy for ACh-SWCNT and RAV-SWCNT interactions is obtained when the radius of nanotube in the range 3.86 Å < rc < 4.07 Å as shown in Figures 4-6.
For the three proposed configurations, we note that the both antiviral compounds, ACh-SWCNT and RAV-SWCNT, are repulsive and unstable when
and
Å, respectively, and the (9, 2) SWCNT of radius
Å is the most favorable nanotube followed by the condition where
, 3.775, 4.615, 3.590, 3.523, 5.523, 6.102, 7.551 and 3.391 Å, respectively. For all interactions, ACh-SWCNT (connected-spheres), RAV-SWCNT (spheroidal) and RAV-SWCNT (cylindrical), are with minimum energies of approximately −0.664, −1.059 and 1.204 kcal/mol, respectively. Furthermore, we can noticeably see that the magnitude of the minimum energy for RAV-SWCNT (an ellipsoid structure) interaction is slightly smaller than that of RAV-SWCNT (perfect cylinder). We also note that the perfect cylinder (RAV) has the
Table 1. The Lennard-Jones constants (ε: Bond length, σ: Non-bond distance and ζ: Non-bond energy) (single bond: sb, double bond: db) [24] [27] [28] .
Table 2. Parameters for carbon nanotubes, ACh and RAV molecules.
Table 3. Numerical values of the significant constants (A and B) involved in this model.
Figure 4. Interaction energy (E) arising from the encapsulation of two antiviral compounds (ACh and RAV) inside SWCNTs with varian radii
a)
b)
c)
d)
Å.
Figure 5. Interaction energy (E) arising from the encapsulation of two antiviral compounds (ACh and RAV) inside SWCNTs with varian radii
a)
b)
c)
d)
Å.
Figure 6. Interaction energy (E) arising from the encapsulation of two antiviral compounds (ACh and RAV) inside SWCNTs with varian radii
a)
b)
c)
d)
Å.
maximum binding energy because of its volume 502.804 Å3 being larger than that of an ellipsoid structure (RAV) which is 167.601 Å3. This means that smaller size of an ellipsoid ends requires smaller size of nanotube to accommodate the spheroidal shell (RAV) compared with that of cylindrical structure (RAV) despite having similar dimensions. Moreover, we observe that our results consistently agree with the most recent research findings, for example, Dresselhaus et al. [1] predict that the (5, 5) CNT of
could be the most significant and smallest effective physical nanotube, and ACh drug can be carried with SWCNT of radius in the range of 4 Å ≤ rc ≤ 8 Å (8 Å < diameter = 2rc < 16 Å) delivered to the target and infected cells [23] .
4. Conclusion
In this study, the Lennard-Jones potential and continuum approach are adopted to evaluate the minimum energy for each configuration. The proposed model obtained mathematically by representing each molecule using the rectangular coordinate
as a reference system. Through investigation, we find that the SWCNT plays a significant role by increasing the effectiveness of the antiviral compounds against the growth and symptoms of the AD. The SWCNT is a selective tool because of its distinct properties, such as high conductivity and low solubility in aqueous media. It can be concluded that the RAV antiviral compound is more effective against the AD growth, and both antiviral compounds, ACh and RAV, would not be accepted when
Å. For all possible configurations, we note that the lowest minimum energy obtained when
Å. Our results are in very good agreement with Yang’s work who has shown that the ACh antiviral compound is successfully carried out and conjugated with SWCNTs with variant radii
[23] .
Acknowledgements
The author acknowledges financial support from the Research and Consultation Centre (RCC) at the University of Business and Technology.
Funding
This project funded by the Research and Consultation Centre (RCC) at the University of Business and Technology.
Ethical Approval
This article does not contain any studies with animals performed by any of authors.