Under investigation in this paper are the coupled higher-order nonlinear Schrodinger equations with variable coefficients, which represent the propagation of femtosecond soliton pulses comprising of two fields with the left and right polarization in the inhomogeneous optical fiber media. Infinitely-many conservation laws are obtained based on the Lax pair. Via the Hirota method and symbolic computation, bilinear forms, bilinear Backlund transformations, one-and two-soliton-like solutions are also derived. With different coefficients, bell-shaped, periodic-changing, quadratic-varying, exponential-decreasing and exponential-increasing soliton-like profiles are seen, to describe the propagation and interactions of the femtosecond soliton pulses. Head-on and overtaking elastic interactions are shown, which are decided by the directions of the velocities. We also get the bound states with periodic attraction and repulsion between two solitons. (C) 2013 Elsevier B.V. All rights reserved.