Crystal structure and Hirshfeld surface analysis of 3-(cyclo-propyl-meth-oxy)-4-(di-fluoro-meth-oxy)-N-(pyridin-2-ylmeth-yl)benzamide

The title compound, C18H18F2N2O3, crystallizes with two independent mol-ecules (A and B) in the asymmetric unit. They differ essentially in the orientation of the pyridine ring with respect to the benzene ring; these two rings are inclined to each other by 53.3 (2)° in mol-ecule A and by 72.9 (2)° in mol-ecule B. The 3-(cyclo-propyl-meth-oxy) side chain has an extended conformation in both mol-ecules. The two mol-ecules are linked by a pair of C-H⋯O hydrogen bonds and two C-H⋯π inter-actions, forming an A-B unit. In the crystal, this unit is linked by N-H⋯O hydrogen bonds, forming a zigzag -A-B-A-B- chain along [001]. The chains are linked by C-H⋯N and C-H⋯F hydrogen bonds to form layers parallel to the ac plane. Finally, the layers are linked by a third C-H⋯π inter-action, forming a three-dimensional structure. The major contributions to the Hirshfeld surface are those due to H⋯H contacts (39.7%), followed by F⋯H/H⋯F contacts (19.2%).

Crystal structures of disordered Z-type helices

The X-ray crystal structures of a decamer sequence d(CGCGTACGCG)₂ and a tetradecamer sequence d(CGCGCGTACGCGCG)₂ are presented here. Both sequences are alternating pyrimidine-purine repeat sequences and they form disordered, pseudo-continuous left handed Z-type helices. They demonstrate interesting variants of the 'bundles of columns of helices' mode of packing.

Protein-small molecule docking with receptor flexibility in iMOLSDOCK

We have earlier reported the iMOLSDOCK technique to perform 'induced-fit' peptide-protein docking. iMOLSDOCK uses the mutually orthogonal Latin squares (MOLSs) technique to sample the conformation and the docking pose of the small molecule ligand and also the flexible residues of the receptor protein, and arrive at the optimum pose and conformation. In this paper we report the extension carried out in iMOLSDOCK to dock nonpeptide small molecule ligands to receptor proteins. We have benchmarked and validated iMOLSDOCK with a dataset of 34 protein-ligand complexes as well as with Astex Diverse dataset, with nonpeptide small molecules as ligands. We have also compared iMOLSDOCK with other flexible receptor docking tools GOLD v5.2.1 and AutoDock Vina. The results obtained show that the method works better than these two algorithms, though it consumes more computer time. The source code and binary of MOLS 2.0 (under a GNU Lesser General Public License) are freely available for download at https://sourceforge.net/projects/mols2-0/files/ .

[BCIgEPRED-a Dual-Layer Approach for Predicting Linear IgE Epitopes]

Allergy is a common health problem worldwide, especially food allergy. Since B cell epitopes that are recognized by the IgE antibodies act as antigenic determinants for allergy, they play a vital role in diagnostics. Hence, knowledge of an IgE binding epitope in a protein is of particular interest for identifying aller-genic proteins. Though IgE epitopes maybe conformational or linear, identification of the later is useful especially in food allergens that undergo processing or digestion. Very few computational tools are available for the prediction of linear IgE epitopes. Here we report a prediction system that predicts the exact linear IgE epitope. Since our earlier study on linear B cell epitope prediction demonstrated the effectiveness of using an exact epitope dataset (in contrast to epitope containing region datasets), the dataset in this study uses only experimentally verified exact IgE, IgG, IgM and IgA epitopes. Models for Support Vector Machine (SVM) and Random Forest (RF) were constructed adopting Dipeptide Deviation from the Expected mean (DDE) feature vector. Extensive validation procedures including five-fold cross validation and two different independent dataset tests have been performed to validate the proposed method, which achieved a balanced accuracy ranging from 74 to 78% with area under receiver operator curve greater than 0.8. Performance of the proposed method was observed to be better (accuracy difference of 16-28%) in comparison to the existing available method. The proposed method is developed as a standalone tool that could be used for predicting IgE epitopes as well as to be incorporated into any allergen prediction toolhttps://github.com/brsaran/BCIgePred.

