The effect of pressure on the crystal structure of hexagonal L-cystine

The crystal structure of hexagonal L-cystine has been determined at room temperature at pressures between 0.4 and 3.7 GPa; unit-cell dimensions were measured up to 6.4 GPa. The structure of this phase consists of molecules in their zwitterionic form, and crystallizes in the hexagonal space group P6₁22. The structure consists of hydrogen-bonded layers which are strongly reminiscent of those seen in α-glycine, and consist of (16) hydrogen-bonded ring motifs. These layers are connected on one side by the disulfide bridges within the cystine molecules, and on the other by NH⋯O hydrogen bonds to other glycine-like layers. The most compressible unit-cell dimension, and the direction of greatest strain in the structure, is along the c-axis, and application of pressure pushes the layers closer together. The compression occurs approximately equally in the regions of the interlayer hydrogen bonds and the disulfide bridges; in the latter, changes in the C—S—S—C torsion angles allow the cystine molecules to act like springs. The effects of pressure can be interpreted in terms of closing-up of voids in the structure, and this leads to (i) a lessening of the N—C—C—O and C—S—S—C torsional angles, (ii) shortening of the N—H⋯O hydrogen bonds by 0.10–0.60 Å and (iii) a further shortening of an already short S⋯S contact from 3.444 (4) Å to 3.264 (4) Å.

1. Introduction

The application of pressure is becoming an ever-more popular technique with which to analyse the nature of intermolecular interactions [for example, Oswald, Allan, Motherwell & Parsons (2005), Oswald, Allan, Day, Motherwell & Parsons (2005), Wunschel et al. (2003), Boldyreva (2004a,b), Allan et al. (1999, 2001)]. Application of pressure (up to ∼10 GPa, approximately 10⁵ atm) to non-cubic organic crystals commonly results in the anisotropic closure of voids in the structure along with an increase in the number and strength of weak interactions, such as CH⋯O hydrogen bonds. In some cases the discovery of new phases of materials is observed with increasing pressure (Moggach et al., 2005). In α-amino acids the compression of numerous weak CH⋯O interactions may be extremely important, as they have a supporting role to medium-strength hydrogen bonds under ambient pressure conditions, e.g. NH⋯O (Desiraju & Steiner, 1999; Derewenda et al., 1995). The compression of such interactions has been discussed in other studies, such as the compression of L-serine-I (Moggach et al., 2005).

Most molecular systems that have been studied previously are relatively simple with small unit cells with volumes ranging from a few hundred up to a few thousand ų. In this study the effect of pressure on hexagonal L-cystine (I) was studied on station 9.8 at the CCLRC Daresbury Laboratory (Cernik et al., 1997) as part of our on-going research of the effect of pressure on amino acids and related compounds. Although L-cystine is a moderately small organic molecule containing 14 non-H atoms, it crystallizes in a high symmetry space group (P6₁22), with an unusually long axis (c = 55.9 Å). Attempts to study the compound at high pressure in the home laboratory were frustrated by difficulties in resolving data along c*. The low divergence of synchrotron radiation, however, enabled this difficulty to be overcome, and we report here a compression study of L-cystine carried out using the high-pressure facilities on station 9.8 at SRS Daresbury Laboratory. Structural data to 3.7 GPa, and unit-cell dimensions up to 6.4 GPa, are reported.

L-Cystine is a dimer of L-cysteine produced by oxidation of the side-chain thiols to form the so-called disulfide bridge. The disulfide bridge is a structurally important feature in proteins, being the only common covalent linkage to occur within and between polypeptide chains. Naturally occurring disulphide bridges are found in many proteins that are secreted by cells where their presence adds stability to the generally fragile protein structure. Common examples include the defensive immunoglobulin proteins, the hydrolytic enzymes like trypsin and lysozyme, and structural proteins like keratin. Analyses of the conformations adopted by cystine (Thornton, 1981; Richardson, 1981; Morris et al., 1992; Görbitz, 1990) reveals that left- (—C—S—S—C— = χ3 = −90°) and right-handed disulphides are equally distributed in globular proteins. The distributions of the χ1 and χ2 angles are also more or less as expected. χ1 favours the −60° conformation since there is less crowding between the S and the adjacent peptide carbonyl; χ2 tends to adopt a value of around −85°. A more recent analysis, including that of the disulphide bridge environment by Bhattacharaya et al. (2004), confirms and extends these results using data from 1266 polypeptide chains in the Protein Data Bank (Berman et al., 2000) with R-factors ≤ 20% and data resolution ≤ 2 Å. It includes an analysis of the S to main-chain (carbonyl) oxygen distances that determines the average distance to be 3.6 (2) Å, about a quarter of which are interactions with their own carbonyl oxygen.

2. Experimental

2.1. Crystal growth

All starting materials were purchased from Sigma-Aldrich and used as received. L-Cystine (1.48 g, 6.16 mmol) was dissolved in warm ammonium hydroxide solution (ACS reagent, NH₃ 28.0–30.0%, 8 cm³) and crystallized on cooling to ambient temperature (Dahaoui et al., 1999). A high-quality hexagonally shaped crystal of dimensions 0.1 mm × 0.2 mm × 0.2 mm was then loaded into a Merrill–Bassett diamond anvil cell (DAC; Merrill & Bassett, 1974).

2.2. High-pressure crystallography

2.2.1. General procedures

High-pressure experiments were carried out with a Merrill–Bassett DAC equipped with 600 µm culet-cut diamonds and a tungsten gasket. The sample and a chip of ruby (as a pressure calibrant) were loaded into the DAC with a 1:1 mixture of pentane and isopentane as a hydrostatic medium. The ruby fluorescence method was utilized to measure the pressure (Piermarini et al., 1975).

2.2.2. Data collection, reduction and refinement

A sphere of data was collected on a crystal of L-cystine at ambient temperature and pressure in order to provide data for comparison with the high-pressure studies, which were also performed at ambient temperature (see below). Diffraction data were collected on a single crystal of L-cystine on a Bruker SMART APEX diffractometer with graphite-monochromated Mo Kα radiation (λ = 0.71073 Å). These data were integrated using the program SAINT (Bruker-AXS, 2003), while the absorption correction was carried out using the program SADABS (Sheldrick, 2004). Refinement was carried out against |F|² using all data (CRYSTALS; Betteridge et al., 2003) starting from the low-temperature (110 K) coordinates of Dahaoui et al. (1999). The final conventional R-factor was 0.0476 for 898 data. Listings of crystal and refinement data are given in Table 1.¹

Table 1
Crystallographic data for hexagonal L-cystine at ambient temperature and pressure, and at increasing pressures

	Pressure (GPa)
	0	0.4	1.4	2.3	3.7
Crystal data
Chemical formula	C6H12N2O4S2	C6H12N2O4S2	C6H12N2O4S2	C6H12N2O4S2	C6H12N2O4S2
Mr	240.30	240.30	240.30	240.30	240.30
Cell setting, space group	Hexagonal, P6122	Hexagonal, P6122	Hexagonal, P6122	Hexagonal, P6122	Hexagonal, P6122
a, c (Å)	5.4203 (5), 55.980 (12)	5.4030 (1), 55.984 (4)	5.3524 (2), 55.087 (14)	5.3141 (1), 54.291 (3)	5.26630 (10), 53.549 (5)
V (Å^3)	1424.3 (4)	1415.40 (10)	1366.7 (4)	1327.75 (8)	1286.17 (13)
Z	6	6	6	6	6
Dx (Mg m^−3)	1.681	1.691	1.752	1.803	1.861
Radiation type	Mo Kα	Synchrotron	Synchrotron	Synchrotron	Synchrotron
No. of reflections for cell parameters	1792	3156	1652	3365	3341
θ range (°)	9–46	9–54	9–53	9–53	9–51
μ (mm^−1)	0.55	0.56	0.57	0.59	0.61
Temperature (K)	293	293	293	293	293
Crystal form, colour	Hexagonal block, colourless	Hexagonal block, colourless	Hexagonal block, colourless	Hexagonal block, colourless	Hexagonal block, colourless
Crystal size (mm)	0.32 × 0.17 × 0.15	0.20 × 0.20 × 0.10	0.20 × 0.20 × 0.10	0.20 × 0.20 × 0.10	0.20 × 0.20 × 0.10
Data collection
Diffractometer	Bruker APEX	Bruker APEX II	Bruker APEX II	Bruker APEX II	Bruker APEX II
Data collection method	ω	ω	ω	ω	ω
Absorption correction	Multi-scan (based on symmetry-related measurements)	Multi-scan (based on symmetry-related measurements)	Multi-scan (based on symmetry-related measurements)	Multi-scan (based on symmetry-related measurements)	Multi-scan (based on symmetry-related measurements)
Tmin	0.81	0.71	0.62	0.64	0.56
Tmax	0.92	0.95	0.94	0.94	0.94
No. of measured, independent and observed parameters	16236, 965, 898	4719, 539, 490	4645, 492, 464	4229, 492, 438	4229, 467, 439
Criterion for observed reflections	I > 2.00 u(I)	I > 2.00 u(I)	I > 2.00 u(I)	I > 2.00 u(I)	I > 2.00 u(I)
Completeness (%)	99.0	45.4	45.5	45.5	45.9
Rint	0.085	0.047	0.047	0.053	0.063
θmax (°)	26.4	28.5	28.1	28.1	27.9
Range of h, k, l	–6 → h → 6	–7 → h → 7	–6 → h → 6	–7 → h → 6	–6 → h → 6
	–6 → k → 6	–7 → k → 6	–7 → k → 7	–7 → k → 6	–6 → k → 6
	–64 → l → 68	–29 → l → 29	–28 → l → 28	–28 → l → 28	–28 → l → 28
