We have performed molecular dynamics (MD) simulations of multiple copies of unacetylated G5 and G7 and acetylated G5 dendrimers in dimyristoylphosphatidylcholine bilayers with explicit water using the coarse-grained model developed by Marrink et al. (J. Phys. Chem. B 2007, 111, 7812) with the inclusion of long-range electrostatics. When initially clustered together near the bilayer, neutral acetylated dendrimers aggregate, whereas cationic unacetylated dendrimers do not aggregate, but separate from each other, similar to the observations from atomic force microscopy by Mecke et al. (Chem. Phys. Lipids 2004, 132, 3). The bilayers interacting with unacetylated dendrimers of higher concentration are significantly deformed and show pore formation on the positively curved portions, while acetylated dendrimers are unable to form pores. Unacetylated G7 dendrimers bring more water molecules into the pores than do unacetylated G5 dendrimers. These results agree qualitatively with experimental results showing that significant cytoplasmic-protein leakage is produced by unacetylated G7 dendrimers at concentrations as low as 10 nM, but only at a much higher concentration of 400 nM for unacetylated G5 dendrimers (Bioconjugate Chem. 2004, 15, 774). This good qualitative agreement indicates that the effect on pore formation of the concentration and size of large nanoparticles can be studied through coarse-grained MD simulations, provided that long-range electrostatic interactions are included.