All-atom empirical potential for molecular modeling and dynamics studies of proteins

New protein parameters are reported for the all-atom empirical energy function in the CHARMM program. The parameter evaluation was based on a self-consistent approach designed to achieve a balance between the internal (bonding) and interaction (nonbonding) terms of the force field and among the solvent-solvent, solvent-solute, and solute-solute interactions. Optimization of the internal parameters used experimental gas-phase geometries, vibrational spectra, and torsional energy surfaces supplemented with ab initio results. The peptide backbone bonding parameters were optimized with respect to data for N-methylacetamide and the alanine dipeptide. The interaction parameters, particularly the atomic charges, were determined by fitting ab initio interaction energies and geometries of complexes between water and model compounds that represented the backbone and the various side chains. In addition, dipole moments, experimental heats and free energies of vaporization, solvation and sublimation, molecular volumes, and crystal pressures and structures were used in the optimization. The resulting protein parameters were tested by applying them to noncyclic tripeptide crystals, cyclic peptide crystals, and the proteins crambin, bovine pancreatic trypsin inhibitor, and carbonmonoxy myoglobin in vacuo and in crystals. A detailed analysis of the relationship between the alanine dipeptide potential energy surface and calculated protein φ, χ angles was made and used in optimizing the peptide group torsional parameters. The results demonstrate that use of ab initio structural and energetic data by themselves are not sufficient to obtain an adequate backbone representation for peptides and proteins in solution and in crystals. Extensive comparisons between molecular dynamics simulations and experimental data for polypeptides and proteins were performed for both structural and dynamic properties. Energy minimization and dynamics simulations for crystals demonstrate that the latter are needed to obtain meaningful comparisons with experimental crystal structures. The presented parameters, in combination with the previously published CHARMM all-atom parameters for nucleic acids and lipids, provide a consistent set for condensed-phase simulations of a wide variety of molecules of biological interest.

CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields

The widely used CHARMM additive all-atom force field includes parameters for proteins, nucleic acids, lipids, and carbohydrates. In the present article, an extension of the CHARMM force field to drug-like molecules is presented. The resulting CHARMM General Force Field (CGenFF) covers a wide range of chemical groups present in biomolecules and drug-like molecules, including a large number of heterocyclic scaffolds. The parametrization philosophy behind the force field focuses on quality at the expense of transferability, with the implementation concentrating on an extensible force field. Statistics related to the quality of the parametrization with a focus on experimental validation are presented. Additionally, the parametrization procedure, described fully in the present article in the context of the model systems, pyrrolidine, and 3-phenoxymethylpyrrolidine will allow users to readily extend the force field to chemical groups that are not explicitly covered in the force field as well as add functional groups to and link together molecules already available in the force field. CGenFF thus makes it possible to perform "all-CHARMM" simulations on drug-target interactions thereby extending the utility of CHARMM force fields to medicinally relevant systems.

Automation of the CHARMM General Force Field (CGenFF) I: bond perception and atom typing

Molecular mechanics force fields are widely used in computer-aided drug design for the study of drug-like molecules alone or interacting with biological systems. In simulations involving biological macromolecules, the biological part is typically represented by a specialized biomolecular force field, while the drug is represented by a matching general (organic) force field. In order to apply these general force fields to an arbitrary drug-like molecule, functionality for assignment of atom types, parameters, and charges is required. In the present article, which is part I of a series of two, we present the algorithms for bond perception and atom typing for the CHARMM General Force Field (CGenFF). The CGenFF atom typer first associates attributes to the atoms and bonds in a molecule, such as valence, bond order, and ring membership among others. Of note are a number of features that are specifically required for CGenFF. This information is then used by the atom typing routine to assign CGenFF atom types based on a programmable decision tree. This allows for straightforward implementation of CGenFF's complicated atom typing rules and for equally straightforward updating of the atom typing scheme as the force field grows. The presented atom typer was validated by assigning correct atom types on 477 model compounds including in the training set as well as 126 test-set molecules that were constructed to specifically verify its different components. The program may be utilized via an online implementation at https://www.paramchem.org/ .