Refinement
Refinement on	F^2	F^2	F^2	F^2	F^2
R[F^2 > 2σ(F^2)], wR(F^2), S	0.048, 0.105, 1.07	0.035, 0.089, 1.06	0.036, 0.089, 1.05	0.048, 0.112, 1.05	0.093, 0.214, 1.09
No. of reflections	956	513	485	458	451
No. of parameters	64	64	64	64	34
H-atom treatment	Not refined	Not refined	Not refined	Not refined	Not refined
Weighting scheme	Calculated w = 1/[σ^2(F^2) + 0.03 + 3.29P] where P = [max(F o^2,0) + 2F c^2]/3	Calculated w = 1/[σ^2(F^2) + 0.05 + 0.8P] where P = [max(F o^2,0) + 2F c^2]/3	Calculated w = 1/[σ^2(F^2) + 0.05 + 1.06P] where P = [max(F o^2,0) + 2F c^2]/3	Calculated w = 1/[σ^2(F^2) + 0.03 + 3.33P] where P = [max(F o^2,0) + 2F c^2]/3	Calculated w = 1/[σ^2(F^2) + 0.00 + 18.85P] where P = [max(F o^2,0) + 2F c^2]/3
(Δ/σ)max	0.001	<0.0001	<0.0001	<0.0001	<0.0001
Δρmax, Δρmin (e Å^–3)	0.35, −0.35	0.14, −0.19	0.22, −0.17	0.27, −0.21	0.54, −0.60

High-pressure diffraction data were collected with synchrotron radiation on a Bruker APEX II diffractometer at the CCLRC Daresbury Laboratory on station 9.8 (λ = 0.6774 Å); see below. Data were collected in ω-scans in eight settings of 2θ and φ with a frame and step size of 1s and 0.2°, respectively. This data collection strategy was based on that described by Dawson et al. (2004). The data were integrated using the program SAINT using `dynamic masks' to avoid integration of regions of the detector shaded by the body of the pressure cell (Dawson et al., 2004). Additional masks were also added to mask out the three most intense powder rings (100, 002 and 101) from the beryllium backing discs. This additional mask was required because the low divergence of synchrotron radiation gave the beryllium powder pattern a much more textured spotty appearance than if collected using a sealed-tube laboratory X-ray source. Application of this additional mask removed the possibility of unreliable measurements of sample reflections which coincided with the beryllium rings and was shown to improve the merging statistics slightly. Absorption corrections for the DAC and sample were carried out using the programs SHADE (Parsons, 2004) and SADABS, respectively. Data collections were taken in approximately 1.0 GPa steps from 0.4 GPa up to a final pressure of 6.4 GPa.

Refinements of L-cystine were carried out starting from the published low-temperature coordinates (Dahaoui et al., 1999). Refinements were carried out against |F|² using all data (CRYSTALS). Because of the low completeness of the data sets, all 1,2 and 1,3 distances were restrained to the values observed in the ambient pressure structure. Specifically, the applied restraints were as follows: distances (Å) S1—C1 1.818 (20); C3—O2 1.267 (20); C3—O1 1.245 (20); C3—C2 1.538 (20); C2—N1 1.480 (20); C2—C1 1.523 (20); angles (°) O1—C3—O2 125.73 (100); O1—C3—C2 117.25 (100); O2—C3—C2 117.01 (100); N1—C2—C1 112.05 (100); N1—C2—C3 109.67 (100); C1—C2—C3 113.39 (100). H atoms attached to carbon and nitrogen were placed geometrically and not refined. The numbering scheme used is the same as Cambridge Structural Database (CSD) refcode LCYSTI14 (Dahaoui et al., 1999).

All non-H atoms were refined with anisotropic displacement parameters to 2.3 GPa. At 3.7 GPa only the S atoms were refined anisotropically, while all other non-H atoms were refined with isotropic displacement parameters. This parameterization strategy was dictated by the quality of the data collected: reflections broadened with increasing pressure, making indexing and integration increasingly difficult. Acceptable refinements could be obtained for the final 5.0 and 6.4 GPa data sets; unit-cell dimensions at these pressures were a = 5.21970 (10), c = 52.793 (6) Å and a = 5.1852 (3), c = 52.172 (7) Å, respectively. The deterioration of the data quality was the result of increasing mosaic spread with increasing pressure, a frequent problem with crystallographic pressure studies on soft materials. Listings of crystal and refinement data are given in Table 1.

2.2.3. High-pressure data collection on station 9.8, SRS, Daresbury Laboratory

High-pressure single-crystal diffraction experiments may be carried out at SRS, Daresbury Laboratory on stations 9.8 and 16.2SMX. Both facilities are equipped with Bruker D8 three-circle fixed-χ goniometers with APEXII detectors. On most currently available models of home-laboratory X-ray diffractometer it is necessary to re-machine the collimator and beam stop in order that full φ-rotations may be performed. This is not necessary on a synchrotron beamline because the fixed distance constraints of tube-to-collimator housing are absent. In `high-pressure mode' at the SRS the collimator is retracted towards the X-ray source by 30 mm into the space where the tube housing and optics would usually sit. This increases the collimator-to-sample distance to allow full 360° rotation in φ; rotations in ω are restricted only by the position of the detector. The distance from sample to beam stop a custom beam was set at 70 mm. Alignment procedures follow Dawson et al. (2004). We have found that the ease and simplicity of set-up combined with small beam divergence and high brightness allows high pressure to be a routine option for users of the facilities at the SRS.

2.2.4. Software for structure analysis

Crystal structures were visualized using the programs CAMERON (Watkin et al., 1993), DIAMOND (Crystal Impact, 2004), MERCURY (Bruno et al., 2002) and XP (Sheldrick, 1997). Analyses were carried out using PLATON (Spek, 2004), as incorporated in the WIN-GX suite (Farrugia, 1999). Searches of the CSD (Allen, 2002; Allen & Motherwell, 2002) were performed using the program CONQUEST and version 5.26 of the database with updates up to November 2004.

Topological calculations of void distributions were carried out using TOPOS (Blatov et al., 2000; Blatov & Shevchenko, 2003). Considerable simplification of the void distributions can be gained by clustering; voids were therefore clustered using what the program calls the `floating' method with the `size' parameter specified as 1.0.

Equation-of-state calculations were carried out using EOSFIT (Angel, 2002).

3. Results

3.1. The structure of hexagonal L-cystine at ambient pressure

L-Cystine crystallizes in two polymorphic forms: a tetragonal phase (P4₁) that was characterized by Chaney & Steinrauf (1974) and a hexagonal phase (P6₁22) first investigated by Oughton & Harrison (1959) but more recently in a charge density study by Dahaoui et al. (1999). We have re-determined the structure at room temperature (Fig. 1) in order to be able to compare all structures obtained here under similar conditions of temperature. Though we quote our own room-temperature structural parameters, our interpretation of the structure owes much to the paper of Dahaoui et al. (1999), which, in addition to describing the structure of cystine at high resolution, also contains a database survey of geometries and interactions of crystal structures containing —CSSC— moieties.

Figure 1 The molecular structure of L-cystine at ambient temperature and pressure showing atom labelling. Ellipsoids are drawn at the 30% probability level, and H atoms are shown as spheres of arbitrary radius.

