Drug-drug co-crystals: Temperature-dependent proton mobility in the molecular complex of isoniazid with 4-aminosalicylic acid

Two drug-drug co-crystals of the anti-tuberculosis drugs isoniazid (INH), pyrazinamide (PYR) and 4-aminosalicylic acid (PAS) are reported. The first is the 1 : 1 molecular complex of INH and PAS. The second is the monohydrate of the 1 : 1 complex of PYR and PAS. The crystal structures of both co-crystals are characterized by a number of hydrogen bonded synthons. Hydrogen bonding of the COOH center dot center dot center dot N-pyridine type is found in both cases. In the INH : PAS co-crystal, there are two symmetry independent COOH center dot center dot center dot center dot N-pyridine hydrogen bonds. In one of these, the H-atom is located on the carboxylic group and is indicative of a co-crystal. In the second case, partial proton transfer occurs across the hydrogen bond, and the extent of proton transfer depends on the temperature. This is more indicative of a salt. Drug-drug co-crystals may have some bearing in the treatment of tuberculosis.

Density-matrix renormalization-group study of low-lying excitations of polyacene within a Pariser-Parr-Pople model

We have carried out symmetrized density-matrix renormalization-group calculations to study the nature of excited states of long polyacene oligomers within a Pariser-Parr-Pople Hamiltonian. We have used the C-2 symmetry, the electron-hole symmetry, and the spin parity of the system in our calculations. We find that there is a crossover in the lowest dipole forbidden two-photon state and the lowest dipole allowed excited state with size of the oligomer. In the long system limit, the two-photon state lies below the lowest dipole allowed excited state. The triplet state lies well below the two-photon state and energetically does not correspond to its description as being made up of two triplets. These results are in agreement with the general trends in linear conjugated polymers. However, unlike in linear polyenes wherein the two-photon state is a localized excitation, we find that in polyacenes, the two-photon excitation is spread out over the system. We have doped the systems with a hole and an electron and have calculated the charge excitation gap. Using the charge gap and the optical gap, we estimate the binding energy of the 1(1)B(-) exciton to be 2.09 eV. We have also studied doubly doped polyacenes and find that the bipolaron in these systems, to be composed of two separated polarons, as indicated by the calculated charge-density profile and charge-charge correlation function. We have studied bond orders in various states in order to get an idea of the excited state geometry of the system. We find that the ground state, the triplet state, the dipole allowed state, and the polaron excitations correspond to lengthening of the rung bonds in the interior of the oligomer while the two-photon excitation corresponds to the rung bond lengths having two maxima in the system.

On the importance of backbone loop and peptide flip in the Walker sequence in F-1-ATPase action

The Walker sequence, GXXXXGKT, present in all the six subunits of F-1-ATPase exists in a folded form, known as phosphate-binding loop (P-loop). Analysis of the Ramachandran angles showed only small RMS deviation between the nucleotide-bound and nucleotide-free forms. This indicated a good overlap of the backbone loops. The catalytic beta-subunits (chains D, E and F) showed significant changes in the Ramachandran angles and the side chain torsion angles, but not the structural alpha-subunits (chains A, B and C). Most striking among these are the changes associated with Val160 and Gly161 corresponding to a flip in the peptide unit between them when a nucleotide is bound (chains D or F compared to nucleotide-free chain E). The conformational analysis further revealed a hitherto unnoticed hydrogen bond between amide-N of the flipped Gly161 and terminal phosphate-O of the nucleotide. This assigns a role for this conserved amino acid, otherwise ignored, of making an unusual direct interaction between the peptide backbone of the enzyme protein and the incoming nucleotide substrate. Significance of this interaction is enhanced, as it is limited only to the catalytic subunits, and also likely to involve a mechanical rotation of bonds of the peptide unit. Hopefully this is part of the overall events that link the chemical hydrolysis of ATP with the mechanical rotation of this molecule, now famous as tiny molecular motor.

One-dimensional zinc phosphates with linear chain structure

