source leaprc.protein.ff14SB source leaprc.water.tip3p source leaprc.gaff2 # loading parameters staple molecule loadamberparams link.frcmod # loading chain A with staple coordinates # you may to make bonds in xleap recAL = loadpdb recal.pdb # loading rest of the protein complex recb = loadpdb rec_rest.pdb