linear time series with MA(1) residuals

Dear Friends, I want to run a time series linear regression with the form y(t)=beta*x(t)+e(t), where e(t) follows a MA(1) process. the output looks like y=beta*x+e, and e(t)=a(t)-theta*a(t-1). and don't know how to do it using Matlab commands. the help will be deeply appreciated.