Both the tetragonal and hexagonal phases of L-cystine crystallize with the molecule in its zwitterionic form. The S—S bond distances in the two polymorphs [2.042 (6) and 2.043 (2) Å, respectively] do not differ significantly from the average for such bonds in the CSD [2.039 (2) Å]. The same comment applies to the C—S—S—C torsion angle, which is positive in both forms [69.3 (2)° and 75.18 (5)°, respectively]. The hydrochloride and hydrobromide salts of L-cystine have both been structurally characterized in the anhydrous forms and as dehydrates; interestingly, these all have negative C—S—S—C torsion angles in the range −79° to −83°. The length of S—S bonds in compounds containing C—S—S—C moieties has been shown to be dependent on the torsion angle about the S—S bond: torsion angles which approach 90° result in short S—S bonds because π-interactions are optimized and repulsions between S-based lone pairs are minimized (Dahaoui et al., 1999). The values of χ1 and χ2 (i.e. the torsion angles N1—C2—C1—S1 and C2—C1—S1—S1′) are 56.1 (5)° and 80.1 (3)°, respectively; these values are relatively less common in proteins (see Introduction).

In this study the effect of pressure on the crystal structure of the hexagonal phase was studied. The principal intermolecular interactions take the form of unusually short S⋯S contacts and NH⋯O hydrogen bonds, formed between the ammonium and carboxylate moieties. The NH⋯O hydrogen bonds lead to the formation of layers which lie perpendicular to the c-axis. The layers can be considered to be formed from the intersection of two primary C(5) chain motifs (Bernstein et al., 1995), one via N1—H4⋯O1 interactions, the other via N1—H5⋯O2 interactions (Fig. 2a). The first of these interactions is substantially bifurcated (Jeffrey & Maluszynska, 1982) with donor-to-acceptor distances (i.e. N⋯O) in N1—H4⋯O1/O2 measuring 2.792 (5)/3.147 (5) Å.

Figure 2 (16) ring motifs in (a) hexagonal L-cystine and (b) α-glycine, as viewed perpendicular to the c- and b-axes, respectively; the figures correspond to N⋯O distances in N—H⋯O hydrogen bonds. (c) (16) ring motifs created by interlaced CH⋯O and NH⋯O interactions in cystine as viewed perpendicular to the c-axis. The black dotted lines represent the ring motif as created by NH⋯O hydrogen bonds [cf. (a)], while the orange dotted lines show the interlaced ring motif as created by CH⋯O interactions. In (a) and (c) only the basic α-amino acid unit involved within the ring motif is shown for the sake of clarity. Colour scheme: C, grey; H, white; O, red; N, blue.

The structure of the layers is strongly reminiscent of layers formed perpendicular to the b-axis in α-glycine (Fig. 2b; Boldyreva et al., 2003). In both cases intersection of the C(5) chains leads to formation of a secondary-level (16) ring motif. In the asymmetric unit of LCYSTI14, and all the structures reported here, the chains based on N1—H4⋯O1 and N1—H5⋯O2 interactions are respectively formed along lattice repeats in the [010] and [110] directions, but these directions interchange for successive layers generated by the 6₁ screw axis. Head-to-tail NH⋯OOC motifs are often observed in amino acid structures and are usually associated with a cell dimension of ∼5.5 Å. This is also observed in L-cystine with an a-axis length of 5.4203 (5) Å.

The hydrogen-bonded layers are linked to other such layers on one side by covalent S—S bridges and on the other by NH⋯O hydrogen bonds (Fig. 3). The latter are formed via N1—H6⋯O1 interactions, which connect pairs of molecules in successive layers, forming an R₂²(10) ring motif (Fig. 4).

Figure 3 Hydrogen-bonded glycine-like double layers linked down the c-axis by disulfide bridges. This view is along a, with the c-axis running from left to right. The dashed black lines correspond to N⋯O hydrogen bonds; the dashed orange lines are S⋯S contacts. Sulfur is shown in yellow, but otherwise the colour scheme is as Fig. 2.

Figure 4 R2^2(10) ring motif as viewed perpendicular to the c-axis. (a) Standard ball-and-stick model; (b) space-filling model of (a) on the same scale. Half of each L-cystine molecule has been removed for clarity. Note the apparent lack of free space at the centre of the R2^2(10) ring.

Weak CH⋯O hydrogen bonds are also present at ambient pressure and can play an important role in stabilizing medium-strength hydrogen bonds, e.g. NH⋯O (Desiraju & Steiner, 1999; Derewenda et al., 1995). Under ambient pressure two such interactions are present, both of which occur within the glycine-like layers already described. The first of these, C1—H1⋯O2, runs approximately parallel to N1—H4⋯O1 and bisects the (16) ring motif described above, while the second, C2—H3⋯O1, runs in the same direction as N1—H5⋯O2. If all other hydrogen-bonding interactions are ignored, these weak C—H⋯O interactions combine to produce an (16) ring motif of their own (Fig. 2c).

At ambient pressure, S⋯S intermolecular interactions are formed (Fig. 3); the S⋯S distance [3.444 (4) Å] is much shorter than the sum of the van der Waals radii of the S atoms (3.7 Å). At 110 K this distance is 3.4264 (4) Å (Dahaoui et al., 1999). These interactions run perpendicular to the c-axis direction, and the S⋯S contact is approximately co-linear with the S—S covalent bond. Trends in the geometries of such interactions have been shown to be consistent with donation of a lone pair into an S—S-based σ* molecular orbital (Dahaoui et al., 1999).

3.2. The response of hexagonal L-cystine to pressure up to 3.7 GPa

Hexagonal L-cystine was found to be stable up to 6.4 GPa. However, no structural data could be extracted above 5.0 GPa, and even at this pressure the refined parameters are imprecise. Therefore only structural data to 3.7 GPa were used for comparison with the ambient pressure structure. The S—S distance appears to decrease from 2.043 (2) to 2.038 (5) Å, though this is not statistically significant. Other bond distances and angles were restrained during refinement to their ambient pressure values, and so no conclusions about their variation with pressure should be drawn from these data. Significant changes do occur in the torsion angles, however. The C—S—S—C torsion angle reduces from 75.20 (7)° at ambient pressure to 72.3 (4)° at 3.7 GPa, while the N1—C2—C3—O2 torsion angle (also referred to as NCαCO or ψ in macromolecular crystallography) also reduces from 166.12 (8)° to 164.1 (7)°.

3.2.1. Lattice parameters

The response of the lattice parameters of hexagonal L-cystine to pressure is anisotropic, with the principal component of the strain tensor acting in the c-axes direction. Between ambient pressure and 6.4 GPa the c-axis reduces in length by 6.8%; the reduction in the a-axis length is 4.4%. The volume change in the initial 0.4 GPa range is approximately 1.0%, while from 0.4 to 6.4 GPa the volume changes by approximately 2.8% every 1.0 GPa until the sample became polycrystalline.

Changes in unit-cell dimensions and volume as a function of pressure are plotted in Fig. 5. The strain induced within molecular systems by pressure is usually anisotropic, and in layered structures the greatest direction of strain is often observed to be normal to the layers. Such behaviour has been observed, for example, in the monoclinic and orthorhombic polymorphs of paracetamol (Boldyreva et al., 2000, 2002) and is also observed in L-cystine, where the direction of greatest compressibility is normal to the glycine-like layers formed within the structure.

Figure 5 Variation of lattice parameters (a, c) (Å) and volume (Å^3) of hexagonal L-cystine as a function of pressure (GPa).

The bulk modulus (K₀), refined for a Birch–Murnaghan equation of state (Birch, 1947; Angel, 2004) to second order, is 29.1 (4) GPa. The data set used to calculate this quantity is admittedly rather limited, and the values of V₀, K′ and K′′ were fixed at 1424.3 ų, 4 and −0.1337 GPa⁻¹, respectively; the values of V₀ and K′ did not change significantly from these values on refinement. Molecular solids typically have K₀ < 30 GPa (Angel, 2004), and the value obtained for L-cystine-I is towards the high end of this range. Slebodnick et al. (2004) quote the following K₀ values which are useful for comparison: Ru₃(CO)₁₂ 6.6 GPa, NaCl 25 GPa, quartz 37 GPa, ceramics 50–300 GPa and diamond 440 GPa.

3.2.2. NH⋯O hydrogen bonds

Variation in hydrogen-bonding parameters between ambient and 3.7 GPa are given in Table 2. H atoms were placed in idealized positions during refinement, and so we quantify hydrogen bonding with increasing pressure by using hydrogen bond donor-to-acceptor (N⋯O) distances rather than H⋯O distances.

Table 2
Hydrogen-bonding parameters in hexagonal L-cystine

Distances are in Å, of which only D⋯A distances are given (see Experimental).

	Pressure (GPa)
	0	0.4	1.4	2.3	3.7
