Awards Nomination 20+ Million Readerbase
Indexed In
  • Open J Gate
  • Genamics JournalSeek
  • ResearchBible
  • RefSeek
  • Directory of Research Journal Indexing (DRJI)
  • Hamdard University
  • OCLC- WorldCat
  • Scholarsteer
  • Publons
  • MIAR
  • Euro Pub
  • Google Scholar
Share This Page
Recommended Webinars & Conferences
Journal Flyer
Flyer image

Research Article - (2022) Volume 11, Issue 11

Biostatistical and Mathematical Analysis on Liver Disease during COVID-19 Pandemic
Bin Zhao* and Xia Jiang
Department of Science, Hubei University of Technology, Hubei, China
*Correspondence: Bin Zhao, Department of Science, Hubei University of Technology, Hubei, China, Email:

Received: 19-Sep-2022, Manuscript No. BABCR-22-18260; Editor assigned: 22-Sep-2022, Pre QC No. BABCR-22-18260 (PQ); Reviewed: 07-Oct-2022, QC No. BABCR-22-18260; Revised: 13-Oct-2022, Manuscript No. BABCR-22-18260 (R); Published: 21-Oct-2022, DOI: 10.35248/2161-1009.22.11.465


The liver is the body's biggest organ and is involved in metabolic processes. The liver is involved in both metabolism and dehydration. It helps in the detoxification of harmful chemicals as well as the growth of microorganisms. For systemic removal and preservation of organisms and other creatures, the liver is the most essential organ. As a result, liver injury has significant ramifications. Liver disease is regarded as a serious health issue. Overuse of several medicines, which are occasionally prescribed as part of treatment, can cause organ damage. Hepatotoxicity is caused by other chemical agents, such as those employed in labs (thioacetamide, alcohol, etc.) and industry, as well as natural compounds (such as microcystin). Many over-the-counter medicines aren't prescribed to help people with liver disease, and they can actually harm the liver. As a result, plant-based treatments are utilised to treat liver illness. As a result, several human herbal treatments have been evaluated in experimental animal models for antioxidant and hepatoprotective liver function. Acanthus ilicifolius was used as Traditional Chinese medicine (TCM) and Traditional Indian medicine (TIM). The plants showed many clinical properties. Still, the neurological related functions and disorders are not well explored in this plant. Complex interplay of positive and negative emotions orchestrated by intricately associated neuronal circuits, neurotransmitters coupled with endocrinal influence holds responsible for human behavior, considered as the root of human civilization, is currently facing existential crisis during COVID-19 pandemic. In the present study, an attempt was made to identify the interaction between A. ilicifolius natural compounds and Echinacoside as reference compounds were to study the neurotransmitters functions through biomathematical and computational method. Initially, in silico molecular docking was performed to identify the potent natural compounds against neurological disease. The results show among 8 natural compounds, 26.27-Di (nor)-cholest-5,7, 23-trien-22-ol, 3-methoxymethoxy, Cholest-5-en-3-ol (3, Beta.)-, carbonochloridate, Cholesterol and Echinacoside exhibited maximum interaction with all the target proteins. Especially, Echinacoside exhibited the maximum interaction with (Serotonin) 5-hydroxytryptamine receptor 2A (-17.077), Sodium-dependent serotonin transporter (-15.810) and (Histamine) Histamine H2 receptor (-17.556). These two neurotransmitters act as a major concern related to the mental disorders and neurological functions. The natural compounds may potent inhibitor for neurological disorders.


Liver; Herbs; Neurotransmitters; Acanthus ilicifolius; Echinacoside


The liver is an essential organ that regulates many physiological processes in the body. This contributes to important functions, such as metabolism, excretion, food resistance and energy. The liver is supposed to defend itself from the hazards of medicines and chemicals in addition to performing medicinal functions. The liver plays an important role in many important physiological processes such as glucose homeostasis, essential protein production, lipoprotein and lipid production, bile acid excretion and vitamin accumulation. It is widely reported that foreign substances, hepatotoxins, and chemotherapy impair its function. Damage due to chronic alcohol abuse, bacterial hepatitis, or metabolic disorders, cellular necrosis, fibrosis, increased tissue oxidation, and reduced tissue glutathione levels are all linked to liver injury. By introducing lipid peroxides and other harmful effects, a vast range of hazardous substances harm liver cells. The most common causes of liver disease are viral infections, medications, toxic chemicals, alcohol abuse, and autoimmune problems. A large number of hepatotoxic chemicals destroy liver cells due to lipid peroxidation and other oxidative damage. One of the most serious health issues is liver disease and steroids, antibiotics and antibiotics are recommended as treatments for liver disease. The treatment has many side effects especially when used for long time. There is a universal method of using traditional herbal methods to treat many liver diseases. Natural herbal remedies have been known to be hepatoprotective agents that contain a wide variety of chemicals. Herbal medicines have gained popularity and popularity around the world because to their safety in recent years, security and cost. According to the World Health Organization (WHO), about a quarter of the world’s population uses traditional or herbal medicine for a variety of ailments, including liver disease. Medicinal herbs have been proven to be effective in treating viral and communicable illnesses in several investigations. Allopathic medication can be used to cure hepatitis and liver damage in clinical trials, but reports of numerous adverse effects, such as immunosuppression, bone loss, and renal failure, led other herbal medicine researchers to believe that it had little impact, but without experiencing it. It is quite successful in the treatment of hepatitis.