Automation of the CHARMM General Force Field (CGenFF) II: assignment of bonded parameters and partial atomic charges

Molecular mechanics force fields are widely used in computer-aided drug design for the study of drug candidates interacting with biological systems. In these simulations, the biological part is typically represented by a specialized biomolecular force field, while the drug is represented by a matching general (organic) force field. In order to apply these general force fields to an arbitrary drug-like molecule, functionality for assignment of atom types, parameters, and partial atomic charges is required. In the present article, algorithms for the assignment of parameters and charges for the CHARMM General Force Field (CGenFF) are presented. These algorithms rely on the existing parameters and charges that were determined as part of the parametrization of the force field. Bonded parameters are assigned based on the similarity between the atom types that define said parameters, while charges are determined using an extended bond-charge increment scheme. Charge increments were optimized to reproduce the charges on model compounds that were part of the parametrization of the force field. A "penalty score" is returned for every bonded parameter and charge, allowing the user to quickly and conveniently assess the quality of the force field representation of different parts of the compound of interest. Case studies are presented to clarify the functioning of the algorithms and the significance of their output data.

Development and current status of the CHARMM force field for nucleic acids

The CHARMM27 all-atom force field for nucleic acids represents a highly optimized model for investigations of nucleic acids via empirical force field calculations. The force field satisfactorily treats the A, B, and Z forms of DNA as well as RNA, and it also useful for nucleosides and nucleotides. In addition, it is compatible with the CHARMM force fields for proteins and lipids, allowing for simulation studies of heterogeneous systems.

Molecular dynamics simulation of unsaturated lipid bilayers at low hydration: parameterization and comparison with diffraction studies

A potential energy function for unsaturated hydrocarbons is proposed and is shown to agree well with experiment, using molecular dynamics simulations of a water/octene interface and a dioleoyl phosphatidylcholine (DOPC) bilayer. The simulation results verify most of the assumptions used in interpreting the DOPC experiments, but suggest a few that should be reconsidered. Comparisons with recent results of a simulation of a dipalmitoyl phosphatidylcholine (DPPC) lipid bilayer show that disorder is comparable, even though the temperature, hydration level, and surface area/lipid for DOPC are lower. These observations highlight the dramatic effects of unsaturation on bilayer structure.

CHARMM additive and polarizable force fields for biophysics and computer-aided drug design

BACKGROUND

Molecular Mechanics (MM) is the method of choice for computational studies of biomolecular systems owing to its modest computational cost, which makes it possible to routinely perform molecular dynamics (MD) simulations on chemical systems of biophysical and biomedical relevance.

SCOPE OF REVIEW

As one of the main factors limiting the accuracy of MD results is the empirical force field used, the present paper offers a review of recent developments in the CHARMM additive force field, one of the most popular biomolecular force fields. Additionally, we present a detailed discussion of the CHARMM Drude polarizable force field, anticipating a growth in the importance and utilization of polarizable force fields in the near future. Throughout the discussion emphasis is placed on the force fields' parametrization philosophy and methodology.

MAJOR CONCLUSIONS

Recent improvements in the CHARMM additive force field are mostly related to newly found weaknesses in the previous generation of additive force fields. Beyond the additive approximation is the newly available CHARMM Drude polarizable force field, which allows for MD simulations of up to 1μs on proteins, DNA, lipids and carbohydrates.

GENERAL SIGNIFICANCE

Addressing the limitations ensures the reliability of the new CHARMM36 additive force field for the types of calculations that are presently coming into routine computational reach while the availability of the Drude polarizable force fields offers an inherently more accurate model of the underlying physical forces driving macromolecular structures and dynamics. This article is part of a Special Issue entitled "Recent developments of molecular dynamics".

Intrinsic conformational properties of deoxyribonucleosides: implicated role for cytosine in the equilibrium among the A, B, and Z forms of DNA