N1H4⋯O1^i
N1⋯O1	2.792 (5)	2.777 (2)	2.732 (2)	2.693 (4)	2.643 (7)
N1H4⋯O2^i
N1⋯O2	3.147 (5)	3.119 (2)	3.071 (3)	3.029 (3)	2.970 (7)
N1H5⋯O2^ii
N1⋯O2	2.816 (4)	2.806 (3)	2.769 (3)	2.742 (4)	2.703 (8)
N1H6⋯O1^iii
N1⋯O1	2.842 (4)	2.829 (5)	2.815 (5)	2.781 (8)	2.786 (14)
C1H1⋯O2^i
C1⋯O2	3.259 (6)	3.235 (3)	3.190 (3)	3.153 (5)	3.124 (9)
C2H3⋯O1^iv
C2⋯O1	3.499 (4)	3.480 (2)	3.407 (2)	3.355 (3)	3.285 (6)
Symmetry codes: (i) x, −1 + y, z; (ii) 1 + x, y, z; (iii) 1 + x − y, 2 − y, −z; (iv) −1 + x, −1 + y, z.

The most compressible of the NH⋯O hydrogen bonds was found to be the longer interaction in the bifurcated N1—H4⋯O1/O2 hydrogen bond. Between 0 and 3.7 GPa N1⋯O2 shortens from 3.147 (5) to 2.970 (7) Å at 3.7 GPa. The shorter of the hydrogen bonds, N1—H4⋯O1, also decreases in length from 2.792 (5) Å to 2.643 (7) Å. This head-to-tail hydrogen-bonding interaction is also observed in α-glycine and serine, but in these structures it is less compressible than observed in this study, shortening from 2.7703 (9) to 2.740 (7) Å in α-glycine at 6.2 GPa, and from 2.871 (3) to 2.775 (13) Å in serine at 4.8 GPa.

The other hydrogen bond formed within the glycine-like layers, N1—H5⋯O2, decreases in length from 2.816 (4) Å to 2.690 (18) Å. A search of the CSD reveals that there are no ambient pressure amino acid structures (out of 213 hits) with N—H⋯O interactions shorter than N1⋯O1 at 3.7 GPa, and only three amino acid structures shorter than N1⋯O2 at 3.7 GPa. The shortest interaction at ambient pressure [2.661 (2) Å] was observed in L-arginine L-glutamate trihydrate (DUSMAF; Suresh et al., 1986).

The fourth hydrogen bond, N1—H6⋯O1, which interacts between the layers to form an R₂²(10) ring motif, decreases in length from 2.842 (4) Å to 2.786 (14) Å at 3.7 GPa. The O⋯C contacts formed between the carboxyl groups across the rings are of a similar length at ambient pressure [3.243 (1) Å] to those observed in dipolar C=O⋯C=O contacts (2.92–3.32 Å), though the geometry does not match any of the common types identified by Allen et al. (1998). These contacts also shorten only modestly to 3.200 (14) Å at 3.7 GPa. For comparison, the sheared-parallel C=O⋯C=O contacts formed in acetone at 0 and 1.5 GPa are 3.587 (3) and 3.365 (2) Å in acetone-II (Allan et al., 1999).

This is the least compressible of the head-to-tail NH⋯OOC interactions even though it is formed in the same direction as the most compressible unit-cell dimension, the c-axis.

3.2.3. CH⋯O hydrogen bonds

Within the glycine-like layers the compression of weak C—H⋯O interactions parallels the compression of the N—H⋯O interactions. The more compressible is C2—H3⋯O1, which shortens from C⋯O = 3.499 (4) Å to 3.285 (6) Å. This interaction runs approximately parallel to the N1—H5⋯O2 hydrogen bond. C1—H1⋯O2 (formed in the same direction as N1—H4⋯O1) shortens from 3.259 (6) Å to 3.124 (9) Å.

3.2.4. S⋯S interactions

The structure of L-cystine is notable for the formation of S⋯S contacts which are substantially shorter [3.444 (4) Å] than the sum of the van der Waals radii of the S atoms (3.6 Å; Bondi, 1964). Under pressure, this contact shortens to 3.263 (4) Å at 3.7 GPa. A search of the CSD reveals 19 entries containing a C(sp³)—S—S—C(sp³) fragment in which an S⋯S contact is observed between 2.5 and 3.7 Å. All of these structures were determined at ambient pressure. The distribution is rather uniform in this distance range, but the shortest contact observed so far (refcode KOWMEO; Fujimura et al., 1992) measures 3.31 Å. The application of pressure has therefore caused the S⋯S interaction to reduce in length close to the lower distance limit observed under ambient pressure conditions.

4. Discussion

4.1. Anisotropic compression of L-cystine

In pressure studies of α-glycine and L-serine-I we have ascribed trends in the relative compressibilities of C—H⋯O and N—H⋯O hydrogen bonds to the closing up of voids existing in the structure at ambient pressure. At pressure these voids contract with the shortening of N—H⋯O and C—H⋯O interactions. Similar conclusions on the importance of void closure were drawn in a study of Ru₃(CO)₁₂. Recent work by Blatov & Shevchenko (2003) makes it possible to analyse the sizes and distributions of voids in crystal structures using the Voronoi-Dirichlet formalism.

The environment of an atom in a crystal structure can be visualized using a Voronoi-Dirichlet polyhedron or VDP (Peresypkina & Blatov, 2000a,b). Voronoi-Dirichlet analysis is a method for partitioning space amongst points which occupy that space. A point is separated from a neighbouring point by a plane which bisects the vector between them. This construction is repeated for every pair of points to yield a subdivision of the space into cells which each contain one point. VDP analysis carried out using individual atoms to define the points leads to a molecular VDP, and these correspond to a complete partitioning of space in a crystal structure. Interatomic interactions occur at the points in the structure where faces of VDPs meet, whereas the topological definition of a void is a point at which the vertices of VDPs meet (Blatov & Shevchenko, 2003).

Fig. 6(a) shows the distribution of voids in the crystal structure of L-cystine at ambient pressure. The largest void conglomerate in the structure occurs in the vicinity of the S atoms; this is labelled 1 in Fig. 6(a). These voids straddle the short intermolecular S⋯S contact of 3.444 (4) Å, and the contraction of this distance to 3.263 (4) Å at 3.7 GPa is consistent with the orientation of the void conglomerates in this region of the structure.

Figure 6 (a) The distribution of voids in cystine under ambient conditions, viewed with the c-axis vertical. The numbers 1–4 refer to regions of the structure, as discussed in the text. (b) Voids in the glycine-like (16) layers under ambient conditions. The colour scheme represents the volumes of the VDPs of the voids (in Å^3) as discussed by Blatov & Shevchenko (2003); the largest void is shown in white, voids then progressively decrease in volume from blue to red; precise values are given in the legend.

Fig. 6(b) shows the distribution of voids in the region of the (16) rings formed in the glycine-like layers; this region is labeled 2 in Fig. 6(a). Fig. 7 illustrates the sizes of the voids in cystine at ambient pressure; the close correspondence between the positions of the voids in Fig. 7(a) and the largest void positions in Fig. 6(b) is readily apparent. On compression to 3.7 GPa the sizes of the voids in the (16) ring are reduced, leading to the observed shortening in the CH⋯O and NH⋯O hydrogen bonds which form within these layers. It is notable that the S⋯S interaction is approximately parallel to the N1—H4⋯O1 C(5) chains, and the change of −0.18 Å in the former is reflected in the −0.15 Å change in the latter.

Figure 7 Space-filling plots showing (16) ring-motifs L-cystine at ambient temperature and pressure (a), 2.3 GPa (b) and 3.7 GPa (c), and α-glycine at ambient temperature and pressure (d) and at 2 GPa (e). Note the reduction in the sizes of the voids in the middle of the ring motifs with increasing pressure.

There is a collection of voids between the glycine-like sheets in the region labeled 3 in Fig. 6(a). Voids are less evident within the R₂²(10) ring motifs, which connect pairs of glycine-like layers, than in the (16) rings formed within the layers (see also Fig. 4). This may explain why the N1—H6⋯O1 hydrogen bond is the least compressible of the N—H⋯O interactions in the system. Larger voids do exist elsewhere in region 3 though, as shown in Fig. 6(a). There is a marked contraction of the void between the two O atoms labeled O2 in Fig. 8, with the O⋯O distance reducing from 3.643 (4) to 3.350 (15) Å, between ambient pressure and 3.7 GPa. The reduction of the N1—C2—C3—O2 torsion angle from 166.12 (8)° to 164.1 (7)° may also be associated with the void closure in this region of the structure. Overall the distance between the planes containing the C2 (or C_α) atoms in neighbouring glycine-like layers (labelled A in Fig. 9) reduces from 3.729 (2) Å at ambient pressure to 3.539 (21) Å at 3.7 GPa.