Materials and Methods

Herbal medicine dates back to 2100 BC in ancient China's Xia dynasty, this was written for the first time in 600 BC. There are more than 300 preparations in the traditional Indian medical system for treating jaundice and chronic liver disease. The semi-synthetic properties of these compounds will also be useful in the pharmaceutical sector. The ingredients found in this ingredient provide an excellent therapeutic effect with few side effects and low cost. Antioxidants have also been shown to reduce toxicity. Natural antioxidants are broken down into several compounds as secondary metabolites.

Polyphenols (phenolic acids and flavonoids), terpenoids (carotenoids), and foods containing significant amounts of these chemicals, for example, have an essential role in illness prevention. Plant-derived medicines are utilised in medicine. The current study is based on reported research work on hepatoprotective phytochemicals from different medicinal plants tested for various types of hepatotoxicity [1-17].

Collection and authentication of plant

A. ilicifolius leaves were collected from Cuddalore District, Tamil Nadu, India and authenticated by the Botanical Survey of India, Tamil Nadu Agriculture University, Coimbatore, Tamil Nadu, India. (BSI/SC/5/23/09-10/Tech. 306). A voucher specimen of the plant has been deposited at the Herbarium of Botanical Survey of India.

Preparation of the leaf extracts

The fresh leaves of A. ilicifolius were washed and shade dried at room temperature (28 ± 2°C). The dried leaves were powdered by the electrical blender. 25 gms of A. ilicifolius leaf powder was used for methanol extraction in the Soxhlet apparatus [18]. The solvent was boiled gently at 64°C in a heating mantle until the extraction was done. Then the solvent was evaporated using a rotary vacuum evaporator to yield a viscous dark green residue of methanol leaf extracts.

Identification of phytochemicals

The GC-MS of A. ilicifolius methanolic leaf extract was identified compounds are 26.27-Di (nor)-cholest-5, 7, 23-trien-22-ol, 3-methoxymethoxy (RT=12.31). 9H–purin-6-amine, N, 9-bis (trimethylsilyl)-8-((trimethylsilyl) oxy) (RT=14.09). Cyano colchicines (RT=6.06). 3Beta-methoxy-5-cholesten-19-oic acid (RT=18.46) [19].

Cholest-5-en-3-ol (3, Beta.)-, carbonochloridate (RT=25.978), Cholesterol (RT=27.518), Cholest-5-en-3-ol (3, Beta.)-propionate (RT=28.51) and Echinacoside [20]. Infrared spectroscopy identifies the functional group present in the above-listed compounds, the presence of alcohols and phenols in the O-H region at 3389 cm-1.

In silico studies

Preparation of ligands and standard drug: The hepatoprotective power of the hydrochemical extraction of Aervalanata extract against paracetamol, which induces hepatotoxicity, was discussed. Liver damage can occur with the use of paracetamol up to 3 g/kg. In paracetamol-fed mice, liver enzymes such as ALT, AST, ALP, and serum bilirubin rose substantially (P<0.001) compared to normal controls.

When mice were given Avera lanata extract instead of paracetamol, drug-induced hypersensitivity factors such as liver enzyme (p<0.01), and ALT (p<0.001) changed significantly (p<0.01). Hepatoprotective silymarin (25 mg/kg) given as standard resulted in a substantial decrease in liver enzyme and bilirubin levels (p<0.001) when likened to paracetamol.From the above analysis, it can be assumed that the medicinal plant Aervalanata has been shown to be a liver-protecting herb.

