Conformational Characterization of S-Methyl Dithiocarbazate by Infrared Spectra and Vibrational Assignments

Infrared spectra are recorded for S-methyl dithiocarbazate and its N-deuterated compound in two molecular conformations in the solid state and in solution between 4000 and 30 cm−1. The assignments have been supported from a complete normal coordinate analysis; the conformation sensitive bands of the –CSNHNH2 grouping are discussed. The assignments are compared with those of related molecules to check the internal consistency and to obtain the pattern of the characteristic bands of thiocarbazoyl (–CSNHNH2) group. The magnitudes of the C–N and S–CH3 torsional barriers are estimated from the force constants.

Mutation analysis of TGF-β signaling pathway genes among Finnish patients with primary pulmonary hypertension and hereditary hemorrhagic telangiectasia

Primary pulmonary hypertension (PPH), or according to the recent classification idiopathic pulmonary hypertension (IPAH), is a rare, progressive disease of pulmonary vasculature leading to pulmonary hypertension and right heart failure. Most of the patients are sporadic but in about 6% of cases the disease is familial (FPPH). In 2000 two different groups identified the gene predisposing to PPH. This gene, Bone morphogenetic protein receptor type 2 (BMPR2), encodes a subunit of transforming growth factor β (TGF-β) receptor complex. There is a genetic connection between PPH and hereditary hemorrhagic telangiectasia (HHT), a bleeding disorder characterized by local telangiectasias and sometimes with pulmonary hypertension. In HHT, mutations in ALK1 (activin like kinase type 1) and Endoglin, another members of the TGF-β signaling pathway are found. In this study we identified all of the Finnish PPH patients for the years 1986-1999 using the hospital discharge registries of Finnish university hospitals. During this period we found a total of 59 confirmed PPH patients: 55 sporadic and 4 familial representing 3 different families. In 1999 the prevalence of PPH was 5.8 per million and the annual incidence varied between 0.2-1.3 per million. Among 28 PPH patients studied, heterozygous BMPR2 mutations were found in 12% (3/26) of sporadic patients and in 33% of the PPH families (1/3). All the mutations found were different. Large deletions of BMPR2 were excluded by single-stranded chain polymomorphism analysis. As a candidate gene approach we also studied ALK1, Endoglin, Bone Morphogenetic Receptor Type IA (BMPR1A or ALK3), Mothers Against Decapentaplegic Homolog 4 (SMAD4) and Serotonine Transporter Gene (SLC6A4) using single-strand conformational polymorphism (SSCP) analysis and direct sequencing. Among patients and family members studied, we found two mutations in ALK1 in two unrelated samples. We also identified all the HHT patients treated at the Department of Otorhinolaryngology at Helsinki University Central Hospital between the years of 1990-2005 and 8 of the patients were studied for Endoglin and ALK1 mutations using direct sequencing. A total of seven mutations were found and all the mutations were different. The absence of a founder mutation in the Finnish population in both PPH and HHT was somewhat surprising. This suggests that the mutations of BMPR2, ALK1 and Endoglin are quite young and the older mutations have been lost due to repetitive genetic bottlenecks and/or negative selection. Also, other genes than BMPR2 may be involved in the pathogenesis of PPH. No founder mutations were found in PPH or HHT and thus no simple genetic test is available for diagnostics.

Conformationally Restricted Nucleocyclitols: a Study into their Conformational Preferences and Supramolecular Architecture in the Solid State

With the intent of probing the feasibility of employing annulation as a tactic to engender axial rich conformations in nucleoside analogues, two adenine-derived, ``conformationally restricted'' nucleocylitols, 9 and 10, have been conceptualized as representatives of a hitherto unexplored class of nucleic acid base-cyclitol hybrids. A general synthetic strategy, with an inherent scope for diversification, allowed rapid functionalization of indane and tetralin to furnish 9 and 10 respectively in fair yield. Single-crystal X-ray diffraction analysis revealed that the two nucleocyclitols under study, though homologous, present completely dissimilar modes of molecular packing, marked, in particular, by the nature of involvement of the adenynyl NH2 group in the supramolecular assembly. In addition, the crystal structures of 9 and 10 also exhibit two different conformations of the functionalized cyclohexane ring. Thus, while the six-membered carbocycle in cyclopenta-annulated 9 exists in the expected chair (C) conformation that in cyclohexaannulated 10, which crystallizes as a dihydrate, shows an unusual twist-boat (TB) conformation. From a close analysis of the (HNMR)-H-1 spectroscopic data recorded for 9 and 10 in CD3OD, it was possible to put forth a putative explanation for the uncanny conformational preferences of crystalline 9 and 10.

