Liquid crystalline networks (LCNs) are stimuli-responsive materials formed from polymeric chains cross-linked with rod-like mesogenic segments, which, in the nematic phase, align along a non-polar director. A key characteristic of these nematic systems is the existence of singularities in the director field, known as topological defects or disclinations, and classified by their topological charge. In this talk, we address the open question of modeling mathematically the coupling between mesogens disclination and polymeric network by providing a mathematical framework describing the out-of-plane shape changes of initially flat LCN sheets containing a central topological defect. Adopting a variational approach, we define an energy associated with the deformations consisting of two contributions: an elastic energy term accounting for spatial director variations, and a strain-energy function describing the elastic response of the polymer network. The interplay between nematic elasticity, which seeks to minimize distortions in the director field, variations in the degree of order, with the consequent tendency of monomers in the polymer chains to distribute anisotropically in response to an external stimulus, and mechanical stiffness, which resists deformation, determines the resulting morphology. We analyze the transition to instability of the ground-state flat configuration and characterize the corresponding buckling modes. Building on this framework, we then introduce a continuum mechanical model aimed at describing developing biological systems, accounting for the dynamic and adaptive changes that occur in living tissues over time. Our goal is to understand how topological defects relate to the formation of morphological features in the regenerated organism.