Structural properties of biomolecules are dictated by their intrinsic conformational energetics in combination with environmental contributions. Calculations using high-level ab initio methods on the deoxyribonucleosides have been performed to investigate the influence of base on the intrinsic conformational energetics of nucleosides. Energy minima in the north and south ranges of the deoxyribose pseudorotation surfaces have been located, allowing characterization of the influence of base on the structures and energy differences between those minima. With all bases, chi values associated with the south energy minimum are lower than in canonical B-DNA, while chi values associated with the north energy minimum are close to those in canonical A-DNA. In deoxycytidine, chi adopts an A-DNA conformation in both the north and south energy minima. Energy differences between the A and B conformations of the nucleosides are <0.5 kcal/mol in the present calculations, except with deoxycytidine, where the A form is favored by 2.3 kcal/mol, leading the intrinsic conformational energetics of GC basepairs to favor the A form of DNA by 1.5 kcal/mol as compared with AT pairs. This indicates that the intrinsic conformational properties of cytosine at the nucleoside level contribute to the A form of DNA containing predominately GC-rich sequences. In the context of a B versus Z DNA equilibrium, deoxycytidine favors the Z form over the B form by 1.6 kcal/mol as compared with deoxythymidine, suggesting that the intrinsic conformational properties of cytosine also contribute to GC-rich sequences occurring in Z DNA with a higher frequency than AT-rich sequences. Results show that the east pseudorotation energy barrier involves a decrease in the furanose amplitude and is systematically lower than the inversion barrier, with the energy differences influenced by the base. Energy barriers going from the south (B form) sugar pucker to the east pseudorotation barrier are lower in pyrimidines as compared with purines, indicating that the intrinsic conformational properties associated with base may also influence the sugar pseudorotational population distribution seen in DNA crystal structures and the kinetics of B to A transitions. The present work provides evidence that base composition, in addition to base sequence, can influence DNA conformation.

Active site of human liver aldehyde dehydrogenase

Bromoacetophenone (2-bromo-1-phenylethanone) functions as an affinity reagent for human aldehyde dehydrogenase (EC 1.2.1.3) and has been found specifically to label a unique tryptic peptide in the enzyme. Amino-terminal sequence analysis of the labeled peptide after purification by two different procedures revealed the following sequence: Val-Thr-Leu-Glu-Leu-Gly-Gly-Lys. Radioactivity was found to be associated with the glutamate residue, which was identified as Glu-268 by reference to the known amino acid sequence. This paper constitutes the first identification of an active site of aldehyde dehydrogenase.

Structure, force, and energy of a double-stranded DNA oligonucleotide under tensile loads

The end-to-end stretching of a duplex DNA oligonucleotide has been studied using potential of mean force (PMF) calculations based on molecular dynamics (MD) simulations and atomic force microscopy (AFM) experiments. Near quantitative agreement between the calculations and experiments was obtained for both the extension length and forces associated with strand separation. The PMF calculations show that the oligonucleotide extends without a significant energetic barrier from a length shorter than A-DNA to a length 2.4 times the contour length of B-DNA at which the barrier to strand separation is encountered. Calculated forces associated with the barrier are 0.09 +/- 0.03 nN, based on assumptions concerning tip and thermal-activated barrier crossing contributions to the forces. Direct AFM measurements show the oligonucleotide strands separating at 2.6 +/- 0.8 contour lengths with a force of 0.13 +/- 0.05 nN. Analysis of the energies from the MD simulations during extension reveals compensation between increases in the DNA-self energy and decreases in the DNA-solvent interaction energy, allowing for the barrierless extension of DNA beyond the canonical B form. The barrier to strand separation occurs when unfavorable DNA interstrand repulsion cannot be compensated for by favorable DNA-solvent interactions. The present combination of single molecule theoretical and experimental approaches produces a comprehensive picture of the free energy surface of biological macromolecular structural transitions.

Progress toward chemical accuracy in the computer simulation of condensed phase reactions