The seven major phytochemical compounds considered are: 26.27-Di (nor)-cholest-5, 7, 23-trien-22-ol, 3-methoxymethoxy, 9H –purin-6-amine, N, 9-bis (trimethylsilyl)-8-((trimethylsilyl) oxy), Cyanocolchicines and 3Beta-methoxy-5-cholesten-19-oic acid. Cholest-5-en-3-ol (3, Beta.)-, carbonochloridate, Cholesterol, Cholest-5-en-3-ol (3, Beta.)-, propionate and Echinacoside structures were retrieved from Protein Data Bank (PDB) and ISIS Draw 2.3 software (freeware) ( was used to design the ligands.

Analogs were changed into MOL files and 3D optimization was done by ChemSketch 3D viewer of ACDLABS 8.0.The standard drug was used for comparison. Acetylcholine, Dopamine, GABA, Glutamic acid, Norepinephrine, Serotonin, and Histamine were collected from PubChem.

ADME property of active components

Lipinski rule of five is used to check the Adsorption, Distribution, Metabolism, and Excretion (ADME) orally active drug in humans. This was done by ADME tool. OSIRIS Property Xplorer was used to validate the drug molecule (Active phytochemicals) which has an inhibitory effect on the modeled target protein. Properties such as mutagenicity, irritant, tumorigenic and drug likeliness of the phytochemicals were studied [21,22].

Collection of target proteins for HCC

The lists of target proteins for neurotransmitters were collected through a literature survey and structures were retrieved from the Protein Data Bank (PDB) shown in Table 1.

S. No. Gene name PDB ID Protein name References
1 Acetylcholine 5CXV Muscarinic acetylcholine receptor M1 [23,24]
3UON Muscarinic acetylcholine receptor M2
4U14 Muscarinic acetylcholine receptor M3
5AFH Neuronal acetylcholine receptor subunit alpha-7
2 Dopamine 6CM4 D2 dopamine receptor [25]
3 GABA 4MR7 Gamma-aminobutyric acid type B receptor subunit 1  [26]
4 Glutamic acid 2ZNT Glutamate receptor ionotropic, kainate 1 [27]
5 Norepinephrine 2R4S Beta2 adrenoceptor [28]
6 Serotonin 5TUD 5-hydroxytryptamine receptor 2A [29]
5I71 Sodium-dependent serotonin transporter
7 Histamine 3SN6 Histamine H2 receptor [30]

Table 1: List of neurotransmitters receptors.

Molecular docking of A. ilicifolius phytochemicals against HCC target proteins

The structure was minimized using OPLS-2005 force field with Polack-Ribiere Conjugate Gradient (PRCG) algorithm. The Schrodinger Glide program version 2017 has been used for docking. The best 10 poses and corresponding scores have been evaluated using Glide in single precision mode (GlideSP) for each ligand [23-30]. For each screened ligand, the pose with the lowest Glide SP score has been taken as the input for the Glide calculation in extra precision mode (Glide XP). The docking was carried out with the following non-default settings in Glide SP and Glide XP both [31].

Results and Discussion

The leaves and fruits of Azadirachta indica play a vital function in traditional medicine, and they are a key element of its therapeutic properties. Imbidine, margosic acid, nimbine, and polysaccharides are the active phytochemical components. Medical activities include antiemetic, antineoplastic, emetic, antisyphilitic, emollient, stimulant, aphrodisiac, carminative, antiasthmatic, astringent, antihelminthic, antiseptic, antimalarial and anti-inflammatory.

Azadirachta indica juice was found to decrease paracetamol-induced lipid peroxidation in one research. It is given orally in a dose of 200 mg/kg with freshly squeezed juice. Orally, a dosage of 2 g/kg paracetamol is given. The decrease in sulfhydryl and group hepatocytes was avoided by drinking freshly squeezed juice. The concentration of liver enzymes in the blood increases after taking paracetamol. Enzyme levels in the blood stay constant when Azadirachta indica is administered. Ashfaq et al. discovered a phytochemical (3-deacetyl-3-cinnamoyl-azadirachtin) that has antibacterial action against the HCV NS3 protease, as well as molecular and mimetic mechanisms that can serve as NS3/4A protease inhibitors, in the literature. Gupta and Chaphalkar examined the immunomodulatory and anti-inflammatory properties of A. indica extract in the presence of a range of antibiotics (DPT, hepatitis and rabies). Antimicrobial DPT, HBsAg, and antibiotics showed lower IgG levels, as well as fewer granulocytes, monocytes, and lymphocytes in whole blood, than A. indica. As a consequence, A. indica-containing plants are a strong candidate for commercial vaccines that not only promote, but also sustain, high immune status against HBsAg, DPT, and rabies.

Structure and function of the neurotransmitters by Bio-mathematical model

Bombacaceae are medicinal plants that can come from the tropical regions of the Indian subcontinent. Its treatment has been carried over into traditional health systems such as Ayurveda, Siddha and Unani. The protective function of removing methanol from Bombax cepa L flowers tested against isoniazid and rifampicin produced hepatotoxicity and high levels of albino Wistar. Hepatotoxicity in mice was associated with two antituberculous drugs, isoniazid and rifampicin for 10-21 days intraperitoneally in mice. Three doses of Bombax ceiba L. (150, 300 and 450 mg/kg of active ingredient) were administered to study the protective function of the liver. From the normal range, there was a substantial increase in AST, ALT, ALP, total bilirubin, and a significant decrease in total isoniazid protein and rifampicin levels. Taking different drugs of Bombax Cipa-L can alter the levels of enzymes in the blood. Obviously this plant reduced TBARS levels and increased GSH levels at all doses compared to normal control. The liver tissue of animals treated with Bombax ceiba L. protects against liver damage caused by isoniazid and rifampicin. In addition to the hepatotoxicity of this plant, it is also widely used for its anti-inflammatory, antihypertensive, vascular and antioxidant activity.

Calotropis gigantean ethanol fruit was discovered for its hepatoprotective properties. Ethanol produces liver damage. High levels of AST, ALT, ALP and LPO, as well as low levels of vitamin C, indicate significant liver damage in toxic animals. Treatment before weight gain of 250 to 500 mg/kg of Calotropis gigantean tissue released reduced levels of enzyme markers, lipid peroxidation, and an increase in blood vitamin C. The pharmacological effects of the extract may be due to its antioxidant ability to prevent lipid peroxidation and prevent damage to vitamin C.

Citrus retinol essential oil has been listed for its hepatoprotective properties. Is isoniazid (50 g/kg, po. For 30 days) developed hepatotoxicity. Essential oils (200 mg/kg, p.o.) were administered every 24 hours for 30 days, while a standard received Liv52. At the end of the study, the red blood cell content was estimated. The treated isoniazid showed only a significant increase in ALT, AST and bilirubin and a significant reduction in total protein compared to normal (non-isoniazid-associated) control animals. Animals treated with cross-linked citrus and Liv52 essential oils showed significant reductions in all studies. In other words, the above results suggest that the citrus lattice has a potent hepatoprotective function against liver damage caused by isoniazid.

Hepatotoxic effects of cissaampelos pareira ethanolic root extract at different dosages of 100 mg/kg, 100 mg/kg, 100 mg/kg, and ccl4 in rats were discovered. SGOT, SGPT, ALP, and bilirubin levels were measured to evaluate liver function. Hepatotoxicity can be induced in this case by mixing ccl4 with olive oil (weight 2 ml/kg, i.p.). The hepatotoxicity group had higher levels of biosystems such SGOT, SGPT, ALP, and bilirubin. The ethanol extraction of cissaampelos pareira root has a strong protective effect, lowering the biochemical potential to near-normal levels. Historical research shows control of root extract in ccl4, indicating safety in the liver.

The therapeutic effect of the removal of toxic fruit and sharp seeds of Embeliatsjeriam cottam and liver damage with isoniazid in mice was investigated. One treatment with isoniazid (50 g/kg for 30 days) showed an increase in ALT, AST, bilirubin and significant reduction in total protein compared to normal control. An animal treated with an alcoholic drink showed what showed severe side effects with all the ingredients of the chemical. In vivo peroxidation studies showed that mice from the isoniazid treatment group demonstrated this significant increase in Malondialdehyde (MDA) compared to the high control dose. The use of alcohol and intoxicants reduces MDA levels. There was a decrease in GSH levels, SOD activity and treatment with CAT and isoniazidone over a normal control. Levels of GSH, SOD and CAT increased significantly in alcohol and single drug extraction. The subjects' liver biopsy studies showed normal hepatocytes and normal liver morphology. Liver of killer mice with isoniazid Cavity, sinus congestion, mild inflammation and damage of hepatocytes with central necrosis are present. Both alcohol and embeliatsjeriam-cottam alcohol extracts can reduce hepatotoxicity caused by isoniazid.

Orthosiphon is an important ingredient commonly used in Malaysia to treat various types of liver diseases. Hepatotoxicity and hepatotoxicity of thioacetamide were studied with ethanol supplementation in Sprague Dawley rats. A maximum of 40 samples were taken for the experimental design. The five groups of adult rats were group 1 (control), group 2 thioacetamide as an active control (one hepatotoxic group), and group 3 silymarin as a known medicine (hepatoprotective group) (each comprising eight animals), as well as tiny and micro dosages, belong to groups 4 and 5 (single treatment). All rats were killed after two months of treatment. When compared to normal animals, the hepatotoxic group showed a substantial rise in HCC (ALT, AST, ALP, and bilirubin) and malondialdehyde (MDA) levels, as well as a significant decrease in protein and albumin levels. The high therapeutic effect and ethanol extract of Orthosiphon Stamineus (200 mg/kg) alters liver function enzymes in high doses to be beneficial. These studies indicate that Orthosiphon stamineus is an important Ayurvedic plant with hepatotoxic effects on thioacetamide-induced hepatotoxicity in rats.

The hepatoprotective and antioxidant activity of the methanolic extract of phyllanthus polyphilus against paracetamol-induced hepatotoxicity was discussed. Plant production showed strong hepatoprotective effect and antioxidant activity against paracetamol leading to hepatotoxicity and mice. Liver enzyme and Lipid Peroxidase (LP) were all elevated in animals given 750 mg/kg paracetamol. TC and Total Protein (TP) are reduced, as are Superoxide Dismutase (SOD), catalase, Glutathione Peroxidase (GPx), and Glutathione S-Transferase (GST). Different therapeutic medicines (200 and 300 mg/kg, p.) were used that Phyllanthus polyphyllus has strong liver-in-protective and antioxidant properties.

The protective effect of ethanol extracted from Pyrenacanthastaudtii and Carbon Tetrachloride (CCl4) was induced by hepatotoxicity in rats in experimental studies. Sixteen Wistar rats weighing 100-170 g were divided into four groups, divided into groups I, II, III and IV. The second, third and fourth doses of 5 l/kg of CCl4 weight were injected intraperitoneally while 5 ml/kg of crude oil was injected in the normal control. After 48 hours, Pyrenacanthastaudtii ethanol of groups III and IV (750 mg/kg and weight 1500 mg/kg) was administered, respectively. The mice were then removed 5 days later. Based on the research on drivers, it was shown that two ethanol (750 mg/kg and 1500 mg/kg weight) reduced the levels of liver enzymes which increased significantly with the treatment dose ccl4. CCl4 animals had increased bilirubin levels in mice compared to mice. But two levels of ethanol significantly altered bilirubin levels depending on the type. The release of Pyrenacanthastaudtii has hepatoprotective effects against hepatotoxicity and damage caused by CCl due to the presence of flavonoids and other bioactive substances.

Tecomella undulate

Known for its methanol derivatives of the famous herb, Tecomella leaf damage against alcohol and acetaminophen caused liver damage and excessive albino. Levels of blood enzymes such as AST, ALT, ALP, GGT (gamma-glutamyl transpeptidase) and total bilirubin in the blood associated with the activity of LPO (lipid peroxidase), SOD (superoxide dismutase), CAT catalase), GSH (reductase) GPx A glutathione peroxidase in the liver is treated in this way. Historical changes in liver function were also compared with related controls. The enzymatic component of hepatic oxidative stress markers (SOD, CAT, GSH and GPx) increased significantly (p<0.001) and LPO was significantly reduced in alcohol consumption by 30% by paracetamol. In 30 percent alcohol and paracetamol mice, oral therapy with Tecomella undulate (100 mg/kg, 200 mg/kg, btBt/day) altered all serological and hepatic alterations. The medication silymarin (25 mg/kg bw/day) was also used to compare biochemical data. These findings show that the presence of flavonoids, quinones, and other antioxidants in ticumilla wafer leaves may give liver protection.

