The ratio of the proton to the electron mass was once noticed  as magical
The 1951 ratio can be further adjusted to the most accurate current value by the simple tuning multiplicative factor i.e. 36821-st root of 2
We find here a simple precise analytical formula for the gravitational constant G relating it to the electron charge, electron mass, the vacuum dielectric constant and the speed of light (or magnetic permeability of the vacuum) in power of the fine structure constant as it was only a constant of the Quantum Electrodynamics (QED) ,i.e. relating the gravitational constant to the Planck constant through others which also well exist without the quantum mechanics therefore relating two fundamental constants as not independent through the parameters of the electron and the electromagnetic properties of the vacuum. We asked for the following conditions: 1) The gravitational constant is proportional to the natural strong gravity , gravitational constant for the electron which already can be constructed from the dimensional analysis of units, simply
2) It is proportional as the QED perturbation expansion dimensionless coupling constant i.e. in power of the fine structure constant so it formally looks like the n-th order perturbation correction to the electron-electron Coulomb interaction through the Yukawa exchange mechanism  of the logitudinal photons but is insensitive to the sign of charge. 3) The front rest coefficient is around 1. 4) In similarity to the proton electron mass ratio, it can be expressed as a simple combination of the integer and the irrational numbers and 5) The analytical formula is perfectly matching the best known experimental CODATA value.
2. The Preliminary Newtonian Value
The preliminary formula yields
i.e. gravitational constant is proportional in (4/3) of the 20-th power of the fine structure constant , to the square of the ratio between the electron charge e and its mass divided by times the vacuum dielectric constant or in 21-th power of the fine structure constant to the strong gravitational constant defined here as
which shows immediately the electron-electron gravitational interaction as an effective ultra-small electromagnetic (electrostatic) interaction of some order in and that the gravitational energy corrections to the electron-electron electromagnetic interactions are in the 20-th power of the fine structure constant. The value is that obtained by Newton  by the rough estimate of the Earth density as the 6 times that of water. This is in agreement with  predicting the approximate (when the complex i, is assumed approximately 1 within the complex numbers neighborhood circle is approximately 101 in that sense) “Rydberg” spectrum of the elementary particles with masses from the formula (34) in  as
i.e. electron for
as a Big Bang ultra-high random “Rydberg toss” excitation which is approximately assuming that the elementary particle is a discharged Black Hole which self binds itself orbiting itself on a circular orbit with the Schwarzschild radius without the relativistic mass gain with the speed of light on its own mass with the fractional closed multiple-spiral near-circular orbit “hydrino”  Bohr quantizationi.e. is a fractional Planck particle when only the integer fraction of its Compton wavelength is its Schwarzschild circumference ( times the Schwarzschield radius) of its massi.e.
It also agrees up to the small multiplicative factor with self-consistently solving the exact nonrelativistic Schrödinger equation for the Black Hole to obtain its possible quantum mass spectrum  treated as the self-bound Hydrogen atom while equalizing its Rydberg energy to the rest relativistic energy of the minus half of its effective mass but considering also the “hydrino” fractional quantum numbers i.e. .
The leading order of the Lamb shift  is in the 6-th power of in that sense. The formula immediately predicts the Positronium ground state  gravitational shift (and so the between 1 - 2 levels gap gravitational “Lamb shift” as −3/4 of it) by the effective attractive Coulomb coupling renormalization entering the Rydberg energy in square in the leading order as
and so of that for the Hydrogen and twice of that for the Deuterium ground state respectively. Assuming the somehow perfectly Mössbauer-rigid nucleus in both cases not to change the Kepler reduced mass and eliminate the native isotope effect it leads mathematically to the Trojan or anti-Trojan Wave Packet generation  from the same circular quantum number state in atom Rydberg state of two different isotopes due to the electron-nucleus gravity as big as of 1068 meters radius (1042 times the size of the Universe)i.e. with
. It also takes about 1027 seconds (109 ages of the Universe) (and so only 10−27 of it within one second) to cause one interference fringe shift assuming two counter rotating Trojan Wave Packets around the principal quantum number which belong to two Hydrogen atom of two different isotopes can be interfered between each other with the method of delta pulse ionization described in .
3. The Refined Value
The simple correcting prefactor (inverse of the 38-th root of 2) is further adjusting the value to the best known experimental  CODATA value giving
It may be interpreted as a ratio factor between the gravitational self energy radius and the electromagnetic electron radius when the gravitational self energy deficiency in modulus is assumed to be the quantum fraction of the electrostatic one i.e.
i.e. . This is in inequality agreement with the classical model of the electron with the spin assuming slightly different rotating Gaussian mass and charge distributions inside the electron
predicting the anomalous electrong-factor 2 if only
as it was built from the solid electrified matter but with the locally variable density and non-proportional variable level of electrification. The refined “Rydberg” quantum number of the electron as the excitation itself including the correction is now (Table 1).
We so find (Figure 1)
where G is the gravitational constant,
Figure 1. The picture of the discovered “Theory of Everything” formula expressing the gravitational constant G in terms of the other fundamental constants and the electron parameters as the elementary particle (charge and mass).
Table 1. The table showing the progressive refinement of the value of the strong gravity constant Gs and the gravitation constant G in terms of other fundamental constants.
is the reduced Planck constant, c is the speed of light, is the electron mass and is the fine structure constant i.e. with the perfect agreement with the most recent CODATA  value only if the right side quantities are also the CODATA values. This expresses explicitly that the gravitational interaction between the electrons may be considered as an exotic ultra-weak Coulomb attractive electromagnetic (electrostatic) interaction of the 20-th order in the fine structure constant and therefore is normally neglected in QED as much below the highest 7-th order in normally considered there. It also defines electron strong gravity constant as
The further ultra-fine tuning prefactor placing the value perfectly at the center of the CODATA error field can be found as (hundred sixty thousandth root of 2) i.e.
with the method of power 2 approximants of the number which is close to 1 as for example for
etc. for which the ultra-refined
We find here the “Theory of Everything” formula relating the gravitational constant to the Planck constant through the electron mass, speed of light and the fine structure constant (electric charge and the dielectric permittivity of the vacuum) suggesting that the electron quantum mechanics is a result of the gravitational and dielectric properties of the vacuum or alternatively the gravity of electrons is an ultra-weak QED effect of longitudinal virtual quantum photons interacting with the 1-st scalar component of the electron 4-current and the Yukawa-like interaction between charges, which is insensitive to the sign of the charge. The 4/3 and 1/21/38 (1 over 38-th root of 2) coefficients adjusting the 21-st power order of the so defined strong gravity constant to the first Newtonian 15-teen century value and than the perfect CODATA value are easy to remember as having the 4/3 factor of the sphere volume and 1/38 in the 2 exponent as the calculator LCD display altered value of the approximate fine structure constant 1/137. Yet another interpretation is that the parameters of the electron as the elementary particle like the mass m and charge e are fully determined by the properties of the vacuum and the quantum mechanics i.e. the value of by some Planck particle like mechanism of de Broglie standing wave Black Hole self no-escape for example with fractional angular momentum. Perfect CODATA values of all involved quantities must be used for the G CODATA value agreement.
Author would like to thank W. Tarkowski, the author of  for the valuable discussion regarding his paper used for the partial explanation of the discovered formula formulating the model of the elementary Planck-like particle as the self-bound Black-Hole in the extended 5-dimensional space-time.
 Tarkowski, W. (2004) A Toy Model of the Five-Dimensional Universe with the Cosmological Constant. International Journal of Modern Physics A, 19, 5151-5084.
 Kalinski, M. (1998) Aharonov-Bohm Oscillations in a Hydrogen Atom in a Radiation Field through Electron Self-Interference. Physical Review A, 57, 2239-2242.