A general dynamical theory is presented for the rate constant of weak coupling, nonadiabatic proton-tunneling reactions in solution. The theory incorporates the critical role of the solvent and the vibration of the separation of the heavy particles between which the proton transfers, including their dynamics. The formulation which allows the computation of the quantum rate constant k via classical molecular dynamics simulation techniques is presented, as are a number of approximate analytic results for k in a variety of different important regimes. The frequent appearance of (nearly) classical Arrhenius behavior for k - even though the intrinsic reactive event is quantum proton tunneling - is discussed, together with the solvent and vibrational contributions to the apparent activation energy. In certain weak solvation limits, however, non-Arrhenius behavior for k is found and is related to vibrational Franck-Condon features in the reaction.