In rats, the hepatotoxicity of Ziziphus mauritiana leaf extract was reported to produce hepatotoxicity. Wistar Albino dosages of 100-120 g were given in one of six single doses and the animals were treated for six weeks. Prolonged alcohol intake (40 percent v/v, 1 ml/100 g), increased for six weeks (p 0.05), and liver enzyme and total bilirubin levels might cause hepatotoxicity (TB). In comparison to normal controls, the levels of catalase, glutathione peroxidase, glutathione reductase, and superoxide dismutase all decreased significantly (p 0.05).Meanwhile, therapy with Ziziphus mauritiana 200 mg/kg or silymarin 100 mg/kg resulted in near-normal ALT, AST, ALP, and TB levels that were statistically significant (p 0.05). When Ziziphus mauritiana was given instead of alcohol, excellent antioxidant measures including catalase, glutathione peroxidase, glutathione reductase, and superoxide dismutase were significantly increased. Cell death, fat burning, and other histopathological features of alcoholic liver suggest liver injury. Rehydration treatment with Ziziphus mauritiana or silymarin reduces tissue damage and cell death caused by persistent drinking. Ziziphusmauritiana's protective role is to block antioxidant substances such tannins, saponins, and phenolic compounds, among others.

Typical neuron mathematical equation

A typical neuron is generally classified into three parts namely cell body, dendrites, and axon [32]. The cell body contains the nucleus and associated intracellular structures whereas; dendrites are the extension of the cell body. Axon carries information from the cell body to other cells (receiving cell). Dendrites and axons, both extensions of the cell body, are also referred to as processes.