Comparison of X-ray crystal structures of a tetradecamer sequence d(CCCGGGTACCCGGG)2 at 1.7 Å resolution

We present here a comparison of three different X-ray crystal structures of DNA tetradecamer sequence d(CCCGGGTACCCGGG)₂ all at about 1.7 Å resolution. The sequence was designed as an attempt to form a DNA four-way junction with A-type helical arms. However, in the presence of zinc, magnesium, and in the absence of any metal ion, it does not take up the junction structure, but forms an A-type double helix. This allowed us to study possible conformational changes in the double helix due to the presence of metal ions. Upon addition of the zinc ion, there is a change in the space group from P4₁2₁2 to P4₁. The overall conformation of the duplex remains the same. There are small changes in the interaction of the metal ions with the DNA. In the zinc-bound structure, there are two zinc ions that show direct interaction with the N7 atoms of terminal G₁₃ bases at either end of the molecule. There are small changes in the interhelical contacts. The consequence of these differences is to break some of the symmetry and change the space group.

iMOLSDOCK: Induced-fit docking using mutually orthogonal Latin squares (MOLS)

We have earlier reported the MOLSDOCK technique to perform rigid receptor/flexible ligand docking. The method uses the MOLS method, developed in our laboratory. In this paper we report iMOLSDOCK, the 'flexible receptor' extension we have carried out to the algorithm MOLSDOCK. iMOLSDOCK uses mutually orthogonal Latin squares (MOLS) to sample the conformation and the docking pose of the ligand and also the flexible residues of the receptor protein. The method then uses a variant of the mean field technique to analyze the sample to arrive at the optimum. We have benchmarked and validated iMOLSDOCK with a dataset of 44 peptide-protein complexes with peptides. We have also compared iMOLSDOCK with other flexible receptor docking tools GOLD v5.2.1 and AutoDock Vina. The results obtained show that the method works better than these two algorithms, though it consumes more computer time.

MOLS 2.0: software package for peptide modeling and protein-ligand docking

We previously developed an algorithm to perform conformational searches of proteins and peptides, and to perform the docking of ligands to protein receptors. In order to identify optimal conformations and docked poses, this algorithm uses mutually orthogonal Latin squares (MOLS) to rationally sample the vast conformational (or docking) space, and then analyzes this relatively small sample using a variant of mean field theory. The conformational search part of the algorithm was denoted MOLS 1.0. The docking portion of the algorithm, which allows only "flexible ligand/rigid receptor" docking, was denoted MOLSDOCK. Both are FORTRAN-based command-line-only molecular docking computer programs, though a GUI was developed later for MOLS 1.0. Both the conformational search and the rigid receptor docking parts of the algorithm have been extensively validated. We have now further enhanced the capabilities of the program by incorporating "induced fit" side-chain receptor flexibility for docking peptide ligands. Benchmarking and extensive testing is now being carried out for the flexible receptor portion of the docking. Additionally, to make both the peptide conformational search and docking algorithms (the latter including both flexible ligand/rigid receptor and flexible ligand/flexible receptor techniques) more accessible to the research community, we have developed MOLS 2.0, which incorporates a new Java-based graphical user interface (GUI). Here, we give a detailed description of MOLS 2.0. The source code and binary for MOLS 2.0 are distributed free (under a GNU Lesser General Public License) to the scientific community. They are freely available for download at https://sourceforge.net/projects/mols2-0/files/ .

Restricted mobility of side chains on concave surfaces of solenoid proteins may impart heightened potential for intermolecular interactions

Significant progress has been made in the determination of the protein structures with their number today passing over a hundred thousand structures. The next challenge is the understanding and prediction of protein-protein and protein-ligand interactions. In this work we address this problem by analyzing curved solenoid proteins. Many of these proteins are considered as "hub molecules" for their high potential to interact with many different molecules and to be a scaffold for multisubunit protein machineries. Our analysis of these structures through molecular dynamics simulations reveals that the mobility of the side-chains on the concave surfaces of the solenoids is lower than on the convex ones. This result provides an explanation to the observed preferential binding of the ligands, including small and flexible ligands, to the concave surface of the curved solenoid proteins. The relationship between the landscapes and dynamic properties of the protein surfaces can be further generalized to the other types of protein structures and eventually used in the computer algorithms, allowing prediction of protein-ligand interactions by analysis of protein surfaces.