Conformational Choices for the Stereochemically Constrained gamma-Amino Acid Residue Gabapentin: Theoretical Studies and Correlation With Experimental Results

Gabapentin (1-aminomethylcyclohexaneacetic acid, Gpn) is an achiral, conformationally constrained gamma amino acid residue. A survey of available crystal structures of Gpn peptides reveals that the torsion angles about the C-gamma-C-beta (theta(1)) and C-beta-C-alpha(theta(2)) bonds are overwhelmingly limited to gauche, gauche (g(+)g(+)/g(-)g(-)) conformations. The Gpn residue forms C-7 and C-9 hydrogen bonds in which the donor and acceptor atoms come from the flanking peptide units. In combination with alpha amino acid residues alpha gamma and gamma alpha segments can adopt C-12 hydrogen bonded structures. The conformational choices available to the Gpn residue have been probed using energy calculations, adopting a grid search strategy. Ramachandran phi-psi maps have been constructed for fixed values of theta(1) and theta(2), corresponding to the gauche and trans conformations. The sterically allowed and energetically favorable regions of conformational space have been defined and experimental observations compared. C-7 and C-9 hydrogen bonded conformational families have been identified using a grid search approach in which theta(1) and theta(2) values are varied over a range of +/- 10 degrees about ideal values at 1 degrees intervals. The theoretical analysis together with experimental observations for 59 Gpn residues from 35 crystal structures permits definition of the limited range of conformational possibilities at this gamma amino acid residue. .

Conformational studies on cyclic dipeptides

An analysis of 11 crystal structures of cyclic dipeptides so far reported in the literature is made, with main reference to the internal parameters of these molecules. Preferred conformations of the side chains of cyclic dipeptides with different α-amino acid residues have been studied by classical energy calculations. The possible conformations of the DKP ring are also studied. The significance of the non-bonded interaction in deciding the pathway for conformational change has also been investigated. The agreement between theoretical results and experimental observations is quite good, both with respect to the conformation of these molecules as well as the enthalpy difference as estimated from n.m.r. studies between different conformers.

Allostery and conformational free energy changes in human tryptophanyl-tRNA synthetase from essential dynamics and structure networks

The interdependence of the concept of allostery and enzymatic catalysis, and they being guided by conformational mobility is gaining increased prominence. However, to gain a molecular level understanding of llostery and hence of enzymatic catalysis, it is of utter importance that the networks of amino acids participating in allostery be deciphered. Our lab has been exploring the methods of network analysis combined with molecular dynamics simulations to understand allostery at molecular level. Earlier we had outlined methods to obtain communication paths and then to map the rigid/flexible regions of proteins through network parameters like the shortest correlated paths, cliques, and communities. In this article, we advance the methodology to estimate the conformational populations in terms of cliques/communities formed by interactions including the side-chains and then to compute the ligand-induced population shift. Finally, we obtain the free-energy landscape of the protein in equilibrium, characterizing the free-energy minima accessed by the protein complexes. We have chosen human tryptophanyl-tRNA synthetase (hTrpRS), a protein esponsible for charging tryptophan to its cognate tRNA during protein biosynthesis for this investigation. This is a multidomain protein exhibiting excellent allosteric communication. Our approach has provided valuable structural as well as functional insights into the protein. The methodology adopted here is highly generalized to illuminate the linkage between protein structure networks and conformational mobility involved in the allosteric mechanism in any protein with known structure.

Experimental and Theoretical Charge Density Analysis of Polymorphic Structures: The Case of Coumarin 314 Dye