For the typical neuron, a simple model is as follows:


Where N is the neuron; C is the Cell body; D is the Dendrites; A is the Axon; K is the diffusion coefficient. Factors that relate to the K, are the neuron weight, degree of ionization, neuron space configuration, and the condition that whether it combined with the three distinct parts.

Cell body

Cell body assimilates the synaptic input and transmits the determined message to another cell by the axon. It was found to be responsible for the diversity of biochemical process such as transforming glucose into high-energy compounds to other parts of the neuron, highly active proteins serve as chemical messengers between cells are manufactured and packaged and specialized organelles perform the cell’s function. Chen [33] proposed a mathematical theory for hormonal functions in the cell body were as follows:


Where fi is one hormonal function of the cell body, and fi is made up of three main parts:



Dendrites expand its sensitive receptive surface to the surrounding nervous tissue, reflects the function of the cell and the functional properties can be predicted from the pattern of dendritic branching. The thin branching and treelike forms increase the chance for synaptic connections in the brain [34]. Dendrites in the many neurons present with a special form of synaptic connection called dendritic spines. They are small (1-2 µm), a thorn-like protuberance from the dendrite and are the major anatomical feature of neurons in the human nervous system.

According to the definition of dendritic spines, suppose that the dendritic spines are e = x+y. If for any x (0), y (0) satisfy the condition Image

