X-Pro peptides. A theoretical study of the hydrogen bonded conformations of (alpha-aminoisobutyryl-L -prolyl)n sequences

Intramolecularly hydrogen bonded conformations of (Aib-Pro)n sequences have been analysed theoretically. Both 4-1 (C10 and 3-1 (C7 hydrogen bonded regular structures are shown to be stereochemically feasible. Conformational energies for the helical structures have been estimated using classical potential energy methods. Both C10 and C7 conformations have very similar energies. Pyrrolidine ring puckering has a pronounced effect on the energies, and only Cv-endo puckered Pro residues can be accommodated. The theoretical calculations using spectroscopic data suggest that the recently proposed novel 310 helical conformation for benzyloxycarbonyl(Aib-Pro)4-methyl ester is in solution, is indeed energetically and stereochemically favourable.

Hydrogen bonding with mono- and di-carbonyls

The electronic absorption and i.r. spectroscopic studies are reported for the hydrogen bonding systems involving alcohol and various ketones. It is shown that the hydrogen bonding abilities of ketones are determined by the extent of delocalization of the lone pair electrons in their non-bonding molecular orbitals. Evidence for the formation of very weak intermolecular hydrogen bonds between alcohol and the π-electron part of the dicarbonyls has also been presented from the i.r. studies in the 3400–3700 cm−1 region.

Interaction of Metal Ions with 2'-Deoxyribonucleotides. Crystal and Molecular Structure of a Cobalt(l1) Complex with 2'-Deoxyinosine 5'-Monophosphate

The crystal structure of the cobalt( 11) complex with 2'-deoxyinosine 5'-monophosphate (5'- dlMP), [Co(5'-dlMP) (H,0),]-2H20, has been analysed by X-ray diffraction. The complex crystallizes in the space group P2,2,2, with a = 6.877(3), b = 10.904(2), c = 25.421 (6) A, and Z = 4. The structure was solved by the heavy-atom method and refined to an R value of 0.043 using 1 776 unique reflections. The cobalt ion binds only to the 6-oxopurine base of the nucleotide at the N(7) position, the octahedral co-ordination of the metal being completed by five water oxygens. The phosphate oxygens are involved in hydrogen bonding with the co-ordinated water molecules. The structure is closely similar to that of the corresponding ribonucleotide complex. The nucleotide has the energetically preferred conformation: an anti base, a C(3') -endo sugar pucker, and a gauche-gauche conformation about the C(4')-C( 5') bond. The significance of sugar puckering in the monomeric complexes of general formula [ M (5'-nucleotide) (H20),] is explained in terms of the structural requirements for metal-water-phosphate bridging interactions.

Molecular structure of a cyclic tetrapeptide disulfide. A novel 310 helical conformation with an S-S bridge

The crystal structure of the cyclic peptide disulfide Boc-Cys-Pro-Aib-Cys-NHMe has been determined by X-ray diffraction. The peptide crystallizes in the space group P212121, with A = 8.646(1), B = 18.462(2), C = 19.678(3)Å and Z = 4. The molecules adopt a highly folded compact conformation, stabilized by two intramolecular 4→ 1 hydrogen bonds between the Cys (1) and Pro (2) CO groups and the Cys (4) and methylamide NH groups, respectively. The backbone conformational angles for the peptide lie very close to those expected for a 310 helix. The S-S bridge adopts a right handed twist with a dihedral angle of 82°. The structure illustrates the role of stereochemically constrained residues, in generating novel peptide conformations. Aib, α-aminoisobutyric acid; Z, benzyloxycarbonyl; Boc, t-butyloxycarbonyl; OMe, methyl ester; OBz, benzyl ester; NHMe, N-methylamide; Tosyl, p-toluenesulfonyl.

A vanadate-stimulated NADH oxidase in erythrocyte membrane generates hydrogen peroxide

