The critical role of QM/MM X-ray refinement and accurate tautomer/protomer determination in structure-based drug design

Conventional protein:ligand crystallographic refinement uses stereochemistry restraints coupled with a rudimentary energy functional to ensure the correct geometry of the model of the macromolecule—along with any bound ligand(s)—within the context of the experimental, X-ray density. These methods generally lack explicit terms for electrostatics, polarization, dispersion, hydrogen bonds, and other key interactions, and instead they use pre-determined parameters (e.g. bond lengths, angles, and torsions) to drive structural refinement. In order to address this deficiency and obtain a more complete and ultimately more accurate structure, we have developed an automated approach for macromolecular refinement based on a two layer, QM/MM (ONIOM) scheme as implemented within our DivCon Discovery Suite and "plugged in" to two mainstream crystallographic packages: PHENIX and BUSTER. This implementation is able to use one or more region layer(s), which is(are) characterized using linear-scaling, semi-empirical quantum mechanics, followed by a system layer which includes the balance of the model and which is described using a molecular mechanics functional. In this work, we applied our Phenix/DivCon refinement method—coupled with our XModeScore method for experimental tautomer/protomer state determination—to the characterization of structure sets relevant to structure-based drug design (SBDD). We then use these newly refined structures to show the impact of QM/MM X-ray refined structure on our understanding of function by exploring the influence of these improved structures on protein:ligand binding affinity prediction (and we likewise show how we use post-refinement scoring outliers to inform subsequent X-ray crystallographic efforts). Through this endeavor, we demonstrate a computational chemistry ↔ structural biology (X-ray crystallography) "feedback loop" which has utility in industrial and academic pharmaceutical research as well as other allied fields.

MovableType Software for Fast Free Energy-Based Virtual Screening: Protocol Development, Deployment, Validation, and Assessment

For decades, the complicated energy surfaces found in macromolecular protein:ligand structures, which require large amounts of computational time and resources for energy state sampling, have been an inherent obstacle to fast, routine free energy estimation in industrial drug discovery efforts. Beginning in 2013, the Merz research group addressed this cost with the introduction of a novel sampling methodology termed “Movable Type” (MT). Using numerical integration methods, the MT method reduces the computational expense for energy state sampling by independently calculating each atomic partition function from an initial molecular conformation in order to estimate the molecular free energy using ensembles of the atomic partition functions. In this work, we report a software package, the DivCon Discovery Suite with the MovableType module from QuantumBio Inc., that performs this MT free energy estimation protocol in a fast, fully encapsulated manner. We discuss the computational procedures and improvements to the original work, and we detail the corresponding settings for this software package. Finally, we introduce two validation benchmarks to evaluate the overall robustness of the method against a broad range of protein:ligand structural cases. With these publicly available benchmarks, we show that the method can use a variety of input types and parameters and exhibits comparable predictability whether the method is presented with “expensive” X-ray structures or “inexpensively docked” theoretical models. We also explore some next steps for the method. The MovableType software is available at http://www.quantumbioinc.com/

Accurate macromolecular crystallographic refinement: incorporation of the linear scaling, semiempirical quantum-mechanics program DivCon into the PHENIX refinement package

Macromolecular crystallographic refinement relies on sometimes dubious stereochemical restraints and rudimentary energy functionals to ensure the correct geometry of the model of the macromolecule and any covalently bound ligand(s). The ligand stereochemical restraint file (CIF) requires a priori understanding of the ligand geometry within the active site, and creation of the CIF is often an error-prone process owing to the great variety of potential ligand chemistry and structure. Stereochemical restraints have been replaced with more robust functionals through the integration of the linear-scaling, semiempirical quantum-mechanics (SE-QM) program DivCon with the PHENIX X-ray refinement engine. The PHENIX/DivCon package has been thoroughly validated on a population of 50 protein-ligand Protein Data Bank (PDB) structures with a range of resolutions and chemistry. The PDB structures used for the validation were originally refined utilizing various refinement packages and were published within the past five years. PHENIX/DivCon does not utilize CIF(s), link restraints and other parameters for refinement and hence it does not make as many a priori assumptions about the model. Across the entire population, the method results in reasonable ligand geometries and low ligand strains, even when the original refinement exhibited difficulties, indicating that PHENIX/DivCon is applicable to both single-structure and high-throughput crystallography.

