name ='d-adenosine'
symbol ='DA'

phosphate = Group('na_phosphate')
sugar = Group('desoxyribose')
base = Group('adenine')

bonds = [Bond(sugar.O_5, phosphate.P), Bond(base.N_9, sugar.C_1), ]

chain_links = [phosphate.P, sugar.O_3]

amber_charge = {phosphate.P: 1.1659,
                phosphate.O_1: -0.7761,
                phosphate.O_2: -0.7761,
                sugar.O_5: -0.4954,
                sugar.C_5: -0.0069,
                sugar.H_51: 0.0754,
                sugar.H_52: 0.0754,
                sugar.C_4: 0.1629,
                sugar.H_4: 0.1176,
                sugar.O_4: -0.3691,
                sugar.C_1: 0.0431,
                sugar.H_1: 0.1838,
                base.N_9: -0.0268,
                base.C_8: 0.1607,
                base.H_8: 0.1877,
                base.N_7: -0.6175,
                base.C_5: 0.0725,
                base.C_6: 0.6897,
                base.N_6: -0.9123,
                base.H_61: 0.4167,
                base.H_62: 0.4167,
                base.N_1: -0.7624,
                base.C_2: 0.5716,
                base.H_2: 0.0598,
                base.N_3: -0.7417,
                base.C_4: 0.3800,
                sugar.C_3: 0.0713,
                sugar.H_3: 0.0985,
                sugar.C_2: -0.0854,
                sugar.H_21: 0.0718,
                sugar.H_22: 0.0718,
                sugar.O_3: -0.5232,
                }
