Simulation of large-scale nonlinear dynamical systems on hardware with a high resemblance to their mathematical equivalents has been always a challenge in engineering. This paper presents a novel currentinput current-output circuit supporting a systematic synthesis procedure of log-domain circuits capable of computing bilateral dynamical systems with considerably low power consumption and acceptable precision. Here, the application of the method is demonstrated by synthesizing four different case studies: 1) a relatively complex two-dimensional (2-D) nonlinear neuron model, 2) a chaotic 3-D nonlinear dynamical system Lorenz attractor having arbitrary solutions for certain parameters, 3) a 2-D nonlinear Hopf oscillator including bistability phenomenon sensitive to initial values and 4) three small neurosynaptic networks comprising three FHN neuron models variously coupled with excitatory and inhibitory synapses. The validity of our approach is verified by nominal and Monte Carlo simulated results with realistic process parameters from the commercially available AMS 0.35 µm technology. The resulting continuous-time, continuous-value and low-power circuits exhibit various bifurcation phenomena, nominal time-domain responses in good agreement with their mathematical counterparts and fairly acceptable process variation results (less than 5% STD).