Increased Peptide Contacts Govern High Affinity Binding of a Modified TCR Whilst Maintaining a Native pMHC Docking Mode

Natural T cell receptors (TCRs) generally bind to their cognate pMHC molecules with weak affinity and fast kinetics, limiting their use as therapeutic agents. Using phage display, we have engineered a high affinity version of the A6 wild-type TCR (A6wt), specific for the human leukocyte antigen (HLA-A^∗0201) complexed with human T cell lymphotropic virus type 1_11–19 peptide (A2-Tax). Mutations in just 4 residues in the CDR3β loop region of the A6wt TCR were selected that improved binding to A2-Tax by nearly 1000-fold. Biophysical measurements of this mutant TCR (A6c134) demonstrated that the enhanced binding was derived through favorable enthalpy and a slower off-rate. The structure of the free A6c134 TCR and the A6c134/A2-Tax complex revealed a native binding mode, similar to the A6wt/A2-Tax complex. However, concordant with the more favorable binding enthalpy, the A6c134 TCR made increased contacts with the Tax peptide compared with the A6wt/A2-Tax complex, demonstrating a peptide-focused mechanism for the enhanced affinity that directly involved the mutated residues in the A6c134 TCR CDR3β loop. This peptide-focused enhanced TCR binding may represent an important approach for developing antigen specific high affinity TCR reagents for use in T cell based therapies.

Loss of T cell antigen recognition arising from changes in peptide and major histocompatibility complex protein flexibility: implications for vaccine design

Modification of the primary anchor positions of antigenic peptides to improve binding to major histocompatibility complex (MHC) proteins is a commonly used strategy for engineering peptide-based vaccine candidates. However, such peptide modifications do not always improve antigenicity, complicating efforts to design effective vaccines for cancer and infectious disease. Here we investigated the MART-1(27-35) tumor antigen, for which anchor modification (replacement of the position two alanine with leucine) dramatically reduces or ablates antigenicity with a wide range of T cell clones despite significantly improving peptide binding to MHC. We found that anchor modification in the MART-1(27-35) antigen enhances the flexibility of both the peptide and the HLA-A*0201 molecule. Although the resulting entropic effects contribute to the improved binding of the peptide to MHC, they also negatively impact T cell receptor binding to the peptide·MHC complex. These results help explain how the "anchor-fixing" strategy fails to improve antigenicity in this case, and more generally, may be relevant for understanding the high specificity characteristic of the T cell repertoire. In addition to impacting vaccine design, modulation of peptide and MHC flexibility through changes to antigenic peptides may present an evolutionary strategy for the escape of pathogens from immune destruction.

TCRs used in cancer gene therapy cross-react with MART-1/Melan-A tumor antigens via distinct mechanisms

T cells engineered to express TCRs specific for tumor Ags can drive cancer regression. The first TCRs used in cancer gene therapy, DMF4 and DMF5, recognize two structurally distinct peptide epitopes of the melanoma-associated MART-1/Melan-A protein, both presented by the class I MHC protein HLA-A*0201. To help understand the mechanisms of TCR cross-reactivity and provide a foundation for the further development of immunotherapy, we determined the crystallographic structures of DMF4 and DMF5 in complex with both of the MART-1/Melan-A epitopes. The two TCRs use different mechanisms to accommodate the two ligands. Although DMF4 binds the two with a different orientation, altering its position over the peptide/MHC, DMF5 binds them both identically. The simpler mode of cross-reactivity by DMF5 is associated with higher affinity toward both ligands, consistent with the superior functional avidity of DMF5. More generally, the observation of two diverging mechanisms of cross-reactivity with the same Ags and the finding that TCR-binding orientation can be determined by peptide alone extend our understanding of the mechanisms underlying TCR cross-reactivity.

Probing dimerization and structural flexibility of mammalian lipoxygenases by small-angle X-ray scattering

