A modular treatment of molecular traffic through the active site of cholinesterase.

We present a model for the molecular traffic of ligands, substrates, and products through the active site of cholinesterases (ChEs). First, we describe a common treatment of the diffusion to a buried active site of cationic and neutral species. We then explain the specificity of ChEs for cationic ligands and substrates by introducing two additional components to this common treatment. The first module is a surface trap for cationic species at the entrance to the active-site gorge that operates through local, short-range electrostatic interactions and is independent of ionic strength. The second module is an ionic-strength-dependent steering mechanism generated by long-range electrostatic interactions arising from the overall distribution of charges in ChEs. Our calculations show that diffusion of charged ligands relative to neutral isosteric analogs is enhanced approximately 10-fold by the surface trap, while electrostatic steering contributes only a 1.5- to 2-fold rate enhancement at physiological salt concentration. We model clearance of cationic products from the active-site gorge as analogous to the escape of a particle from a one-dimensional well in the presence of a linear electrostatic potential. We evaluate the potential inside the gorge and provide evidence that while contributing to the steering of cationic species toward the active site, it does not appreciably retard their clearance. This optimal fine-tuning of global and local electrostatic interactions endows ChEs with maximum catalytic efficiency and specificity for a positively charged substrate, while at the same time not hindering clearance of the positively charged products.

The Macromolecular Route to Chiral Amplification.

Cooperative phenomena, described by one-dimensional statistical physical methods, are observed between the enantiomeric characteristics of monomeric materials and the polymers they produce. The effect of minute energies associated with this amplified chirality, although currently not interpretable, can be easily measured. Nonlinear relationships between enantiomeric excess or enantiomeric content and polymer properties may offer the possibility of developing chiral catalysts and chiral chromatographic materials in which the burden of large enantiomeric excess or content may be considerably alleviated. New approaches to information and sensor technology may become possible.

A model of prebiotic replication: survival of the fittest versus extinction of the unfittest.

A model of prebiotic replication discussed by Szathmary and Maynard Smith (1997) is refined in accordance with Lifson's (1997) theory on the crucial stages in the origin of animate matter. The refined model accounts for two processes associated with replication: (1) Replicators always decompose at some rate. (2) Replicators deplete their reactants at a rate equal to the rate of their own replication. Consequently, replicators and their reactants are linked by a non-linear feedback process that keeps them within limits and leads them toward a steady state. The model suggests that: (a) All competing replicators, including those which replicate at sub-exponential "parabolic" and super-exponential "hyperbolic" rates, are subject to natural selection. (b) Survival/extinction is determined by positive/negative net-replication irrespective of the mechanism of replication. (c) Being fit/unfit is the consequence of survival/extinction rather than its cause. In other words, natural selection and survival of the fit is the outcome of continual extinction of the unfit.

On the crucial stages in the origin of animate matter.

Theories of the origin of life have proposed hypotheses to link inanimate to animate matter. The theory proposed here derived the crucial stages in the origin of animate matter directly from the basic properties of inanimate matter. It asked what were the general characteristics of the link, rather than what might have been its chemical details. Life and its origin are shown to be one continuous physicochemical process of replication, random variation, and natural selection. Since life exists here and now, animate properties must have been initiated in the past somewhere. According to the theory, life originated from an as yet unknown elementary autocatalyst which occurred spontaneously, then replicated autocatalytically. As it multiplied to macroscopic abundance, its replicas gradually exhausted their reactants. Random chemical drift initiated diversity among autocatalysts. Diversity led to competition. Competition and depletion of reactants slowed down the rates of net replication of the autocatalysts. Some reached negative rates and became extinct, while those which stayed positive "survived." Thus chemical natural selection appeared, the first step in the transition from inanimate to animate matter. It initiated the first animate property, fitness, i.e., the capacity to adapt to the environment and to survive. As the environment was depleted of reactants, it was enriched with sequels-namely, with decomposition products and all other products which accompany autocatalysis. The changing environment exerted a selective pressure on autocatalysts to replace dwindling reactants by accumulating sequels. Sequels that were incorporated into the autocatalytic process became internal components of complex autocatalytic systems. Primitive forms of metabolism and organization were thus initiated. They evolved further by the same mechanism to ever higher levels of complexity, such as homochirality (handedness) and membranal enclosure. Subsequent evolution by the same mechanism generated cellular metabolism, cell division, information carriers, and a genetic code. Theories of self-organization without natural selection are refuted.

