The plastic behavior of the insensitive energetic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) is investigated through molecular dynamics simulations.A recent method, built to follow any prescribed deformation path, is used to apply directional shear and compressive deformations to a TATB single crystal, leading to the tridimensional characterization of its nucleation von Mises stresses Ď v (θ, Ď), where θ and Ď are the two angles (latitude and longitude, respectively) that define the loading direction. Furthermore, the local computation of the deformation gradient tensor helps to identify the mechanisms of the irreversible deformation. Two main types of plasticity mechanisms have been identified for the TATB single crystal: first, molecular dynamics simulations predict the existence of dislocations with an unusual local through-plane dilatancy process. Various slip systems among four different non-basal planes have been identified, namely ( 101), (101), (0 11) and (011) planes. Secondly, every deformation containing a basal-plane compressive component involves buckling deformation. A deformation path allowing a perfect twinning of the TATB triclinic cell has been found. This structure has been verified through molecular dynamics (MD) simulations. In order to understand the buckling mechanism, the TATB single crystal behavior under compression along its basal plane is studied in detail.