On the basis of adaptive control methods, we can give the following equations and the stability of neural networks with the dendritic structure:


Where Image

If t →∞, then || e ||→0 and system (1) and system (2) achieve modifiable structures indicating the possibility that the pineal gland, a primary source of dendrites.

If we put (3) and system (1) to the system (2), then the following error equations can be obtained between the groups for some fractional differential equations and the central nervous system:


Where Image are the parameter estimation errors.

Next, according to (4), we design the adaptive update law for each parameter estimation error:


Where 0 1,(5,6,7,8,9) < β < = i i and (5) are obtained by eliminating the dendritic potentials from the underlying compartmental model or cable equations [35].


Axon is the excitable membrane that extends to the region of synaptic contact and generates or propagate the action potential. Generally, cells contain one axon but there may be off branches or collaterals to transmit the action potential to the brain. The distinctive length of the axon is the action potential and there occurs a Turing-like instability condition as a precursor for pattern formation in a spatially organized network.

According to Image and (5), we can get the parameters of the adaptive control law:


According to (4) and (5), we get the total error of the system:


Where ImageImage

Then we consider Eq. (7), and expand the formula, we obtain:


Setting P = E9 . Then we obtain the following result:

P = E9

AP+PAT = A+AT = -Q



It is easy to see that Image is a semi-positive definite matrix. Then, the state variable of (7) Imageis asymptotically stable, that is, Image approach zero asymptotically with time. Therefore er we achieve the neuronal functions in the endoplasmic reticulum by a number of the adaptive robust set of fractional differential equations anti-synchronization indicating how the dispersion relation depends on the spatial distribution of the axon-dendritic weights with respect to both network and dendritic coordinates [36].

These primary afferent axons come in different diameters and can be divided into different groups based on their size. Here, in order of decreasing size, are the different nerve fiber groups: A-alpha (13-20 µm), A-beta (6-12 µm), A-delta (1-5 µm) and C-nerve fibers (0.2-1.5 µm). A-alpha, A-beta, and A-delta nerve fibers are insulated with myelin. C-nerve fibers are unmyelinated. The thickness of the nerve fiber is correlated to the speed with which information travels in it-the thicker the nerve fiber, the faster information travels in it.

The model is described by a hyperbolic system of equations


Where Image and 0 <∈=1. Here pi(x,t) is the thickness of the nerve fiber in one of n nerve fiber groups, and x is the sizeof the nerve fiber.



Where λm is determined by the boundary conditions at x = 0 and is a weighted average of the velocities vi (vi can be positive or negative). It is easy to prove that

Image Where Qm(s,t) is the bounded solution [37].

Overall, Cell body, dendrites, and axon are the three main parts of a neuron. We first consider the existence of positive solutions to parabolic nonlinear genetic equations of the form


Where Image The dendrite receives the signal from other neurons; then the signal is computed at the synapse and transmitted to the cell body.

And c* represents the Legendre transform of a function Image


for z∈Rd Here, the bounded domain of Rd is Ω including dendrites, cell body with a nucleus, axon; the outward unit normal to ∂Ω is v . b : R→R is a monotone non-decreasing function; V :Ω→ R is a potential function; Image is a convex function; f is a non-negative real-valued function, and μ0 :Ω→R is a measurable function.

The rest is Image

If the signal into the cell body exceeds the holding threshold, the cell will fire and send the signal down to other neurons through axon [38].

Function and transportation of neurotransmitters

Neurotransmitters are the signaling molecules in neurons that play a vital role in transmitting neural signals through specific receptors, cytomembranes, and postsynaptic membranes. The end of axons of nerve cell secrete neurotransmitters (chemical agents), diffuse and transmit a signal to adjoining cells like muscle cells, neurons, and glands across the synaptic gap by altering its electrical state.