Human lipoxygenases (LOXs) and their metabolites have a great impact on human homeostasis and are of interest for targeted drug design. This goal requires detailed knowledge of their structures and an understanding of structure-function relationship. At the moment, there are two complete crystal structures for mammalian LOX [rabbit 12/15LOX (r-12/15LOX) and human 5LOX (h-5LOX)] and a fragment of human 12LOX. The low-resolution structures in solution for various LOX isoforms have brought about controversial results. Here we explored the behavior of r-12/15LOX in aqueous solution under different conditions (salt and pH) by small-angle X-ray scattering (SAXS) and compared it with human platelet-type 12S-LOX (hp-12LOX) and h-5LOX. Thermodynamic calculations concerning the stability of molecular assemblies, thermal motion analysis [TLSMD (translation, libration, and screw rotation motion detection based on crystallographic temperature factor B(j))], and results of SAXS analyses brought about the following conclusions: (i) in contrast to its crystal structure, r-12/15LOX functions as a monomer that dominates in solution; (ii) it dimerizes at higher protein concentrations in the presence of salt and with increasing degree of motional freedom of the N-terminal PLAT domain, as suggested by the Y98,614→R double mutant; (iii) in aqueous solutions, hp-12LOX is stable as a dimer, in contrast to h-5LOX and r-12/15LOX, which are monomeric; and (iv) all three mammalian isozymes show a high level of flexibility not only for the PLAT domain but also for other subdomains of the catalytic part in TLS (translation, libration, and screw rotation) analysis and hp-12LOX in SAXS.

Conformational melding permits a conserved binding geometry in TCR recognition of foreign and self molecular mimics

Molecular mimicry between foreign and self Ags is a mechanism of TCR cross-reactivity and is thought to contribute to the development of autoimmunity. The αβ TCR A6 recognizes the foreign Ag Tax from the human T cell leukemia virus-1 when presented by the class I MHC HLA-A2. In a possible link with the autoimmune disease human T cell leukemia virus-1-associated myelopathy/tropical spastic paraparesis, A6 also recognizes a self peptide from the neuronal protein HuD in the context of HLA-A2. We found in our study that the complexes of the HuD and Tax epitopes with HLA-A2 are close but imperfect structural mimics and that in contrast with other recent structures of TCRs with self Ags, A6 engages the HuD Ag with the same traditional binding mode used to engage Tax. Although peptide and MHC conformational changes are needed for recognition of HuD but not Tax and the difference of a single hydroxyl triggers an altered TCR loop conformation, TCR affinity toward HuD is still within the range believed to result in negative selection. Probing further, we found that the HuD-HLA-A2 complex is only weakly stable. Overall, these findings help clarify how molecular mimicry can drive self/nonself cross-reactivity and illustrate how low peptide-MHC stability can permit the survival of T cells expressing self-reactive TCRs that nonetheless bind with a traditional binding mode.

Structures of native and affinity-enhanced WT1 epitopes bound to HLA-A*0201: implications for WT1-based cancer therapeutics

Presentation of peptides by class I or class II major histocompatibility complex (MHC) molecules is required for the initiation and propagation of a T cell-mediated immune response. Peptides from the Wilms Tumor 1 transcription factor (WT1), upregulated in many hematopoetic and solid tumors, can be recognized by T cells and numerous efforts are underway to engineer WT1-based cancer vaccines. Here we determined the structures of the class I MHC molecule HLA-A*0201 bound to the native 126-134 epitope of the WT1 peptide and a recently described variant (R1Y) with improved MHC binding. The R1Y variant, a potential vaccine candidate, alters the positions of MHC charged side chains near the peptide N-terminus and significantly reduces the peptide/MHC electrostatic surface potential. These alterations indicate that the R1Y variant is an imperfect mimic of the native WT1 peptide, and suggest caution in its use as a therapeutic vaccine. Stability measurements revealed how the R1Y substitution enhances MHC binding affinity, and together with the structures suggest a strategy for engineering WT1 variants with improved MHC binding that retain the structural features of the native peptide/MHC complex.

Structures of MART-126/27-35 Peptide/HLA-A2 complexes reveal a remarkable disconnect between antigen structural homology and T cell recognition