Experimental charge density distributions in two known conformational polymorphs (orange and yellow) of coumarin 314 dye are analyzed based on multipole modeling of X-ray diffraction data collected at 100 K. The experimental results are compared with the charge densities derived from multipole modeling of theoretical structure factors obtained from periodic quantum calculation with density functional theory (DFT) method and B3LYP/6-31G(d,p) level of theory. The presence of disorder at the carbonyl oxygen atom of ethoxycarbonyl group in the yellow form, which was not identified earlier, is addressed here. The investigationof intermolecular interactions, based on Hirshfeld surface analysis and topological properties via quantum theory of atoms in molecule and total electrostatic interaction energies, revealed significant differences between the polymorphs. The differences of electrostatic nature in these two polymorphic forms were unveiled via construction of three-dimensional deformation electrostatic potential maps plotted over the molecular surfaces. The lattice energies evaluated from ab initio calculations on the two polymorphic forms indicate that the yellow form is likely to be the most favorable thermodynamically. The dipole moments derived from experimental and theoretical charge densities and also from Lorentz tensor approach are compared with the single-molecule dipole moments. In each case, the differences of dipole moments between the polymorphs are identified.

Conformational characteristics of asparaginyl residues in proteins

Backbone conformations at 1064 asparaginyl residues in 123 non-homologous, high-resolution X-ray structures of proteins were analysed. Asn adopts conformations in left-handed x-helical region and other partially allowed regions in the Ramachandran map more readily than any other non-glycyl residue. Asn conformational clusters in the (phi,psi) regions of left-handed alpha-helix, right-handed alpha-helix and extended (beta) strands were investigated in detail for their occurrence in various secondary structures, especially in beta-turn regions. Preferences were observed for Asn conformations in different positions in various beta-turn types, including the first and fourth positions of the turn. Asparaginyl residues with extended conformations are found to occur frequently in irregular regions, although they are expected to occur predominantly in extended strands or in the third position of type II beta-turns. Asn conformations at the N-cap positions of helices strongly prefer extended conformation than alpha(L), which seems to be characteristic of non-glycyl residues at that position. In the linkers connecting two extended strands and those connecting an alpha-helix and an extended strand, Asn with alpha(L) or alpha(R) conformation is more favoured than Asn with the beta-conformation. Analysis of Asn-Asn doublets and Asn-X-Asn triplets permitted identification of conformational families in such sequences. Results of this investigation provide useful hints in modelling Asn-rich regions in proteins such as malaria parasite coat protein. (C) Munksgaard 1994.

Stereochemical analysis of the antigenic tip of the V3 loop peptide of HIV-1 gp120

A novel multiple turn conformation has been observed for a segment GPGRAFY in the crystal structure of a complex of HIV-1 gp120 V3 loop peptide with the Fab fragment of a neutralizing antibody [Ghiara ct al. (1994) Science 264, 82-85]. A structural motif has been defined for the peptide segment, employing idealized backbone conformations characterized by ranges of virtual C-alpha torsion angles and bond angles. A search of 122 high-resolution protein crystal structures has permitted identification of 24 examples of similar structural motifs. Two major conformational families have been identified, which differ primarily in the conformation at residue 3. The observed conformation at residue 3 in family 1 is left-handed helical (alpha(L)) and that in family 2 is right-handed helical (alpha(R)). Of the 10 examples in family 1, 9 examples have Gly residues at position 3. Of the 12 examples in family 2, 7 examples have Asn/Asp at position 3. Computer modeling of the V3 loop tip sequence using the two backbone conformational families as starting points leads to minimum-energy conformations in which antigenically important side-chains occupy similar spatial arrangements. This stereochemical analysis of the V3 loop tip sequence suggests a rational basis for the design of synthetic analog peptides for use as viral antagonists or synthetic antigens. (C) Munksgaard 1995.

Stereochemical Analysis of Higher α,α-Dialkylglycine Containing Peptides. Characterization of Local Helical Conformations at Dipropylglycine Residues and Observation of a Novel Hydrated Multiple β-Turn Structure in Crystals of a Glycine Rich Peptide