Oxidation of NADH by rat erythrocyte plasma membrane was stimulated by about 50-fold on addition of decavanadate, but not other forms of vanadate like orthovanadate, metavanadate aad vanadyl sulphate. The vanadate-stimulated activity was observed only in phosphate buffer while other buffers like Tris, acetate, borate and Hepes were ineffective. Oxygen was consumed during the oxidation of NADH and the products were found to be NAD+ and hydrogen peroxide. The reaction had a stoichiometry of one mole of oxygen consumption and one mole of H2O2 production for every mole of NADH that was oxidized. Superoxide dismutase and manganous inhibited the activity indicating the involvement of superoxide anions. Electron spin resonance in the presence of a spin trap, 5, 5prime-dimethyl pyrroline N-oxide, indicated the presence of superoxide radicals. Electron spin resonance studies also showed the appearance of VIV species by reduction of VV of decavanadate indicating thereby participation of vanadate in the redox reaction. Under the conditions of the assay, vanadate did not stimulate lipid peroxidation in erythrocyte membranes. Extracts from lipid-free preparations of the erythrocyte membrane showed full activity. This ruled out the possibility of oxygen uptake through lipid peroxidation. The vanadate-stimulated NADH oxidation activity could be partially solubilized by treating erythrocyte membranes either with Triton X-100 or sodium cholate. Partially purified enzyme obtained by extraction with cholate and fractionation by ammonium sulphate and DEAE-Sephadex was found to be unstable.

Water-ethanol mixtures by molecular dynamics and x-ray Compton scattering

Water-ethanol mixtures are commonly used in industry and house holds. However, quite surprisingly their molecular-level structure is still not completely understood. In particular, there is evidence that the local intermolecular geometries depend significantly on the concentration. The aim of this study was to gain information on the molecular-level structures of water-ethanol mixtures by two computational methods. The methods are classical molecular dynamics (MD), where the movement of molecules can be studied, and x-ray Compton scattering, in which the scattering cross section is sensitive to the electron momentum density. Firstly, the water-ethanol mixtures were studied with MD simulations, with the mixture concentration ranging from 0 to 100%. For the simulations well-established force fields were used for the water and ethanol molecules (TIP4P and OPLS-AA, respectively). Moreover, two models were used for ethanol, rigid and non-rigid. In the rigid model the intramolecular bond lengths are fixed, whereas in the non-rigid model the lengths are determined by harmonic potentials. Secondly, mixtures with three different concentrations employing both ethanol models were studied by calculating the experimentally observable x-ray quantity, the Compton profile. In the MD simulations a slight underestimation in the density was observed as compared to experiment. Furthermore, a positive excess of hydrogen bonding with water molecules and a negative one with ethanol was quantified. Also, the mixture was found more structured when the ethanol concentration was higher. Negligible differences in the results were found between the two ethanol models. In contrast, in the Compton scattering results a notable difference between the ethanol models was observed. For the rigid model the Compton profiles were similar for all the concentrations, but for the non-rigid model they were distinct. This leads to two possibilities of how the mixing occurs. Either the mixing is similar in all concentrations (as suggested by the rigid model) or the mixing changes for different concentrations (as suggested by the non-rigid model). Either way, this study shows that the choice of the force field is essential in the microscopic structure formation in the MD simulations. When the sources of uncertainty in the calculated Compton profiles were analyzed, it was found that more statistics needs to be collected to reduce the statistical uncertainty in the final results. The obtained Compton scattering results can be considered somewhat preliminary, but clearly indicative of the behaviour of the water-ethanol mixtures when the force field is modified. The next step is to collect more statistics and compare the results with experimental data to decide which ethanol model describes the mixture better. This way, valuable information on the microscopic structure of water-ethanol mixtures can be found. In addition, information on the force fields in the MD simulations and on the ability of the MD simulations to reproduce the microscopic structure of binary liquids is obtained.

Vibration spectra of ferroelectric alkali trihydrogen selenites and the nature of the hydrogen bond potential in these crystals

Raman spectra of the ferroelectric LiH3 (SeO3)2 and NaH3(SeO3)2 and the anti-ferroelectric KH3 (SeO3)2 have been recorded at room temperature using a He-Ne and also an Ar-ion laser source. The infrared absorption spectra of these crystals and their deuterated analogues have been recorded in the region 400–4000 cm−1 both below and above the Curie temperature. From an analysis of the spectrum in the region 400–900 cm−1 it is concluded that (i) in LiH3 (SeO3)2 the protons are ordered in an asymmetric double minimum potential with a low barrier and the spectrum can be interpreted in terms of HSeO3− and H2SeO3 vibrations, (ii) in NaH3 (SeO3)2 all three protons occupy a single minimum potential at room temperature and below the transition temperature the groups HSeO3− and H2SeO3 are present, (iii) the proton at the inversion centre in KH3(SeO3)2 is in a broad troughed potential well and the low temperature spectrum is more likely to be due to H3SeO3+ and SeO32− species. This deviation of the spectrum from that of the previous two crystals is attributed to the difference in H-bond scheme and hence the absence of any cooperative motion of protons in this crystal.

