Reactivity of chars and carbons: New insights through molecular modeling

A new model proposed for the gasification of chars and carbons incorporates features of the turbostratic nanoscale structure that exists in such materials. The model also considers the effect of initial surface chemistry and different reactivities perpendicular to the edges and to the faces of the underlying crystallite planes comprising the turbostratic structure. It may be more realistic than earlier models based on pore or grain structure idealizations when the carbon contains large amounts of crystallite matter. Shrinkage of the carbon particles in the chemically controlled regime is also possible due to the random complete gasification of crystallitic planes. This mechanism can explain observations in the literature of particle size reduction. Based on the model predictions, both initial surface chemistry and the number of stacked planes in the crystallites strongly influence the reactivity and particle shrinkage. Its test results agree well with literature data on the air-oxidation of Spherocarb and show that it accurately predicts the variation of particle size with conversion. Model parameters are determined entirely from rate measurements.

Rates of electronic energy transfer in conformationally flexible bichromophoric macrocyclic complexes: a combined experimental and molecular modeling study

Electronic energy transfer (EET) rate constants between a naphthalene donor and anthracene acceptor in [ZnL4a](ClO4)(2) and [ZnL4b](ClO4)(2) were determined by time-resolved fluorescence where L-4a and L-4b are the trans and cis isomers of 6-((anthracen-9-yl-methyl)amino)-6,13-dimethyl-13-((naphthalen-1-yl-methyl)amino)-1,4,8,11-tetraazacyclotetradecane, respectively. These isomers differ in the relative disposition of the appended chromophores with respect to the macrocyclic plane. The trans isomer has an energy transfer rate constant (k(EET)) of 8.7 x 10(8) s(-1), whereas that of the cis isomer is significantly faster (2.3 x 10(9) s(-1)). Molecular modeling was used to determine the likely distribution of conformations in CH3CN solution for these complexes in an attempt to identify any distance or orientation dependency that may account for the differing rate constants observed. The calculated conformational distributions together with analysis by H-1 NMR for the [ZnL4a](2+) trans complex in the common trans-III N-based isomer gave a calculated Forster rate constant close to that observed experimentally. For the [ZnL4b](2+) cis complex, the experimentally determined rate constant may be attributed to a combination of trans-Ill and trans-I N-based isomeric forms of the complex in solution.

Identification and molecular modeling of a novel, plant-like, human purple acid phosphatase

Purple acid phosphatases are a family of binuclear metallohydrolases that have been identified in plants, animals and fungi. Only one isoform of similar to 35 kDa has been isolated from animals, where it is associated with bone resorption and microbial killing through its phosphatase activity, and hydroxyl radical production, respectively. Using the sensitive PSI-BLAST search method, sequences representing new purple acid phosphatase-like proteins have been identified in mammals, insects and nematodes. These new putative isoforms are closely related to the similar to 55 kDa purple acid phosphatase characterized from plants. Secondary structure prediction of the new human isoform further confirms its similarity to a purple acid phosphatase from the red kidney bean. A structural model for the human enzyme was constructed based on the red kidney bean purple acid phosphatase structure. This model shows that the catalytic centre observed in other purple acid phosphatases is also present in this new isoform. These observations suggest that the sequences identified in this study represent a novel subfamily of plant-like purple acid phosphatases in animals and humans. (c) 2006 Elsevier B.V. All rights reserved.

Solid-state NMR structure determination of melittin in a lipid environment

Solid-state C-13 NMR spectroscopy was used to investigate the three-dimensional structure of melittin as lyophilized powder and in ditetradecylphosphatidylcholine (DTPC) membranes. The distance between specifically labeled carbons in analogs [1-C-13]Gly3-[2-C-13]Ala4, [1-C-13]Gly3-[2-C-13]Leu6, [1-C-13]Leu13-[2-C-13]Ala15, [2-C-13]Leu13-[1-C-13]Ala15, and [1-C-13]Leu13-[2-C-13]Leu16 was measured by rotational resonance. As expected, the internuclear distances measured in [1-C-13]Gly3-[2-C-13]Ala4 and [1-C-13]Gly3-[2-C-13]Leu6 were consistent with alpha -helical structure in the N-terminus irrespective of environment. The Internuclear distances measured in [1-C-13]Leu13-[2-C-13]Ala15, [2-C-13]Leu13-[1-C-13]Ala15, and [1-C-13]Leu13-[2-C-13]Leu16 revealed, via molecular modeling, some dependence upon environment for conformation in the region of the bend in helical structure induced by Pro14. A slightly larger interhelical angle between the N- and C-terminal helices was indicated for peptide in dry or hydrated gel state DTPC (139 degrees -145 degrees) than in lyophilized powder (121 degrees -139 degrees) or crystals (129 degrees). The angle, however, is not as great as deduced for melittin in aligned bilayers of DTPC in the liquid-crystalline state (similar to 160 degrees) (R. Smith, F. Separovic, T. J. Milne, A. Whittaker, F. M. Bennett, B. A. Cornell, and A. Makriyannis, 1994, J. Mol, Biol 241:456-466). The study illustrates the utility of rotational resonance in determining local structure within peptide-lipid complexes.

