Team:Tec-Chihuahua/Model

Erwinions








From the very moment we made this our project, we were repeatedly asked the same question: “how can you make sure enough AMPs get to the bee larvae?” This is a valid and crucial concern. After all, one of the key elements to our project’s success is that the antimicrobial peptides reach the bee larvae with a high enough peptide concentration to inhibit pathogenic bacteria. The method we chose was PLGA nanoencapsulation.

Therefore, in order to measure the final amount of AMPs that will be delivered to each larva, we model the amount of peptide released from the nanoparticles in their path through the insides of the nurse bee and into the jelly fed to larvae. The model we developed will help us to determine the initial dosage of AMPs necessary to reach their final destination in desirable concentrations.


Theoretical Foundation

Given previous research that confirms antibiotics supplemented to nurses in food get to the brood1, the nanoparticles must travel through the nurse bee body along with the artificial food. Therefore, is important to examine the alimentary canal as well as the hemolymph system, since this is the path we postulate the nanoparticles will travel.

The nanoparticles must travel through the nurse bee body firstly along with the artificial food, therefore, is important to examine the alimentary canal, the blood and the organs of circulation of the bee since these are the main systems in the pathway from PLGA nanoparticles to the larvae.2

The cellular layer segregates juices and enzymes and the resulting products of digestion will be absorbed by the peritrophic membrane and travel through the stomach walls to the hemolymph. Therefore nanoparticles will go to the blood or hemolymph system. Once they are in the hemolymph they will circulate through the body by the single blood vessel which is a long slender tube that extends from the back abdomen through thorax and ends at the head.2

The abdomen segment of the vessel tube is the heart and by rhythmic pulsations, it makes the blood circulate. The blood with the nanoparticles will travel to the aorta and will irrigate all the structures of the head, including the hypopharyngeal and mandibular glands. The hypopharyngeal synthesize the royal jelly that will mix with the secretions of the mandibular glands and the PLGA nanoparticles. The mix will then be delivered to the larvae through the excretory ducts that are the normal pathway for which the royal jelly is delivered to the larvae.2


The polylactic-co-glycolic acid (PLGA) is a biodegradable and biocompatible polymer that is used for the fabrication of drug delivery devices. PLGA could be used to transport proteins, peptides and macromolecules such as DNA and RNA. PLGA is ideal for drug delivery because its degradation can be used to control the drug’s release profile.

While zero order kinetics would be the most desirable case, this is currently unachivable. There are three possible non-zero ordered behaviours for polymer degradation. Polymer degradation/erosion (A) is dependent on the degradation of the polymer coating. Diffusion (B) is currently the most common release mechanism, wherein release is dependent on the dissolved drug’s concentration. Finally, solvent penetration or device swelling (C) is dependent on the permeability of the device.

The release behavior of PLGA is linked to the degradation which occurs by hydrolysis. The degradation profile is through bulk degradation as a consequence of the water penetration into the particle and by the polymer degradation. The process can be divided into two phases:

  1. The initial burst where there is a quantity of drug released at time zero related to the drug type, drug concentration, and polymer hydrophobicity.
  2. The progressive drug liberation through the weakened PLGA layer. The solvent hydrolyzes the polymer and creates a channel for where the drug can be released by diffusion and erosion until polymer degradation is complete.

The factors that influence the PLGA degradation are the following:3

  1. Polymer composition determines the hydrophilicity and the rate of degradation of the particles. A high percentage of glycolic acid in the composition increases weight loss.
  2. Polymer coating, which can result in a different reaction, e.g. Chitosan coating makes for sturdier nanoparticles that begin degradation in more acidic environments.
  3. The average molecular weight is related to the polymer chain size. Heavier molecular weights correspond to lower degradation rates.
  4. The drug or molecule type in the PLGA matrix modifies the degradation mechanism and the rate of degradation. It is clear that the chemical properties of the drug affect the release profile. However, there is not a clear correlation between the drug chemistry and the behavior.
  5. The size and shape also impacts matrix degradation. A higher surface area in comparison to volume shows higher matrix degradation making drug liberation faster.
  6. Environment pH may also affect the rate degradation of PLGA. Strongly alkaline and acidic media leads to a faster degradation whilst less acidic and neutral media exhibits slower acceleration on the rate degradation.
  7. Also, another important parameter is the amount of drug load. PLGA loaded with high drug concentrations increase the initial burst release.
  8. Finally, in biological systems, the degradation may be influenced by enzymatic activity. However, there are not conclusive results in how the enzymes play a role in the degradation process.