External and internal electrostatic potentials of cholinesterase models.

The electrostatic potentials for the three-dimensional structures of cholinesterases from various species were calculated, using the Delphi algorithm, on the basis of the Poisson-Boltzmann equation. We used structures for Torpedo californica and mouse acetylcholinesterase, and built homology models of the human, Bungarus fasciatus, and Drosophila melanogaster acetylcholinesterases and human butyrylcholinesterase. All these structures reveal a negative external surface potential, in the area around the entrance to the active-site gorge, that becomes more negative as the rim of the gorge is approached. Moreover, in all cases, the potential becomes increasingly more negative along the central axis running down the gorge, and is largest at the base of the gorge, near the active site. Ten key acidic residues conserved in the sequence alignments of AChE from various species, both in the surface area near the entrance of the active-site gorge and at its base, appear to be primarily responsible for these potentials. The potentials are highly correlated among the structures examined, down to sequence identities as low as 35%. This indicates that they are a conserved property of the cholinesterase family, could serve to attract the positively charged substrate into and down the gorge to the active site, and may play other roles important for cholinesterase function.

A helical polymer with a cooperative response to chiral information.

Polyisocyanates, long studied as theoretical models for wormlike chains in dilute solution and liquid crystals, differ from their biological helical analogs in the absence of a pre-determined helical sense. These polymers have an unusual sensitivity to chiral effects that arises from a structure in which alternating right- and left-handed long helical blocks are separated by infrequent and mobile helical reversals. Statistical thermodynamic methods yield an exact description of the polymer and the cooperative nature of its chiral properties. Minute energies that favor one of the helical senses drive easily measurable conformational changes, even though such energies may be extremely difficult to calculate from structural theory. In addition, the chiral nature of the polymer can be used to test theoretical ideas concerned with cholesteric liquid crystals, one of which solves the problem of assigning the helical sense.

Chemical selection, diversity, teleonomy and the second law of thermodynamics. Reflections on Eigen's theory of self-organization of matter.

Two fundamental properties of animate matter, specific complexity and purposeful organization (teleonomy), are traced to their origin, applying Eigen's theory of self-organization of matter. Template-replicating copolymers possess the three dynamic properties that are essential for prebiotic evolution: autocatalysis, diversification and selection. By autocatalysis, even a single microscopic molecule replicates exponentially to macroscopic quantities. By diversification, it extends to a divergent distribution of such molecules. By selection, the distribution converges to a 'quasi-species' that possesses properties like 'survival' and 'adaptation' to its environment. These are teleonomic properties that evolved from a nonteleonomic distribution by selection. Alternating divergent and convergent courses of chemical evolution lead such distributions to ever-growing complexity, including mutual catalytic interactions between the template-replicating copolymers and their chemical environment. Thus, chemical evolution may have started from even a single step, a de novo synthesis of a template-replicating copolymer, and arrived at a primordial living cell, just as biological evolution has started from a primordial cell and arrived at the biological world of today.

Enniatin B and valinomycin as ion carriers: an empirical force field analysis.

The alkali-ion binding properties of two natural depsipeptide ion carriers, enniatin B (EnB) and valinomycin (VM), are examined and compared by the empirical force field method. While VM has been shown to bind preferentially K+, Rb+, and Cs+ over Na+ in most solvents, EnB is considerably less specific. We find that EnB forms two kinds of complexes, internal and external. In internal complexes, the ion binds to all six carbonyl oxygens, while in external ones, only three oxygens, preferentially those of the D-hydroxy-isovaleryl residues, are bound. The size of the internal cavity is best suited for Na+, while K+ and Rb+ squeeze in asymmetrically by distorting the molecule, and Cs+ not at all. External binding is much less specific. Since internal complexes possess much higher strain energies than external ones, the latter may be at least as stable as the former, even in fairly non-polar solvents. VM is calculated to bind only internally, and with much less strain energy than EnB. The size of its internal cavity is well suited for binding the ions K+, Rb+, and Cs+, but is too big for Na+. The difference between the binding energies of Na+ and K+ is much smaller than that between the corresponding hydration enthalpies, thus explaining the binding preference for the latter ion.