diff --git a/openff/interchange/components/mdconfig.py b/openff/interchange/components/mdconfig.py
index b584623f..98f5752a 100644
--- a/openff/interchange/components/mdconfig.py
+++ b/openff/interchange/components/mdconfig.py
@@ -207,6 +207,7 @@ def write_lammps_input(
         self,
         interchange: "Interchange",
         input_file: str = "run.in",
+        data_file: str = "out.lmp",
     ) -> None:
         """Write a LAMMPS input file for running single-point energies."""
         # TODO: Get constrained angles
@@ -338,14 +339,19 @@ def _get_coeffs_of_constrained_bonds_and_angles(
             vdw_cutoff = round(self.vdw_cutoff.m_as(unit.angstrom), 4)
             coul_cutoff = round(self.coul_cutoff.m_as(unit.angstrom), 4)
 
-            if self.coul_method == _PME:
-                lmp.write(f"pair_style lj/cut/coul/long {vdw_cutoff} {coul_cutoff}\n")
-            elif self.coul_method == "cutoff":
-                lmp.write(f"pair_style lj/cut/coul/cut {vdw_cutoff} {coul_cutoff}\n")
+            if self.vdw_method == "pme" and self.coul_method == "pme":
+                lmp.write("pair_style lj/long/coul/long")
+            elif self.vdw_method == "cutoff":
+                if self.coul_method == _PME:
+                    lmp.write("pair_style lj/cut/coul/long")
+                elif self.coul_method == "cutoff":
+                    lmp.write("pair_style lj/cut/coul/cut")
             else:
                 raise UnsupportedExportError(
-                    f"Unsupported electrostatics method {self.coul_method}",
+                    "Unsupported nonbonded methods found: "
+                    f"{self.coul_method=}, {self.vdw_method=}",
                 )
+            lmp.write(f" {vdw_cutoff} {coul_cutoff}\n")
 
             if self.mixing_rule == "lorentz-berthelot":
                 lmp.write("pair_modify mix arithmetic tail yes\n\n")
@@ -355,8 +361,7 @@ def _get_coeffs_of_constrained_bonds_and_angles(
                 raise UnsupportedExportError(
                     f"Mixing rule {self.mixing_rule} not supported",
                 )
-
-            lmp.write("read_data out.lmp\n\n")
+            lmp.write(f"read_data {data_file}\n\n")
             lmp.write(
                 "thermo_style custom ebond eangle edihed eimp epair evdwl ecoul elong etail pe\n\n",
             )