Interactions of Mn2+ with a non-self-complementary Z-type DNA duplex

Crystal structures of the hexanucleotide d(CACGCG)·d(CGCGTG) were determined in two crystal lattices when different concentrations of the counterion Mn2+ were used in crystallization. The availability of Mn2+ during the crystallization process appears to play an important role in inducing different crystal packings that lead to crystals belonging to the two space groups P2(1) and P6(5). Analysis of the molecular interactions of Mn2+ with the Z-form duplexes shows direct coordination to the purine residues G and A.

Crystallographic and spectroscopic studies of d(CCGGTACCGG)

The decanucleotide sequence d(CCGGTACCGG) crystallizes as a four-way junction at low cobalt ion concentrations (i.e., 1 mM). When the cobalt concentration in the crystallization solution is increased to 5 mM, the sequence crystallizes as resolved B-DNA duplexes. Gel retardation studies of the decamer show both a faint slow-moving band and a much thicker fast-moving band at low cobalt ion concentrations, and only the intense fast-moving band at higher ion concentration. Circular dichroism (CD) spectroscopy of the decamer indicates a structural transition as the cobalt ion concentration in the solution is increased, probably from B-type to A-type DNA. These studies revealed that the oligomer sequence has several conformations and structures accessible to it, in a manner dependent on sequence, ion concentration, and DNA concentration. [Supplementary materials are available for this article. Go to the publisher's online edition of Nucleosides, Nucleotides & Nucleic Acids for the following free supplemental resources(s): Supplementary Figures 1, 2, and 3.].

Structure of d(CGGGTACCCG)4 as a four-way Holliday junction

The crystal structure of the decamer sequence d(CGGGTACCCG)(4) as a four-way Holliday junction has been determined at 2.35 Å resolution. The sequence was designed in order to understand the principles that govern the relationship between sequence and branching structure. It crystallized as a four-way junction structure with an overall geometry similar to those of previously determined Holliday junction structures.

Age, co-morbidity and poor mobility: no evidence of predicting in-patient death and acute hospital length of stay in the oldest old

BACKGROUND

The oldest old (aged over 90 years) are the fastest growing section of the UK population. Limited data exist regarding the effect of age, location, co-morbidity and physical performance status on outcome of acute illness in this age group.

METHODS

We performed a prospective study in people aged ≥ 90 years using hospital audit data in three hospitals in England and Scotland. We examined the characteristics of those admitted over three consecutive calendar months and calculated risk ratios of death and prolonged length of acute hospital stay (>7 days).

RESULTS

A total of 419 patients were included in this study (68% female, median age 93 years). There were similarities in presentation and diagnoses, but patients in Scotland (n = 164) were more likely to be admitted from sheltered housing or nursing homes than those in England (n = 255). Patients in England were significantly less likely to be able to mobilize < 10 m (41 vs. 34%, P < 0.001) but had lower prevalence of hypertension (40 vs. 55%, P = 0.02), ischaemic heart disease (30% vs. 45%, P = 0.02) and fewer prescribed medications (median 2 vs. 3, P < 0.001). Mortality was similar for the England and Scotland centres (P = 0.98). Previously recognized risk factors for death following hospital admission and length of stay e.g. older age, higher number of co-morbidities and poor mobility were not predictive in this study.

CONCLUSION

The 'oldest old' should not be considered as a homogenous group and findings from single-centre studies involving this age group may not be generalizable. We found no conclusive evidence that patient-related factors predict outcome in this age group in acute medical admission settings.

The sequence d(CGGCGGCCGC) self-assembles into a two dimensional rhombic DNA lattice

We report here the crystal structure of the partially self-complementary decameric sequence d(CGGCGGCCGC), which self assembles to form a four-way junction with sticky ends. Each junction binds to four others through Watson-Crick base pairing at the sticky ends to form a rhombic structure. The rhombuses bind to each other and form two dimensional tiles. The tiles stack to form the crystal. The crystal diffracted in the space group P1 to a resolution of 2.5Å. The junction has the anti-parallel stacked-X conformation like other junction structures, though the formation of the rhombic net noticeably alters the details of the junction geometry.

The structure of a full turn of an A-DNA duplex d(CGCGGGTACCCGCG)₂

We report the 2.6Å resolution crystal structure of the tetra-decamer d(CGCGGGTACCCGCG) in the tetragonal space group P4₃. This sequence contains the KpnI restriction site GGTACC in the centre which is flanked by alternating 'CG' sequences, and has a 'TA' step at the centre. These are features could favour the left-handed Z type helix. Despite this, overall the molecule has the A form. This is the first tetra-decamer crystallized in the A-DNA conformation, i.e. more than one full turn of the A helix. The crystallographic asymmetric unit consists of one tetra-decamer duplex. The helical twist and slide, as well as the base pair-base pair stacking interactions show alternations at the alternating pyrimidine-purine and purine-pyrimidine base steps. This variation is reminiscent of the dinucleotide repeat in left-handed Z-DNA helices. The crystal packing is unlike other A-DNA crystal structures, with each helix having a large number of contacts of many different types with symmetry-related neighbours.

Protein-ligand docking using mutually orthogonal Latin squares (MOLSDOCK)

The theoretical prediction of the association of a flexible ligand with a protein receptor requires efficient sampling of the conformational space of the ligand. Several docking methodologies are currently available. We have proposed a docking technique that performs well at low computational cost. The method uses mutually orthogonal Latin squares to efficiently sample the docking space. A variant of the mean field technique is used to analyze this sample to arrive at the optimum. The method has been previously applied to search through both the conformational space of a peptide as well its docking space. Here we extend this method to simultaneously identify both the low energy conformation as well as a high scoring docking mode for the small organic ligand molecules. Application of the method to 45 protein-ligand complexes, in which the number of rotatable torsions varies from 2 to 19, and comparisons with AutoDock 4.0, showed that the method works well.

The structure of d(CACACG).d(CGTGTG)

The crystal structure of d(CACACG).d(CGTGTG) was solved to a resolution of 2.05 A in space group P2(1). The duplex assumes the left-handed Z-DNA structure. The presence of two A.T base pairs in the hexamer does not greatly affect the conformation. The most significant changes compared with the regular structure of Z-DNA are in the values of twist in the central portion of the helix. This variation, as well as others in the values of roll, inclination etc., follow the pattern observed previously in the structure of d(CGCACG).d(CGTGCG).

Hexammineruthenium(III) ion interactions with Z-DNA

The hexamer duplex d(CGCGCA).d(TGCGCG) was crystallized with hexammineruthenium(III) ions in an orthorhombic space group; the crystals diffracted to 1.54 A resolution. Strong ion interactions with the adenine base induce a tautomeric shift from the amino to the imino form. Consequently, the A.T base pairing is disrupted. This structural study may be relevant to metal toxicity.

Exploring conformational space using a mean field technique with MOLS sampling

The computational identification of all the low energy structures of a peptide given only its sequence is not an easy task even for small peptides,due to the multiple-minima problem and combinatorial explosion. We have developed an algorithm, called the MOLS technique,that addresses this problem, and have applied it to a number of different aspects of the study of peptide and protein structure. Conformational studies of oligopeptides, including loop sequences in proteins have been carried out using this technique. In general the calculations identified all the folds determined by previous studies,and in addition picked up other energetically favorable structures. The method was also used to map the energy surface of the peptides. In another application, we have combined the MOLS technique, using it to generate a library of low energy structures of an oligopeptide, with a genetic algorithm to predict protein structures. The method has also been applied to explore the conformational space of loops in protein structures.Further, it has been applied to the problem of docking a ligand in its receptor site, with encouraging results.

Exploring the conformational space of protein loops using a mean field technique with MOLS sampling

We have recently developed a computational technique that uses mutually orthogonal Latin square sampling to explore the conformational space of oligopeptides in an exhaustive manner. In this article, we report its use to analyze the conformational spaces of 120 protein loop sequences in proteins, culled from the PDB, having the length ranging from 5 to 10 residues. The force field used did not have any information regarding the sequences or structures that flanked the loop. The results of the analyses show that the native structure of the loop, as found in the PDB falls at one of the low energy points in the conformational landscape of the sequences. Thus, a large portion of the structural determinants of the loop may be considered intrinsic to the sequence, regardless of either adjacent sequences or structures, or the interactions that the atoms of the loop make with other residues in the protein or in neighboring proteins.

Effect of DNA structural flexibility on promoter strength--molecular dynamics studies of E. coli promoter sequences

To study possible correlations between promoter activity and the structural flexibility of the DNA helix, we have carried out unrestrained molecular dynamics simulations of the -10 consensus region sequence and five variants forming the -10 region of various Escherichia coli promoter sequences. Analyses of the trajectories obtained from the simulations show that the consensus sequence has a pattern of two structurally flexible nucleotide steps sandwiched between two stiff steps. In the other sequences, this pattern varies in consonance with the change in the sequence. The variations in the patterns show correlation with the promoter strength.

Protein structure prediction using mutually orthogonal Latin squares and a genetic algorithm

We combine a new, extremely fast technique to generate a library of low energy structures of an oligopeptide (by using mutually orthogonal Latin squares to sample its conformational space) with a genetic algorithm to predict protein structures. The protein sequence is divided into oligopeptides, and a structure library is generated for each. These libraries are used in a newly defined mutation operator that, together with variation, crossover, and diversity operators, is used in a modified genetic algorithm to make the prediction. Application to five small proteins has yielded near native structures.

Principal component analysis of DNA oligonucleotide structural data

The microstructure of a DNA helix is characterized by several base pair and base step parameters such as twist, rise, roll, propeller twist, etc., in addition to conformational parameters such as the backbone and the glycosidic torsion angles. Among these only a few, which are independent of all others and of each other, may be used to precisely characterize the helix. The problem however is to identify these independent parameters. We have used principal component analysis to identify a relatively small set of independent parameters, with which to characterize each DNA helix. We show that these principal components clearly discriminate between A and B DNA helical types. The calculations further suggest that the microstructure of a DNA helix is better characterized using dinucleotides.

Structure of d(TGCGCG).d(CGCGCA) in two crystal forms: effect of sequence and crystal packing in Z-DNA

The sequence d(TGCGCG).d(CGCGCA) crystallized in two crystal forms, orthorhombic and hexagonal, in the presence of cobalt hexammine chloride, a known inducer of the left-handed Z-form of DNA. The crystal structures have been solved and refined at 1.71 A resolution in space group P2(1)2(1)2(1) and 2.0 A resolution in space group P6(5). The orthorhombic structure contains one Z-DNA hexamer duplex, while the hexagonal structure contains two hexamer duplexes in the structure. Of the latter, one is situated on a crystallographic sixfold screw axis, leading to disorder. This paper reports the effects of sequence and crystal packing on the structure of Z-type DNA. The structures lend additional support to the authors' earlier conclusion that a stretch of four C.G base pairs is sufficient to nucleate and define the regular model of the left-handed helix based on the structure of d(CGCGCG)2.

Space-group degeneracy in the packing of a non-selfcomplementary Z-DNA hexamer

The X-ray diffraction pattern of the crystals of the non-selfcomplementary hexadeoxyribonucleotide d(CGCACG).d(CGTGCG) can be indexed in four different space groups: (i) P6(5) and P2(1), with cell parameters a = 17.75 (1), b = 17.76 (1), c = 42.77 (3) A, alpha = 90, beta = 90, gamma = 120 degrees, and (ii) P2(1)2(1)2(1) and C2, with cell parameters a = 17.75 (1), b = 30.74 (2), c = 42.77 (3) A, alpha = 90, beta = 90, gamma = 90 degrees. While the R(merge) for the equivalent reflections in the different space groups indicates that P2(1) is the correct choice in the present case, it is demonstrated that the near degeneracy of the space groups arises out of the fact that the DNA molecule is nearly cylindrical. A perfect cylinder would show perfect degeneracy.

Cobalt hexammine induced tautomeric shift in Z-DNA: the structure of d(CGCGCA)*d(TGCGCG) in two crystal forms

We report here the crystal structure of the DNA hexamer duplex d(CGCGCA)*d(TGCGCG) at 1.71 A resolution. The crystals, in orthorhombic space group, were grown in the presence of cobalt hexammine, a known inducer of the left-handed Z form of DNA. The interaction of this ion with the DNA helix results in a change of the adenine base from the common amino tautomeric form to the imino tautomer. Consequently the A:T base pair is disrupted from the normal Watson-Crick base pairing to a 'wobble' like base pairing. This change is accommodated easily within the helix, and the helical parameters are those expected for Z-DNA. When the cobalt hexammine concentration is decreased slightly in the crystallization conditions, the duplex crystallizes in a different, hexagonal space group, with two hexamer duplexes in the asymmetric unit. One of these is situated on a crystallographic 6-fold screw axis, leading to disorder. The tautomeric shift is not observed in this space group. We show that the change in inter-helix interactions that lead to the two different space groups probably arise from the small decrease in ion concentration, and consequently disordered positions for the ion.

Conformational studies on enkephalins using the MOLS technique

Conformational studies of two linear enkephalin molecules, Met-enkephalin and Leu-enkephalin, have been carried using the mutually orthogonal Latin squares (MOLS) technique with the ECEPP/3 force field. This technique was developed recently in our laboratory to perform an unbiased search of the conformational space of peptides and to locate low energy conformations. The present study identified all the folds predicted by other studies, and in addition picked up other energetically favorable structures. The results suggest that the peptide backbone exists as a mixture of folded and unfolded forms (approximately 50% each). The study also provides information on the distribution of the low energy conformations that we have classified on the basis of structural motifs, backbone hydrogen-bonding patterns, and root mean square deviations in atomic positions.

An application of experimental design using mutually orthogonal Latin squares in conformational studies of peptides

We address the question-can we use experimental design methods to investigate peptide conformation and identify conformational parameters that may contribute more significantly to the potential energy than others? We used mutually orthogonal Latin square design to sample the conformational space of peptides and analysed the samples using analysis of variance. We examined the equality of the effect of the torsion angles on the conformational potential energy. The results showed that different torsion angles contributed differently to the conformational energy. We are able to identify those parameters that may have to be more carefully considered in conformational studies of peptides.

Correlations between nucleotide frequencies and amino acid composition in 115 bacterial species

We studied the correlations between amino acid composition and mononucleotide and dinucleotide frequencies in 115 bacterial genomes of varying G+C content. Observed amino acid frequencies were compared with those expected from the actual mononucleotide and dinucleotide frequencies. Both mononucleotide and dinucleotide frequencies correlate well with the amino acid frequency, with dinucleotide frequencies doing so better. Despite the strong correlations, some of the observed amino acid frequencies, in particular for Arg, Val, Asp, Glu, Ser, and Cys, were consistently different from predicted values in all genomes. We suggest that this variation from predicted values is a consequence of selection pressure at the level of amino acids, while the close correspondence to the predictions in residues such as Thr, Phe, Lys, and Asn arises only from mutation and selection pressure at the level of the nucleic acid sequences.

Enhanced sampling of the molecular potential energy surface using mutually orthogonal latin squares: application to peptide structures

The computational identification of the optimal three-dimensional fold of even a small peptide chain from its sequence, without reference to other known structures, is a complex problem. There have been several attempts at solving this by sampling the potential energy surface of the molecule in a systematic manner. Here we present a new method to carry out the sampling, and to identify low energy conformers of the molecule. The method uses mutually orthogonal Latin squares to select (of the order of) n(2) points from the multidimensional conformation space of size m(n), where n is the number of dimensions (i.e., the number of conformational variables), and m specifies the fineness of the search grid. The sampling is accomplished by first calculating the value of the potential energy function at each one of the selected points. This is followed by analysis of these values of the potential energy to obtain the optimal value for each of the n-variables separately. We show that the set of the n-optimal values obtained in this manner specifies a low energy conformation of the molecule. Repeated application of the method identifies other low energy structures. The computational complexity of this algorithm scales as the fourth power of the size of the molecule. We applied this method to several small peptides, such as the neuropeptide enkephalin, and could identify a set of low energy conformations for each. Many of the structures identified by this method have also been previously identified and characterized by experiment and theory. We also compared the best structures obtained for the tripeptide (Ala)(3) by the present method, with those obtained by an exhaustive grid search, and showed that the algorithm is successful in identifying all the low energy conformers of this molecule.

Structure of d(TGCGCA)2 at 293 K: comparison of the effects of sequence and temperature

The crystal structure of a hexameric DNA fragment with the sequence d(TGCGCA)(2) has been solved and refined at 293 K at a resolution of 1.64 A. The molecule adopts a left-handed Z-type helical conformation which is common in alternating pyrimidine-purine sequences. The presence of A.T base pairs at the two terminals does not perturb the structure to any great degree. However, several sequence-specific microstructural changes are noticeable. The structure of the identical sequence determined at 120 K involving somewhat different crystallization conditions has been reported previously [Harper et al. (1998), Acta Cryst. D54, 1273-1284]. A comparison of the present structure with that at low temperature and with that of d(CGCGCG)(2) shows that the effect of the change in sequence is greater than the combined effect of changes in temperature and environment.

Crystal structure at 1.63 A resolution of the native form of porcine beta-trypsin: revealing an acetate ion binding site and functional water network

The active center of a serine protease is the catalytic triad composed of His-57, Ser-195 and Asp-102. The existing crystal structure data on serine proteases have not fully answered a number of fundamental questions relating to the catalytic activity of serine proteases. The new high resolution native porcine beta-trypsin (BPT) structure is aimed at extending the knowledge on the conformation of the active site and the ordered water structure within and around the active site. The crystal structure of BPT has been determined at 1.63 A resolution. An acetate ion bound at the active site of a trypsin molecule by both classical hydrogen bonds and C-HellipsisO hydrogen bonds has been identified for the first time. A large network of water molecules extending from the recognition amino acid Asp-184 to the entry of the active site has been observed in the BPT structure. A detailed comparison with inhibitor complexes and autolysates indicates that the sulfate ion and the acetate ion bind at the same site of the trypsin molecule. The Ser-195 Cbeta-Ogamma-His-57 Nepsilon angle in the catalytic triad of BPT is intermediate between the corresponding values of the complex and native structure due to acetate ion binding. The network of waters from the recognition amino acid to the active site entry is probably the first ever complete picture of functional waters around the active site. Structural comparisons show that the functional waters involved in the binding of small molecule inhibitors and protease inhibitors are distinctly different.

Structure of d(CACGCG).d(CGCGTG) in crystals grown in the presence of ruthenium III hexammine chloride

Hexammine ions are strong inducers of the transition from the B-form to the left-handed Z-form in DNA. Here the structure of d(CACGCG). d(CGCGTG) obtained from crystals grown from a drop containing [Ru(NH3)6]Cl3 is reported. The structure is clearly characterized as Z-DNA. When compared with the structure of d(CACGCG).d(CGCGTG)/MgCl2 and that of d(CGCGCG)2, subtle differences are seen, most noticeably in the water structure. Since stable well diffracting crystals grow easily in the presence of [Ru(NH3)6]Cl3 and since this ion is not visible in the electron density it is concluded that the ion plays a non-specific role in stabilizing Z-DNA.

Sequence-dependent microheterogeneity of Z-DNA: the crystal and molecular structures of d(CACGCG).d(CGCGTG) and d(CGCACG).d(CGTGCG)

We have solved and refined the crystal structures of the first two non-self-complementary hexadeoxyribonucleotide duplexes with Watson-Crick base pairs, namely d(CACGCG).d(CGCGTG) and d(CGCACG).d(CGTGCG). Both the hexamers crystallize in the left-handed Z-DNA conformation. The packing of the molecules is similar in the two crystals: the hexamers are stacked on top of each other to form columns of infinite length, which are arranged in a close-packed hexagonal pattern. In spite of the similar packing, crystals of the first duplex belong to the space group P2(1)2(1)2(1), isomorphous to previous Z-DNA hexamers, while crystals of the other duplex belong to the space group P2(1), which has not been observed before in this context. The molecular structures of the two duplexes are also different. The first is very similar in conformation to the other Z-DNA hexamers and to the idealized ZI model. The other duplex shows large differences in the stacking and the relative disposition of its bases, leading to its possible classification as a novel Z-DNA variant. We conclude that the differences in the structure under similar environmental conditions point to a sequence-dependent plasticity in the DNA molecule which is visible even under close-packed conditions.