The peptide Boc-Gly-Dpg-Gly-Gly-Dpg-Gly-NHMe (1) has been synthesized to examine the conformational preferences of Dpg residues in the context of a poor helix promoting sequence. Single crystals of 1 were obtained in the space group P21/c with a = 13.716(2) Å, b = 12.960(2) Å, c = 22.266(4) Å, and β = 98.05(1)°; R = 6.3% for 3660 data with |Fo| > 4σ. The molecular conformation in crystals revealed that the Gly(1)-Dpg(2) segment adopts φ, ψ values distorted from those expected for an ideal type II‘ β-turn (φGly(1) = +72.0°, ψGly(1) = −166.0°; φDpg(2) = −54.0°, ψDpg(2) = −46.0°) with an inserted water molecule between Boc-CO and Gly(3)NH. The Gly(3)-Gly(4) segment adopts φ, ψ values which lie broadly in the right handed helical region (φGly(3) = −78.0°, ψGly(3) = −9.0°; φGly(4) = −80.0°, ψGly(4) = −18.0°). There is a chiral reversal at Dpg(5) which takes up φ, ψ values in the left handed helical region. The Dpg(5)-Gly(6) segment closely resembles an ideal type I‘ β-turn (φDpg(5) = +56.0°, ψDpg(5) = +32.0°; φGly(6) = +85.0°, ψGly(6) = −3.0°). Molecules of both chiral senses are found in the centrosymmetric crystal. The C-terminus forms a hydrated Schellman motif, with water insertion into the potential 6 → 1 hydrogen bond between Gly(1)CO and Gly(6)NH. NMR studies in CDCl3 suggest substantial retention of the multiple turn conformation observed in crystals. In solution the observed NOEs support local helical conformation at the two Dpg residues.

Mechanism of lectin-cell interactions: thermodynamic and kinetic analysis of the binding of winged bean acidic lectin (WBA II) to human erythrocytes

In order to identify the forces involved in the binding and to understand the mechanism involved, equilibrium and kinetic studies were performed on the binding of the winged bean acidic lectin to human erythrocytes. The magnitudes of delta S and delta H were positive and negative respectively, an observation differing markedly from the lectin-simple sugar interactions where delta S and delta H are generally negative. Analysis of the sign and magnitudes of these values indicate that ionic and hydrogen bonded interactions prevail over hydrophobic interactions resulting in net -ve delta H (-37.12 kJ.mol-1) and +ve delta S (14.4 J.mole-1 K-1 at 20 degrees C), thereby suggesting that this entropy driven reaction also reflects conformational changes in the lectin and/or the receptor. Presence of two kinds of receptors for WBA II on erythrocytes, as observed by equilibrium studies, is consistent with the biexponential dissociation rate constants (at 20 degrees C K1 = 1.67 x 10(-3) M-1 sec-1 and K2 = 11.1 x 10(-3) M-1 sec-1). These two rate constants differed by an order of magnitude accounting for the difference in the association constants of the two receptors of WBA II. However, the association process remains monoexponential suggesting no observable difference in the association rates of the lectin molecule with both the receptors, under the experimental conditions studied. The thermodynamic parameters calculated from kinetic data correlate well with those observed by equilibrium. A two-step binding mechanism is proposed based on the kinetic parameters for WBA II-receptor interaction

Comparative analysis of thermophilic and mesophilic proteins using Protein Energy Networks

Background: Thermophilic proteins sustain themselves and function at higher temperatures. Despite their structural and functional similarities with their mesophilic homologues, they show enhanced stability. Various comparative studies at genomic, protein sequence and structure levels, and experimental works highlight the different factors and dominant interacting forces contributing to this increased stability. Methods: In this comparative structure based study, we have used interaction energies between amino acids, to generate structure networks called as Protein Energy Networks (PENs). These PENs are used to compute network, sub-graph, and node specific parameters. These parameters are then compared between the thermophile-mesophile homologues. Results: The results show an increased number of clusters and low energy cliques in thermophiles as the main contributing factors for their enhanced stability. Further more, we see an increase in the number of hubs in thermophiles. We also observe no community of electrostatic cliques forming in PENs. Conclusion: In this study we were able to take an energy based network approach, to identify the factors responsible for enhanced stability of thermophiles, by comparative analysis. We were able to point out that the sub-graph parameters are the prominent contributing factors. The thermophiles have a better-packed hydrophobic core. We have also discussed how thermophiles, although increasing stability through higher connectivity retains conformational flexibility, from a cliques and communities perspective.

