There is a strongly geometric characterization of those algebras which arise as $C^\infty(M)$ for $M$ compact and orientable, recently proved by Connes, see 
<a href="http://arxiv.org/abs/0810.2088">here</a>. This has come up on MO before, e.g. in Joel Fine's answer to this question:

https://mathoverflow.net/questions/5344/algebraic-description-of-compact-smooth-manifolds

Like the proofs of most major theorems in differential topology, Connes's approach invokes (a) Riemannian metrics, and (b) hard analysis. Spectral geometry is not my area, so this will be an amateurish explanation... 

If $M$ is a compact smooth manifold, $C^\infty(M)$ is represented faithfully on the Hilbert space $H$ of $L^2$ sections of any hermitian vector bundle $S$. If $D$ is a first-order differential operator acting (unboundedly) on $H$ then we can recover the projective $C^\infty(M)$-module of sections $C^\infty(M;S)$ within $H$ as $\bigcap_{k>0}{dom(D^k)}$ (these domains will actually by Sobolev spaces, I believe). The algebraic counterpart of being first-order is that $[[D,f],g]=0$ for any $f,g \in C^\infty(M)$. 

$D$ has particularly nice properties when it's elliptic. There's no canonical elliptic operator over a smooth manifold until one chooses a Riemannian metric; there's then the signature operator $D=d+d^\ast$ acting on the complexified differential forms. This is an example of a Dirac operator (its square is a Laplacian -  this is a condition on the symbol of the operator). As such, it's formally self-adjoint, Fredholm, and its (real) spectrum has known growth rate depending on $\dim(M)$. 

Connes (see Theorem 11.4) shows that a commutative $\mathbb{R}$-algebra $A$ arises as $C^\infty(M)$ for a smooth manifold structure on the space $M$ Gelfand-dual to $A$ provided that it's part of a "spectral triple" $(A,H,D)$ of the right kind. This means that $A$ should act on a Hilbert space $H$ carrying an unbounded symmetric operator $D$ satisfying various properties. I've hinted at some of these; the most sophisticated property is an "orientation" condition invoking a Hochschild cycle $c\in Z_{\dim M}(A,A)$. This cycle is something like a volume form, and from its components Connes rebuilds local charts.