Molecular Mechanisms and Targets of Cyclic Guanosine Monophosphate (cGMP) in Vascular Smooth Muscles

Molecular mechanisms and targets of cyclic guanosine monophosphate (cGMP) accounting for vascular smooth muscles (VSM) contractility are reviewed. Mathematical models of five published mechanisms are presented, and four novel mechanisms are proposed. cGMP, which is primarily produced by the nitric oxide (NO) dependent soluble guanylate cyclase (sGC), activates cGMP-dependent protein kinase (PKG). The NO/cGMP/PKG signaling pathway targets are the mechanisms that regulate cytosolic calcium ([Ca2+]i) signaling and those implicated in the Ca2+-desensitization of the contractile apparatus. In addition to previous mathematical models of cGMP-mediated molecular mechanisms targeting [Ca2+]i regulation, such as large-conductance Ca2+-activated K+ channels (BKCa), Ca2+-dependent Cl− channels (ClCa), Na+/Ca2+ exchanger (NCX), Na+/K+/Cl− cotransport (NKCC), and Na+/K+-ATPase (NKA), other four novel mechanisms are proposed here based on the existing but perhaps overlooked experimental results. These are the effects of cGMP on the sarco−/endo- plasmic reticulum Ca2+-ATPase (SERCA), the plasma membrane Ca2+-ATPase (PMCA), the inositol 1,4,5-trisphosphate (IP3) receptor channels type 1 (IP3R1), and on the myosin light chain phosphatase (MLCP), which is implicated in the Ca2+-desensitization. Different modeling approaches are presented and discussed, and novel model descriptions are proposed

MLC-kinase/phosphatase control of Ca[sup]2+ signal transduction in airway smooth muscles

In airway smooth muscles, kinase/phosphatase-dependent phosphorylation and dephosphorylation of the myosin light chain (MLC) have been revealed by many authors as important steps in calcium ▫$(Ca^{2+})$▫ signalling pathway from the variation of ▫$Ca^{2+}$▫ concentration in cytosol to the force development. Here, a theoretical analysis of the control action of MLC-kinase (MLCK) and MLC-phosphatase (MLCP) in ▫$Ca^{2+}$▫ signalling is presented and related to the general control principles of these enzymes, which were previously studied by Reinhart Heinrich and his co-workers. The kinetic scheme of the mathematical model considers interactions among ▫$Ca^{2+}$▫, calmodulin (CaM) and MLCK and the well-known 4-state actomyosin latch bridge model, whereby a link between them is accomplished by the conservation relation of all species of MLCK. The mathematical model predicts the magnitude and velocity of isometric force in smooth muscles upon transient biphasic ▫$Ca^{2+}$▫ signal. The properties of signal transduction in the system such as the signalling time, signal duration and signal amplitude, which are reflected in the properties of force developed, are studied by the principles of the metabolic control theory. The analysis of our model predictions confirms as shown by Reinhart Heinrich and his co-workers that MLCK controls the amplitude of signal more than its duration, whereas MLCP controls both. Finally, the simulations of elevated total content of MLCK, a typical feature of bronchial muscles of asthmatic subjects and spontaneously hypertensive rats as well as potentiation of MLCP catalytic activity, are carried out and are discussed in view of an increase in the force magnitude

Self-organization of enzyme-catalyzed reactions studied by the maximum entropy production principle

The self-organization of open reaction systems is closely related to specific mechanisms that allow the export of internally generated entropy from systems to their environment. According to the second law of thermodynamics, systems with effective entropy export to the environment are better internally organized. Therefore, they are in thermodynamic states with low entropy. In this context, we study how self-organization in enzymatic reactions depends on their kinetic reaction mechanisms. Enzymatic reactions in an open system are considered to operate in a non-equilibrium steady state, which is achieved by satisfying the principle of maximum entropy production (MEPP). The latter is a general theoretical framework for our theoretical analysis. Detailed theoretical studies and comparisons of the linear irreversible kinetic schemes of an enzyme reaction in two and three states are performed. In both cases, in the optimal and statistically most probable thermodynamic steady state, a diffusion-limited flux is predicted by MEPP. Several thermodynamic quantities and enzymatic kinetic parameters, such as the entropy production rate, the Shannon information entropy, reaction stability, sensitivity, and specificity constants, are predicted. Our results show that the optimal enzyme performance may strongly depend on the number of reaction steps when linear reaction mechanisms are considered. Simple reaction mechanisms with a smaller number of intermediate reaction steps could be better organized internally and could allow fast and stable catalysis. These could be features of the evolutionary mechanisms of highly specialized enzymes

Theoretical and experimental investigation of calcium-contraction coupling in airway smooth muscle

We investigated theoretically and experimentally the ▫$Ca^{2+}$▫-contraction couplingin rat tracheal smooth muscle. ▫$[Ca^{2+}]_i$▫, isometric contraction and myosin light chain (MLC) phosphorylation were measured in response to 1 mM carbachol. Theoretical modeling consisted in coupling a model of ▫$Ca^{2+}-dependent$▫ MLC kinase (MLCK) activation with a four-state model of smooth muscle contractile apparatus. Stimulation resulted in a short-time contraction obtained within 1 min, followed by a long-time contraction up to the maximal force obtained in 30 min. ML-7 and Wortmannin (MLCK inhibitors) abolished the contraction. Chelerythrine (PKC inhibitor) did not change the short-time, but reduced the long-time contraction. ▫$[Ca^{2+}]_i$▫ responses of isolated myocytes recorded during the first 90 s consisted in a fast peak, followed by a plateau phase and, in 28 % of the cells, superimposed ▫$Ca^{2+}$▫ oscillations. MLC phosphorylation was maximal at 5 s and then decreased, whereas isometric contraction followed a Hill-shaped curve. The model properlypredicts the time course of MLC phosphorylation and force of the short-time response. With oscillating ▫$Ca^{2+}$▫ signal, the predicted force does not oscillate. According to the model, the amplitude of the plateau and the frequency of oscillations encode for the amplitude of force, whereas the peak encodes for force velocity. The long-time phase of the contraction, associated with a second increase in MLC phosphorylation, may be explained, at least partially, by MLC phosphatase (MLCP) inhibition, possibly via PKC inhibition