Comparative analysis of thermophilic and mesophilic proteins using Protein Energy Networks

Background: Thermophilic proteins sustain themselves and function at higher temperatures. Despite their structural and functional similarities with their mesophilic homologues, they show enhanced stability. Various comparative studies at genomic, protein sequence and structure levels, and experimental works highlight the different factors and dominant interacting forces contributing to this increased stability. Methods: In this comparative structure based study, we have used interaction energies between amino acids, to generate structure networks called as Protein Energy Networks (PENs). These PENs are used to compute network, sub-graph, and node specific parameters. These parameters are then compared between the thermophile-mesophile homologues. Results: The results show an increased number of clusters and low energy cliques in thermophiles as the main contributing factors for their enhanced stability. Further more, we see an increase in the number of hubs in thermophiles. We also observe no community of electrostatic cliques forming in PENs. Conclusion: In this study we were able to take an energy based network approach, to identify the factors responsible for enhanced stability of thermophiles, by comparative analysis. We were able to point out that the sub-graph parameters are the prominent contributing factors. The thermophiles have a better-packed hydrophobic core. We have also discussed how thermophiles, although increasing stability through higher connectivity retains conformational flexibility, from a cliques and communities perspective.

Analysis of short interproton distances in proline peptides as a guide in the interpretation of nuclear overhauser effects

The conformational dependence of interproton distances in model proline peptides has been investigated in order to facilitate interpretation of the results of Nuclear Overhauser Effect (NOE) studies on such peptides. For this purpose two model systems, namely, Ac-Pro-NHMe and Ac-Pro-X-NHMe have been chosen and used. In the former, short interproton distances detectable in NOE experiments permit a clear distinction between conformations with Pro ψ = -300 (helical region) and those in which ψ is around 1200 (polyproline region). For the latter, the variation of distances between the protons of methyl amide and the Pro ring have been studied by superimposing on the Ramachandran map in the (φ3, ψ3) plane. The results show that β-turns and non-β-turn conformations can be readily distinguished from NOE data and such long range NOEs should be detectable for specific non-β-turn conformations. NOEs involving Cβ and Cγ protons are particularly sensitive to the state of pyrrolidine ring puckering.

Elucidation of the conformational free energy landscape in H.pylori LuxS and its implications to catalysis

Background: One of the major challenges in understanding enzyme catalysis is to identify the different conformations and their populations at detailed molecular level in response to ligand binding/environment. A detail description of the ligand induced conformational changes provides meaningful insights into the mechanism of action of enzymes and thus its function. Results: In this study, we have explored the ligand induced conformational changes in H. pylori LuxS and the associated mechanistic features. LuxS, a dimeric protein, produces the precursor (4,5-dihydroxy-2,3-pentanedione) for autoinducer-2 production which is a signalling molecule for bacterial quorum sensing. We have performed molecular dynamics simulations on H. pylori LuxS in its various ligand bound forms and analyzed the simulation trajectories using various techniques including the structure network analysis, free energy evaluation and water dynamics at the active site. The results bring out the mechanistic details such as co operativity and asymmetry between the two subunits, subtle changes in the conformation as a response to the binding of active and inactive forms of ligands and the population distribution of different conformations in equilibrium. These investigations have enabled us to probe the free energy landscape and identify the corresponding conformations in terms of network parameters. In addition, we have also elucidated the variations in the dynamics of water co-ordination to the Zn2+ ion in LuxS and its relation to the rigidity at the active sites. Conclusions: In this article, we provide details of a novel method for the identification of conformational changes in the different ligand bound states of the protein, evaluation of ligand-induced free energy changes and the biological relevance of our results in the context of LuxS structure-function. The methodology outlined here is highly generalized to illuminate the linkage between structure and function in any protein of known structure.