Ultrafast detection and control of molecular dynamics

Many elementary chemical and physical processes such as the breaking of a chemical bond or the vibrational motion of atoms within a molecule take place on a femtosecond (fs = 10−15 s) or picosecond (ps = 10−12 s) time scale. It is now possible to monitor these events as a function of time with temporal resolution well below 100 fs. This capability is based on the pump-probe technique where one optical pulse triggers a reaction and a second delayed optical pulse probes the changes that ensue. To illustrate this capability, the dynamics of ligand motion within a protein are presented. Moving beyond casual observation of a reaction to active control of its outcome requires additional experimental and theoretical effort. To illustrate the concept of control, the effect of optical pulse duration on the vibrational dynamics of a tri-atomic molecule are discussed. The experimental and theoretical resources currently available are poised to make the dream of reaction control a reality for certain molecular systems.

Ionophores and receptors using cation-π interactions: Collarenes

Cation-π interactions are important forces in molecular recognition by biological receptors, enzyme catalysis, and crystal engineering. We have harnessed these interactions in designing molecular systems with circular arrangement of benzene units that are capable of acting as ionophores and models for biological receptors. [n]Collarenes are promising candidates with high selectivity for a specific cation, depending on n, because of their structural rigidity and well-defined cavity size. The interaction energies of [n]collarenes with cations have been evaluated by using ab initio calculations. The selectivity of these [n]collarenes in aqueous solution was revealed by using statistical perturbation theory in conjunction with Monte Carlo and molecular dynamics simulations. It has been observed that in [n]collarenes the ratio of the interaction energies of a cation with it and the cation with the basic building unit (benzene) can be correlated to its ion selectivity. We find that collarenes are excellent and efficient ionophores that bind cations through cation-π interactions. [6]Collarene is found to be a selective host for Li+ and Mg2+, [8]collarene for K+ and Sr2+, and [10]collarene for Cs+ and Ba2+. This finding indicates that [10]collarene and [8]collarene could be used for effective separation of highly radioactive isotopes, 137Cs and 90Sr, which are major constituents of nuclear wastes. More interestingly, collarenes of larger cavity size can be useful in capturing organic cations. [12]Collarene exhibits a pronounced affinity for tetramethylammonium cation and acetylcholine, which implies that it could serve as a model for acetylcholinestrase. Thus, collarenes can prove to be novel and effective ionophores/model-receptors capable of heralding a new direction in molecular recognition and host-guest chemistry.

Quantitative modeling of stochastic systems in molecular biology by using stochastic Petri nets

An integrated understanding of molecular and developmental biology must consider the large number of molecular species involved and the low concentrations of many species in vivo. Quantitative stochastic models of molecular interaction networks can be expressed as stochastic Petri nets (SPNs), a mathematical formalism developed in computer science. Existing software can be used to define molecular interaction networks as SPNs and solve such models for the probability distributions of molecular species. This approach allows biologists to focus on the content of models and their interpretation, rather than their implementation. The standardized format of SPNs also facilitates the replication, extension, and transfer of models between researchers. A simple chemical system is presented to demonstrate the link between stochastic models of molecular interactions and SPNs. The approach is illustrated with examples of models of genetic and biochemical phenomena where the UltraSAN package is used to present results from numerical analysis and the outcome of simulations.

Molecular Interaction Map of the Mammalian Cell Cycle Control and DNA Repair Systems

Eventually to understand the integrated function of the cell cycle regulatory network, we must organize the known interactions in the form of a diagram, map, and/or database. A diagram convention was designed capable of unambiguous representation of networks containing multiprotein complexes, protein modifications, and enzymes that are substrates of other enzymes. To facilitate linkage to a database, each molecular species is symbolically represented only once in each diagram. Molecular species can be located on the map by means of indexed grid coordinates. Each interaction is referenced to an annotation list where pertinent information and references can be found. Parts of the network are grouped into functional subsystems. The map shows how multiprotein complexes could assemble and function at gene promoter sites and at sites of DNA damage. It also portrays the richness of connections between the p53-Mdm2 subsystem and other parts of the network.

Molecular keys to speciation: DNA polymorphism and the control of genetic exchange in enterobacteria

Speciation involves the establishment of genetic barriers between closely related organisms. The extent of genetic recombination is a key determinant and a measure of genetic isolation. The results reported here reveal that genetic barriers can be established, eliminated, or modified by manipulating two systems which control genetic recombination, SOS and mismatch repair. The extent of genetic isolation between enterobacteria is a simple mathematical function of DNA sequence divergence. The function does not depend on hybrid DNA stability, but rather on the number of blocks of sequences identical in the two mating partners and sufficiently large to allow the initiation of recombination. Further, there is no obvious discontinuity in the function that could be used to define a level of divergence for distinguishing species.

On the controllability of distributed systems