Small structural changes in peptides presented by major histocompatibility complex (MHC) molecules often result in large changes in immunogenicity, supporting the notion that T cell receptors are exquisitely sensitive to antigen structure. Yet there are striking examples of TCR recognition of structurally dissimilar ligands. The resulting unpredictability of how T cells will respond to different or modified antigens impacts both our understanding of the physical bases for TCR specificity as well as efforts to engineer peptides for immunomodulation. In cancer immunotherapy, epitopes and variants derived from the MART-1/Melan-A protein are widely used as clinical vaccines. Two overlapping epitopes spanning amino acid residues 26 through 35 are of particular interest: numerous clinical studies have been performed using variants of the MART-1 26-35 decamer, although only the 27-35 nonamer has been found on the surface of targeted melanoma cells. Here, we show that the 26-35 and 27-35 peptides adopt strikingly different conformations when bound to HLA-A2. Nevertheless, clonally distinct MART-1(26/27-35)-reactive T cells show broad cross-reactivity towards these ligands. Simultaneously, however, many of the cross-reactive T cells remain unable to recognize anchor-modified variants with very subtle structural differences. These dichotomous observations challenge our thinking about how structural information on unligated peptide/MHC complexes should be best used when addressing questions of TCR specificity. Our findings also indicate that caution is warranted in the design of immunotherapeutics based on the MART-1 26/27-35 epitopes, as neither cross-reactivity nor selectivity is predictable based on the analysis of the structures alone.

T cell receptor recognition via cooperative conformational plasticity

Although T cell receptor cross-reactivity is a fundamental property of the immune system and is implicated in numerous autoimmune pathologies, the molecular mechanisms by which T cell receptors can recognize and respond to diverse ligands are incompletely understood. In the current study we examined the response of the human T cell lymphotropic virus-1 (HTLV-1) Tax-specific T cell receptor (TCR) A6 to a panel of structurally distinct haptens coupled to the Tax 11-19 peptide with a lysine substitution at position 5 (Tax5K, LLFG[K-hapten]PVYV). The A6 TCR could cross-reactively recognize one of these haptenated peptides, Tax-5K-4-(3-Indolyl)-butyric acid (IBA), presented by HLA-A*0201. The crystal structures of Tax5K-IBA/HLA-A2 free and in complex with A6 reveal that binding is mediated by a mechanism of cooperative conformational plasticity involving conformational changes on both sides of the protein-protein interface, including the TCR complementarity determining region (CDR) loops, Valpha/Vbeta domain orientation, and the hapten-modified peptide. Our findings illustrate the complex role that protein dynamics can play in TCR cross-reactivity and highlight that T cell receptor recognition of ligand can be achieved through diverse and complex molecular mechanisms that can occur simultaneously in the interface, not limited to molecular mimicry and CDR loop shifts.

Effect of crystal freezing and small-molecule binding on internal cavity size in a large protein: X-ray and docking studies of lipoxygenase at ambient and low temperature at 2.0 A resolution

Flash-freezing is a technique that is commonly used nowadays to collect diffraction data for X-ray structural analysis. It can affect both the crystal and molecular structure and the molecule's surface, as well as the internal cavities. X-ray structural data often serve as a template for the protein receptor in docking calculations. Thus, the size and shape of the binding site determines which small molecules could be found as potential ligands in silico, especially during high-throughput rigid docking. Data were analyzed for wild soybean lipoxygenase-3 (MW 97 kDa) at 293 and 93 K and compared with the results from studies of its molecular complexes with known inhibitors, structures published by others for a derivative of the same enzyme (98 K) or a topologically close isozyme lipoxygenase-1 (at ambient temperature and 100 K). Analysis of these data allows the following conclusions. (i) Very small changes in the relative orientation of the molecules in the crystal can cause major changes in the crystal reciprocal lattice. (ii) The volume of the internal cavities can ;shrink' by several percent upon freezing even when the unit-cell and the protein molecular volume show changes of only 1-2%. (iii) Using a receptor structure determined based on cryogenic data as a target for computational screening requires flexible docking to enable the expansion of the binding-site cavity and sampling of the alternative conformations of the crucial residues.

Unraveling a hotspot for TCR recognition on HLA-A2: evidence against the existence of peptide-independent TCR binding determinants