Figure 8 Voids formed between glycine-like layers at ambient pressure (a and b) and at 3.7 GPa (c). Note the reduction in the sizes of the voids in (b) and (c) with increasing pressure.

Figure 9 Ball-and-stick models of L-cystine at (a) ambient pressure and (b) 3.7 GPa as viewed down (). The red planes are through the Cα atoms in each glycine-like layer running perpendicular to the c-axes of each plot. The overall compression of L-cystine caused via the `twisting' of the molecules can be seen in compression between the layers (A), and between pairs of layers (B). The distance between layers A and B reduces from 3.729 (2) to 3.539 (21) Å, and 5.597 (2) to 5.386 (22) Å from ambient to 3.7 GPa, respectively. The distance between the first and last adjacent pairs of layers for both the ambient and 3.7 GPa data sets is shown.

The final conglomerate of voids (4) occurs in the region formed by the NH₃—CH—CH₂—S moieties. The C—S—S—C torsion angle reduces from 75.27 (18)° at ambient pressure to 72.3 (4)° at 3.7 GPa. Since the S—S bond is approximately perpendicular to the c-axis the conformational change reduces the height of the molecule in the c-direction, and the distance between the planes containing the C_α atoms either side of the disulfide bridges (labelled B in Fig. 9) reduces from 5.597 (2) Å at ambient pressure to 5.386 (22) Å at 3.7 GPa. The twisting action of the L-cystine molecules could be described as a molecular `spring', whose action allows the compression of the structure along the c-axis direction.

4.2. Comparison of compression in α-glycine and L-cystine

We have previously determined the crystal structure of α-glycine at 2.0 GPa (Dawson et al., 2005), and in this paper we report the crystal structure of L-cystine at 2.3 GPa. Given the similarity in the layer structures of L-cystine and α-glycine it is interesting to compare the effect of pressures around 2 GPa on this part of the two crystal structures. Fig. 7 compares the voids in the two structures as a function of pressure.

α-Glycine is the more compressible of the two structures and at 2 GPa its unit-cell volume is 90% of its ambient pressure value; the corresponding figure of cystine is 93%. Comparison of Figs. 7(a) and 7(d) suggests that the voids in the centres of the (16) rings are slightly larger in α-glycine than in cystine, and this observation can be quantified by comparison of the respective areas of the ac and ab unit-cell faces in the two structures.² At ambient pressure these areas are 25.904 (2) Ų for glycine and 25.444 (3) Ų for cystine. At 2 GPa the areas are remarkably similar: 24.429 (8) Ų for glycine and 24.446 (1) Ų for cystine.

The overall similarity of the structures of the glycine-like layers is reflected in the similarity of the mean values of the principal components of the strain tensor which reside in the layers in both structures, viz. −0.03 for glycine and −0.02 for cystine. That the value for glycine is slightly higher than for cystine is presumably ascribable to the more efficient packing in cystine at ambient pressure.

Compression in the ab plane of cystine is constrained by the crystal symmetry to be isotropic, and this is reflected in the similar responses of the N1—H4⋯O1 and N1—H5⋯O2 hydrogen-bond distances to pressure (Table 2), which shorten by 0.10 and 0.08 Å, respectively. In α-glycine [which is monoclinic with the b-direction perpendicular to the view shown in Figs. 7(d) and 7(e)] the strain is markedly anisotropic. The largest strain component (−0.06) acts to close up the voids in the structure by forming a CH⋯O interaction as shown in Figs. 7(d) and 7(e). The N—H⋯O hydrogen bond which forms parallel to this C—H⋯O interaction shortens by 0.10 Å, whereas the other barely changes in length at all.

5. Conclusions

We have described the effect of pressure on the structure of hexagonal L-cystine at room temperature up to 3.7 GPa, with data on unit-cell dimensions reported up to 6.4 GPa. The structure can be described as being composed of pairs of glycine-like layers which lie in the ab plane. Each layer is composed of head-to-tail NH⋯OOC hydrogen bonds which combine to produce an (16) ring motif. Another NH⋯OOC hydrogen-bonding interaction then links the layers together in pairs via a R₂²(10) ring motif perpendicular to the c-axis. Paired layers then repeat down the c-axis separated by disulfide bridges about the 6₁ screw axes.

Under pressure, the structure undergoes anisotropic compression with the direction of greatest strain acting in the c-axis direction. This direction is normal to the glycine-like layers within the structure, which up to 3.7 GPa are pushed closer together along the c-axis direction. This compression occurs with the closing up of voids which are distributed along the c-direction (in regions 3 and 4 in Fig. 6). As this closure proceeds, two intramolecular torsion angles decrease, C1—S1—S1′—C1′ and N1—C2—C3—O2. Compression in the ab planes is constrained to be isotropic by symmetry and acts to decrease an already short S⋯S interaction and to close up voids within the (16) rings with a shortening of both N—H⋯O and C—H⋯O distances.

Supporting information

S1. Refinement

034_ALERT_1_C No Flack Parameter Given. Z. GT. Si, NonCentro. ?

The absolute structure was known prior to data collection.

127_ALERT_1_C Implicit Hall Symbol Inconsistent with Explicit P 61 2 (0 0 5)

This symbol was pasted from the PLATON Chk file.

199_ALERT_1_C Check the Reported _cell_measurement_temperature 293 K 200_ALERT_1_C Check the Reported _diffrn_ambient_temperature. 293 K

Confirm: This determination was performed at room temperature.

340_ALERT_3_C Low Bond Precision on C—C bonds (x 1000) Ang ··· 6 790_ALERT_4_C Centre of Gravity not Within Unit Cell: Resd. # 1 C6 H12 N2 O4 S2

No action taken.

910_ALERT_3_C # Missing FCF Reflections Below TH(Min) ······.. 1 911_ALERT_3_C # Missing FCF Refl. Between TH(Min) & STH/L=0.6 6 912_ALERT_3_C # Missing FCF Reflections Above STH/L=0.6 ······ 3

No action. Coverage stats are below:

================================================================================ Resolution & Completeness Statistics (Cumulative) ================================================================================ Theta sin(th)/Lambda Complete Expected Measured Missing ——————————————————————————– 20.82 0.500 0.982 397 390 7 23.01 0.550 0.986 500 493 7 25.24 0.600 0.989 650 643 7 ———————————————————— ACTA Min. Res. —- 26.36 0.625 0.986 717 707 10

Computing details

For all compounds, data collection: SMART (Siemens, 1993); cell refinement: SAINT (Siemens ,1995); data reduction: SAINT (Siemens ,1995). Program(s) used to solve structure: SIR92 (Altomare et al., 1994) for cystin00; SHELXS 86 (Sheldrick, 1986) for cystin04; USER DEFINED STRUCTURE SOLUTION for cystin14, cystin23, cystin37. For all compounds, program(s) used to refine structure: CRYSTALS (Betteridge et al., 2003); molecular graphics: CAMERON (Watkin et al., 1996); software used to prepare material for publication: CRYSTALS (Betteridge et al., 2003).

Figures

(cystin00) L-(-)-3,3'-Dithiobis(2-aminopropanoic acid)

Crystal data

C6H12N2O4S2	Dx = 1.681 Mg m^−3
Mr = 240.30	Mo Kα radiation, λ = 0.71073 Å
Hexagonal, P6122	Cell parameters from 1792 reflections
Hall symbol: P 61 2 (0 0 5)	θ = 9–46°
a = 5.4203 (5) Å	µ = 0.55 mm^−1
c = 55.980 (12) Å	T = 293 K
V = 1424.3 (4) Å^3	Hexagonal block, colourless
Z = 6	0.32 × 0.17 × 0.15 mm
F(000) = 756

Data collection

Bruker SMART APEX diffractometer	898 reflections with I > 2.00u(I)
Graphite monochromator	Rint = 0.085
ω scans	θmax = 26.4°, θmin = 4.3°
Absorption correction: multi-scan SADABS (Siemens, 1996)	h = −6→6
Tmin = 0.81, Tmax = 0.92	k = −6→6
16236 measured reflections	l = −64→68
965 independent reflections

Refinement

