The 1912 Mürefte and 1999 Izmit 7.4 earthquakes are the last devastating events of the western and eastern sections of the Marmara region, respectively. The center of the Sea of Marmara, the seismic gap along the North Anatolian Fault (NAF) between the rupture zones of these two earthquakes, is prone to creating one or more large earthquakes. The main objective of our study is to determine 3D dynamic earthquake rupture scenarios, considering non-planar and heterogeneous stress distribution in the Sea of Marmara. Recent studies show that some segments of the North Anatolian Fault (NAF) beneath the Marmara are partially creeping. In this study, it is the first time that we attempt to generate realistic earthquake scenarios by putting constrains on initial stress of the fault using regional stress from earthquake focal mechanisms, in addition to stress release during past earthquakes and strain accumulation during interseismic period using geodetical measurements on slip-rate and locking depth at various segments along the NAF beneath the Sea of Marmara. In order to constrain the regional stress we combine stress tensor inversion analyses from our studies, which contain six cluster of earthquakes, and other previous local studies in the Sea of Marmara. We use 3D Finite Element Method (PyLith) for dynamic earthquake rupture simulations and tetragonal mesh for better smoothing at the fault bends, which allows us to implement nonplanar fault geometry and initial stress heterogeneity using slip-weakening friction law. We place constraints on initial shear stress from geodetic and seismic studies of locking depth and interseismic strain accumulation. We consider around 100 scenarios and calculate slip distribution, rupture velocity and moment magnitude in addition to slip-rate and traction on the fault surface, and displacement and velocity on the ground surface. We find that for the most of the scenarios possible earthquake magnitude does not exceed Mw7.2. In addition, in none of the possible scenarios we obtain super-shear rupture velocity. We find that depending on the location of the initiation point, asperities in the partially creeping segments and loaded initial stress, the rupture may not extend into the Prince’s Island Segment.