Assessment of the ability to model proteins with leucine-rich repeats in light of the latest structural information

The three-dimensional structures of leucine-rich repeat (LRR) -containing proteins from five different families were previously predicted based on the crystal structure of the ribonuclease inhibitor. using an approach that combined homology-based modeling, structure-based sequence alignment of LRRs, and several rational assumptions. The structural models have been produced based on very limited sequence similarity, which, in general. cannot yield trustworthy predictions. Recently, the protein structures from three of these five families have been determined. In this report we estimate the quality of the modeling approach by comparing the models with the experimentally determined structures. The comparison suggests that the general architecture, curvature, interior/exterior orientations of side chains. and backbone conformation of the LRR structures can be predicted correctly. On the other hand. the analysis revealed that, in some cases. it is difficult to predict correctly the twist of the overall super-helical structure. Taking into consideration the conclusions from these comparisons, we identified a new family of bacterial LRR proteins and present its structural model. The reliability of the LRR protein modeling suggests that it would be informative to apply similar modeling approaches to other classes of solenoid proteins.

Characterization of the retinoid orphan-related receptor-alpha coactivator binding interface: A structural basis for ligand-independent transcription

The retinoid orphan-related receptor-alpha (RORalpha) is a member of the ROR subfamily of orphan receptors and acts as a constitutive activator of transcription in the absence of exogenous ligands. To understand the basis of this activity, we constructed a homology model of Rill using the closely related TRbeta as a template. Molecular modeling suggested that bulky hydrophobic side chains occupy the RORa ligand cavity leaving a small but distinct cavity that may be involved in receptor stabilization. This model was subject to docking simulation with a receptor-interacting peptide from the steroid receptor coactivator, GR-interacting protein-1, which delineated a coactivator binding surface consisting of the signature motif spanning helices 3-5 and helix 12 [activation function 2 (AF2)]. Probing this surface with scanning alanine mutagenesis showed structural and functional equivalence between homologous residues of RORalpha and TRbeta. This was surprising (given that Rill is a ligand-independent activator, whereas TRbeta has an absolute requirement for ligand) and prompted us to use molecular modeling to identify differences between Rill and TRbeta in the way that the All helix interacts with the rest of the receptor. Modeling highlighted a nonconserved amino acid in helix 11 of RORa (Phe491) and a short-length of 3.10 helix at the N terminus of AF2 which we suggest i) ensures that AF2 is locked permanently in the holoconformation described for other liganded receptors and thus 2) enables ligand-independent recruitment of coactivators. Consistent with this, mutation of RORa Phe491 to either methionine or alanine (methionine is the homologous residue in TRbeta), reduced and ablated transcriptional activation and recruitment of coactivators, respectively. Furthermore, we were able to reconstitute transcriptional activity for both a deletion mutant of Ill lacking All and Phe491 Met, by overexpression of a GAL-AF2 fusion protein, demonstrating ligand-independent recruitment of AF2 and a role for Phe491 in recruiting AF2.

Structural basis and prediction of substrate specificity in protein serine/threonine kinases

The large number of protein kinases makes it impractical to determine their specificities and substrates experimentally. Using the available crystal structures, molecular modeling, and sequence analyses of kinases and substrates, we developed a set of rules governing the binding of a heptapeptide substrate motif (surrounding the phosphorylation site) to the kinase and implemented these rules in a web-interfaced program for automated prediction of optimal substrate peptides, taking only the amino acid sequence of a protein kinase as input. We show the utility of the method by analyzing yeast cell cycle control and DNA damage checkpoint pathways. Our method is the only available predictive method generally applicable for identifying possible substrate proteins for protein serine/threonine kinases and helps in silico construction of signaling pathways. The accuracy of prediction is comparable to the accuracy of data from systematic large-scale experimental approaches.