We describe a procedure for the generation of chemically accurate computer-simulation models to study chemical reactions in the condensed phase. The process involves (i) the use of a coupled semiempirical quantum and classical molecular mechanics method to represent solutes and solvent, respectively; (ii) the optimization of semiempirical quantum mechanics (QM) parameters to produce a computationally efficient and chemically accurate QM model; (iii) the calibration of a quantum/classical microsolvation model using ab initio quantum theory; and (iv) the use of statistical mechanical principles and methods to simulate, on massively parallel computers, the thermodynamic properties of chemical reactions in aqueous solution. The utility of this process is demonstrated by the calculation of the enthalpy of reaction in vacuum and free energy change in aqueous solution for a proton transfer involving methanol, methoxide, imidazole, and imidazolium, which are functional groups involved with proton transfers in many biochemical systems. An optimized semiempirical QM model is produced, which results in the calculation of heats of formation of the above chemical species to within 1.0 kcal/mol (1 kcal = 4.18 kJ) of experimental values. The use of the calibrated QM and microsolvation QM/MM (molecular mechanics) models for the simulation of a proton transfer in aqueous solution gives a calculated free energy that is within 1.0 kcal/mol (12.2 calculated vs. 12.8 experimental) of a value estimated from experimental pKa values of the reacting species.

Ab initio conformational analysis of nucleic acid components: intrinsic energetic contributions to nucleic acid structure and dynamics

In recent years, the use of high-level ab initio calculations has allowed for the intrinsic conformational properties of nucleic acid building blocks to be revisited. This has provided new insights into the intrinsic conformational energetics of these compounds and its relationship to nucleic acids structure and dynamics. In this article we review recent developments and present new results. New data include comparison of various levels of theory on conformational properties of nucleic acid building blocks, calculations on the abasic sugar, known to occur in vivo in DNA, on the TA conformation of DNA observed in the complex with the TATA box binding protein, and on inosine. Tests of the Hartree-Fock (HF), second-order Møller-Plesset (MP2), and Density Functional Theory/Becke3, Lee, Yang and Par (DFT/B3LYP) levels of theory show the overall shape of backbone torsional energy profiles (for gamma, epsilon, and chi) to be similar for the different levels, though some systematic differences are identified between the MP2 and DFT/B3LYP profiles. The east pseudorotation energy barrier in deoxyribonucleosides is also sensitive to the level of theory, with the HF and DFT/B3LYP east barriers being significantly lower (approximately 2.5 kcal/mol) than the MP2 counterpart (approximately 4.0 kcal/mol). Additional calculations at various levels of theory suggest that the east barrier in deoxyribonucleosides is between 3.0 and 4.0 kcal/mol. In the abasic sugar, the west pseudorotation energy barrier is found to be slightly lower than the east barrier and the south pucker is favored more than in standard nucleosides. Results on the TA conformation suggest that, at the nucleoside level, this conformation is significantly destabilized relative to the global energy minimum, or relative to the A- and B-DNA conformations. Deoxyribocytosine would destabilize the TA conformation more than other bases relative to the A-DNA conformation, but not relative to the B-DNA conformation.

Targeting of an interrupted polypurine:polypyrimidine sequence in mammalian cells by a triplex-forming oligonucleotide containing a novel base analogue

The DNA triple helix consists of a third strand of nucleic acid lying in the major groove of an intact DNA duplex. The most stable triplexes form on polypurine:polypyrimidine sequences, and pyrimidine interruptions in the purine strand are destabilizing. Sequence stringency is imparted by specific Hoogsteen hydrogen bonds between third strand bases and the purine bases in the duplex. Appropriate base and sugar modifications of triple helix-forming oligonucleotides (TFOs) confer chromosome targeting activity in living cells. However, broad utilization of TFOs as gene targeting reagents in mammalian cells has been limited by the requirement for homopurine target sequences. Although there have been a number of base analogues described that appear to be promising as candidates for triplex target expansion, none has been examined in a biological system. We have employed a postsynthetic strategy to prepare a collection of TFOs with base analogues at a defined position. Following assessment of affinity for a triplex target with a single C:G inversion, TFOs with a second generation of analogues were synthesized. One of these, TFO-5a, with 2'-OMe-guanidinylethyl-5-methylcytosine at the position corresponding to the C:G interruption in the target sequence, was further modified to confer bioactivity. The activity of this TFO, linked to psoralen, was measured in a mammalian cell line that was engineered by directed sequence conversion to carry a triplex target with a single C:G interruption. TFO-5a was active against this target and inactive against the corresponding target with an uninterrupted polypurine:polypyrimidine sequence.