T cell receptor (TCR) recognition of peptide takes place in the context of the major histocompatibility complex (MHC) molecule, which accounts for approximately two-thirds of the peptide/MHC buried surface. Using the class I MHC HLA-A2 and a large panel of mutants, we have previously shown that surface mutations that disrupt TCR recognition vary with the identity of the peptide. The single exception is Lys66 on the HLA-A2 alpha1 helix, which when mutated to alanine disrupts recognition for 93% of over 250 different T cell clones or lines, independent of which peptide is bound. Thus, Lys66 could serve as a peptide-independent TCR binding determinant. Here, we have examined the role of Lys66 in TCR recognition of HLA-A2 in detail. The structure of a peptide/HLA-A2 molecule with the K66A mutation indicates that although the mutation induces no major structural changes, it results in the exposure of a negatively charged glutamate (Glu63) underneath Lys66. Concurrent replacement of Glu63 with glutamine restores TCR binding and function for T cells specific for five different peptides presented by HLA-A2. Thus, the positive charge on Lys66 does not serve to guide all TCRs onto the HLA-A2 molecule in a manner required for productive signaling. Furthermore, electrostatic calculations indicate that Lys66 does not contribute to the stability of two TCR-peptide/HLA-A2 complexes. Our findings are consistent with the notion that each TCR arrives at a unique solution of how to bind a peptide/MHC, most strongly influenced by the chemical and structural features of the bound peptide. This would not rule out an intrinsic affinity of TCRs for MHC molecules achieved through multiple weak interactions, but for HLA-A2 the collective mutational data place limits on the role of any single MHC amino acid side-chain in driving TCR binding in a peptide-independent fashion.

Aspirin blocks binding of photosensitizer SnET2 into human serum albumin: implications for photodynamic therapy

Tin etiopurpurin dichloride (SnET2) is one of the photosensitizers under investigation to be used in photodynamic therapy of prostate cancer. The drug is delivered intravenously, transported in vivo by liposomes and plasma proteins and localized within the prostate. SnET2 exists in two tautomeric forms (I - closed ring, II - open ring) with I converting spontaneously into the more energetically stable form II at physiological pH. Up to approximately 50% of the drug can be carried by serum albumin, although this association can increase photo-bleaching and diminish the drug efficiency. Molecular modeling and force field calculations indicate that Sudlow Site I in human serum albumin (HSA) is the most probable binding site for both forms of SnET2, with the porphyrin moiety nestling between domains IIA and IB, and the esterolytic side group oriented toward domain IIIA of HSA. Other drugs, including aspirin, bind to the same part of HSA. SnET2 does not bind to HSA when pre-incubated with aspirin, which confirms that its place of binding to this protein must be located near Lys199. This observation could be exploited to improve photo-efficiency of SnET2 by finding drugs that could compete with the photosensitizer for binding into Sudlow Site I of HSA.

Increased immunogenicity of an anchor-modified tumor-associated antigen is due to the enhanced stability of the peptide/MHC complex: implications for vaccine design

The use of "anchor-fixed" altered peptide ligands is of considerable interest in the development of therapeutic vaccines for cancer and infectious diseases, but the mechanism by which successful altered peptide ligands elicit enhanced immunity is unclear. In this study, we have determined the crystallographic structure of a major tumor rejection Ag, gp100(209-217), in complex with the HLA-A*0201 (HLA-A2) molecule, as well as the structure of a modified version of the peptide which substitutes methionine for threonine at position 2 (T2M; gp100(209-2M)). The T2M-modified peptide, which is more immunogenic in vitro and in vivo, binds HLA-A2 with a approximately 9-fold greater affinity and has a approximately 7-fold slower dissociation rate at physiological temperature. Within the limit of the crystallographic data, the T2M substitution does not alter the structure of the peptide/HLA-A2 complex. Consistent with this finding, in peripheral blood from 95 human subjects, we were unable to identify higher frequencies of T cells specific for either the native or modified peptide. These data strongly support the conclusion that the greater immunogenicity of the gp100(209-2M) peptide is due to the enhanced stability of the peptide/MHC complex, validating the anchor-fixing approach for generating therapeutic vaccine candidates. Thermodynamic data suggest that the enhanced stability of the T2M-modified peptide/HLA-A2 complex is attributable to the increased hydrophobicity of the modified peptide, but the gain due to hydrophobicity is offset considerably by the loss of a hydrogen bond made by the native peptide to the HLA-A2 molecule. Our findings have broad implications for the optimization of current vaccine-design strategies.

Lipoxygenase interactions with natural flavonoid, quercetin, reveal a complex with protocatechuic acid in its X-ray structure at 2.1 A resolution