Enzyme molecular choreography : studies on soluble inorganic pyrophosphatases

Inorganic pyrophosphatases (PPases, EC 3.6.1.1) hydrolyse pyrophosphate in a reaction that provides the thermodynamic 'push' for many reactions in the cell, including DNA and protein synthesis. Soluble PPases can be classified into two families that differ completely in both sequence and structure. While Family I PPases are found in all kingdoms, family II PPases occur only in certain prokaryotes. The enzyme from baker's yeast (Saccharomyces cerevisiae) is very well characterised both kinetically and structurally, but the exact mechanism has remained elusive. The enzyme uses divalent cations as cofactors; in vivo the metal is magnesium. Two metals are permanently bound to the enzyme, while two come with the substrate. The reaction cycle involves the activation of the nucleophilic oxygen and allows different pathways for product release. In this thesis I have solved the crystal structures of wild type yeast PPase and seven active site variants in the presence of the native cofactor magnesium. These structures explain the effects of the mutations and have allowed me to describe each intermediate along the catalytic pathway with a structure. Although establishing the ʻchoreographyʼ of the heavy atoms is an important step in understanding the mechanism, hydrogen atoms are crucial for the mechanism. The most unambiguous method to determine the positions of these hydrogen atoms is neutron crystallography. In order to determine the neutron structure of yeast PPase I perdeuterated the enzyme and grew large crystals of it. Since the crystals were not stable at ambient temperature, a cooling device was developed to allow neutron data collection. In order to investigate the structural changes during the reaction in real time by time-resolved crystallography a photolysable substrate precursor is needed. I synthesised a candidate molecule and characterised its photolysis kinetics, but unfortunately it is hydrolysed by both yeast and Thermotoga maritima PPases. The mechanism of Family II PPases is subtly different from Family I. The native metal cofactor is manganese instead of magnesium, but the metal activation is more complex because the metal ions that arrive with the substrate are magnesium different from those permanently bound to the enzyme. I determined the crystal structures of wild type Bacillus subtilis PPase with the inhibitor imidodiphosphate and an inactive H98Q variant with the substrate pyrophosphate. These structures revealed a new trimetal site that activates the nucleophile. I also determined that the metal ion sites were partially occupied by manganese and iron using anomalous X- ray scattering.

Molecular Modelling of Estrogen-Producing Enzymes CYP450 Aromatase and 17beta-Hydroxysteroid Dehydrogenase Type 1