New insights into the structure of abasic DNA from molecular dynamics simulations

Abasic (AP) sites constitute a common form of DNA damage, arising from the spontaneous or enzymatic breakage of the N-glycosyl bond and the loss of a nucleotide base. To examine the effects of such damage on DNA structure, especially in the vicinity of the abasic sugar, four 1.5 ns molecular dynamics simulations of double-helical DNA dodecamers with and without a single abasic (tetrahydrofuran, X) lesion in a 5'-d(CXT) context have been performed and analyzed. The results indicate that the abasic site does not maintain a hole or gap in the DNA, but instead perturbs the canonical structure and induces additional flexibility close to the abasic site. In the apurinic simulations (i.e., when a pyrimidine is opposite the AP site), the abasic sugar flipped in and out of the minor groove, and the gap was water filled, except during the occurrence of a novel non-Watson-Crick C-T base pair across the abasic site. The apyrimidinic gap was not penetrated by water until the abasic sugar flipped out and remained extrahelical. Both AP helices showed kinks of 20-30 degrees at the abasic site. The Watson-Crick hydrogen bonds are more transient throughout the DNA double helices containing an abasic site. The abasic sugar displayed an unusually broad range of sugar puckers centered around the northern pucker. The increased motion of the bases and backbone near the abasic site appear to correlate with sequence-dependent helical stability. The data indicate that abasic DNA contorts more easily and in specific ways relative to unmodified DNA, an aspect likely to be important in abasic site recognition and hydrolysis.

(-)-3 beta-Substituted ecgonine methyl esters as inhibitors for cocaine binding and dopamine uptake

Ten 3 beta-ecgonine analogues were synthesized and characterized by 1H and 13C NMR, MS, and elemental analysis. The compounds were synthesized as (-)-stereoisomers from (-)-cocaine. These compounds were assessed for their ability to inhibit [3H]cocaine binding to rat striatal tissue and to inhibit [3H]DA uptake into rat striatal synaptosomes. In this series of compounds, the length of the spacer between the aryl group and the tropane skeleton ranged from 1 to 4 bond distances, and conformational flexibility of the linkage and orientation of the aryl ring system were controlled by various types of linkages. The most potent of the analogues was methyl-(1R-2-exo-3-exo)-8-methyl-3-(beta-styrenyl)-8-azabicyclo[3. 2.1] octane-2-carboxylate. One of the less potent compounds was found to inhibit [3H]cocaine binding and [3H]DA uptake with significantly different IC50 values, in contrast to 14 other 3 beta-substituted analogues. Molecular modeling and CoMFA analysis were used to obtain a rigorous structure-function relationship for the studied compounds. The results showed that the potencies of these 3 beta-substituted ecgonine methyl esters were dominated by steric effects and were acutely sensitive to the distance between the aryl ring and the tropane skeleton and to the orientation of the aryl ring system relative to the tropane skeleton. The current study provides a clearer picture of the shape and size of the putative hydrophobic binding pocket for the 3 beta substituent at the cocaine receptor as well as emphasizing the importance of a drug's free energy of solvation in obtaining structure-activity relationships.

Molecular dynamics simulations of ribonuclease T1: analysis of the effect of solvent on the structure, fluctuations, and active site of the free enzyme