Refinement on F^2	Primary atom site location: structure-invariant direct methods
Least-squares matrix: full	Hydrogen site location: inferred from neighbouring sites
R[F^2 > 2σ(F^2)] = 0.048	H-atom parameters not refined
wR(F^2) = 0.105	w = 1/[σ^2(F^2) + 0.03 + 3.29P] where P = (max(Fo^2,0) + 2Fc^2)/3
S = 1.07	(Δ/σ)max = 0.001
956 reflections	Δρmax = 0.35 e Å^−3
64 parameters	Δρmin = −0.35 e Å^−3
0 restraints

Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Ų)

	x	y	z	Uiso*/Ueq
S1	0.2003 (2)	0.0314 (2)	−0.253692 (16)	0.0352
C1	0.0929 (8)	−0.2640 (8)	−0.27400 (6)	0.0317
C2	0.0528 (8)	−0.2125 (7)	−0.29990 (6)	0.0246
C3	−0.2003 (8)	−0.1667 (8)	−0.30422 (5)	0.0265
N1	0.3156 (6)	0.0347 (7)	−0.30997 (5)	0.0272
O1	−0.1505 (5)	0.0747 (6)	−0.31055 (4)	0.0299
O2	−0.4416 (5)	−0.3775 (6)	−0.30124 (5)	0.0395
H1	0.2418	−0.3214	−0.2734	0.0401*
H2	−0.0922	−0.4246	−0.2681	0.0401*
H3	0.0165	−0.3857	−0.3091	0.0308*
H6	0.2854	0.0601	−0.3254	0.0360*
H4	0.4613	0.0001	−0.3090	0.0360*
H5	0.3583	0.1928	−0.3016	0.0360*

Atomic displacement parameters (Ų)

	U^11	U^22	U^33	U^12	U^13	U^23
S1	0.0297 (5)	0.0469 (6)	0.0255 (5)	0.0164 (5)	0.0000 (4)	−0.0015 (4)
C1	0.035 (2)	0.036 (2)	0.029 (2)	0.0218 (18)	0.0003 (16)	0.0021 (15)
C2	0.0269 (19)	0.0251 (17)	0.0250 (17)	0.0155 (15)	0.0003 (15)	−0.0016 (14)
C3	0.0289 (18)	0.0319 (19)	0.0210 (17)	0.0169 (16)	0.0007 (14)	−0.0030 (14)
N1	0.0267 (15)	0.0384 (19)	0.0247 (16)	0.0225 (15)	0.0002 (12)	0.0018 (13)
O1	0.0286 (14)	0.0353 (14)	0.0344 (14)	0.0224 (12)	0.0020 (11)	0.0057 (11)
O2	0.0286 (15)	0.0329 (16)	0.0566 (18)	0.0151 (12)	0.0036 (13)	0.0006 (13)

Geometric parameters (Å, º)

S1—S1^i	2.043 (2)	C2—H3	1.000
S1—C1	1.806 (4)	C3—O1	1.248 (4)
C1—C2	1.513 (5)	C3—O2	1.245 (4)
C1—H1	1.000	N1—H6	0.900
C1—H2	1.000	N1—H4	0.900
C2—C3	1.530 (5)	N1—H5	0.900
C2—N1	1.495 (5)
S1^i—S1—C1	104.23 (13)	C3—C2—H3	107.4
S1—C1—C2	116.9 (3)	N1—C2—H3	107.4
S1—C1—H1	107.6	C2—C3—O1	118.3 (3)
C2—C1—H1	107.6	C2—C3—O2	116.5 (3)
S1—C1—H2	107.6	O1—C3—O2	125.3 (3)
C2—C1—H2	107.6	C2—N1—H6	109.5
H1—C1—H2	109.5	C2—N1—H4	109.5
C1—C2—C3	113.9 (3)	H6—N1—H4	109.5
C1—C2—N1	111.5 (3)	C2—N1—H5	109.5
C3—C2—N1	108.9 (3)	H6—N1—H5	109.5
C1—C2—H3	107.4	H4—N1—H5	109.5

Symmetry code: (i) −x+y, y, −z−1/2.

(cystin04) L-(-)-3,3'-Dithiobis(2-aminopropanoic acid)

Crystal data

C6H12N2O4S2	Dx = 1.691 Mg m^−3
Mr = 240.30	Synchrotron radiation, λ = 0.6774 Å
Hexagonal, P6122	Cell parameters from 3156 reflections
a = 5.4030 (1) Å	θ = 9–54°
c = 55.984 (4) Å	µ = 0.56 mm^−1
V = 1415.40 (10) Å^3	T = 293 K
Z = 6	Hexagonal block, colourless
F(000) = 756	0.20 × 0.20 × 0.10 mm

Data collection

Bruker APEX-II diffractometer	539 independent reflections
Radiation source: synchrotron	490 reflections with I > 2.00u(I)
Silicon (111) monochromator	Rint = 0.047
ω scans	θmax = 28.5°, θmin = 4.5°
Absorption correction: multi-scan SADABS (Siemens, 1996)	h = −7→7
Tmin = 0.71, Tmax = 0.95	k = −7→6
4719 measured reflections	l = −29→29

Refinement

Refinement on F^2	Primary atom site location: structure-invariant direct methods
Least-squares matrix: full	Hydrogen site location: inferred from neighbouring sites
R[F^2 > 2σ(F^2)] = 0.035	H-atom parameters not refined
wR(F^2) = 0.089	w = 1/[σ^2(F^2) + 0.05 + 0.8P] where P = (max(Fo^2,0) + 2Fc^2)/3
S = 1.06	(Δ/σ)max = 0.000153
513 reflections	Δρmax = 0.14 e Å^−3
64 parameters	Δρmin = −0.19 e Å^−3
67 restraints

Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Ų)

	x	y	z	Uiso*/Ueq
S1	0.03254 (16)	0.83088 (15)	0.07982 (3)	0.0300
O1	0.0774 (3)	1.2248 (3)	0.02287 (8)	0.0284
C3	−0.1663 (3)	1.0335 (3)	0.02879 (6)	0.0244
O2	−0.3798 (3)	1.0657 (3)	0.03203 (7)	0.0347
C2	−0.2104 (2)	0.7332 (2)	0.03313 (4)	0.0226
C1	−0.2643 (4)	0.6421 (4)	0.05923 (4)	0.0267
N1	0.0337 (4)	0.7156 (3)	0.02312 (6)	0.0231
H2	−0.4249	0.6684	0.0650	0.0303*
H1	−0.3233	0.4349	0.0599	0.0303*
H3	−0.3844	0.5958	0.0240	0.0270*
H6	0.0589	0.7709	0.0077	0.0282*
H5	0.1929	0.8309	0.0315	0.0282*
H4	−0.0022	0.5343	0.0241	0.0282*

Atomic displacement parameters (Ų)

	U^11	U^22	U^33	U^12	U^13	U^23
S1	0.0399 (4)	0.0373 (4)	0.021 (2)	0.0258 (3)	−0.0006 (6)	−0.0013 (5)
O1	0.0298 (10)	0.0175 (8)	0.037 (6)	0.0111 (7)	0.0058 (13)	0.0056 (11)
C3	0.0253 (11)	0.0182 (10)	0.032 (7)	0.0124 (9)	0.0004 (16)	0.0022 (17)
O2	0.0293 (9)	0.0297 (9)	0.051 (6)	0.0191 (8)	0.0011 (13)	−0.0012 (13)
C2	0.0236 (11)	0.0164 (9)	0.027 (4)	0.0097 (9)	0.0019 (16)	0.0011 (16)
C1	0.0307 (12)	0.0195 (10)	0.026 (3)	0.0094 (9)	0.0056 (14)	−0.0012 (16)
N1	0.0274 (11)	0.0148 (9)	0.028 (5)	0.0114 (8)	0.0040 (15)	0.0003 (13)

Geometric parameters (Å, º)

S1—S1^i	2.0417 (15)	C2—H3	1.000
S1—C1	1.8181 (7)	C1—H2	1.000
O1—C3	1.2456 (9)	C1—H1	1.000
C3—O2	1.2628 (9)	N1—H6	0.900
C3—C2	1.5366 (8)	N1—H5	0.900
C2—C1	1.523 (1)	N1—H4	0.900
C2—N1	1.4793 (9)
S1^i—S1—C1	103.88 (10)	C2—C1—H2	107.6
O1—C3—O2	125.68 (5)	S1—C1—H2	107.6
O1—C3—C2	117.38 (5)	C2—C1—H1	107.6
O2—C3—C2	116.93 (5)	S1—C1—H1	107.6
C3—C2—C1	113.37 (5)	H2—C1—H1	109.5
C3—C2—N1	109.61 (4)	C2—N1—H6	109.5
C1—C2—N1	112.04 (5)	C2—N1—H5	109.5
C3—C2—H3	107.2	H6—N1—H5	109.5
C1—C2—H3	107.2	C2—N1—H4	109.5
N1—C2—H3	107.2	H6—N1—H4	109.5
C2—C1—S1	117.03 (15)	H5—N1—H4	109.5