Active site mutations and substrate inhibition in human sulfotransferase 1A1 and 1A3

Human SULT1A1 is primarily responsible for sulfonation of xenobiotics, including the activation of promutagens, and it has been implicated in several forms of cancer. Human SULT1A3 has been shown to be the major sulfotransferase that sulfonates dopamine. These two enzymes shares 93% amino acid sequence identity and have distinct but overlapping substrate preferences. The resolution of the crystal structures of these two enzymes has enabled us to elucidate the mechanisms controlling their substrate preferences and inhibition. The presence of two p-nitrophenol (pNP) molecules in the crystal structure of SULT1A1 was postulated to explain cooperativity at low and inhibition at high substrate concentrations, respectively. In SULT1A1, substrate inhibition occurs with pNP as the substrate but not with dopamine. For SULT1A3, substrate inhibition is found for dopamine but not with pNP. We investigated how substrate inhibition occurs in these two enzymes using molecular modeling, site-directed mutagenesis, and kinetic analysis. The results show that residue Phe-247 of SULT1A1, which interacts with both p-nitrophenol molecules in the active site, is important for substrate inhibition. Mutation of phenylalanine to leucine at this position in SULT1A1 results in substrate inhibition by dopamine. We also propose, based on modeling and kinetic studies, that substrate inhibition by dopamine in SULT1A3 is caused by binding of two dopamine molecules in the active site. © 2004 by The American Society for Biochemistry and Molecular Biology, Inc.

Human sulfotransferases and their role in chemical metabolism

Sulfonation is an important reaction in the metabolism of numerous xenobiotics, drugs, and endogenous compounds. A supergene family of enzymes called sulfotransferases (SULTs) catalyze this reaction. In most cases, the addition of a sulfonate moiety to a compound increases its water solubility and decreases its biological activity. However, many of these enzymes are also capable of bioactivating procarcinogens to reactive electrophiles. In humans three SULT families, SULT1, SULT2, and SULT4, have been identified that contain at least thirteen distinct members. SULTs have a wide tissue distribution and act as a major detoxification enzyme system in adult and the developing human fetus. Nine crystal structures of human cytosolic SULTs have now been determined, and together with site-directed mutagenesis experiments and molecular modeling, we are now beginning to understand the factors that govern distinct but overlapping substrate specificities. These studies have also provided insight into the enzyme kinetics and inhibition characteristics of these enzymes. The regulation of human SULTs remains as one of the least explored areas of research in the field, though there have been some recent advances on the molecular transcription mechanism controlling the individual SULT promoters. Interindividual variation in sulfonation capacity may be important in determining an individual's response to xenobiotics, and recent studies have begun to suggest roles for SULT polymorphism in disease susceptibility. This review aims to provide a summary of our present understanding of the function of human cytosolic sulfotransferases.

Comparison of the binding of 3-fluoromethyl-7-sulfonyl-1,2,3,4-tetrahydroisoquinolines with their isosteric sulfonamides to the active site of phenylethanolamine N-methyltransferase

3-Fluoromethyl-7-(N-substituted aminosulfonyl)-1,2,3,4-tetrahydroisoquinolines (14, 16, and 18-22) are highly potent and selective inhibitors of phenylethanolamine N-methyltransferase (PNMT). Molecular modeling studies with 3-fluoromethyl-7-(N-alkyl aminosulfonyl)-1,2,3,4-tetrahydroisoquinolines, such as 16, suggested that the sulfonamide -NH-could form a hydrogen bond with the side chain of Lys57. However, SAR studies and analysis of the crystal structure of human PNMT (hPNMT) in complex with 7 indicated that the sulfonamide oxygens, and not the sulfonamide -NH-, formed favorable interactions with the enzyme. Thus, we hypothesized that replacement of the sulfonamide -NH-with a methylene group could result in compounds that would retain potency at PNMT and that would have increased lipophilicity, thus increasing the likelihood they will cross the blood brain barrier. A series of 3-fluoromethyl-7-sulfonyl-1,2,3,4-tetrahydroisoquinolines (23-30) were synthesized and evaluated for their PNMT inhibitory potency and affinity for the R2-adrenoceptor. A comparison of these compounds with their isosteric sulfonamides (14, 16, and 18-22) showed that the sulfones were more lipophilic but less potent than their corresponding sulfonamides. Sulfone 24 (hPNMT K-i = 1.3 mu M) is the most potent compound in this series and is quite selective for PNMT versus the R2-adrenoceptor, but 24 is less potent than the corresponding sulfonamide, 16 (hPNMT K-i = 0.13 mu M). We also report the crystal structure of hPNMT in complex with sulfonamide 15, from which a potential hydrogen bond acceptor within the hPNMT active site has been identified, the main chain carbonyl oxygen of Asn39. The interaction of this residue with the sulfonamide -NH-is likely responsible for much of the enhanced inhibitory potency of the sulfonamides versus the sulfones.