Molecular dynamics simulations were performed on ribonuclease T1 (RNase T1; EC 3.1.27.3) to determine a structure for the free enzyme. Simulations starting with the X-ray coordinates for the 2'GMP-RNase T1 complex were done in vacuo and with an 18-A water ball around the active site using stochastic boundary conditions to understand the influence of water on both the structure and fluctuations of the enzyme. Removal of 2'GMP caused structural changes in the loop regions, including those directly interacting with the bound inhibitor in the crystal structure, while regions of secondary structure were less affected. The presence of solvent in the simulation damped the structural changes observed, which may be related to the use of full charges in both simulations. Fluctuations were also affected by the water, which generally increased both at the surface and in the interior of the protein. The active site in vacuo collapsed upon itself, forming a number of protein-protein hydrogen bonds leading to larger structural changes and lowered fluctuations while the presence of water kept the active site open, minimized structural changes, and increased fluctuations. Such fluctuations in the active site may be important for the binding of inhibitors or substrates to the enzyme. Lastly, results from the water simulation allow the prediction of a motion for a hypothetical tryptophan at position 45, which can ultimately be tested experimentally via time-resolved fluorescence using a site-specific mutant of the enzyme.

Protein dynamics. A time-resolved fluorescence, energetic and molecular dynamics study of ribonuclease T1

Studies using time-resolved fluorescence depolarization were performed on the internal motion of Trp 59 of ribonuclease T1 (EC 3.1.27.3) in the free enzyme, 2'-GMP-enzyme complex and 3'-GMP-enzyme complex. The Trp 59 motion was also studied in the free enzyme using molecular dynamics simulations. Energetic analysis of activation barriers to the Trp 59 motion was performed using both the transition state theory and Kramers' theory. The activation parameters showed a dependence on solvent viscosity indicating the transition state approach in aqueous solution to be inadequate. When taking solvent viscosity contributions into account agreement between the transition state and Kramers' theories was obtained. The results indicate the three enzyme forms to have different conformations with the free enzyme and 3'-GMP-enzyme complex being similar. Comparison of the experimental and theoretical results showed a good agreement on the Trp 59 motion in the free enzyme. Trp 59 appears to vibrate rapidly, with a relaxation time of the order of 1 ps, within free space in the protein matrix and to have a slower motion, with a relaxation time of the order of 100 ps, which is related to breathing of the surrounding protein matrix. Molecular dynamics results indicate high mobility in regions of the enzyme involved in the interaction with the guanine base of the inhibitor or substrate while much lower mobility occurred in residues involved in the catalytic mechanism of ribonuclease T1.

Human aldehyde dehydrogenase: kinetic identification of the isozyme for which biogenic aldehydes and acetaldehyde compete

Michaelis constants and maximal velocities for phenylacetaldehyde (a metabolite of phenylethylamine), 3,4-dihydroxyphenylacetaldehyde (a metabolite of dopamine), 5-hydroxyindole acetaldehyde (a metabolite of serotonin), and 3,4-dihydroxyphenylglycolaldehyde (a metabolite of epinephrine and norepinephrine) have been determined for both cytoplasmic (E1) and mitochondrial (E2) isozymes of human liver aldehyde dehydrogenase (EC 1.2.1.3). Kinetic constants with biogenic aldehydes have never been previously determined for individual homogeneous isozymes of aldehyde dehydrogenase from any species. Mathematical treatment of these constants suggests that competition with acetaldehyde during alcohol metabolism would severely inhibit dehydrogenation of biogenic aldehydes with the mitochondrial and not the cytoplasmic isozyme of human liver aldehyde dehydrogenase.

Human mitochondrial aldehyde dehydrogenase inhibition by diethyldithiocarbamic acid methanethiol mixed disulfide: a derivative of disulfiram

A derivative and possible physiological metabolite of disulfiram, diethyldithiocarbamic acid methanethiol mixed disulfide, is shown here for the first time to inactivate the mitochondrial low-Km isozyme of human aldehyde dehydrogenase (EC 1.2.1.3). By comparing inactivating effects of diethyldithiocarbamic acid mixed disulfides with thiols of increasing chain length evidence is provided that steric hindrance is the reason for lack of inhibition of the mitochondrial enzyme by disulfiram in vitro. Since methanethiol is a normal metabolite [(1983) Annu. Rev. Biochem. 52, 187-222] the results also suggest a mechanism by which aldehyde dehydrogenase is inhibited by disulfiram and diethyldithiocarbamic acid in vivo.