Three one-dimensional zinc phosphates, [C5N2H14][Zn(HPO4)2], I, [C10N4H26][Zn(HPO4)2].2H2O II, and [C4N2H6]2[Zn(HPO4)], III, have been prepared employing hydro/solvothermal methods in the presence of organic amines. While I and II consist of linear chains of corner-shared four-membered rings, III is a polymeric wire where the amine molecule is directly bonded to the metal center. The wire, as well as the chain in these structures, are held together by hydrogen bond interactions involving the amine and the framework oxygens. The polymeric zinc phosphate with wire-like architecture, III, is only the second example of such architecture. Crystal data: I, monoclinic, P21/c (no. 14), a=8.603(2), b=13.529(2), c=10.880(1) Å, β=94.9(1)°, V=1261.6(1) Å3, Z=4, ρcalc.=1.893 gcm−3, μ(MoKα)=2.234 mm−1, R1=0.032, wR2=0.086, [1532 observed reflections with I>2σ(I)], II, orthorhombic, Pbca (no. 61), a=8.393(1), b=15.286(1), c=22.659(1) Å, V=2906.9(2) Å3, Z=8, ρcalc.=1.794 gcm−3, μ(MoKα)=1.957 mm−1, R1=0.055, wR2=0.11, [1565 observed reflections with I>2σ(I) and III, monoclinic, P21/c (no. 14), a=8.241(1), b=13.750(2), c=10.572(1) Å, β=90.9(1)°, V=1197.7(2) Å3, Z=4, ρcalc.=1.805 gcm−3, μ(MoKα)=2.197 mm−1, R1=0.036, wR2=0.10, [1423 observed reflections with I>2σ(I)].

Revisiting the Mechanism of the Triosephosphate Isomerase Reaction: The Role of the Fully Conserved Glutamic Acid 97 Residue

An analysis of 503 available triosephosphate isomerase sequences revealed nine fully conserved residues. Of these, four residues-K12, H95, E97 and E165-are capable of proton transfer and are all arrayed around the dihydroxyacetone phosphate substrate in the three-dimensional structure. Specific roles have been assigned to the residues K12, H95 and E165, but the nature of the involvement of E97 has not been established. Kinetic and structural characterization is reported for the E97Q and E97D mutants of Plasmodium falciparum triosephosphate isomerase (Pf TIM). A 4000-fold reduction in k(cat) is observed for E97Q, whereas the E97D mutant shows a 100-fold reduction. The control mutant, E165A, which lacks the key catalytic base, shows an approximately 9000-fold drop in activity. The integrity of the overall fold and stability of the dimeric structure have been demonstrated by biophysical studies. Crystal structures of E97Q and E97D mutants have been determined at 2.0 angstrom resolution. In the case of the isosteric replacement of glutamic acid by glutamine in the E97Q mutant a large conformational change for the critical K12 side chain is observed, corresponding to a trans-to-gauche transition about the C gamma-C delta (chi(3)) bond. In the E97D mutant, the K12 side chain maintains the wild-type orientation, but the hydrogen bond between K12 and D97 is lost. The results are interpreted as a direct role for E97 in the catalytic proton transfer cycle. The proposed mechanism eliminates the need to invoke the formation of the energetically unfavourable imidazolate anion at H95, a key feature of the classical mechanism.

Conformational studies of the phosphazenyl side chain in cyclophosphazenes. Part 3. Crystal and molecular structure of N3P3Cl4(NEt2)(NPPh3)

The crystal structure of N3P3Cl4(NEt2)(NPPh3) has been determined. The crystals are orthorhombic, space group Pbca, with a= 8.208(1), b= 21.890(1), c= 31.722(2)Å, Z= 8, and m.p. = 146.5 °C. The structure was solved by direct methods and refined to a final R value of 0.045 for 2 025 independent reflections. The analysis reveals significant variations in the ring P–N bond lengths. The two nitrogenous substituents, NPPh3 and NEt2, reside on the same phosphorus atom. The latter, NEt2, has an almost exact type II conformation (the plane NC2 almost perpendicular to the local NPN plane)(the first observed for a dialkylamino-group in cyclophosphazenes), the former, NPPh3, deviates from type II towards type III (in type III the plane Pring–N–Pexo makes an angle of ca. 45° with the local N–P–N ring plane). The present structure is compared with others of triphenylphosphazenyl-cyclophosphazenes and the conformation of the NPPh3 substituent and its electron supply in the ground and perturbed states are discussed.

Microwave spectrum and structure of C(6)H(5)CCH center dot center dot center dot H(2)S complex

Rotational spectra of C(6)H(5)CCH center dot center dot center dot H(2)S, C(6)H(5)CCH center dot center dot center dot H(2)(34)S, C(6)H(5)CCH center dot center dot center dot HDS, C(6)H(5)CCH center dot center dot center dot D(2)S and C(6) H(5)CCD center dot center dot center dot H(2)S complexes have been observed using a pulsed nozzle Fourier transform microwave spectrometer. The observed spectrum is consistent with a structure in which hydrogen sulfide is located over the phenyl ring pi cloud and the distance between the centers of masses of the two monomers is 3.74 +/- 0.01 angstrom. In the complex, the H(2)S unit is shifted from the phenyl ring center towards the acetylene group. The vibrationally averaged structure has an effective Cs symmetry. Ab initio calculations were performed at MP2/aug-cc-pVDZ level of theory to locate the possible geometries of the complex. The calculations reveal the experimentally observed structure to be more stable than a coplanar arrangement of the monomers, which was observed for the C(6)H(5)CCH center dot center dot center dot H(2)O complex. Atoms in molecule theoretical analysis shows the presence of S-H center dot center dot center dot pi hydrogen bond. For the parent isotopologue, each transition frequency was found to split into two resulting from an interchange of the equivalent hydrogens of H(2)S unit in the complex. (C) 2011 Elsevier Inc. All rights reserved.

2-{[{2-Hydroxy-3-[2-methyl-5-(propan-2-yl)phenoxy]propyl}(pyridin -2-ylmethyl)amino]methyl}phenol

In the title racemic compound, C(26)H(32)N(2)O(3), an intramolecular O-H center dot center dot center dot N hydrogen bond is formed between the phenolic OH group and the tertiary amine N atom. Another O-H center dot center dot center dot N hydrogen bond that is formed between the OH group and the pyridine N atom links the molecules into a polymeric chain extending along the a axis. The structure is further stabilized by intramolecular and intermolecular C-H center dot center dot center dot O interactions.

Dynamic and thermodynamic anomalies of water at low temperatures: from bulk water to reverse micelles and DNA hydration layer

Liquid water is known to exhibit remarkable thermodynamic and dynamic anomalies, ranging from solvation properties in supercritical state to an apparent divergence of the linear response functions at a low temperature. Anomalies in various dynamic properties of water have also been observed in the hydration layer of proteins, DNA grooves and inside the nanocavity, such as reverse micelles and nanotubes. Here we report studies on the molecular origin of these anomalies in supercooled water, in the grooves of DNA double helix and reverse micelles. The anomalies have been discussed in terms of growing correlation length and intermittent population fluctuation of 4- and 5-coordinated species. We establish correlation between thermodynamic response functions and mean squared species number fluctuation. Lifetime analysis of 4- and 5-coordinated species reveals interesting differences between the role of the two species in supercooled and constrained water. The nature and manifestations of the apparent and much discussed liquid-liquid transition under confinement are found to be markedly different from that in the bulk. We find an interesting `faster than bulk' relaxation in reverse micelles which we attribute to frustration effects created by competition between the correlations imposed by surface interactions and that imposed by hydrogen bond network of water.

Phenylboronic acids in crystal engineering: Utility of the energetically unfavorable syn,syn-conformation in co-crystal design

Phenylboronic acids can exist, in principle, in three different conformers (syn,syn; syn,anti and anti,anti) with distinct energy profiles. In their native state, these compounds prefer the energetically favored syn, anti-conformation. In molecular complexes, however, the functionality exhibits conformational diversity. In this paper we report a series of co-crystals, with N-donor compounds, prepared by a design strategy involving the synthons based on the syn, syn-conformation of the boronic acid functionality. For this purpose, we employed compounds with the 1,2-diazo fragment (alprazolam, 1H-tetrazole, acetazolamide and benzotriazole), 1,10-phenanthroline and 2,2'-bipyridine for the co-crystallization experiments. However, our study shows that the mere presence of the 1,2-diazo fragment in the coformer does not guarantee the successful formation of co-crystals with a syn, syn-conformation of the boronic acid. [GRAPHICS] The -B(OH)(2) fragment makes unsymmetrical O-H center dot center dot center dot N heterosynthons with alprazolam (ALP) and 1,10-phenanthroline (PHEN). In the co-crystals of phenylboronic acids with 1H-tetrazole (TETR) and 2,2'-bipyridine (BPY), the symmetrical boronic acid dimer is the major synthon. In the BPY complex, boronic acid forms linear chains and the pyridine compound interacts with the lateral OH of boronic acid dimers that acts as a connector, thus forming a ladder structure. In the TETR complex, each heterocycle interacts with three boronic acids. While two boronic acids interact using the phenolic group, the third molecule generates O-H center dot center dot center dot N hydrogen bonds using the extra OH group, of -B(OH)(2) fragment, left after the dimer formation. Thus, although molecules were selected retrosynthetically with the 1,2-diazo fragment or with nearby hetero-atoms to induce co-crystal formation using the syn,syn-orientation of the -B(OH)(2) functionality, co-crystal formation is in fact selective and is probably driven by energy factors. Acetazolamide (ACET) contains self-complementary functional groups and hence creates stable homosynthons. Phenylboronic acids being weak competitors fail to perturb the homosynthons and hence the components crystallize separately. Therefore, besides the availability of possible hydrogen bond acceptors in the required position and orientation, the ability of the phenyl-boronic acid to perturb the existing interactions is also a prerequisite to form co-crystals. This is illustrated in the table below. In the case of ALP, PHEN and BPY, the native structures are stabilized by weak interactions and may be influenced by the boronic acid fragment. Thus phenylboronic acids can attain co-crystals with those compounds, wherein the cyclic O-H center dot center dot center dot N hydrogen bonds are stronger than the individual homo-interactions. This can lower the lattice energy of the molecular complex as compared with the individual crystals. [GRAPHICS] Phenylboronic acids show some selectivity in the formation of co-crystals with N-heterocycles. The differences in solubility of the components fall short to provide a possible reason for the selective formation of co-crystals only with certain compounds. These compounds, being weak acids, do not follow the Delta pK(a) analysis and hence fail to provide any conclusive observation. Theoretical results show that of the three conformers possible, the syn,anti conformer is the most stable. The relative stabilities of the three conformers syn,anti,syn,syn and anti,anti are 0.0, 2.18 and 3.14 kcal/mol, respectively. The theoretical calculations corroborate the fact that only energetically favorable synthons can induce the formation of heterosynthons, as in ALP and PHEN complexes. From a theoretical and structural analysis it is seen that phenylboronic acids will form interactions with those molecules wherein the heterocyclic and acidic fragments can interrupt the homosynthons. However, the energy profile is shallow and can be perturbed easily by the presence of competing functional groups (such as OH and COOH) in the vicinity. [GRAPHICS] .

NMRof peptides

Developments and applications of NMR spectroscopy especially with biomolecules has taken big strides over the decades. This review gives a brief overview of peptide analysis by NMR as carried out in the author’s laboratory. A brief introduction to peptide biomolecules and NMR useful parameters are discussed in the beginning. This is followed by diagnostics features observed in NMR for identification of secondary structures. It further goes on to show how a three dimensional structure could be obtained by all-important NOE and hydrogen bond information. Use of heteronuclear experiments, which could be done at natural abundance is also highlighted in getting more details of peptide structures.Applications using Solid state NMR at natural abundance in connecting peptide solution and x-ray structures is demonstrated with couple of examples.

Conformation-Changing Aggregation in Hydroxyacetone: A Combined Low-Temperature FTIR, Jet, and Crystallographic Study

Aggregation in hydroxyacetone (HA) is studied using low-temperature FTIR, supersonic jet expansion, and X-ray crystallographic (in situ cryocrystallization) techniques. Along with quantum chemical methods (MP2 and DFT), the experiments unravel the conformational preferences of HA upon aggregation to dinners and oligomers. The O-H center dot center dot center dot O=C intramolecular hydrogen bond present in the gas-phase monomer partially opens upon aggregation in supersonic expansions, giving rise to intermolecular cooperatively enhanced O-H center dot center dot center dot O-H hydrogen bonds in competition with isolated O-H center dot center dot center dot O=C hydrogen bonds. On the other hand, low-temperature IR studies on the neat solid and X-ray crystallographic data reveal that HA undergoes profound conformational changes upon crystallization, with the HOCC dihedral angle changing from similar to 0 degrees in the gas phase to similar to 180 degrees in the crystalline phase, hence giving rise to a completely new conformation. These conclusions are supported by theoretical calculations performed on the geometry derived from the crystalline phase.

From anomalies in neat liquid to structure, dynamics and function in the biological world

Water brings its remarkable thermodynamic and dynamic anomalies in the pure liquid state to biological world where water molecules face a multitude of additional interactions that frustrate its hydrogen bond network. Yet the water molecules participate and control enormous number of biological processes in manners which are yet to be understood at a molecular level. We discuss thermodynamics, structure, dynamics and properties of water around proteins and DNA, along with those in reverse micelles. We discuss the roles of water in enzyme kinetics, in drug-DNA intercalation and in kinetic-proof reading ( the theory of lack of errors in biosynthesis). We also discuss how water may play an important role in the natural selection of biomolecules. (C) 2011 Elsevier B. V. All rights reserved.

Hydration dynamics of protein molecules in aqueous solution: Unity among diversity

Dielectric dispersion and NMRD experiments have revealed that a significant fraction of water molecules in the hydration shell of various proteins do not exhibit any slowing down of dynamics. This is usually attributed to the presence of the hydrophobic residues (HBR) on the surface, although HBRs alone cannot account for the large amplitude of the fast component. Solvation dynamics experiments and also computer simulation studies, on the other hand, repeatedly observed the presence of a non-negligible slow component. Here we show, by considering three well-known proteins (lysozyme, myoglobin and adelynate kinase), that the fast component arises partly from the response of those water molecules that are hydrogen bonded with the backbone oxygen (BBO) atoms. These are structurally and energetically less stable than those with the side chain oxygen (SCO) atoms. In addition, the electrostatic interaction energy distribution (EIED) of individual water molecules (hydrogen bonded to SCO) with side chain oxygen atoms shows a surprising two peak character with the lower energy peak almost coincident with the energy distribution of water hydrogen bonded to backbone oxygen atoms (BBO). This two peak contribution appears to be quite general as we find it for lysozyme, myoglobin and adenylate kinase (ADK). The sharp peak of EIED at small energy (at less than 2 k(B)T) for the BBO atoms, together with the first peak of EIED of SCO and the HBRs on the protein surface, explain why a large fraction (similar to 80%) of water in the protein hydration layer remains almost as mobile as bulk water Significant slowness arises only from the hydrogen bonds that populate the second peak of EIED at larger energy (at about 4 k(B)T). Thus, if we consider hydrogen bond interaction alone, only 15-20% of water molecules in the protein hydration layer can exhibit slow dynamics, resulting in an average relaxation time of about 5-10 ps. The latter estimate assumes a time constant of 20-100 ps for the slow component. Interestingly, relaxation of water molecules hydrogen bonded to back bone oxygen exhibit an initial component faster than the bulk, suggesting that hydrogen bonding of these water molecules remains frustrated. This explanation of the heterogeneous and non-exponential dynamics of water in the hydration layer is quantitatively consistent with all the available experimental results, and provides unification among diverse features.

Relaxation and jump dynamics of water at the mica interface

The orientational relaxation dynamics of water confined between mica surfaces is investigated using molecular dynamics simulations. The study illustrates the wide heterogeneity that exists in the dynamics of water adjacent to a strongly hydrophilic surface such as mica. Analysis of the survival probabilities in different layers is carried out by normalizing the corresponding relaxation times with bulk water layers of similar thickness. A 10-fold increase in the survival times is observed for water directly in contact with the mica surface and a non-monotonic variation in the survival times is observed moving away from the mica surface to the bulk-like interior. The orientational relaxation time is highest for water in the contact layer, decreasing monotonically away from the surface. In all cases the ratio of the relaxation times of the 1st and 2nd rank Legendre polynomials of the HH bond vector is found to lie between 1.5 and 1.9 indicating that the reorientational relaxation in the different water layers is governed by jump dynamics. The orientational dynamics of water in the contact layer is particularly novel and is found to undergo distinct two-dimensional hydrogen bond jump reorientational dynamics with an average waiting time of 4.97 ps. The waiting time distribution is found to possess a long tail extending beyond 15 ps. Unlike previously observed jump dynamics in bulk water and other surfaces, jump events in the mica contact layer occur between hydrogen bonds formed by the water molecule and acceptor oxygens on the mica surface. Despite slowing down of the water orientational relaxation near the surface, life-times of water in the hydration shell of the K ion are comparable to that observed in bulk salt solutions. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4717710]