Comparative modeling and genomics for galactokinase (Gal1p) enzyme

The Gal1p (Galactokinase) protein is known for regulation of D-galactose metabolism. It catalyzes the formation of galactose -1-phosphate from alpha - D-galactose, which is an important step in galactose catabolism. The knowledge of Gal1p protein structure, its protein interacting partners and enumeration of functional site residues will provide great insight in understanding the functional role of Gal1p. These studies are lacking in case of the Gal11p kinase enzyme. Structure of this enzyme has already been determined in S. cerevisiae, however, no structural information for this protein is available for K. lactis and E. coli. We used the homology modeling based approach to model the structures of Gal1p for K. lactis and E. coli. Furthermore, functional residues were predicted for these Gal1 proteins and the strength of interaction between Gal1p and other Gal proteins was determined by protein-protein interaction studies via patchdock software. The interaction studies revealed that the affinity for Gal1p for other Gal proteins varies in different organisms. Sequence and structural based comparison of Gal1p kinase enzyme showed that the orthologs in K.lactis and S. cervisiae are more similar to each other as compared to the ortholog in E. coli. These studies carried out by us will help in better understanding of the galactose metabolism. Our sequence and structure comparison studies revealed that Human Gal1p shows more homology for Gal1p protein of E. coli. The above studies may be applied to Human Gal1p, where it can help in gaining useful insight into Galactosemia disease.     

Open interface and large quaternary structure movements in 3D domain swapped proteins: insights from molecular dynamics simulations of the C-terminal swapped dimer of ribonuclease A

Bovine pancreatic ribonuclease (RNase A) forms two three-dimensional (3D) domain swapped dimers. Crystallographic investigations have revealed that these dimers display completely different quaternary structures: one dimer (N-dimer), which presents the swapping of the N-terminal helix, is characterized by a compact structure, whereas the other (C-dimer), which is stabilized by the exchange of the C-terminal end, shows a rather loose assembly of the two subunits. The dynamic properties of monomeric RNase A and of the N-dimer have been extensively characterized. Here, we report a molecular dynamics investigation carried out on the C-dimer. This computational experiment indicates that the quaternary structure of the C-dimer undergoes large fluctuations. These motions do not perturb the proper folding of the two subunits, which retain the dynamic properties of RNase A and the N-dimer. Indeed, the individual subunits of the C-dimer display the breathing motion of the beta-sheet structure, which is important for the enzymatic activity of pancreatic-like ribonucleases. In contrast to what has been observed for the N-dimer, the breathing motion of the two subunits of the C-dimer is not coupled. This finding suggests that the intersubunit communications in a 3D domain swapped dimer strongly rely on the extent of the interchain interface. Furthermore, the observation that the C-dimer is endowed with a high intrinsic flexibility holds interesting implications for the specific properties of 3D domain swapped dimers. Indeed, a survey of the quaternary structures of the other 3D domain swapped dimers shows that large variations are often observed when the structural determinations are conducted in different experimental conditions. The 3D domain swapping phenomenon coupled with the high flexibility of the quaternary structure may be relevant for protein-protein recognition, and in particular for the pathological aggregations.     

Structure of the house dust mite allergen Der f 2: implications for function and molecular basis of IgE cross-reactivity

The X-ray structure of the group 2 major allergen from Dermatophagoides farinae (Der f 2) was determined to 1.83 A resolution. The overall Der f 2 structure comprises a single domain of immunoglobulin fold with two anti-parallel beta-sheets. A large hydrophobic cavity is formed in the interior of Der f 2. Structural comparisons to distantly related proteins suggest a role in lipid binding. Immunoglobulin E (IgE) cross-reactivity between group 2 house dust mite major allergens can be explained by conserved surface areas representing IgE binding epitopes.     

Molecular dynamics simulations of glycosyltransferase LgtC

Molecular dynamics simulations have been performed on fully solvated alpha-(1-->4)-galactosyltransferase LgtC from Neisseria meningitidis with and without the donor substrate UDP-Gal and in the presence of the manganese ion. The analysis of the trajectories revealed a limited movement in the loop X (residues 75-80) and a larger conformational change in the loop Y (residues 246-251) in the simulation, when UDP-Gal was not present. In this case, the loops X and Y open by almost 10A, exposing the active site to the solvent. The 'hinge region' responsible for the opening is composed of residues 246-247. We have also analyzed the behavior of the manganese ion in the simulations. The coordination number is 6 when UDP-Gal is present and it increases to 7 when it is absent. In the latter case, three water molecules become coordinated to the ion. In both cases, the coordination is very stable implying that the manganese ion is tightly bound in the active site of the enzyme even if UDP-Gal is not present. Further analysis of the structural water molecules location confirmed that the mobility of water molecules in the active site and the accessibility of this site for solvent are higher in the absence of the substrate.     

Activation of the Ghrelin Receptor is Described by a Privileged Collective Motion: A Model for Constitutive and Agonist-induced Activation of a Sub-class A G-Protein Coupled Receptor (GPCR)

