MMonCa: An Object Kinetic Monte Carlo simulator for damage irradiation evolution and defect diffusion

In this work, we introduce the Object Kinetic Monte Carlo (OKMC) simulator MMonCa and simulate the defect evolution in three different materials. We start by explaining the theory of OKMC and showing some details of how such theory is implemented by creating generic structures and algorithms in the objects that we want to simulate. Then we successfully reproduce simulated results for defect evolution in iron, silicon and tungsten using our simulator and compare with available experimental data and similar simulations. The comparisons validate MMonCa showing that it is powerful and flexible enough to be customized and used to study the damage evolution of defects in a wide range of solid materials.

Definition of amide protection factors for early kinetic intermediates in protein folding

Hydrogen–deuterium exchange experiments have been used previously to investigate the structures of well defined states of a given protein. These include the native state, the unfolded state, and any intermediates that can be stably populated at equilibrium. More recently, the hydrogen–deuterium exchange technique has been applied in kinetic labeling experiments to probe the structures of transiently formed intermediates on the kinetic folding pathway of a given protein. From these equilibrium and nonequilibrium studies, protection factors are usually obtained. These protection factors are defined as the ratio of the rate of exchange of a given backbone amide when it is in a fully solvent-exposed state (usually obtained from model peptides) to the rate of exchange of that amide in some state of the protein or in some intermediate on the folding pathway of the protein. This definition is straightforward for the case of equilibrium studies; however, it is less clear-cut for the case of transient kinetic intermediates. To clarify the concept for the case of burst-phase intermediates, we have introduced and mathematically defined two different types of protection factors: one is Pstruc, which is more related to the structure of the intermediate, and the other is Papp, which is more related to the stability of the intermediate. Kinetic hydrogen–deuterium exchange data from disulfide-intact ribonuclease A and from cytochrome c are discussed to explain the use and implications of these two definitions.

An optimal Mg2+ concentration for kinetic folding of the Tetrahymena ribozyme

Divalent metal ions, such as Mg2+, are generally required for tertiary structure formation in RNA. Although the role of Mg2+ binding in RNA-folding equilibria has been studied extensively, little is known about the role of Mg2+ in RNA-folding kinetics. In this paper, we explore the effect of Mg2+ on the rate-limiting step in the kinetic folding pathway of the Tetrahymena ribozyme. Analysis of these data reveals the presence of a Mg2+-stabilized kinetic trap that slows folding at higher Mg2+ concentrations. Thus, the Tetrahymena ribozyme folds with an optimal rate at 2 mM Mg2+, just above the concentration required for stable structure formation. These results suggest that thermodynamic and kinetic folding of RNA are cooptimized at a Mg2+ concentration that is sufficient to stabilize the folded form but low enough to avoid kinetic traps and misfolding.

The kinetic mechanism of myosin V

Myosin V is an unconventional myosin proposed to be processive on actin filaments, analogous to kinesin on a microtubule [Mehta, A. D., et al. (1999) Nature (London) 400, 590–593]. To ascertain the unique properties of myosin V that permit processivity, we undertook a detailed kinetic analysis of the myosin V motor. We expressed a truncated, single-headed myosin V construct that bound a single light chain to study its innate kinetics, free from constraints imposed by other regions of the molecule. The data demonstrate that unlike any previously characterized myosin a single-headed myosin V spends most of its kinetic cycle (>70%) strongly bound to actin in the presence of ATP. This kinetic tuning is accomplished by increasing several of the rates preceding strong binding to actin and concomitantly prolonging the duration of the strongly bound state by slowing the rate of ADP release. The net result is a myosin unlike any previously characterized, in that ADP release is the rate-limiting step for the actin-activated ATPase cycle. Thus, because of a number of kinetic adaptations, myosin V is tuned for processive movement on actin and will be capable of transporting cargo at lower motor densities than any other characterized myosin.

Kinetic stability as a mechanism for protease longevity

The folding of the extracellular serine protease, α-lytic protease (αLP; EC 3.4.21.12) reveals a novel mechanism for stability that appears to lead to a longer functional lifetime for the protease. For αLP, stability is based not on thermodynamics, but on kinetics. Whereas this has required the coevolution of a pro region to facilitate folding, the result has been the optimization of native-state properties independent of their consequences on thermodynamic stability. Structural and mutational data lead to a model for catalysis of folding in which the pro region binds to a conserved β-hairpin in the αLP C-terminal domain, stabilizing the folding transition state and the native state. The pro region is then proteolytically degraded, leaving the active αLP trapped in a metastable conformation. This metastability appears to be a consequence of pressure to evolve properties of the native state, including a large, highly cooperative barrier to unfolding, and extreme rigidity, that reduce susceptibility to proteolytic degradation. In a test of survival under highly proteolytic conditions, homologous mammalian proteases that have not evolved kinetic stability are much more rapidly degraded than αLP. Kinetic stability as a means to longevity is likely to be a mechanism conserved among the majority of extracellular bacterial pro-proteases and may emerge as a general strategy for intracellular eukaryotic proteases subject to harsh conditions as well.

Kinetic proofreading models for cell signaling predict ways to escape kinetic proofreading