To “control” a system is to make it behave (hopefully) according to our “wishes,” in a way compatible with safety and ethics, at the least possible cost. The systems considered here are distributed—i.e., governed (modeled) by partial differential equations (PDEs) of evolution. Our “wish” is to drive the system in a given time, by an adequate choice of the controls, from a given initial state to a final given state, which is the target. If this can be achieved (respectively, if we can reach any “neighborhood” of the target) the system, with the controls at our disposal, is exactly (respectively, approximately) controllable. A very general (and fuzzy) idea is that the more a system is “unstable” (chaotic, turbulent) the “simplest,” or the “cheapest,” it is to achieve exact or approximate controllability. When the PDEs are the Navier–Stokes equations, it leads to conjectures, which are presented and explained. Recent results, reported in this expository paper, essentially prove the conjectures in two space dimensions. In three space dimensions, a large number of new questions arise, some new results support (without proving) the conjectures, such as generic controllability and cases of decrease of cost of control when the instability increases. Short comments are made on models arising in climatology, thermoelasticity, non-Newtonian fluids, and molecular chemistry. The Introduction of the paper and the first part of all sections are not technical. Many open questions are mentioned in the text.

Biosynthetic control of molecular weight in the polymerization of the octasaccharide subunits of succinoglycan, a symbiotically important exopolysaccharide of Rhizobium meliloti

Succinoglycan, a symbiotically important exopolysaccharide of Rhizobium meliloti, is composed of polymerized octasaccharide subunits, each of which consists of one galactose and seven glucoses with succinyl, acetyl, and pyruvyl modifications. Production of specific low molecular weight forms of R. meliloti exported and surface polysaccharides, including succinoglycan, appears to be important for nodule invasion. In a previous study of the roles of the various exo gene products in succinoglycan biosynthesis, exoP, exoQ, and exoT mutants were found to synthesize undecaprenol-linked fully modified succinoglycan octasaccharide subunits, suggesting possible roles for their gene products in polymerization or transport. Using improved techniques for analyzing succinoglycan biosynthesis by these mutants, we have obtained evidence indicating that R. meliloti has genetically separable systems for the synthesis of high molecular weight succinoglycan and the synthesis of a specific class of low molecular weight oligosaccharides consisting of dimers and trimers of the octasaccharide subunit. Models to account for our unexpected findings are discussed. Possible roles for the ExoP, ExoQ, and ExoT proteins are compared and contrasted with roles that have been suggested on the basis of homologies to key proteins involved in the biosynthesis of O-antigens and of certain exported or capsular cell surface polysaccharides.

Kinetic equilibrium of forces and molecular events in muscle contraction

In biomolecular systems, the mechanical transfer of free energy occurs with both high efficiency and high speed. It is shown here that such a transfer can be achieved only if the participating free-energy-storing elements exhibit opposing relationships between their content of free energy and the force they exert in the transfer direction. A kinetic equilibrium of forces (KEF) results, in which the transfer of free energy is mediated essentially by thermal molecular motion. On the basis of present evidence, KEF is used as a guiding principle in developing a mechanical model of the crossbridge cycle in muscle contraction. The model allows the basic features of molecular events to be visualized in terms of plausible structures. Real understanding of the process will require identification of the elements that perform the functions described here. Besides chemomechanical energy transduction, KEF may have a role in other biomolecular processes in which free energy is transferred mechanically over large distances.

Molecular characterization and assembly of the needle complex of the Salmonella typhimurium type III protein secretion system

Many bacterial pathogens of plants and animals have evolved a specialized protein-secretion system termed type III to deliver bacterial proteins into host cells. These proteins stimulate or interfere with host cellular functions for the pathogen's benefit. The Salmonella typhimurium pathogenicity island 1 encodes one of these systems that mediates this bacterium's ability to enter nonphagocytic cells. Several components of this type III secretion system are organized in a supramolecular structure termed the needle complex. This structure is made of discrete substructures including a base that spans both membranes and a needle-like projection that extends outward from the bacterial surface. We demonstrate here that the type III secretion export apparatus is required for the assembly of the needle substructure but is dispensable for the assembly of the base. We show that the length of the needle segment is determined by the type III secretion associated protein InvJ. We report that InvG, PrgH, and PrgK constitute the base and that PrgI is the main component of the needle of the type III secretion complex. PrgI homologs are present in type III secretion systems from bacteria pathogenic for animals but are absent from bacteria pathogenic for plants. We hypothesize that the needle component may establish the specificity of type III secretion systems in delivering proteins into either plant or animal cells.

Collecting and harvesting biological data: the GPCRDB and NucleaRDB information systems

The amount of genomic and proteomic data that is entered each day into databases and the experimental literature is outstripping the ability of experimental scientists to keep pace. While generic databases derived from automated curation efforts are useful, most biological scientists tend to focus on a class or family of molecules and their biological impact. Consequently, there is a need for molecular class-specific or other specialized databases. Such databases collect and organize data around a single topic or class of molecules. If curated well, such systems are extremely useful as they allow experimental scientists to obtain a large portion of the available data most relevant to their needs from a single source. We are involved in the development of two such databases with substantial pharmacological relevance. These are the GPCRDB and NucleaRDB information systems, which collect and disseminate data related to G protein-coupled receptors and intra-nuclear hormone receptors, respectively. The GPCRDB was a pilot project aimed at building a generic molecular class-specific database capable of dealing with highly heterogeneous data. A first version of the GPCRDB project has been completed and it is routinely used by thousands of scientists. The NucleaRDB was started recently as an application of the concept for the generalization of this technology. The GPCRDB is available via the WWW at http://www.gpcr.org/7tm/ and the NucleaRDB at http://www.receptors.org/NR/.

