VMD has a psfgen utility that can do what you want but all of the bonds, angles and dihedrals are defined by the residues within the topology files. What you are asking to do is chemically inconsistent: -CH2-OH => =C=O , for example. These two fragments will not behave in the same fashion dynamically. So, what is your objective? Are you thinking to produce some sort of coarse-grained model with new atom types? OH and CH2, for example? In that case, go read about coarse-grained simulations, for which there are many examples.