We have investigated the potential of sedimentation velocity analytical ultracentrifugation for the measurement of the second virial coefficients of proteins, with the goal of developing a method that allows efficient screening of different solvent conditions. This may be useful for the study of protein crystallization. Macromolecular concentration distributions were modeled using the Lamm equation with the approximation of linear concentration dependencies of the diffusion constant, D = D(o) (1 + k(D)c), and the reciprocal sedimentation coefficient s = s(o)/(1 + k(s)c). We have studied model distributions for their information content with respect to the particle and its non-ideal behavior, developed a strategy for their analysis by direct boundary modeling, and applied it to data from sedimentation velocity experiments on halophilic malate dehydrogenase in complex aqueous solvents containing sodium chloride and 2-methyl-2,4-pentanediol, including conditions near phase separation. Using global modeling for three sets of data obtained at three different protein concentrations, very good estimates for k(s) and s degrees and also for D degrees and the buoyant molar mass were obtained. It was also possible to obtain good estimates for k(D) and the second virial coefficients. Modeling of sedimentation velocity profiles with the non-ideal Lamm equation appears as a good technique to investigate weak inter-particle interactions in complex solvents and also to extrapolate the ideal behavior of the particle.