PUFA metabolites have a profound effect on inflammatory diseases and cancer progression. Blocking their production by inhibiting PUFA metabolizing enzymes (dioxygenases: cyclooxygenases and LOXs) might be a successful way to control and relieve such problems, if we learn to better understand their actions at a molecular level. Compounds with strong antioxidative and free radical scavenging properties, such as polyphenols, could be effective in blocking PUFA activities, and natural flavonoids possess such qualities. Quercetin belongs to the group of natural catecholic compounds and is known as a potent, competitive inhibitor of LOX. Structural analysis reveals that quercetin entrapped within LOX undergoes degradation, and the resulting compound has been identified by X-ray analysis as protocatechuic acid (3,4-dihydroxybenzoic acid) positioned near the iron site. Its C3-OH group points toward His523, C4-OH forms a hydrogen bond with O=C from the enzyme's C-terminus, and the carboxylic group is incorporated into the hydrogen bonding network of the active-site neighborhood via Gln514. This unexpected result, together with our previous observations concerning other polyphenols, yields new evidence about the metabolism of natural flavonoids. These compounds might be vulnerable to the co-oxidase activity of LOX, leading to enzyme-stimulated oxidative degradation, which results in an inhibitor of a lower molecular weight.

Soybean lipoxygenase-3 in complex with 4-nitrocatechol

4-Nitrocatechol (4NC) is a known inhibitor of lipoxygenase. This work presents the X-ray structure of soybean lipoxygenase-3 in complex with 4NC refined at 2.15 A resolution. The X-ray analysis shows 4NC near iron with partial occupancy, blocking access to Fe but not covalently bound to it. The two hydroxyl groups interact with the C-terminus (4-OH) and His523 ND1 (3-OH). The residues surrounding the iron cofactor, His518*, His523, His709, Ile857* COO(-) and water, form a trigonal bipyramid with the residues marked with asterisks in the axial positions. The water bound to iron and the presence of the inhibitor seem to be responsible for the rearrangements and changes in the geometry of the ligand distribution and confirm the displacement of His518 from iron coordination. A description of the catechol binding contributes to the understanding of lipoxygenase inhibition and the participation of its co-oxidative activity in the utilization of natural flavonoids.

Indole- and carbazole-substituted pyridinium iodide salts: a rare case of conformational isomerism in crystals

A series of indole- and carbazole-substituted pyridinium iodide salts has been synthesized and characterized. X-ray analysis revealed that the iodide salt of the indole-substituted cation (E)-4-(1H-indol-3-ylvinyl)-N-methylpyridinium (IMPE(+)), C(16)H(15)N(2)(+) x I(-), (I), has two polymorphic modifications, (Ia) and (Ib), and a hemihydrate structure, C(16)H(15)N(2)(+) x I(-) x 0.5H(2)O, (II). Until now, only one crystal modification has been identified for the (E)-4-(9-ethyl-9H-carbazol-3-ylvinyl)-N-methylpyridinium (ECMPE(+)) iodide salt, C(22)H(21)N(2)(+) x I(-), (III). Crystals of (Ia) and (Ib) comprise stacks of antiparallel cations with iodide anions located in the channels between the stacks. Due to the presence of the water molecules, the packing in (II) is quite different to that found in (Ia) and (Ib), and positional disorder involving a statistical superposition of two rotamers of IMPE(+), with different orientations of the indole fragment, was found. Crystals of (III) contain two independent ECMPE(+) rotamers with different orientations of their carbazole substituents. The cations are packed in stacks, with the iodide anions located in the channels between the stacks. In (III), the iodide was found to be disordered over two sites, with occupancies of 0.83 and 0.17.

Conformational polymorphs of 22-cyano-N-methyl-5-phenylpent-2-en-4-ynamide

Although the two polymorphic modifications, (I) and (II), of the title compound, C(13)H(10)N(2)O, crystallize in the same space group (P2(1)/c), their asymmetric units have Z' values of 1 and 2, respectively. These are conformational polymorphs, since the molecules in phases (I) and (II) adopt different rotations of the phenyl ring with respect the central 2-cyanocarboxyaminoprop-2-enyl fragment. Calculations of crystal packing using Cerius(2) [Molecular Simulations (1999). 9685 Scranton Road, San Diego, CA 92121, USA] have shown that (I) is more stable than (II), by 1.3 kcal mol(-1) for the crystallographically determined structures and by 1.56 kcal mol(-1) for the optimized structures (1 kcal mol(-1) = 4.184 kJ mol(-1)). This difference is mainly attributed to the different strengths of the hydrogen bonding in the two forms.