Correlation of loss of activity of human aldehyde dehydrogenase with reaction of bromoacetophenone with glutamic acid-268 and cysteine-302 residues. Partial-sites reactivity of aldehyde dehydrogenase

Bromoacetophenone (2-bromo-1-phenylethanone) has been characterized as an affinity reagent for human aldehyde dehydrogenase (EC 1.2.1.3) [MacKerell, MacWright & Pietruszko (1986) Biochemistry 25, 5182-5189], and has been shown to react specifically with the Glu-268 residue [Abriola, Fields, Stein, MacKerell & Pietruszko (1987) Biochemistry 26, 5679-5684] with an apparent inactivation stoichiometry of two molecules of bromoacetophenone per molecule of enzyme. The specificity of bromoacetophenone for reaction with Glu-268, however, is not absolute, owing to the extreme reactivity of this reagent. When bromo[14C]acetophenone was used to label the human cytoplasmic E1 isoenzyme radioactively and tryptic fragmentation was carried out, peptides besides that containing Glu-268 were found to have reacted with reagent. These peptides were purified by h.p.l.c. and analysed by sequencing and scintillation counting to quantify radioactive label in the material from each cycle of sequencing. Reaction of bromoacetophenone with the aldehyde dehydrogenase molecule during enzyme activity loss occurs with two residues, Glu-268 and Cys-302. The activity loss, however, appears to be proportional to incorporation of label at Glu-268. The large part of incorporation of label at Cys-302 occurs after the activity loss is essentially complete. With both Glu-268 and Cys-302, however, the incorporation of label stops after one molecule of bromoacetophenone has reacted with each residue. Reaction with other residues continues after activity loss is complete.

Reevaluation of stereoelectronic contributions to the conformational properties of the phosphodiester and N3'-phosphoramidate moieties of nucleic acids

The anomeric effect in the phosphodiester backbone of nucleic acids is a stereoelectronic effect that has conventionally been linked to interactions between lone pairs on the O(ester) atoms and P-O(ester) antibonding orbitals. The present study demonstrates that the anomeric effect in the phosphodiester backbone is significantly more complex than portrayed by this description. The presence of multiple lone pairs and antibonding orbitals around the phosphorus atom leads to additional contributions to the anomeric effect, especially involving the anionic oxygen lone pairs. On the basis of the structural changes and Natural Bond Orbital analysis it is shown that a complex balance between stereoelectronic effects involving both the ester and anionic oxygen lone pairs governs the conformational properties of the phosphodiester backbone. The N3'-phosphoramidate DNA backbone differs from the phosphodiester backbone due to the N3'-H moiety having only a single lone pair instead of the two lone pairs present on the O3' atom substituted. The present study uses N3'-phosphoramidate as a control to understand the changes in stereoelectronic effects as a result of changes in the structure and conformation. Two previously uncharacterized properties of the N3'-phosphoramidate backbone are also observed and explained through the complex balance of the postulated electronic delocalizations. The first observation is that the N3'-H moiety in N3'-phosphoramidate is a flexible moiety that can change the orientation of its hydrogen through inversion without a significant energetic penalty in both the gas phase and the aqueous phase. The second observation is that the stabilization of the C3'-endo conformation in N3'-phosphoramidate is primarily due to aqueous solvation rather than intrinsic gas-phase effects involving the reduced electronegativity of the 3'-substituent.

Molecular modeling and dynamics of neuropeptide Y

