Visual Analysis of Biomolecular Surfaces

Surface models of biomolecules have become crucially important for the study and understanding of interaction between biomolecules and their environment. We argue for the need for a detailed understanding of biomolecular surfaces by describing several applications in computational and structural biology. We review methods used to model, represent, characterize, and visualize biomolecular surfaces focusing on the role that geometry and topology play in identifying features on the surface. These methods enable the development of efficient computational and visualization tools for studying the function of biomolecules.

MolBridge: a program for identifying nonbonded interactions in small molecules and biomolecular structures

Identification and analysis of nonbonded interactions within a molecule and with the surrounding molecules are an essential part of structural studies, given the importance of these interactions in defining the structure and function of any supramolecular entity. MolBridge is an easy to use algorithm based purely on geometric criteria that can identify all possible nonbonded interactions, such as hydrogen bond, halogen bond, cation-pi, pi-pi and van der Waals, in small molecules as well as biomolecules. The user can either upload three-dimensional coordinate files or enter the molecular ID corresponding to the relevant database. The program is available in a standalone form and as an interactive web server with Jmol and JME incorporated into it. The program is freely downloadable and the web server version is also available at http://nucleix.mbu.iisc.ernet.in/molbridge/index.php.

Hierarchical sampling for metastable conformers determines biomolecular recognition: the case of malectin and diglucosylated N-glycan interactions

Structural information over the entire course of binding interactions based on the analyses of energy landscapes is described, which provides a framework to understand the events involved during biomolecular recognition. Conformational dynamics of malectin's exquisite selectivity for diglucosylated N-glycan (Dig-N-glycan), a highly flexible oligosaccharide comprising of numerous dihedral torsion angles, are described as an example. For this purpose, a novel approach based on hierarchical sampling for acquiring metastable molecular conformations constituting low-energy minima for understanding the structural features involved in a biologic recognition is proposed. For this purpose, four variants of principal component analysis were employed recursively in both Cartesian space and dihedral angles space that are characterized by free energy landscapes to select the most stable conformational substates. Subsequently, k-means clustering algorithm was implemented for geometric separation of the major native state to acquire a final ensemble of metastable conformers. A comparison of malectin complexes was then performed to characterize their conformational properties. Analyses of stereochemical metrics and other concerted binding events revealed surface complementarity, cooperative and bidentate hydrogen bonds, water-mediated hydrogen bonds, carbohydrate-aromatic interactions including CH-pi and stacking interactions involved in this recognition. Additionally, a striking structural transition from loop to beta-strands in malectin CRD upon specific binding to Dig-N-glycan is observed. The interplay of the above-mentioned binding events in malectin and Dig-N-glycan supports an extended conformational selection model as the underlying binding mechanism.

Structural and optical investigation of semiconductor CdSe/CdS core-shell quantum dot thin films

Highly luminescent CdSe/CdS core-shell nanocrystals have been assembled on indium tin oxide (ITO) coated glass substrates using a wet synthesis route. The physical properties of the quantum dots (QD) have been investigated using X-ray diffraction, transmission electron microscopy and optical absorption spectroscopy techniques. These quantum dots showed a strong enhancement in the near band edge absorption. The in situ luminescence behavior has been interpreted in the light of the quantum confinement effect and induced strain in the core-shell structure.

Molecule Modeling of Human Pentameric alpha(7) Neuronal Nicotinic Acetylcholine Receptor and Its Interaction with its Agonist and Competitive Antagonist