Three homology models of the human ghrelin receptor (GHS-R1a) have been generated from the available X-ray structures of rhodopsin (RHO model), opsin (OPS model) and beta-2 adrenergic receptor (B2 model). The latter was used as a starting point for combined molecular dynamics simulation (MDS) and full atom normal modes analysis (NMA). A low-frequency normal mode (mode 16) perfectly reproduced the intracellular motions observed between B2 and RHO models; in the opposite direction along the same mode, the generated structures are closer to the OPS model, suggesting a direct link with GHS-R1a activation. This was in agreement with motions of the seven transmembranous segments, increase of the solvent accessibility of the 140-ERY-142 sequence, and flip of the Trp276 (C WLP) residue, some features related to GPCRs activation. According to our model, His280 was proposed to stabilize Trp276 in the active state; this was verified by site-directed mutagenesis and biochemical characterization of the resulting H280A and H280S mutants, which were fully functional but sharing an important decrease of their basal activities. Docking performed with short ghrelin derivatives Gly-Ser-Ser _[octa]-Phe-NH ₂ and Gly-Ser-Ser _[octa]-Phe-Leu-NH ₂ allowed the identification of a robust position of these peptides in the active site of the receptor. This model was refined by MDS and validated by docking experiments performed on a set of 55 ghrelin receptor ligands based on the 1,2,4- triazole scaffold. Finally, NMA performed on the obtained peptide-receptor complex suggested stabilization of the Trp276 residue and of the whole receptor in the active state, preventing the motion observed along mode 16 computed for the unbound receptor. Our results show that NMA offers a powerful approach to study the conformational diversity and the activation mechanism of GPCRs.     

Chaperone-like activity of α-cyclodextrin via hydrophobic nanocavity to protect native structure of ADH

The chaperone action of α-cyclodextrin (α-CyD), based on providing beneficial microenvironment of hydrophobic nanocavity to form molecular complex with alcohol dehydrogenase (ADH) was examined by experimental and computational techniques. The results of UV-vis and dynamic light scattering (DLS) indicated that the chaperone-like activity of α-CyD depends on molecular complex formation between α-CyD and ADH, which caused to decrease the amount and size of polymerized molecules. Computational calculations of molecular dynamic (MD) simulations and blind docking (BD) demonstrated that α-CyD acts as an artificial chaperone because of its high affinity to the region of ADH’s two chains interface. The hydrophobic nanocavity of α-CyD has the ability to form inclusion complex due to the presence of phenyl ring of aromatic phenylalanine (Phe) residue in the dimeric intersection area. Delocalization of ADH subunits, which causes the exposure of Phe110, takes part in the enzyme polymerization and has proven to be beneficial for aggregation inhibition and solubility enhancement within the host α-CyD-nanocavity.     

Structure of translation initiation factor 1 from Mycobacterium tuberculosis and inferred binding to the 30S ribosomal subunit

The crystal structure of the free form of IF1 from Mycobacterium tuberculosis has been determined at 1.47 Å resolution. The structure adopts the expected OB fold and matches the high structural conservation among IF1 orthologues. In order to further explore the function of Mtb-IF1, we built a model of its interaction with the 30S ribosomal subunit based on the crystal structure of the complex from Thermus thermophilus. The model suggests that several functionally important side chain residues undergo large movements while the rest of the protein in complex shows only very limited conformational change as compared to its form in solution.     

The regulatory factor SipA is a highly stable β-II class protein with a SH3 fold

The small regulator SipA, interacts with the ATP-binding domain of non-bleaching sensor histidine kinase (NblS), the most conserved histidine kinase in cyanobacteria. NblS regulates photosynthesis and acclimation to a variety of environmental conditions. We show here that SipA is a highly stable protein in a wide pH range, with a thermal denaturation midpoint of 345 K. Circular dichroism and 1D ¹H NMR spectroscopies, as well as modelling, suggest that SipA is a β-II class protein, with short strands followed by turns and long random-coil polypeptide patches, matching the SH3 fold. The experimentally determined m-value and the heat capacity change upon thermal unfolding (ΔC_p) closely agreed with the corresponding theoretical values predicted from the structural model, further supporting its accuracy.     

Disulfide linkage in the coiled-coil domain of subunit H of A1AO ATP synthase from Methanocaldococcus jannaschii and the NMR structure of the C-terminal segment H85-104

The C-terminal residues 98-104 are important for structure stability of subunit H of A₁A_O ATP synthases as well as its interaction with subunit A. Here we determined the structure of the segment H_85-104 of H from Methanocaldococcus jannaschii, showing a helix between residues Lys90 to Glu100 and flexible tails at both ends. The helix-helix arrangement in the C-terminus was investigated by exchange of hydrophobic residues to single cysteine in mutants of the entire subunit H (H_I93C, H_L96C and H_L98C). Together with the surface charge distribution of H_85-104, these results shine light into the A-H assembly of this enzyme.     

A model for co-translational translocation: ribosome-regulated nascent polypeptide translocation at the protein-conducting channel

The protein-conducting channel (PCC) must allow both the translocation of soluble polypeptide regions across, and the lateral partitioning of hydrophobic transmembrane helices (TMHs) into, the membrane. We have analyzed existing structures of ribosomes and ribosome-PCC complexes and observe conformational changes suggesting that the ribosome may sense and orient the nascent polypeptide and also facilitate conformational changes in the PCC, subsequently directing the nascent polypeptide into the appropriate PCC-mediated translocation mode. The PCC is predicted to be able to accommodate one central, consolidated channel or two segregated pores with different lipid accessibilities, which may enable the lipid-mediated partitioning of a TMH from one pore, while the other, aqueous, pore allows translocation of a hydrophilic polypeptide segment. Our hypothesis suggests a plausible mechanism for the transitioning of the PCC between different configurations.