Modeling

The model is based on the diffusion from the antimicrobial peptide in the interior of the PLGA nanoparticle to the exterior of it. Since our peptides will be delivered through PLGA 50:50 nanoencapsulation, we decided to use a simple Fickian diffusion model suggested by previous research4, 5, 6.
The diffusion equation can be seen here:

Where:

Symbol Variable Unit
Mt/M∞ Drug released in time t in relation to the total amount of nanoencapsulated drug Unitless, from 0 to 1
T Time Seconds
D Diffusion coefficient, dependent on the drug-polymer interations Centimeter squared per second
R Radius Centimeters
⍺ The initial burst, which is a normal behavior of PLGA-drug release systems where quantity of initial drug is released independently of time. A fraction, in relation to the total amount of drug


We adapted this equation to our needs, adding a new variable ẟ that represents the remainder of dosage still within the nanocapsule such that ẟ at a given time is equal to the subtraction of the amount of peptide released at that time from the full load.

Since our capsules will travel through two stages, inside an adult bee and within royal jelly, we also established a delivery time (td) and expressed the diffusion equations in relation to this travel path.

Time (s) Diffusion coefficient (cm2/s) Radius (cm) Significant upper bound (n at which the summation term is smaller than 0.00001)
t=3,024,000s D=5x10-19 R=2.2x10-5 1
R=0.00013 4
D=4x10-16 R=2.2x10-5 15
R=0.00013 62
t=0s Any value 317


Experiment

The limitation that calculating D provides is clear. This coefficient is an interaction between the drug and its polymer coating, along with the medium. To overcome this, we chose to experimentally evaluate the drug’s release mechanism and pick the best-fit line. We decided upon apidaecin as the model drug. We did this under two conditions: one simulating hemolymph and another for jelly. For this, we created two buffer solution at these pH ranges, and kept them at the desired temperature.


Conditions Simulated pH Temperature
Bee Hemolymph 8.35,6 32 - 35 ºC 9
Worker/Royal Jelly 3.5 - 4.57, 3.98 32 - 35 ºC 9, 30 - 33 ºC 10

With help from our friends at Tec-Monterrey, we had apidaecin nanoencapsulated in PLGA 50:50 (51 kDa) coated with Chitosan and sent back to us, where we could periodically test the amount released over a period of time in the buffer solution.


The best fit line for each case let’s us figure out the diffusion coefficient for these conditions.


Condition PLGA - Apidaecin diffusion coefficient
Bee Hemolymph X
Worker/Royal Jelly X

Limitations to our approach

There are several challenges regarding this experimental approach. Here we present them in an upfront manner, and explain what is to be done:

  1. Since the initial burst is also determined by drug concentration, once we find the diffusion coefficient, this will not be enough to successfully predict another concentration’s behavior but rather provide an approximation. Further experimentation is required.
  2. This model for PLGA diffusion has been shown to lose prediction accuracy beginning on the 10th day. This limitation is insignificant to the scope of our project; as shown by our collaboration with XXXXX, bee metabolism takes only X days and larvae are also most susceptible 12-36 hours after hatching. Not only this, but the brood need only 9-11 days to have their cells capped and start becoming pupae, which rarely die from Foulbrood (since it can no longer infect them).13
  3. The viscosity of jelly might result in a reduced diffusion capacity, which must be taken into account for the amount of remaining peptide.
  4. Bees have other pH values in their guts, head, and glands. However, these all fall in between the pH values chosen for hemolymph and royal jelly.14
  5. We underwent experimentation only in apidaecin. We’d need to repeat this process for the remaining AMPs or their combination.

Results & Discussion

Matlab Coding

Here we present the codes used for our modeling, to be used as confirmation or reference (values provided are examples):


clear all;
close all;
clc;
D=5*10^-19;
t=0;
R=2.2*10^-5;
x=exp(-(D*t*(pi^2))/(R^2));
z=0;
for n=1:1:inf

zi=(x^(n^2))/(n^2);
if zi<0.00001
disp(n);
break;
end;
z=z+zi