Symmetry code: (i) x, x−y+2, −z+1/6.

Hydrogen-bond geometry (Å, º)

D—H···A	D—H	H···A	D···A	D—H···A
N1—H6···O1^ii	0.900	2.050	2.829	144.2
N1—H5···O2^iii	0.900	2.003	2.806	147.9
N1—H4···O1^iv	0.900	1.925	2.777	157.4

Symmetry codes: (ii) x−y+1, −y+2, −z; (iii) x+1, y, z; (iv) x, y−1, z.

(cystin14) L-(-)-3,3'-Dithiobis(2-aminopropanoic acid)

Crystal data

C6H12N2O4S2	Dx = 1.752 Mg m^−3
Mr = 240.30	Synchrotron radiation, λ = 0.6774 Å
Hexagonal, P6122	Cell parameters from 1652 reflections
Hall symbol: P 61 2 (0 0 5)	θ = 9–53°
a = 5.3524 (2) Å	µ = 0.58 mm^−1
c = 55.087 (14) Å	T = 293 K
V = 1366.7 (4) Å^3	Hexagonal block, colourless
Z = 6	0.20 × 0.20 × 0.10 mm
F(000) = 756

Data collection

Bruker APEX-II diffractometer	492 independent reflections
Radiation source: synchrotron	464 reflections with I > 2.00u(I)
Silicon (111) monochromator	Rint = 0.047
ω scans	θmax = 28.1°, θmin = 4.5°
Absorption correction: multi-scan SADABS (Siemens, 1996)	h = −6→6
Tmin = 0.62, Tmax = 0.94	k = −7→7
4645 measured reflections	l = −28→28

Refinement

Refinement on F^2	Primary atom site location: structure-invariant direct methods
Least-squares matrix: full	Hydrogen site location: inferred from neighbouring sites
R[F^2 > 2σ(F^2)] = 0.036	H-atom parameters not refined
wR(F^2) = 0.089	w = 1/[σ^2(F^2) + 0.05 + 1.06P] where P = (max(Fo^2,0) + 2Fc^2)/3
S = 1.06	(Δ/σ)max = 0.000260
485 reflections	Δρmax = 0.22 e Å^−3
64 parameters	Δρmin = −0.17 e Å^−3
67 restraints

Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Ų)

	x	y	z	Uiso*/Ueq
S1	0.03587 (16)	0.83005 (16)	0.08014 (4)	0.0274
O1	0.0838 (3)	1.2266 (3)	0.02353 (8)	0.0229
C3	−0.1649 (3)	1.0326 (3)	0.02870 (7)	0.0216
O2	−0.3803 (3)	1.0660 (4)	0.03180 (7)	0.0292
C2	−0.2109 (3)	0.7291 (2)	0.03308 (4)	0.0204
C1	−0.2669 (4)	0.6364 (5)	0.05957 (4)	0.0251
N1	0.0359 (4)	0.7113 (3)	0.02303 (6)	0.0227
H2	−0.4287	0.6634	0.0654	0.0286*
H1	−0.3265	0.4273	0.0603	0.0286*
H3	−0.3860	0.5910	0.0237	0.0242*
H6	0.0624	0.7675	0.0074	0.0272*
H5	0.1961	0.8273	0.0316	0.0272*
H4	−0.0006	0.5282	0.0240	0.0272*

Atomic displacement parameters (Ų)

	U^11	U^22	U^33	U^12	U^13	U^23
S1	0.0348 (4)	0.0322 (4)	0.023 (2)	0.0224 (3)	−0.0015 (6)	−0.0017 (5)
O1	0.0259 (10)	0.0164 (8)	0.026 (7)	0.0102 (7)	0.0051 (14)	0.0042 (12)
C3	0.0215 (11)	0.0160 (11)	0.028 (8)	0.0103 (9)	−0.0008 (17)	−0.0012 (18)
O2	0.0255 (10)	0.0264 (10)	0.041 (7)	0.0171 (8)	−0.0008 (14)	−0.0019 (14)
C2	0.0188 (11)	0.0150 (10)	0.027 (4)	0.0079 (9)	0.0015 (17)	0.0016 (17)
C1	0.0291 (13)	0.0167 (11)	0.026 (3)	0.0085 (10)	0.0039 (15)	−0.0006 (16)
N1	0.0227 (11)	0.0119 (9)	0.033 (6)	0.0085 (9)	0.0045 (15)	0.0009 (14)

Geometric parameters (Å, º)

S1—S1^i	2.0418 (15)	C2—H3	1.000
S1—C1	1.8178 (7)	C1—H2	1.000
O1—C3	1.2447 (9)	C1—H1	1.000
C3—O2	1.2633 (9)	N1—H6	0.900
C3—C2	1.5353 (8)	N1—H5	0.900
C2—C1	1.5224 (10)	N1—H4	0.900
C2—N1	1.4784 (9)
S1^i—S1—C1	103.29 (11)	C2—C1—H2	107.8
O1—C3—O2	125.72 (5)	S1—C1—H2	107.8
O1—C3—C2	117.21 (5)	C2—C1—H1	107.8
O2—C3—C2	116.98 (5)	S1—C1—H1	107.8
C3—C2—C1	113.35 (5)	H2—C1—H1	109.5
C3—C2—N1	109.50 (4)	C2—N1—H6	109.5
C1—C2—N1	112.01 (5)	C2—N1—H5	109.5
C3—C2—H3	107.2	H6—N1—H5	109.5
C1—C2—H3	107.2	C2—N1—H4	109.5
N1—C2—H3	107.2	H6—N1—H4	109.5
C2—C1—S1	115.95 (16)	H5—N1—H4	109.5

Symmetry code: (i) x, x−y+2, −z+1/6.

Hydrogen-bond geometry (Å, º)

D—H···A	D—H	H···A	D···A	D—H···A
N1—H6···O1^ii	0.900	2.035	2.815	144.1
N1—H5···O2^iii	0.900	1.969	2.769	147.3
N1—H4···O1^iv	0.900	1.882	2.732	156.8

Symmetry codes: (ii) x−y+1, −y+2, −z; (iii) x+1, y, z; (iv) x, y−1, z.

(cystin23) L-(-)-3,3'-Dithiobis(2-aminopropanoic acid)

Crystal data

C6H12N2O4S2	Dx = 1.803 Mg m^−3
Mr = 240.30	Synchrotron radiation, λ = 0.6774 Å
Hexagonal, P6122	Cell parameters from 3365 reflections
a = 5.3141 (1) Å	θ = 9–53°
c = 54.291 (3) Å	µ = 0.59 mm^−1
V = 1327.75 (8) Å^3	T = 293 K
Z = 6	Hexagonal block, colourless
F(000) = 756	0.20 × 0.20 × 0.10 mm

Data collection

Bruker APEX-II diffractometer	492 independent reflections
Radiation source: synchrotron	438 reflections with I > 2.00u(I)
Silicon (111) monochromator	Rint = 0.053
ω scans	θmax = 28.1°, θmin = 4.5°
Absorption correction: multi-scan SADABS (Siemens, 1996)	h = −7→6
Tmin = 0.64, Tmax = 0.94	k = −7→6
4229 measured reflections	l = −28→28

Refinement

Refinement on F^2	Primary atom site location: structure-invariant direct methods
Least-squares matrix: full	Hydrogen site location: inferred from neighbouring sites
R[F^2 > 2σ(F^2)] = 0.048	H-atom parameters not refined
wR(F^2) = 0.112	w = 1/[σ^2(F^2) + 0.03 + 3.33P] where P = (max(Fo^2,0) + 2Fc^2)/3
S = 1.05	(Δ/σ)max = 0.000120
458 reflections	Δρmax = 0.27 e Å^−3
64 parameters	Δρmin = −0.21 e Å^−3
67 restraints

Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Ų)

	x	y	z	Uiso*/Ueq
S1	0.0379 (2)	0.8292 (2)	0.08036 (6)	0.0316
O1	0.0899 (5)	1.2288 (4)	0.02371 (12)	0.0218
C3	−0.1619 (4)	1.0333 (4)	0.02859 (10)	0.0256
O2	−0.3791 (5)	1.0675 (5)	0.03148 (11)	0.0331
C2	−0.2098 (4)	0.7269 (3)	0.03295 (7)	0.0247
C1	−0.2695 (6)	0.6317 (6)	0.05979 (6)	0.0291
N1	0.0390 (6)	0.7081 (5)	0.02292 (10)	0.0255
H2	−0.4323	0.6592	0.0657	0.0330*
H1	−0.3298	0.4210	0.0606	0.0330*
H3	−0.3855	0.5889	0.0233	0.0285*
H6	0.0676	0.7658	0.0071	0.0295*
H5	0.1997	0.8237	0.0317	0.0295*
H4	0.0011	0.5232	0.0238	0.0295*