Selenoprotein oxidoreductase with specificity for thioredoxin and glutathione systems

Thioredoxin (Trx) and glutathione (GSH) systems are considered to be two major redox systems in animal cells. They are reduced by NADPH via Trx reductase (TR) or oxidized GSH (GSSG) reductase and further supply electrons for deoxyribonucleotide synthesis, antioxidant defense, and redox regulation of signal transduction, transcription, cell growth, and apoptosis. We cloned and characterized a pyridine nucleotide disulfide oxidoreductase, Trx and GSSG reductase (TGR), that exhibits specificity for both redox systems. This enzyme contains a selenocysteine residue encoded by the TGA codon. TGR can reduce Trx, GSSG, and a GSH-linked disulfide in in vitro assays. This unusual substrate specificity is achieved by an evolutionary conserved fusion of the TR and glutaredoxin domains. These observations, together with the biochemical probing and molecular modeling of the TGR structure, suggest a mechanism whereby the C-terminal selenotetrapeptide serves a role of a protein-linked GSSG and shuttles electrons from the disulfide center within the TR domain to either the glutaredoxin domain or Trx.

Symmetry, stability, and dynamics of multidomain and multicomponent protein systems

Symmetry is commonly observed in many biological systems. Here we discuss representative examples of the role of symmetry in structural molecular biology. Point group symmetries are observed in many protein oligomers whose three-dimensional atomic structures have been elucidated by x-ray crystallography. Approximate symmetry also occurs in multidomain proteins. Symmetry often confers stability on the molecular system and results in economical usage of basic components to build the macromolecular structure. Symmetry is also associated with cooperativity. Mild perturbation from perfect symmetry may be essential in some systems for dynamic functions.

Molecular Cloning of Vacuolar H+-Pyrophosphatase and Its Developmental Expression in Growing Hypocotyl of Mung Bean1

Vacuolar proton-translocating inorganic pyrophosphatase and H+-ATPase acidify the vacuoles and power the vacuolar secondary active transport systems in plants. Developmental changes in the transcription of the pyrophosphatase in growing hypocotyls of mung bean (Vigna radiata) were investigated. The cDNA clone for the mung bean enzyme contains an uninterrupted open reading frame of 2298 bp, coding for a polypeptide of 766 amino acids. Hypocotyls were divided into elongating and mature regions. RNA analysis revealed that the transcript level of the pyrophosphatase was high in the elongating region of the 3-d-old hypocotyl but was extremely low in the mature region of the 5-d-old hypocotyl. The level of transcript of the 68-kD subunit of H+-ATPase also decreased after cell maturation. In the elongating region, the proton-pumping activity of pyrophosphatase on the basis of membrane protein was 3 times higher than that of H+-ATPase. After cell maturation, the pyrophosphatase activity decreased to 30% of that in the elongating region. The decline in the pyrophosphatase activity was in parallel with a decrease in the enzyme protein content. These findings indicate that the level of the pyrophosphatase, a main vacuolar proton pump in growing cells, is negatively regulated after cell maturation at the transcriptional level.

Unifying theory of hypoxia tolerance: molecular/metabolic defense and rescue mechanisms for surviving oxygen lack.

We develop a unifying theory of hypoxia tolerance based on information from two cell level models (brain cortical cells and isolated hepatocytes) from the highly anoxia tolerant aquatic turtle and from other more hypoxia sensitive systems. We propose that the response of hypoxia tolerant systems to oxygen lack occurs in two phases (defense and rescue). The first lines of defense against hypoxia include a balanced suppression of ATP-demand and ATP-supply pathways; this regulation stabilizes (adenylates) at new steady-state levels even while ATP turnover rates greatly decline. The ATP demands of ion pumping are down-regulated by generalized "channel" arrest in hepatocytes and by "spike" arrest in neurons. Hypoxic ATP demands of protein synthesis are down-regulated probably by translational arrest. In hypoxia sensitive cells this translational arrest seems irreversible, but hypoxia-tolerant systems activate "rescue" mechanisms if the period of oxygen lack is extended by preferentially regulating the expression of several proteins. In these cells, a cascade of processes underpinning hypoxia rescue and defense begins with an oxygen sensor (a heme protein) and a signal-transduction pathway, which leads to significant gene-based metabolic reprogramming-the rescue process-with maintained down-regulation of energy-demand and energy-supply pathways in metabolism throughout the hypoxic period. This recent work begins to clarify how normoxic maintenance ATP turnover rates can be drastically (10-fold) down-regulated to a new hypometabolic steady state, which is prerequisite for surviving prolonged hypoxia or anoxia. The implications of these developments are extensive in biology and medicine.