end;






clear all;
clc;
close all;
R=4.5*10^-5; %in cm
D=3*10^-18; %in cm2 / s
t=0:3600:3456000; %in s
Z=0;
a=0.3;
for n=1:1:317

Z=Z+(1/(n^2))*exp(-((n^2)*(pi^2)*D.*t/(R^2)))

end;
y=1-(6/(pi^2)).*Z+a; %fraction
plot(t,y);
xticks([432000 864000 1.296*10^6 1.728*10^6 2.16*10^6 2.592*10^6 3.024*10^6 3.456*10^6]);
xticklabels({'5','10','15','20','25','30','35','40'});
xlabel('Time (days)');
ylabel('Fraction released');




  1. Elzen, P. J., Westervelt, D., Causey, D., Ellis, J., Hepburn, H. R., & Neumann, P. (2002). Method of Application of Tylosin, an Antibiotic for American Foulbrood Control, with Effects on Small Hive Beetle (Coleoptera: Nitidulidae) Populations. Journal of Economic Entomology, 95(6), 1119–1122. doi:10.1603/0022-0493-95.6.1119
  2. Grout, R. (1949). The Hive and the Honey Bee. Hamilton, IL: Dandant & Sons. p.504-509.
  3. Makadia, H. K., & Siegel, S. J. (2011). Poly Lactic-co-Glycolic Acid (PLGA) as Biodegradable Controlled Drug Delivery Carrier. Polymers, 3(3), 1377–1397.https://doi.org/10.3390/polym3031377
  4. Budhian, A., Siegel, S. J., & Winey, K. I. (2008). Controlling the in vitro release profiles for a system of haloperidol-loaded PLGA nanoparticles. International Journal of Pharmaceutics, 346(1–2), 151–159. https://doi.org/10.1016/j.ijpharm.2007.06.011
  5. Siepmann, J., Elkharraz, K., Siepmann, F., & Klose, D. (2005). How Autocatalysis Accelerates Drug Release from PLGA-Based Microparticles: A Quantitative Treatment. Biomacromolecules, 6(4), 2312–2319. doi:10.1021/bm050228k
  6. Ford Versypt, A. N., Pack, D. W., & Braatz, R. D. (2013). Mathematical modeling of drug delivery from autocatalytically degradable PLGA microspheres — A review. Journal of Controlled Release, 165(1), 29–37. https://doi.org/10.1016/j.jconrel.2012.10.015
  7. Rojas-López, C. A., Padilla-Salas, V. R., López-Cortés, G. I., Ordaz-Briseño, R., & Islas-Rodríguez, A. E. (2017). Sistema inmune innato de las abejas. Revista De Ciencias Ambientales Y Recursos Naturales, 3(7), 25-32.
  8. Nation, J. L., Sr. (2015). 15.4.5 Hemolymph pH and hemplymph buffers. In Insect Physiology and Biochemistry (3rd ed., pp. 424-425). CRC Press
  9. Chen, C. (1995). Changes in protein components and storage stability of Royal Jelly under various conditions. Food Chemistry, 54(2), 195–200. https://doi.org/10.1016/0308-8146(95)00031-d
  10. Sanhueza-Rojas, O.H. (2016). Caracterización química multifactorial de miel en relación a la infección por nosema ceranae en abejas, actividad antimicrobiana y origen geográfico del producto. Thesis. Universidad de Chile.
  11. Tautz J (2008) The Buzz about Bees. Springer-Verlag, Berlin Heidelberg.
  12. Heinrich, B. (1979). Bumblebee economics. Cambridge: Harvard University Press
  13. Genersch, E. (2010). American Foulbrood in honeybees and its causative agent, Paenibacillus larvae. Journal of Invertebrate Pathology, 103, S10–S19. doi:10.1016/j.jip.2009.06.015
  14. Herbert, E. W., Jr., & Shimanuki, H. (1983). Effect of diet pH on the consumption, brood rearing, and pH of worker jelly produced by caged honey bees. Apidologie, 14(3), 191-196. Retrieved from https://hal.archives-ouvertes.fr/hal-00890597.



Follow us on
social media:


Instagram Facebook Twitter Snapchat Gmail

We thank our
sponsors:


Thermo Neb IDT Snap QTEc ITESM