A combination of molecular modeling and molecular dynamics (MD) is used to determine a theoretical structure for neuropeptide Y (NPY). Starting with the X-ray structure for avian pancreatic polypeptide (APP), the substituted amino acids were mutated, the side chains oriented to local potential energy minima, and the entire structure minimized and subjected to an MD simulation. Comparison of the resulting NPY structure with APP X-ray and MD results showed secondary structural elements to be maintained and RMS fluctuations to be similar, although differences in both were observed. The approach presented offers a means to study the structure-function relationships of NPY and other similar polypeptides when combined with pharmacological measurements.

Molecular dynamics simulations of ribonuclease T1: comparison of the free enzyme and the 2' GMP-enzyme complex

Molecular dynamics simulations were performed on free RNase T1 and the 2'GMP-RNase T1 complex in vacuum and with water in the active site along with crystallographically identified waters, allowing analysis of both active site and overall structural and dynamics changes due to the presence of 2'GMP. Differences in the active site include a closing in the presence of 2'GMP, which is accompanied by a decrease in mobility of active site residues. The functional relevance of the active site fluctuations is discussed. 2'GMP alters the motion of Tyr-45, suggesting a role for that residue in providing a hydrophobic environment for the protein-nucleic acid interactions responsible for the specificity of RNase T1. The presence of 2'GMP causes a structural change of the C-terminus of the alpha-helix, indicating the transmission of structural changes from the active site through the protein matrix. Overall fluctuations of both the free and 2'GMP enzyme forms are in good agreement with X-ray temperature factors. The motion of Trp-59 is influenced by 2'GMP, indicating differences in enzyme dynamics away from the active site, with the calculated changes following those previously seen in time-resolved fluorescence experiments.

Analysis of structure-function relationships of neuropeptide Y using molecular dynamics simulations and pharmacological activity and binding measurements

Studies on the structure-function relationship of neuropeptide Y (NPY) were undertaken using a combination of in vacuo molecular dynamics (MD) simulations and pharmacological receptor binding and biological activity measurements. Following a conformational search of NPY from which a theoretical structure was determined, a study of the structural and dynamic changes in the region of amino acids 25-36 was performed in a variety of NPY fragments and in the NPY free acid. Results revealed an increased structural change as the fragment size was decreased. Also, the mobility appears to be lowest in the full NPY vs the NPY fragments. Pharmacological measurements showed a decreased receptor binding and biological activity as fragment size decreased. Combination of the two approaches suggests a model where conformational maintenance and low configurational entropy of the 25-36 region of NPY favors both receptor binding and biological activity. Furthermore, the possibility of two receptor interaction modes is suggested. Analysis of the NPY structure suggests the direct importance of the amidated C-terminus, Gln34 and His26, an indirect importance of the Tyr1 sidechain as well as the potential importance of an apparent electric 'dipole' in NPY for receptor binding and biological activity.

Bromoacetophenone as an affinity reagent for human liver aldehyde dehydrogenase

Human liver aldehyde dehydrogenase isozymes E1 and E2 (EC 1.2.1.3) are both completely and irreversibly inactivated by bromoacetophenone (2-bromo-1-phenylethanone). Steady-state kinetics with both acetophenone and chloroacetophenone indicated interaction with the same enzyme form as the aldehyde substrate. Saturation kinetics with chloroacetophenone and bromoacetophenone indicated interaction at a specific site on the enzyme surface and gave a dissociation constant similar to that from steady-state kinetics, suggesting that the same processes were being observed by both methods and that the active site may be involved. Protection against inactivation was afforded by chloral and NAD together. Stoichiometry of inactivation showed the first 2 equiv per tetramer to abolish the majority of catalytic activity; 4 equiv inactivated both isozymes with complete loss of esterase, NAD-stimulated esterase, and dehydrogenase activities. Peptide mapping of enzyme modified with [carbonyl-14C]bromoacetophenone of CNBr digests (E1) and tryptic digests (E1 and E2) showed one peptide to be preferentially labeled. The above results together with the similarity of bromoacetophenone to the substrate benzaldehyde suggest bromoacetophenone may react with a residue in the active site of aldehyde dehydrogenase. Amino acid analysis of the labeled E1 tryptic fragment indicated reaction with a different peptide from that with which iodoacetamide reacts.