1. Introduction
Designing complex fluorophore-conjugated metal nanostructure (FMNS) constructs is an important component in the world of bio-imaging [
1]. These constructs allow us to use the plasmonic qualities of metal nanostructures along with the qualities of fluorophores (such as radiative fluorescence lifetime (FLT) and fluorescence intensity (FI)), and allow performing multi-modal bio-imaging, such as FLT imaging microscopy (FLIM) [
2,
3], diffusion reflection (DR) [
4,
5], and X-ray computerized tomography (CT) scans [
6,
7,
8], using the same contrast agent. The proximity between plasmonic nanostructures and dye molecules creates special nearfield effects that change the fluorescence intensity, lifetime, and quantum yield of the dye molecule [
9], called metal-enhanced fluorescence (MEF) and quenching. The radiative FLT and quantum yield (QY) are affected by the distance and orientation with respect to the nanostructure (NS) [
10]. Therefore, optimizing the geometry of the FMNS will improve the effectiveness of the MEF mechanism, resulting in better components for the measurements mentioned above.
Genetic algorithms (GAs) are an optimization method based on a natural selection process that acts like biological evolution [
11]. The algorithm is iterative and modifies a population of individual solutions for the optimization problem, allowing the evaluation of a vast number of possibilities. We propose a method to design FMNS constructs using a GA. GAs are being used in other fields of engineering such as radio frequency (RF) antennas, protein structure engineering, solar cells, and more [
12,
13,
14]. Without a computed search tool, such as a GA, it would be extremely hard to examine most of the possibilities for orientation, distance, and NSs shape.
2. Genetic Algorithm
A GA is an iterative method for finding a global extremum to a problem [
11,
15,
16]. It originated in the 1970s like a search algorithm and is lately being used as an optimization method in different fields of engineering like RF, computational biology, plasmonics, and more [
12,
14]. A GA models the natural process of “survival of the fittest,” starting with an initial “population” and using evaluation, death, reproduction, and mutation sub-processes until it converges to the “best valued” population. The (n + 1)th generation is composed of the fitter population members of the nth generation and their “children” (created from combinations of qualities from both parents and a random change that models mutation).
Each generation consists of a population of M members. In nature, a chromosome is a compact package of genes (actually it is a compact state of DNA molecules, containing the genes). An organism in nature has a constant number of chromosomes, typical for its species. In our work, a nanostructure will be modeled by a single array of parameters, analog to a specimen with one chromosome. Each specimen (population member) has a set of parameters that determine its fitness (the parameters are referred to as “genes”). Every member of the population is evaluated according to a value that is set according to the member’s “genes” through a fitness function. After the evaluation process, the top-valued members go through a “reproduction” process, creating the next generation [
13], which involves adding random changes in some of the new members (mutations) and creates the next generation. The GA then selects and builds the next generation based on the high ranked members.
In our realization of the algorithm, specimens with a better fitness value have a greater chance of reproduction, while the diversity of the new population is maintained by the crossover operator (the method of gene exchange between pairs of specimens in the reproduction stage) and the mutation operator, which creates a new specimen by changing an existing one. The mathematical choice of these operators is the heart of the genetic algorithm implementation. There is more than one way of representing the population members and choosing the different GA operators. Yet, two main principles remain in each of GA’s representations: maintaining the population’s diversity in order to keep a wide solution range and promoting “fitter” solutions (specimens).
As mentioned earlier, the crossover operator has a profound influence over the efficiency of the optimization using a GA since it may create solutions of the search space that are not accessible to either of the two “parents.” This creation of new solutions that allow for a combination of traits from each “parent” further optimizes each new population, making a GA an efficient search algorithm, mainly because of its traits as a cooperative process (the “interaction” between two different specimens) [
17]. MATLAB software (1994-2019 The MathWorks, Inc., Natick, MA, USA) offers a solid genetic algorithm realization through its toolboxes.
3. Physical Model
The theory for energy transfer from fluorophores to plasmonic nanostructures is moderately complex, and similar equations have been derived from classical and quantum mechanical considerations. We will describe only the final equations. Readers interested in the physical basis of energy transfer from fluorophores to plasmonic nanostructures are referred to the excellent review by Clegg [
18].
Mathematically, the link between the geometrical parameters (coded as the genome) and physical properties that indicate the quality of fluorescence are the fluorescence lifetime (FLT), enhancement rate, and quenching rate. The model is linked to the GA through the cost parameter. The cost parameter is calculated in a way that takes into account the fluorescence parameters.
The optimization process uses calculations of the scattering cross-section for the nanorod and absorption cross-section for each fluorophore that surrounds the rod. In the heart of the optimization process stands the evaluation function, which checks the value of possible solutions. Every solution (e.g., specimen) is composed of a variation of orientation angles and distance between the fluorophore and the gold nanorod (GNR), as shown in
Figure 1, where both shapes were approximated as oval to simplify the amount of information used by the metallic nanoparticle boundary element method (MNPBEM) solver throughout the GA’s iterations. This will affect the specimen’s orientation factor
(also known as the orientation factor):
where
is the angle between the fluorophore and the plasmon dipole moments;
is the angle between the line that connects the fluorophore’s and the GNR’s centers and the plasmon dipole;
is the angle between the line that connects the fluorophore’s and the GNR’s centers and the fluorophore dipole [
19].
In the calculations mentioned above, we assumed that the distance is independent of the orientation (zero order dependence approximation).
First, we calculate the quenching and enhancement as if the fluorophore is right above the upper edge of the GNR with all angles equal to 0. Then, we multiply our result with an orientation factor: for quenching and factor for enhancement.
Note that the multiplication by the cosine-squared factor is because an antenna or a dipole radiation field is maximal along the direction of the antenna (or dipole) and it is reduced as a cosine of the orientation angle such that the radiation field is zero in perpendicular to the antenna (or dipole). For intensity, the factor is cosine squared. The rate of transfer from fluorophores to plasmonic nanostructures separated by a distance
r is dependent on
κ squared. The
κ-squared term is a factor describing the relative orientation in space of the transition dipoles of the fluorophore and the nanostructures. More about the orientation factor
κ can be find in Haas et al [
19]. There, one may find the full description of the effect of the orientation of the donor and acceptor on the probability of energy transfer involving electronic transitions of mixed polarizations.
The evaluation function computes the radiative enhancement through an expression, which takes into account the quenching and the MEF effects. The quenching effect is the transfer of energy from the fluorophore to the quencher (in this case, the metallic nanostructure’s plasmons), which causes a decrease in fluorescence. In our case, the fluorophores are fixed around the GNR so that the quenching mechanism should describe a static and distant interaction. A good model for that kind of interaction is the use of a fluorescence-resonance energy transfer (FRET), in which the donor is the fluorophore and the acceptor is the GNR [
2], where the quenching energy loss rate is described by Equation (2) and the enhancement energy transfer rate is described by Equation (3). The energy transfer is influenced by the spectral overlap between the GNR absorption and the fluorophore’s emission. This overlap is quantified using an overlap integral that takes into account a molar absorption of the GNR (GNR’s cross section normalized by the total energy, concentration, and cuvette size), and the emission line shape of the fluorophore, described by Equation (2):
where
QR is the quenching rate,
is the nanorod (NR) concentration in the solution (assumed to be = 1 nM), and
is the container (cuvette) size, which equals 1 cm.
is the overlapping integral that represents a normalized evaluation standard for the spectral overlap between the lineshapes of the GNR’s absorption plasmon resonance (
) and the fluorophore’s photon emission (
).
Radiative enhancement was considered using an equation that describes the enhancement of the fluorophore near a gold nanorod (modeled by an ellipsoid) [
20]. The GNR material, height, and width affect its absorption and scattering spectra and change the overlap integral. We assumed that the GNR was excited by an electromagnetic field that corresponds to its longitudinal plasmon resonance [
21]. We multiplied the normalized enhancement expression by 1/τ
o to get the enhancement rate (ER):
After calculating the enhancement rate and the quenching rate of a conjugated structure (a “specimen” in a generation of the GA), we examined two types of optimization processes, where a minimum cost led to a “better” specimen. The first was optimizing the differential enhancement factor so that high is desirable and the GA’s cost parameter The second process was based on the assumption that when enhancement conditions are met, the FLT will become equal to (the natural lifetime) such that maximal FLT is an indication of maximal enhancement. This is an indirect optimization since we use FLT as a way to optimize the enhancement effect. According to the latter, the cost , with being the lifetime.
The differential enhancement factor (DEF) is set by the difference between the enhancement rate and quenching rate, as described in Equation (4):
Enhancement efficiency (TrE) and quenching efficiency (TrQ) were determined in order to quantify the effects of quenching and metal enhancement (as shown in Equations (5) and (6)).
The FLT is modeled as a double exponential, while the relative effect of the metal interaction lifetime is a function of the efficiencies and the natural lifetime . Our tool evaluated the level of spectral overlap. We experimented on different nanostructures and different dye molecules and examined each of its possible combinations.
A scheme for our FLT optimization process is presented in
Figure 2.
If the level of overlap is larger or equal to 40%, we presumed that metal enhancement is the dominant effect. We used different expressions for in order to have a function that fits the expected behavior. If the level of overlap was smaller than 40%, we presumed that quenching was the dominant effect and the weight of the quenching related element was a function of the quenching efficiency (TrQ) only. However, note that the base assumption that quenching is dominant due to spectral overlap being low is not accurate.
4. Validation Procedure
Each generation’s population can be described as a “gene-pool,” which is an array containing all of the current generation’s “genetic” information. Population size is the size of the gene-pool vector, given by Equation (7):
A data structure representing a fluorophore–GNR (FGNR) structure contains information about the number of fluorophores, their distance from the GNR, and their angular position with respect to the GNR long axis (e.g., one specimen). Each specimen gets evaluated according to its fluorescence level or QY. The specimens that are evaluated with the highest “grades” are selected as parents and their “genetic data” is transferred to their offspring. The procedure repeats itself until it reaches the highest radiative emission coefficient (
or maximum number of generations. This process is illustrated by
Figure 3.
The GA optimization was done using a MATLAB solver that receives an evaluation (target) function that is created by the user. The structure of the genome, as well as other optimization parameters, are determined by the main script before the optimization function activates the GA solver. The optimization function also calculates the magnitude of spectral overlap, depending on the nanostructure’s shape and material, and the scattering and absorption spectra of the dye molecule. An indication whether the spectral overlap is above the threshold is then passed to the evaluation function for further calculations of fluorescence parameters and for the estimation of the “cost” parameter, which is needed for the evaluation of each specimen (conjugated structure).
The evaluation function is the link between the physical model’s equations and the GA. The physical equations provide the function quantitative data for parameters like enhancement rate, FLT, and Forster distance. After the calculation of several fluorescence parameters, the evaluation function calculates and returns the GA’s cost parameter for a given specimen. The genome of each specimen is transferred to the evaluation function for the calculation of the cost parameter for this specimen, allowing the GA to evaluate its fitness.
Genes 1–5 are analogous to a genotype, where each of them has a role in expressing the physical separation distance. The separation distance is a geometric average of the genes, which will be used after its calculation in our fluorescence equations. Genes 6–8 contain information regarding the orientation of the fluorophore’s absorption dipole (as seen in
Figure 4). The angles obtained from genes 6–8 are used for the calculation of the orientation factor
.
The cost parameter is calculated differently according to the fluorescence property we wish to optimize. We tested two variations, optimizing , via the cost function of , and optimizing FLT, via the cost function of .
After an optimal conjugated structure is chosen, our tool copies its genetic data and modifies only the genes related to the separation distance in order to view the fluorescence parameters of every possible solution within a given range of distances.
MNPBEM was used as our spectral analysis tool. We compared our results to spectrums of the same shapes, as calculated in El Sayed’s work [
22]. Spectral analysis is needed for the calculation of the spectral overlap ratio and integral, which are important for the optimization process. The comparison between a spectrum of a metal nanostructure and the absorption spectrum of fluorescein can be seen in
Figure 5.
6. Conclusions and Discussions
The evaluation of field intensity over wavelength has shown good correlation with the calculation of lateral intensity dispersion as seen by the brightest intensity graphs matching the peak of the
vs. wavelength graph. We inspected some lateral simulations of
. These simulations were performed at a different wavelength in every simulation, some with excitation polarized longitudinally (a–c) and some with transverse polarization (d–f). The highest
values matched the resonance seen in the spectral simulations above. We can observe these results in
Figure 9.
The optimization processes we initiated has shown us that our tool provides the expected results optimization-wise and it will surely be beneficial for the future design of contrast agents. However, the quality of its results depends greatly on the validity of the physical model that is being used in its evaluation function. In our FLT optimization, we used FLT as an indicator for optimal MEF interactions and we conclude that a model in which the different energy transfer rates are calculated directly, without the use of circumvention through optimization of FLT, may provide valid results. Nevertheless, our optimization tool has provided a proof of concept for the benefit of machine learning algorithms and may become much more useful with more accurate physical equations. Under the assumptions we made, for a specific FLT and enhancement factor models and assuming that the spectral overlap is the main cause of whether MEF is dominant or quenching is dominant, we could deduce that a silver nanorod would be a good conjugate to a green dye molecule, gold nanostructures would be better for infra-red dye molecules, and silver NPs would work best with purple and UV molecules.
The physical model that was used does not express the physics accurately enough, thus affecting the quality of our optimization. A more accurate model needs to be introduced in order to allow our optimization tool to provide us with an engineered solution. Choosing FLT as an indicator of maximal enhancement did not provide us with the information we needed. Using the enhancement rate directly will most likely provide better solutions for optimized intensity-based contrast agents.
As GA optimization effectiveness depends greatly on the structure of the genome, population size, mutation rate, and more, these parameters need to be fine-tuned in order to optimize the process effectively. GAs will most likely provide an acceptable solution, though not necessarily the superior one. GA-based tools, like the one we have created, can be fast, simple, and best suited for multi-variables optimization problems.