The nicotinic Acetylcholine Receptor (nAChR) is the major class of neurotransmitter receptors that is involved in many neurodegenerative conditions such as schizophrenia, Alzheimer's and Parkinson's diseases. The N-terminal region or Ligand Binding Domain (LBD) of nAChR is located at pre- and post-synaptic nervous system, which mediates synaptic transmission. nAChR acts as the drug target for agonist and competitive antagonist molecules that modulate signal transmission at the nerve terminals. Based on Acetylcholine Binding Protein (AChBP) from Lymnea stagnalis as the structural template, the homology modeling approach was carried out to build three dimensional model of the N-terminal region of human alpha(7)nAChR. This theoretical model is an assembly of five alpha(7) subunits with 5 fold axis symmetry, constituting a channel, with the binding picket present at the interface region of the subunits. alpha-netlrotoxin is a potent nAChR competitive antagonist that readily blocks the channel resulting in paralysis. The molecular interaction of alpha-Bungarotoxin, a long chain alpha-neurotoxin from (Bungarus multicinctus) and human alpha(7)nAChR seas studied. Agonists such as acetylcholine, nicotine, which are used in it diverse array of biological activities, such as enhancements of cognitive performances, were also docked with the theoretical model of human alpha(7)nAChR. These docked complexes were analyzed further for identifying the crucial residues involved i interaction. These results provide the details of interaction of agonists and competitive antagonists with three dimensional model of the N-terminal region of human alpha(7)nAChR and thereby point to the design of novel lead compounds.

Allostery in tRNA Synthetases Elucidated from MD Simulations and Protein Structure Networks

tRNA synthetases (aaRS) are enzymes crucial in the translation of genetic code. The enzyme accylates the acceptor stem of tRNA by the congnate amino acid bound at the active site, when the anti-codon is recognized by the anti-codon site of aaRS. In a typical aaRS, the distance between the anti-codon region and the amino accylation site is approximately 70 Å. We have investigated this allosteric phenomenon at molecular level by MD simulations followed by the analysis of protein structure networks (PSN) of non-covalent interactions. Specifically, we have generated conformational ensembles by performing MD simulations on different liganded states of methionyl tRNA synthetase (MetRS) from Escherichia coli and tryptophenyl tRNA synthetase (TrpRS) from Human. The correlated residues during the MD simulations are identified by cross correlation maps. We have identified the amino acids connecting the correlated residues by the shortest path between the two selected members of the PSN. The frequencies of paths have been evaluated from the MD snapshots[1]. The conformational populations in different liganded states of the protein have been beautifully captured in terms of network parameters such as hubs, cliques and communities[2]. These parameters have been associated with the rigidity and plasticity of the protein conformations and can be associated with free energy landscape. A comparison of allosteric communication in MetRS and TrpRS [3] elucidated in this study highlights diverse means adopted by different enzymes to perform a similar function. The computational method described for these two enzymes can be applied to the investigation of allostery in other systems.

Structural Polymorphism of the Nucleosomal DNA and Implications for Protein Binding

A nucleosome forms a basic unit of the chromosome structure. A biologically relevant question is how much of the nucleosomal conformational space is accessible to protein-free DNA, and what proportion of the nucleosomal conformations are induced by bound histones. To investigate this, we have analysed high resolution xray crystal structure datasets of DNA in protein-free as well as protein-bound forms, and compared the dinucleotide step parameters for the two datasets with those for high resolution nucleosome structures. Our analysis shows that most of the dinucleotide step parameter values for the nucleosome structures lie within the range accessible to protein-free DNA, indirectly indicating that the histone core plays more of a stabilizing role. The nucleosome structures are observed to assume smooth and nearly planar curvature, implying that ‘normal’ B-DNA like parameters can give rise to a curved geometry at the gross structural level. Different nucleosome

Effect of configuration of the inhibitors on the mode of binding to the enzyme, thermolysin

Preferred conformations of the competitive inhibitors glycyl-L-phenylalanine and glycyl-D-phenylalanine and their mode of binding to thermolysin have been studied. The difference in configuration is shown to affect significantly the mode of binding to thermolysin. Gly-D-Phe prefers to enter the active site in the global minimum conformation whereas Gly-L-Phe may enter in a higher energy conformation. Moreover, D-enantiomer is shown to have a better fit than the L-counterpart in the active site.