Suppose that the neurotransmitters signals function Image whereImage

If Image

Then the nonlinear genetic equations (9) reduces to:


Where Image represents the transportation of neural signals by the timeImage and position x the summation of extracellular neurotransmitter concentration is:


Through the above equations (10) and (11), we can find out that Methionine Enkephalin (ME), Leucine enkephalin (LE), Dopamine (DA) are able to diffuse freely into both 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine membranes and are guided by the aromatic residues Tyr and Phe. Only a limited number of these neurotransmitters are allowed to penetrate into the membrane, which suggests an intrinsic mechanism by which the membrane is protected from being destroyed by excessive inserted neurotransmitters [39].

Transportation of GABA

GABA, a universal nonprotein amino acid, and functions variedly in different organisms (plants, fungi, and bacteria) and mammalian tissues. It acts as an inhibitory neurotransmitter and helps the neurons to recover from the worry, anxiety, and fretfulness [40].

The corresponding system of equations reduces to


The major predictions from the above equations are as follows: (1) Uptake of GABA is totally sodium‐dependent. (2) Although plots of 1/v versus 1/[Na]2 are nonlinear, the coupling ratio for transport (Na/GABA) is 2. (3) For transport to take place, the order of combination with a carrier must be Na, Na, GABA. (4) Maximal velocity will occur only at infinite Na and GABA concentrations. (5) There is a sigmoidal relationship between apparent maximal velocity (Va) and [Na]. (6) Kt, the [GABA] that gives a velocity equal to Va/2, rises and then falls as [Na] is increased from zero, passing through a maximum at 33.52 mM [Na]. (7) The relationship between initial velocity and [Na] is sigmoidal. (8) Jm, the rate of uptake with infinite [Na], is hyperbolically related to [GABA]; Jm approaches Vmax as [GABA] becomes very large. (9) KNa, the [Na] giving a velocity equal to Jm/2, declines rapidly from 10-7 M to 10-5 M GABA, but is essentially constant at 10-4 M and above. (10) One GABA molecule is translocated per carrier molecule.

Transportation of serotonin

Kogofsky [41] and other contributors cited that serotonin was one of the major neurotransmitters responsible for many biological processes like appetite, mood disorders, sleep, digestion, depression and generalized well-being

The corresponding system of equations reduces to



The above equations of serotonin transporter provide a novel genetic and behavioral primate model to study the molecular, neurodevelopmental, and psychopharmacological mechanisms that underlie genetic variation-associated complex behaviors, with specific implications for the understanding of normal and abnormal serotonin actions and the development of personalized pharmacological treatments for psychiatric disorders.

Transportation of acetylcholine

Acetylcholine was found mostly in neuromuscular junctions and are catalyzed by the acetylcholinesterase enzyme. It is responsible for learning, voluntary movement, sleep and memory and too much of it lead to depression and dementia in case of low level in the hippocampus region [42].

The corresponding system of equations reduces to



It is easy to see that the first intron of the ChAT gene encompasses the open reading frame encoding another protein, Vesicular Acetylcholine Transporter (VAChT), which is responsible for the transportation of acetylcholine from the cytoplasm into the synaptic vesicles.

Transportation of dopamine

Dopamine the inhibitory and excitatory neurotransmitter play the main role in the regulation of reward circuitry and pleasure centers and a dynamic brain chemical for memory and motor skills [43].

The corresponding system of equations reduces to



The major predictions from the above equations involving the effects of Dopamine Transporter (DAT) overexpression in MN-9D cells on the transportation of Dopamine (DA) are as follows: some individuals may be simultaneously more responsive to the effects of environmental adversity and enrichment (i.e., differential susceptibility).

Transportation of epinephrine

Epinephrine, otherwise called as adrenaline, a hormone responsible for its metabolism. It plays a key role in mental focus, attention, arousal, cognition, inhibits insulin excretion and elevates the number of fatty acids in the blood [44].

The corresponding system of equations reduces to



In the above equations, we established a dynamic mathematical model for detection of diabetes in blood with the help of parameters as epinephrine. In addition to this, we also incorporated a new parameter in the existing model i.e. beta cells which has a great impact on the insulin.

Transportation of glutamate

The exciting glutamate neurotransmitter required for memory and learning. Low level of glutamate results in poor brain activity and tiredness and high-level cause death to the neurons in the brain [45]. The corresponding system of equations reduces


The above equations of glutamate transporters tell us that the control of glutamate concentrations is critical to the normal functioning of the central nervous system, and how glutamate transporters regulate glutamate concentrations to maintain dynamic signaling mechanisms between neurons.

Transportation of histamine

Histamine plays a major role in allergic reactions, affect emotions and behavior, control the sleep-wake cycle and promote the release of epinephrine and norepinephrine



