In this communication, we represent a self-consistent systematic optimization procedure for the development of optimally tuned (OT) range-separated hybrid (RSH) functionals from first principles. This is an offshoot of our recent work, which employed a purely numerical approach for efficient computation of exact exchange contribution in the conventional global hybrid functionals through a range-separated (RS) technique. We make use of the size-dependency based ansatz i.e., RS parameter, γ, is a functional of density, ρ(r), of which not much is known. To be consistent with this ansatz, a novel procedure is presented that relates the characteristic length of a given system (where ρ(r) exponentially decays to zero) with γ self-consistently via a simple mathematical constraint. In practice, γ OT is obtained through an optimization of total energy as follows: γ OT ≡ opt γ E tot,γ . It is found that the parameter γ OT , estimated as above can show better performance in predicting properties (especially from frontier orbital energies) than conventional respective RSH functionals, of a given system. We have examined the nature of highest fractionally occupied orbital from exact piece-wise linearity behavior, which reveals that this approach is sufficient to maintain this condition. A careful statistical analysis then illustrates the viability and suitability of the current approach. All the calculations are done in a Cartesian-grid based pseudopotential (G)KS-DFT framework.