The de Rham complex on a smooth manifold $M$ is a complex of sheaves of $\mathbb{R}$-modules, which I will denote by $$\Omega^\bullet \xrightarrow{d} \Omega^{\bullet+1}.$$
A connection (not assumed to be flat!) on a vector bundle $E \to M$ with sheaf of sections $\mathcal{E}$ can be given by a covariant derivative, which is an $\mathbb{R}$-linear sheaf homomorphism:
$$\mathcal{E} \xrightarrow{\nabla} \mathcal{E} \otimes_{\mathcal{C}^\infty} \Omega^1 .$$
Is there an algebraic (in terms of category theory, sheaves and/or homological algebra) way to characterize the construction of the exterior covariant derivative $$\mathcal{E} \otimes_{\mathcal{C}^\infty} \Omega^\bullet \xrightarrow{d^\nabla} \mathcal{E} \otimes_{\mathcal{C}^\infty}\Omega^{\bullet+1}$$ from this data?
Naively, we might just stick stick zero sheaves on either side of the covariant derivative to get a complex, and then tensor with the de Rham complex (over the constant sheaf $\mathbb{R}$). But this is obviously wrong because we get a complex, and the connection is not assumed to be flat. Even in the flat case it is easy to see by looking at the individual sheaves that this approach gives the wrong complex.