Mathematically complete nucleotide and protein sequence searching using Ssearch.

In this unit a protocol is described for predicting the structure of simple transmembrane a-helical bundles. The protocol is based on a global molecular dynamics search (GMDS) of the configuration space of the helical bundle, yielding several candidate structures. The correct structure among these candidates is selected using information from silent amino acid substitutions, employing the premise that only the correct structure must (by definition) accept all of the silent amino acid substitutions.

Mechanistic implications of the cysteine-nicotinamide adduct in aldehyde dehydrogenase based on quantum mechanical/molecular mechanical simulations.

Recent computer simulations of the cysteine nucleophilic attack on propanal in human mitochondrial aldehyde dehydrogenase (ALDH2) yielded an unexpected result: the chemically reasonable formation of a dead-end cysteine-cofactor adduct when NAD+ was in the "hydride transfer" position. More recently, this adduct found independent crystallographic support in work on formyltetrahydrofolate dehydrogenase, work which further found evidence of the same adduct on re-examination of deposited electron densities of ALDH2. Although the experimental data showed that this adduct was reversible, several mechanistic questions arise from the fact that it forms at all.

A quantum chemical study of the mechanism of action of Vitamin K epoxide reductase (VKOR) II. Transition states.

A reaction path including transition states is generated for the Silverman mechanism [R.B. Silverman, Chemical model studies for the mechanism of Vitamin K epoxide reductase, J.

A quantum chemical study of the mechanism of action of Vitamin K carboxylase (VKC) III. Intermediates and transition states.

A reaction path including transition states is generated for the Dowd mechanism [P. Dowd, R. Hershlne, S.

Molecular recognition of aldehydes by aldehyde dehydrogenase and mechanism of nucleophile activation.

Experimental structural data on the state of substrates bound to class 3 Aldehyde Dehydrogenases (ALDH3A1) is currently unknown. We have utilized molecular mechanics (MM) simulations, in conjunction with new force field parameters for aldehydes, to study the atomic details of benzaldehyde binding to ALDH3A1. Our results indicate that while the nucleophilic Cys243 must be in the neutral state to form what are commonly called near-attack conformers (NACs), these structures do not correlate with increased complexation energy calculated with the MM-Generalized Born Molecular Volume (GBMV) method.

PM3-compatible zinc parameters optimized for metalloenzyme active sites.

Recent studies have shown that semiempirical methods (e.g., PM3 and AM1) for zinc-containing compounds are unreliable for modeling structures containing zinc ions with ligand environments similar to those observed in zinc metalloenzymes.

Initial catalytic events in class 3 aldehyde dehydrogenase: MM and QM/MM simulations.

A novel enzyme mechanism has been predicted by computer simulations for formation of the thiohemiacetal intermediate in the rat ALDH3A1 enzyme. We used molecular mechanics simulations to study the atomic details of substrate binding and quantum mechanical/molecular mechanical methods to study the Cys-243 thiolate attack on benzaldehyde (BA) substrate. BA was found to produce more reactive conformers when aligned for formation of the tetrahedral thiohemiacetal in the R-configuration.

Strategies for multiple sequence alignment.

We present an overview of multiple sequence alignments to outline the practical consequences for the choices among different techniques and parameters. We begin with a discussion of the scoring methods for quantifying the quality of a multiple sequence alignment, followed by a discussion of the algorithms implemented within a variety of multiple sequence alignment programs. We also discuss additional alignment details such as gap penalty and distance metrics.

The emerging science of very early detection of disease outbreaks.

A surge of development of new public health surveillance systems designed to provide more timely detection of outbreaks suggests that public health has a new requirement: extreme timeliness of detection. The authors review previous work relevant to measuring timeliness and to defining timeliness requirements. Using signal detection theory and decision theory, the authors identify strategies to improve timeliness of detection and position ongoing system development within that framework.

Strategies for searching sequence databases.

We provide a detailed overview of the choices inherent in performing a sequence database search, including the choice of algorithm, substitution matrix and gap model. Each of these choices has implications that can be described as restrictions on the underlying model of sequence evolution, the expected degree of divergence between the query sequence and the database sequences (if one uses an evolutionary based matrix), as well as the sensitivity and selectivity of the search. We conclude with a series of recommendations for researchers performing these searches based on our experience and literature studies.

Classification of the environment of protein residues.

We have studied the classification of the environment of residues within protein structures. Eisenberg's original idea created environmental categories to discriminate between similar residues [Bowie et al., Science (1991), 253, 164-170].

The first solvation shell of magnesium and calcium ions in a model nucleic acid environment: an ab initio study.

The interaction of organophosphate anions with divalent metal ions is central to many biological catalytic events. While experimental structural studies can give insight into the likely geometries that can be adopted, quantum mechanics allows for a more complete exploration of the competing forms. Ab initio quantum mechanical calculations have been performed on a series of complexes comprised of dimethyl phosphate, a divalent metal ion (either Mg(II) or Ca(II)) and water of hydration.

Is Ca(II) ion binding to prothrombin fragment 1 intrinsically cooperative, or is the cooperative binding accounted for by self-association?

The recent suggestion that the apparent cooperativity seen in the binding of Ca(II) ions to prothrombin fragment 1 is due to protein aggregation is evaluated. Since (1) we find that the Ca(II) ion binding is not dependent upon protein concentration, (2) the analytical expression for the equilibrium constant of the aggregation model is unrealistically large when evaluated at realistic Ca(II) ion concentrations, and (3) a very simple allosteric cooperative binding model (Monod) can be shown to fit the experimental data, we conclude that the aggregation explanation for the apparent cooperativity in the Ca(II) ion binding by prothrombin fragment 1 is not correct.

The first solvation shell of magnesium ion in a model protein environment with formate, water, and X-NH3, H2S, imidazole, formaldehyde, and chloride as ligands: an Ab initio study.

The first coordination shell of an Mg(II) ion in a model protein environment is studied. Complexes containing a model carboxylate, an Mg(II) ion, various ligands (NH3, H2S, imidazole, and formaldehyde) and water of hydration about the divalent metal ion were geometry optimized. We find that for complexes with the same coordination number, the unidentate carboxylate-Mg(II) ion is greater than 10 kcal mol-1 more stable than the bidentate orientation.

A simple approach for the distribution of computationally intense tasks in an heterogeneous environment: distribution of the MDPP image-processing package.

A method has been developed to link the display ability of a high-resolution graphics workstation with the computational power of a local mainframe or a remote supercomputer via an electronic data network. The method allows this link to be established in a manner largely transparent to the user. The application of the method is illustrated by our successful distribution of the computationally intensive portions of an imaging program (MDPP) from a small VAX workstation to a VAX mainframe and Cray Y-MP8/832 using a simple message-passing technique.

Calcium ion binding to human and bovine factor X.

Human and bovine factor X contain 11 and 12 glutamyl residues, respectively, within the first 40 amino terminal residues that are post-translationally modified to gamma-carboxyglutamyl (Gla) residues. We have measured calcium ion binding to human factor X by equilibrium dialysis. This is the first examination of calcium ion binding to human factor X.

Divalent metal ion mediated interaction of proteins with negatively charged membranes. A model study employing molecular mechanics.

Previous molecular mechanics calculations on the effect of Ca(II) and Mg(II) ions on the conformation of the 18-23 loop of bovine prothrombin [Maynard et al. 1988, Int. J.

Generation of potential structures for the G-domain of chloroplast EF-Tu using comparative molecular modeling.

Comparative molecular modeling has been used to generate several possible structures for the G-domain of chloroplast elongation factor Tu (EF-Tu(chl)) based on the crystallographic data of the homologous E. coli protein. EF-Tu(chl) contains a 10 amino acid insertion not present in the E.

The role of hydrated divalent metal ions in the bridging of two anionic groups. An ab initio quantum chemical and molecular mechanics study of dimethyl phosphate and formate bridged by calcium and magnesium ions.

Ab initio quantum chemical (Gaussian82) and molecular mechanics (AMBER2.0) computational techniques are employed to investigate the interaction of two anions (formate an dimethylphosphate) and a central divalent metal cation (magnesium or calcium). These systems are models for the essential GDP binding unit of the G-proteins (e.

Salt or ion bridges in biological systems: a study employing quantum and molecular mechanics.

Equilibrium geometries and binding energies of model "salt" or "ion" bridge systems have been computed by ab initio quantum chemistry techniques (GAUSSIAN82) and by empirical force field techniques (AMBER2.0). Formate and dimethyl phosphate served as anions in the model compounds while interacting with several organic cations, including methyl ammonium, methyl guanidinium, and divalent metal ion (either Mg2+ or Ca2+) without and with an additional chloride; and a divalent metal ion (either Mg2+ or Ca2+), chloride, and four water molecules of hydration about the metal ion.

A theoretical study of the effect of methylation or ethylation at O6-guanine in the structure and energy of DNA double strands.

Quantum and molecular mechanical calculations were employed to examine the effect on binding energies and structure of methylation and ethylation at O6-guanine in double-stranded DNA. Ab initio quantum chemical calculations (STO-3G, 3-21G) were initially used to pseudo-optimize the structure of the 9-methyl derivative of O6-methylguanine. The distal orientation for the O6-methyl group was found to be lower in energy than the proximal orientation.

Chemical modification of bovine prothrombin fragment 1 in the presence of Tb3+ ions. Sequence studies on 3-gamma-MGlu-fragment.

Chemical modification of the gamma-carboxyglutamyl (Gla) residues of bovine prothrombin fragment 1 using the formaldehyde-morpholine method in the presence of 100 Kappm Tb3+ ions at pH 5.0 provided a modified protein containing 3 gamma-methyleneglutamyl residues (gamma-MGlu) and 7 Gla residues (bovine 3-gamma-MGlu-fragment 1). The modified protein bound the same number of Ca2+ ions as the native protein (six to seven), exhibited 28Mg2+-binding properties identical to native fragment 1 (five Mg2+ ions bound), exhibited the metal ion-promoted quenching of the intrinsic fluorescence in a manner similar to the native protein, but did not bind to phosphatidylserine (PS)/phosphatidylcholine (PC) vesicles in the presence of Ca2+ ions.

Effect of calcium (II) and magnesium (II) ions on the 18-23 gamma-carboxyglutamic acid containing cyclic peptide loop of bovine prothrombin. An AMBER molecular mechanics study.

The effect of calcium (II) and magnesium (II) ions on the conformation of the 18-23 cyclic peptide loop of bovine prothrombin are investigated by the molecular mechanics program AMBER (Assisted Model Building with Energy Refinement). The work is an extension of an earlier paper (Eastman et al., Int.

Studies on Ca(II) binding to gamma-carboxyglutamic acid. Use of thermal decarboxylation to probe metal ion/gamma-carboxyglutamic acid interactions.

The thermal decarboxylation of N-benzyloxycarbonyl-L-gamma-carboxyglutamic acid alpha-methyl ester [Z)-L-Gla-OMe) has been studied. In the presence of increasing amounts of calcium or magnesium ions, lyophilized powders of (Z)-L-Gla-OMe exhibit a corresponding increase in thermal stability. Both magnesium and calcium form relatively tight, thermally stable complexes with (Z)-L-Gla-OMe at high metal ion concentrations.

Relative affinity of Ca(II) and Mg(II) ions for human and bovine prothrombin and fragment 1.

Equilibrium dialysis results are presented for Ca(II) and Mg(II) ion binding to human and bovine prothrombin and fragment 1. Ca(II) ions bind cooperatively, Mg(II) does not.