Electronic energy-transfer rate constants for geometrical isomers of a bichromophoric macrocyclic complex

The rate of electronic energy transfer (EET) between a naphthalene donor and an anthracene acceptor in [ZnL3]-(ClO4)(2) and [ZnL4](ClO4)(2) was determined by time-resolved fluorescence measurements, where L 3 and L 4 are the geometrical isomers of 6-[(anthracen-9-ylmethyl)amino]-trans-6,13-dimethyl-1,4,8,11-tetraazacyclotetradecane-13-amine (L-2), substituted with either a naphthalen-1-ylmethyl or naphthalen-2-ylmethyl donor, respectively. The energy transfer rate constant, k(EET), was determined to be (0.92 +/- 0.02) x 10(9) s(-1) for the naphthalen-1-ylmethyl-substituted isomer, while that for the naphthalen-2-ylmethyl-substituted isomer is somewhat faster, with k(EET) = (1.31 +/- 0.01) x 10(9) s(-1). The solid-state structure of [(ZnLCl)-Cl-3]ClO4 has been determined, and using molecular modeling calculations, the likely distributions of solution conformations in CH3CN have been evaluated for both complexes. The calculated conformational distributions in the common trans-III N-based isomeric form gave Forster EET rate constants that account for the differences observed and are in excellent agreement with the experimental values. It is shown that the full range of conformers must be considered to accurately reproduce the observed EET kinetics.

Fast, scalable master equation solution algorithms. IV. Lanczos iteration with diffusion approximation preconditioned iterative inversion

In this paper we propose a second linearly scalable method for solving large master equations arising in the context of gas-phase reactive systems. The new method is based on the well-known shift-invert Lanczos iteration using the GMRES iteration preconditioned using the diffusion approximation to the master equation to provide the inverse of the master equation matrix. In this way we avoid the cubic scaling of traditional master equation solution methods while maintaining the speed of a partial spectral decomposition. The method is tested using a master equation modeling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long-lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Modeling molecular transport in slit pores

We examine the transport of methane in microporous carbon by performing equilibrium and nonequilibrium molecular dynamics simulations over a range of pore sizes, densities, and temperatures. We interpret these simulation results using two models of the transport process. At low densities, we consider a molecular flow model, in which intermolecular interactions are neglected, and find excellent agreement between transport diffusion coefficients determined from simulation, and those predicted by the model. Simulation results indicate that the model can be applied up to fluid densities of the order to 0.1-1 nm(-3). Above these densities, we consider a slip flow model, combining hydrodynamic theory with a slip condition at the solid-fluid interface. As the diffusion coefficient at low densities can be accurately determined by the molecular flow model, we also consider a model where the slip condition is supplied by the molecular flow model. We find that both density-dependent models provide a useful means of estimating the transport coefficient that compares well with simulation. (C) 2004 American Institute of Physics.

Modeling the molecular weight distribution of block copolymer formation in a reversible addition-fragmentation chain transfer mediated living radical polymerization

Block copolymers have become an integral part of the preparation of complex architectures through self-assembly. The use of reversible addition-fragmentation chain transfer (RAFT) allows blocks ranging from functional to nonfunctional polymers to be made with predictable molecular weight distributions. This article models block formation by varying many of the kinetic parameters. The simulations provide insight into the overall polydispersities (PDIs) that will be obtained when the chain-transfer constants in the main equilibrium steps are varied from 100 to 0.5. When the first dormant block [polymer-S-C(Z)=S] has a PDI of 1 and the second propagating radical has a low reactivity to the RAFT moiety, the overall PDI will be greater than 1 and dependent on the weight fraction of each block. When the first block has a PDI of 2 and the second propagating radical has a low reactivity to the RAFT moiety, the PDI will decrease to around 1.5 because of random coupling of two broad distributions. It is also shown how we can in principle use only one RAFT agent to obtain block copolymers with any desired molecular weight distribution. We can accomplish this by maintaining the monomer concentration at a constant level in the reactor over the course of the reaction. (c) 2005 Wiley Periodicals, Inc.

Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.