In the context of cell signaling, kinetic proofreading was introduced to explain how cells can discriminate among ligands based on a kinetic parameter, the ligand-receptor dissociation rate constant. In the kinetic proofreading model of cell signaling, responses occur only when a bound receptor undergoes a complete series of modifications. If the ligand dissociates prematurely, the receptor returns to its basal state and signaling is frustrated. We extend the model to deal with systems where aggregation of receptors is essential to signal transduction, and present a version of the model for systems where signaling depends on an extrinsic kinase. We also investigate the kinetics of signaling molecules, “messengers,” that are generated by aggregated receptors but do not remain associated with the receptor complex. We show that the extended model predicts modes of signaling that exhibit kinetic discrimination for some range of parameters but for other parameter values show little or no discrimination and thus escape kinetic proofreading. We compare model predictions with experimental data.

Kinetic coupling between electron and proton transfer in cytochrome c oxidase: simultaneous measurements of conductance and absorbance changes.

Bovine heart cytochrome c oxidase is an electron-current driven proton pump. To investigate the mechanism by which this pump operates it is important to study individual electron- and proton-transfer reactions in the enzyme, and key reactions in which they are kinetically and thermodynamically coupled. In this work, we have simultaneously measured absorbance changes associated with electron-transfer reactions and conductance changes associated with protonation reactions following pulsed illumination of the photolabile complex of partly reduced bovine cytochrome c oxidase and carbon monoxide. Following CO dissociation, several kinetic phases in the absorbance changes were observed with time constants ranging from approximately 3 microseconds to several milliseconds, reflecting internal electron-transfer reactions within the enzyme. The data show that the rate of one of these electron-transfer reactions, from cytochrome a3 to a on a millisecond time scale, is controlled by a proton-transfer reaction. These results are discussed in terms of a model in which cytochrome a3 interacts electrostatically with a protonatable group, L, in the vicinity of the binuclear center, in equilibrium with the bulk through a proton-conducting pathway, which determines the rate of proton transfer (and indirectly also of electron transfer). The interaction energy of cytochrome a3 with L was determined independently from the pH dependence of the extent of the millisecond-electron transfer and the number of protons released, as determined from the conductance measurements. The magnitude of the interaction energy, 70 meV (1 eV = 1.602 x 10(-19) J), is consistent with a distance of 5-10 A between cytochrome a3 and L. Based on the recently determined high-resolution x-ray structures of bovine and a bacterial cytochrome c oxidase, possible candidates for L and a physiological role for L are discussed.

Kinetic analysis of peptide loading onto HLA-DR molecules mediated by HLA-DM.

The nonclassical major histocompatibility complex class II molecule HLA-DM (DM) has recently been shown to play a central role in the class II-associated antigen presentation pathway: DM releases invariant chain-derived CLIP peptides (class II-associated invariant chain protein peptide) from HLA-DR (DR) molecules and thereby facilitates loading with antigenic peptides. Some observations have led to the suggestion that DM acts in a catalytic manner, but so far direct proof is missing. Here, we investigated in vitro the kinetics of exchange of endogenously bound CLIP for various peptides on DR1 and DR2a molecules: we found that in the presence of DM the peptide loading process follows Michaelis-Menten kinetics with turnover numbers of 3-12 DR molecules per minute per DM molecule, and with KM values of 500-1000 nM. In addition, surface plasmon resonance measurements showed that DM interacts efficiently with DR-CLIP complexes but only weakly with DR-peptide complexes isolated from DM-positive cells. Taken together, our data provide evidence that DM functions as an enzyme-like catalyst of peptide exchange and favors the generation of long-lived DR-peptide complexes that are no longer substrates for DM.

Well posedness of an integrodifferential kinetic model of Fokker-Planck type for angiogenesis.

Tumor induced angiogenesis processes including the effect of stochastic motion and branching of blood vessels can be described coupling a (nonlocal in time) integrodifferential kinetic equation of Fokker–Planck type with a diffusion equation for the tumor induced ingiogenic factor. The chemotactic force field depends on the flux of blood vessels through the angiogenic factor. We develop an existence and uniqueness theory for this system under natural assumptions on the initial data. The proof combines the construction of fundamental solutions for associated linearized problems with comparison principles, sharp estimates of the velocity integrals and compactness results for this type of kinetic and parabolic operators

Kinetic study of the decompositions involved in the thermal degradation of commercial azodicarbonamide

The transitions and reactions involved in the thermal treatment of several commercial azodicarbonamides (ADC) in an inert atmosphere have been studied by dynamic thermogravimetry analysis (TGA), mass spectrometry and Fourier transform infrared (FTIR) spectroscopy. A pseudo-mechanistic model, involving several competitive and non-competitive reactions, has been suggested and applied to the correlation of the weight loss data. The model applied is capable of accurately representing the different processes involved, and can be of great interest in the understanding and quantification of such phenomena, including the simulation of the instantaneous amount of gases evolved in a foaming process. In addition, a brief discussion on the methodology related to the mathematical modeling of TGA data is presented, taking into account the complex thermal behaviour of the ADC.