The modes of binding methyl-alpha (and beta)-D-glucopyranosides and some of their derivatives to concanavalin A--a theoretical approach

The probable modes of binding of Methyl--alpha (and beta)-D-glucopyranosides and some of their derivatives to concanavalin A have been proposed from theoretical studies. Theory predicts that beta-MeGlcP can bind to ConA in three different modes whereas alpha-MeGlcP can bind only in one mode. beta-MeGlcP in its most favourable mode of binding differs from alpha-MeGlcP in its alignment in the active-site of the lectin where it binds in a flipped or inverted orientation. Methyl substitution at the C-2 atom of the alpha-MeGlcP does not significantly affect the possible orientations of the sugar in the active-site of the lectin. Methyl substitution at C-3 or C-4, however, affects the allowed orientations drastically leading to the poor inhibiting power of Methyl-3-O-methyl-alpha-D-glucopyranoside and the inactivity of Methyl-4-O-methyl-alpha-D-glycopyranoside. These studies suggest that the increased activity of the alpha-MeGlcP over beta-MeGlcP may be due to the possibility of formation of better hydrogen bonds and to hydrophobic interactions rather than to steric factors as suggested by earlier workers. These models explain the available NMR and other binding studies.

A detailed study of Li-DNA fibres at various salt concentrations reveals a non-helical B-DNA and a possible similarity of solution and solid state structures

A thorough investigation of salt concentration dependence of lithium DNA fibres is made using X-ray diffraction. While for low salt the C-form pattern is obtained, crystalline B-type diffraction patterns result on increasing the salt concentration. The salt content in the gel (from which fibres are drawn) is estimated by equilibrium dialysis using the Donnan equilibrium principle. The salt range giving the best crystalline B pattern is determined. It is found that in this range meridional reflections occur on the fourth and sixth layer lines. In addition, the tenth layer meridian is absent at a particular salt concentration. These results strongly suggest the presence of non-helical features in the DNA molecule. Preliminary analysis of the diffraction patterns indicates a structural variability within the B-form itself. Further, the possibility of the structural parameters of DNA being similar in solid state and in solution is discussed.

CSSI-PRO: a method for secondary structure type editing, assignment and estimation in proteins using linear combination of backbone chemical shifts

Estimation of secondary structure in polypeptides is important for studying their structure, folding and dynamics. In NMR spectroscopy, such information is generally obtained after sequence specific resonance assignments are completed. We present here a new methodology for assignment of secondary structure type to spin systems in proteins directly from NMR spectra, without prior knowledge of resonance assignments. The methodology, named Combination of Shifts for Secondary Structure Identification in Proteins (CSSI-PRO), involves detection of specific linear combination of backbone H-1(alpha) and C-13' chemical shifts in a two-dimensional (2D) NMR experiment based on G-matrix Fourier transform (GFT) NMR spectroscopy. Such linear combinations of shifts facilitate editing of residues belonging to alpha-helical/beta-strand regions into distinct spectral regions nearly independent of the amino acid type, thereby allowing the estimation of overall secondary structure content of the protein. Comparison of the predicted secondary structure content with those estimated based on their respective 3D structures and/or the method of Chemical Shift Index for 237 proteins gives a correlation of more than 90% and an overall rmsd of 7.0%, which is comparable to other biophysical techniques used for structural characterization of proteins. Taken together, this methodology has a wide range of applications in NMR spectroscopy such as rapid protein structure determination, monitoring conformational changes in protein-folding/ligand-binding studies and automated resonance assignment.

Infrared and Raman spectra of 2,4-dithiobiuret and its normal vibrations

The infrared spectra of 2,4-dithiobiuret(DTB), N-deuterated dithiobiuret(DTB-d5) and the laser Raman spectrum of DTB are reported. Normal coordinate treatments of DTB and DTB-d5 have been carried out to aid the assignment of the vibrational frequencies. A trans—cis conformation is favoured for DTB molecule in the solid state.