Atomic displacement parameters (Ų)

	U^11	U^22	U^33	U^12	U^13	U^23
S1	0.0322 (6)	0.0292 (6)	0.041 (3)	0.0209 (5)	−0.0013 (9)	−0.0005 (8)
O1	0.0224 (14)	0.0147 (12)	0.028 (9)	0.0086 (11)	0.007 (2)	0.0037 (17)
C3	0.0187 (16)	0.0163 (16)	0.043 (10)	0.0099 (13)	0.001 (3)	−0.001 (3)
O2	0.0225 (14)	0.0244 (15)	0.058 (10)	0.0157 (12)	0.000 (2)	0.001 (2)
C2	0.0157 (16)	0.0127 (14)	0.043 (4)	0.0049 (13)	0.002 (2)	0.001 (2)
C1	0.0240 (19)	0.0153 (17)	0.043 (4)	0.0063 (15)	0.0057 (19)	0.002 (2)
N1	0.0183 (16)	0.0112 (14)	0.044 (7)	0.0054 (13)	0.003 (2)	−0.001 (2)

Geometric parameters (Å, º)

S1—S1^i	2.043 (2)	C2—H3	1.000
S1—C1	1.8177 (7)	C1—H2	1.000
O1—C3	1.2446 (9)	C1—H1	1.000
C3—O2	1.2651 (9)	N1—H6	0.900
C3—C2	1.5359 (8)	N1—H5	0.900
C2—C1	1.5226 (9)	N1—H4	0.900
C2—N1	1.4791 (9)
S1^i—S1—C1	102.90 (16)	C2—C1—H2	108.1
O1—C3—O2	125.73 (5)	S1—C1—H2	108.1
O1—C3—C2	117.19 (5)	C2—C1—H1	108.1
O2—C3—C2	116.98 (5)	S1—C1—H1	108.1
C3—C2—C1	113.36 (5)	H2—C1—H1	109.5
C3—C2—N1	109.56 (5)	C2—N1—H6	109.5
C1—C2—N1	112.03 (5)	C2—N1—H5	109.5
C3—C2—H3	107.2	H6—N1—H5	109.5
C1—C2—H3	107.2	C2—N1—H4	109.5
N1—C2—H3	107.2	H6—N1—H4	109.5
C2—C1—S1	114.7 (2)	H5—N1—H4	109.5

Symmetry code: (i) x, x−y+2, −z+1/6.

Hydrogen-bond geometry (Å, º)

D—H···A	D—H	H···A	D···A	D—H···A
N1—H6···O1^ii	0.900	2.007	2.780	143.1
N1—H5···O2^iii	0.900	1.946	2.742	146.5
N1—H4···O1^iv	0.900	1.846	2.693	155.8

Symmetry codes: (ii) x−y+1, −y+2, −z; (iii) x+1, y, z; (iv) x, y−1, z.

(cystin37) L-(-)-3,3'-Dithiobis(2-aminopropanoic acid)

Crystal data

C6H12N2O4S2	Dx = 1.861 Mg m^−3
Mr = 240.30	Synchrotron radiation, λ = 0.6774 Å
Hexagonal, P6122	Cell parameters from 3341 reflections
a = 5.2663 (1) Å	θ = 9–51°
c = 53.549 (5) Å	µ = 0.61 mm^−1
V = 1286.17 (13) Å^3	T = 293 K
Z = 6	Hexagonal block, colourless
F(000) = 756	0.20 × 0.20 × 0.10 mm

Data collection

Bruker APEX-II diffractometer	467 independent reflections
Radiation source: synchrotron	439 reflections with I > 2.00u(I)
Silicon (111) monochromator	Rint = 0.063
ω scans	θmax = 27.9°, θmin = 4.5°
Absorption correction: multi-scan SADABS (Siemens, 1996)	h = −6→6
Tmin = 0.56, Tmax = 0.94	k = −6→6
4229 measured reflections	l = −28→28

Refinement

Refinement on F^2	Primary atom site location: structure-invariant direct methods
Least-squares matrix: full	Hydrogen site location: inferred from neighbouring sites
R[F^2 > 2σ(F^2)] = 0.093	H-atom parameters not refined
wR(F^2) = 0.214	w = 1/[σ^2(F^2) + 0.00 + 18.85P] where P = (max(Fo^2,0) + 2Fc^2)/3
S = 1.09	(Δ/σ)max = 0.000065
451 reflections	Δρmax = 0.54 e Å^−3
34 parameters	Δρmin = −0.60 e Å^−3
31 restraints

Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Ų)

	x	y	z	Uiso*/Ueq
S1	0.0419 (5)	0.8295 (5)	0.08057 (11)	0.0255
C1	−0.2720 (12)	0.6261 (13)	0.06016 (12)	0.0209 (14)*
C2	−0.2060 (8)	0.7249 (6)	0.03304 (13)	0.0180 (12)*
C3	−0.1546 (8)	1.0357 (8)	0.0287 (2)	0.0181 (13)*
N1	0.0457 (13)	0.7045 (10)	0.0232 (2)	0.0179 (14)*
O1	0.1011 (9)	1.2328 (8)	0.0243 (2)	0.0162 (13)*
O2	−0.3749 (9)	1.0704 (10)	0.0307 (2)	0.0201 (14)*
H2	−0.4362	0.6544	0.0661	0.0251*
H1	−0.3332	0.4134	0.0611	0.0251*
H3	−0.3819	0.5876	0.0230	0.0216*
H6	0.0780	0.7645	0.0072	0.0215*
H5	0.2065	0.8191	0.0323	0.0215*
H4	0.0056	0.5173	0.0240	0.0215*

Atomic displacement parameters (Ų)

	U^11	U^22	U^33	U^12	U^13	U^23
S1	0.0300 (12)	0.0262 (13)	0.027 (3)	0.0190 (10)	−0.0057 (15)	−0.0026 (15)

Geometric parameters (Å, º)

S1—S1^i	2.038 (5)	C2—H3	1.000
S1—C1	1.8178 (7)	C3—O1	1.2446 (10)
C1—C2	1.5229 (10)	C3—O2	1.2662 (10)
C1—H2	1.000	N1—H6	0.900
C1—H1	1.000	N1—H5	0.900
C2—C3	1.5373 (9)	N1—H4	0.900
C2—N1	1.4798 (10)
S1^i—S1—C1	102.7 (3)	C3—C2—H3	107.1
S1—C1—C2	112.8 (4)	N1—C2—H3	107.1
S1—C1—H2	108.6	C2—C3—O1	117.23 (6)
C2—C1—H2	108.6	C2—C3—O2	117.00 (5)
S1—C1—H1	108.6	O1—C3—O2	125.72 (6)
C2—C1—H1	108.6	C2—N1—H6	109.5
H2—C1—H1	109.5	C2—N1—H5	109.5
C1—C2—C3	113.38 (5)	H6—N1—H5	109.5
C1—C2—N1	112.04 (5)	C2—N1—H4	109.5
C3—C2—N1	109.64 (5)	H6—N1—H4	109.5
C1—C2—H3	107.1	H5—N1—H4	109.5

Symmetry code: (i) x, x−y+2, −z+1/6.

Hydrogen-bond geometry (Å, º)

D—H···A	D—H	H···A	D···A	D—H···A
N1—H6···O1^ii	0.900	2.017	2.783	142.2
N1—H5···O2^iii	0.900	1.924	2.703	143.9
N1—H4···O1^iv	0.900	1.804	2.643	154.2

Symmetry codes: (ii) x−y+1, −y+2, −z; (iii) x+1, y, z; (iv) x, y−1, z.

 

Footnotes

¹Supplementary data for this paper are available from the IUCr electronic archives (Reference: XD5007 ). Services for accessing these data are described at the back of the journal.

²The unit-cell data for α-glycine at ambient pressure [2.0 GPa] are: monoclinic, a = 5.1047 (3) [4.9669 (9)], b = 11.9720 (14) [11.459 (4)], c = 5.4631 (3) [5.4231 (12)] Å, β = 111.740 (5)° [114.916 (15)°]. All data were measured at room temperature. References: Boldyreva et al. (2003) and Dawson et al. (2005).