Breast cancer is the most common cancer in women in the western countries. Approximately two-thirds of breast cancer tumours are hormone dependent, requiring estrogens to grow. Estrogens are formed in the human body via a multistep route starting from cholesterol. The final steps in the biosynthesis include the CYP450 aromatase enzyme, converting the male hormones androgens (preferred substrate androstenedione ASD) into estrogens(estrone E1), and the 17beta-HSD1 enzyme, converting the biologically less active E1 into the active hormone 17beta-hydroxyestradiol E2. E2 is bound to the nuclear estrogen receptors causing a cascade of biochemical reactions leading to cell proliferation in normal tissue, and to tumour growth in cancer tissue. Aromatase and 17beta-HSD1 are expressed in or near the breast tumour, locally providing the tissue with estrogens. One approach in treating hormone dependent breast tumours is to block the local estrogen production by inhibiting these two enzymes. Aromatase inhibitors are already on the market in treating breast cancer, despite the lack of an experimentally solved structure. The structure of 17beta-HSD1, on the other hand, has been solved, but no commercial drugs have emerged from the drug discovery projects reported in the literature. Computer-assisted molecular modelling is an invaluable tool in modern drug design projects. Modelling techniques can be used to generate a model of the target protein and to design novel inhibitors for them even if the target protein structure is unknown. Molecular modelling has applications in predicting the activities of theoretical inhibitors and in finding possible active inhibitors from a compound database. Inhibitor binding at atomic level can also be studied with molecular modelling. To clarify the interactions between the aromatase enzyme and its substrate and inhibitors, we generated a homology model based on a mammalian CYP450 enzyme, rabbit progesterone 21-hydroxylase CYP2C5. The model was carefully validated using molecular dynamics simulations (MDS) with and without the natural substrate ASD. Binding orientation of the inhibitors was based on the hypothesis that the inhibitors coordinate to the heme iron, and were studied using MDS. The inhibitors were dietary phytoestrogens, which have been shown to reduce the risk for breast cancer. To further validate the model, the interactions of a commercial breast cancer drug were studied with MDS and ligand–protein docking. In the case of 17beta-HSD1, a 3D QSAR model was generated on the basis of MDS of an enzyme complex with active inhibitor and ligand–protein docking, employing a compound library synthesised in our laboratory. Furthermore, four pharmacophore hypotheses with and without a bound substrate or an inhibitor were developed and used in screening a commercial database of drug-like compounds. The homology model of aromatase showed stable behaviour in MDS and was capable of explaining most of the results from mutagenesis studies. We were able to identify the active site residues contributing to the inhibitor binding, and explain differences in coordination geometry corresponding to the inhibitory activity. Interactions between the inhibitors and aromatase were in agreement with the mutagenesis studies reported for aromatase. Simulations of 17beta-HSD1 with inhibitors revealed an inhibitor binding mode with hydrogen bond interactions previously not reported, and a hydrophobic pocket capable of accommodating a bulky side chain. Pharmacophore hypothesis generation, followed by virtual screening, was able to identify several compounds that can be used in lead compound generation. The visualisation of the interaction fields from the QSAR model and the pharmacophores provided us with novel ideas for inhibitor development in our drug discovery project.

Conformation of the LL and LD hairpin bends with internal hydrogen bonds in proteins and peptides

The conformation of three linked peptide units having an internal 4 → 1 type of hydrogen bond has been studied in detail, and the low energy conformations are listed. These conformations all lead to the reversal of the chain direction, and may therefore be called as “hairpin bends” or “U-bends”. Since this bend can occur at the end of two chains hydrogen-bonded in the antiparallel β-conformation, it is also known as the “β-bend”. Two types of conformation are possible when the residues at the second and third Cα atoms are both of type L (the LL bend), while only one type is possible for the LD and the DL bend. The LL bend can also accommodate the sequences LG, GL, GG (G = glycine), while the LD bend can accommodate the sequences LG, GD and GG. The conformations for the sequences DD and DL are exact inverses (or mirror images) of those for the sequences LL and LD, respectively, and have dihedral angles (phi2, ψ2), (phi3, ψ3) of the same magnitudes, but of opposite signs as those for the former types, which are listed, along with the characteristics (length, angle and energy) of the hydrogen bonds. A comparison of the theoretical predictions with experimental data (from X-ray diffraction and NMR studies) on proteins and peptides, show reasonably good agreement. However, a systematic trend is observable in the experimental data, slightly deviating from theory, which indicates that some deformations occur in the shapes of the peptide units forming the bend, differing from that of the standard planar peptide unit.

The 310 Helical Conformation of the Amino Terminal Decapeptide of Suzukacillin. 270 MHz 'H NMR Evidence for Eight Intramolecular Hydrogen Bonds

The amino terminal suzukacillin decapeptide fragment, Boc-Aib-Pro-Val-Aib-Val-Ala-Aib-Ala-Aib-Aitbh-eO Me, two pentapeptides Boc-AibPrc-Val-AibVal-OMe and Boc-Ala-AibAla-AibAibOMe, and the tripeptide Boc-Ala-AibAibOMe have been studied by 270-MHz 'H NMR spectroscopy. By use of solvent dependence of chemical shifts in a CDC13-(CD3),S0 system and temperature dependence of amide NH chemical shifts in (CD3),S0, the intramolecularly hydrogen bonded NH groups in these peptides have been identified. The tripeptide possesses one hydrogen bond, both pentapeptides show evidence for three intramolecular hydrogen bonds, and the decapeptide has eight NH groups participating in hydrogen bonding. An Ala( 1)-Aib(2) @ turn is proposed for the tripeptide. Both pentapeptides favor 310 helical conformations composed of three consecutive B turns. The decapeptide adopts a 310 helical conformation with some flexibility at the Va1(5)-Ala(6) segment. The proposed conformations are consistent with the known stereochemical preferences of Aib residues.