We use τ = t + h instead of the above expression to get



Where Image


In silico molecular docking towards neurodegenerative disorder

ADME/Tox filtering rules such as molecular weight, polar surface area, logP or number of rotatable bonds shown in Table 2. Target proteins were retrieved from PDB (Protein Data Bank). The Schrodinger Glide program version 2017 has been used for docking [46-48] shown in Supplementary Table 1. The hydrogen interactions between ligands and target proteins are shown in Figure 1.

S. No. Ligand LogKaHSA Physical-chemical properties Log P
Molecular H2 H2
weight donors acceptors
1 26.27-Di(nor)- 1.118 414.627 1,000 5,100 5.785
cholest-5,7,23-trien- 22-ol, 3-
2 9H –purin-6-amine, N,9- 0.998 367.672 1,000 4,000 5.157
3 Cyanocolchicines -0.576 424.452 0 9,500 1.902
4 3Beta-methoxy-5- 1.414 430.67 1,000 3,700 7.011
cholesten-19-oic acid
  Cholest-5-en-3-ol (3, Beta.)-, 1.809 383.66 1 1.7 6.916
5 carbonochloridate,
6 Cholesterol 1.843 386.66 1 1.7 6.999
7 Cholest-5-en-3-ol (3, 2.371 442.724 0 2 8.418
Beta.)-, propionate
8 Echinacoside -2.248 786.736 12 28.8 -3.648

Table 2: ADME properties of active phytochemical components.


Figure 1: (A) Interaction between Echinacoside and 5-hydroxytryptamine receptor 2A (B) Interaction between Echinacoside and Sodium-dependent serotonin transporter (C) Interaction between Echinacoside and Histamine H2 receptor.

Totally, eight natural compounds with six standard drugs were docked against eleven target proteins, represent as neurotransmitters. The results show among the eight natural compounds, Echinacoside has shown the highest interaction with 5-hydroxytryptamine receptor 2A, Sodium-dependent serotonin transporter, and Histamine H2 receptor. The receptors are mainly involved in neurogenic disorders in human.


Neurotransmitters are molecules that inhibit or stimulate a postsynaptic cell, which is released into the body by the presynaptic nerve cell to produce a response to a certain stimulus. The development of neurotransmitters and its complex functions are influenced by numerous factors. In this study, some mathematical speculations have been proposed on the basis of structural and functional characteristics of the virtual neuron (especially the physiological phenomena of human beings) with a molecular docking and biomathematical approach to formulating some speculations to the consolidation of the identification of neurotransmitters function. This could pave a way to formulate more mathematical speculations related to the neuron, and finally, these data and approaches will be useful for constructing virtual neuron with the help of biomathematics. The interaction between natural compounds and neurotransmitter studies shows the good interaction with all the compounds. Especially, 26.27-Di (nor)-cholest-5,7,23-trien-22-ol, 3-methoxymethoxy, Cholest-5-en-3-ol (3, Beta.)-, carbonochloridate, Cholesterol and Echinacoside exhibited maximum interaction with all the target proteins. Among the other compounds, echinacoside shows highest interaction with (Serotonin) 5-hydroxytryptamine receptor 2A (-17.077), Sodium-dependent serotonin transporter (-15.810) and (Histamine) Histamine H2 receptor (-17.556).

Serotonin the other major inhibitory neurotransmitter is deemed to be the master neurotransmitter. The imbalance is one of the most often cited contributors to depression and other mood disorders. It is also intimately tied to many biological processes such as sleep, appetite, pain, digestion, and generalized well-being.

Histamine is most commonly known for its role in allergic reactions but it is also involved in neurotransmission and can affect your emotions and behavior as well. Histamine helps control the sleep-wake cycle and promotes the release of epinephrine and norepinephrine.

Serotonin and histamine were the brain monoamines which play a vital role in cognition, emotions, pathophysiology, and treatment of mental disorders. In the current study revealed that the neurotransmitters structure and transportation by mathematical models and in silico molecular docking results strongly shows the Echinacoside is a potent inhibitor in some neurological disorder associated with serotonin and histamine. Further, extend methods adapt to study the mechanism and pathway level interactions between the natural compounds with the disease.

Conflict of Interest

We have no conflict of interests to disclose and the manuscript has been read and approved by all named authors.


This work was supported by the Philosophical and Social Sciences Research Project of Hubei Education Department (19Y049), and the Staring Research Foundation for the Ph.D. of Hubei University of Technology (BSQD2019054), Hubei Province, China.


Citation: Zhao B, Jiang X (2022) Biostatistical and Mathematical Analysis on Liver Disease during COVID-19 Pandemic. Biochem Anal Biochem. 11:465.

Copyright: © 2022 Zhao B, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.