X-Pro peptides. A theoretical study of the hydrogen bonded conformations of (a-aminoisobutyryl-L -prolyl)n sequences

Intramolecularly hydrogen bonded conformations of (Aib-Pro)n sequences have been analysed theoretically. Both 4�1 (C10 and 3�1 (C7 hydrogen bonded regular structures are shown to be stereochemically feasible. Conformational energies for the helical structures have been estimated using classical potential energy methods. Both C10 and C7 conformations have very similar energies. Pyrrolidine ring puckering has a pronounced effect on the energies, and only Cγ-endo puckered Pro residues can be accommodated. The theoretical calculations using spectroscopic data suggest that the recently proposed novel 310 helical conformation for benzyloxycarbonyl(Aib-Pro)4-methyl ester is in solution, is indeed energetically and stereochemically favourable.

Diagnostic Molecular Markers for Phosphine Resistance in U.S. Populations of Tribolium castaneum and Rhyzopertha dominica

Stored product beetles that are resistant to the fumigant pesticide phosphine (hydrogen phosphide) gas have been reported for more than 40 years in many places worldwide. Traditionally, determination of phosphine resistance in stored product beetles is based on a discriminating dose bioassay that can take up to two weeks to evaluate. We developed a diagnostic cleaved amplified polymorphic sequence method, CAPS, to detect individuals with alleles for strong resistance to phosphine in populations of the red flour beetle, Tribolium castaneum, and the lesser grain borer, Rhyzopertha dominica, according to a single nucleotide mutation in the dihydrolipoamide dehydrogenase (DLD) gene. We initially isolated and sequenced the DLD genes from susceptible and strongly resistant populations of both species. The corresponding amino acid sequences were then deduced. A single amino acid mutation in DLD in populations of T.castaneum and R.dominica with strong resistance was identified as P45S in T.castaneum and P49S in R.dominica, both collected from northern Oklahoma, USA. PCR products containing these mutations were digested by the restriction enzymes MboI and BstNI, which revealed presence or absence, respectively of the resistant (R) allele and allowed inference of genotypes with that allele. Seven populations of T.castaneum from Kansas were subjected to discriminating dose bioassays for the weak and strong resistance phenotypes. Application of CAPS to these seven populations confirmed the R allele was in high frequency in the strongly resistant populations, and was absent or at a lower frequency in populations with weak resistance, which suggests that these populations with a low frequency of the R allele have the potential for selection of the strong resistance phenotype. CAPS markers for strong phosphine resistance will help to detect and confirm resistant beetles and can facilitate resistance management actions against a given pest population.

Studies On Hydrogen Bonds: Part IV. Proposed Working Criteria for Assessing Qualitative Strength of Hydrogen Bonds

The data obtained in the earlier parts of this series for the donor and acceptor end parameters of N-H. O and O-H. O hydrogen bonds have been utilised to obtain a qualitative working criterion to classify the hydrogen bonds into three categories: "very good" (VG), "moderately good" (MG) and weak (W). The general distribution curves for all the four parameters are found to be nearly of the Gaussian type. Assuming that the VG hydrogen bonds lie between 0 and ± la, MG hydrogen bonds between ± 1 and ± 2, W hydrogen bonds beyond ± 2 (where is the standard deviation), suitable cut-off limits for classifying the hydrogen bonds in the three categories have been derived. These limits are used to get VG and MG ranges for the four parameters 1 and θ (at the donor end) and ± and ± (at the acceptor end). The qualitative strength of a hydrogen bond is decided by the cumulative application of the criteria to all the four parameters. The criterion has been further applied to some practical examples in conformational studies such as α-helix and can be used for obtaining suitable location of hydrogen atoms to form good hydrogen bonds. An empirical approach to the energy of hydrogen bonds in the three categories has also been presented.

Molecular structures of cytidine-5-diphosphate and cytidine-5-diphospho-choline, and their role in intermediary metabolism

The nucleotide coenzyme cytidine-5-diphospho-choline is highly folded. The CMP-5 parts of the molecules in the crystal structure are strongly linked by metal ligation and hydrogen bonds leaving the phosphoryl-choline residues relatively free. Cytidine-5-diphosphoric acid exists as a zwitterion with N31 protonated. The P−O bond lengths from the anhydride bridging oxygen in the pyrophosphate are significantly different.