Path integral description of polymers using fractional Brownian walks

The statistical properties of fractional Brownian walks are used to construct a path integral representation of the conformations of polymers with different degrees of bond correlation. We specifically derive an expression for the distribution function of the chains’ end‐to‐end distance, and evaluate it by several independent methods, including direct evaluation of the discrete limit of the path integral, decomposition into normal modes, and solution of a partial differential equation. The distribution function is found to be Gaussian in the spatial coordinates of the monomer positions, as in the random walk description of the chain, but the contour variables, which specify the location of the monomer along the chain backbone, now depend on an index h, the degree of correlation of the fractional Brownian walk. The special case of h=1/2 corresponds to the random walk. In constructing the normal mode picture of the chain, we conjecture the existence of a theorem regarding the zeros of the Bessel function.

Structure of Sesbania Mosaic Virus at 4·7 Å Resolution and Partial Sequence of the Coat Protein

Sesbania mosaic virus (SMV) is a plant virus infecting Sesbania grandiflora plants in Andhra Pradesh, India. Amino acid sequence of the tryptic peptides of SMV coat protein were determined using a gas phase sequenator. These sequences showed identical amino acids at 69% of the positions when aligned with the corresponding residues of southern bean mosaic virus (SBMV).Crystals diffracting to better than 3 Å resolution were obtained by precipitating the virus with ammonium sulphate. The crystals belonged to rhombohedral space group R3 with α = 291·4 Å and α = 61·9°. Three-dimensional X-ray diffraction data on these crystals were collected to a resolution of 4·7 Å, using a Siemens-Nicolet area detector system. Self-rotation function studies revealed the icosahedral symmetry of the virus particles, as well as their precise orientation in the unit cell. Cross-rotation function and modelling studies with SBMV showed that it is a valid starting model for SMV structure determination. Low resolution phases computed using a polyalanine model of SBMV were subjected to refinement and extension by real-space electron density averaging and solvent flattening. The final electron density map revealed a polypeptide fold similar to SBMV. The single disulphide bridge of SBMV coat protein is retained in SMV. Four icosahedrally independent cation binding sites have been tentatively identified. Three of these sites, related by a quasi threefold axis, are also found in SBMV. The fourth site is situated on the quasi threefold axis. Aspartic acid residues, which replace Ile218 of SBMV from the quasi threefold-related subunits are suitable ligands to the cation at this site

The variants of the protein toxins abrin and ricin A useful guide to understanding the processing events in the toxin transport

Kinetic data on inhibition of protein synthesis in thymocyte by three abrins and ricin have been obtained. The intrinsic efficiencies of A chains of four toxins to inactivate ribosomes, as analyzed by k1-versus-concentration plots were abrin II, III > ricin > abrin I. The lag times were 90, 66, 75 and 105 min at a 0.0744 nM concentration of each of abrin I, II, III and ricin, respectively. To account for the observed differences in the dose-dependent lag time, functional and structural variables of toxins such as binding efficiency of B chains to receptors and low-pH-induced structural alterations have been analyzed. The association constants obtained by stopped flow studies showed that abrin-I (4.13 × 105 M−1 s−1) association with putative receptor (4-methylumbelliferyl-α-D-galactoside) is nearly two times more often than abrin III (2.6 × 105 M−1 s−1) at 20°C. Equillibrium binding constants of abrin I and II to thymocyte at 37°C were 2.26 × 107 M−1 and 2.8 × 107 M−1 respectively. pH-induced structural alterations as studied by a parallel enhancement in 8-anilino-L-naphthalene sulfonate fluorescence revealed a high degree of qualitative similarity. These results taken with a nearly identical concentration-independent lag time (minimum lag of 41–42 min) indicated that the binding efficiencies and internalization efficiencies of these toxins are the same and that the observed difference in the dose-dependent lag time is causally related to the proposed processing event. The rates of reduction of inter-subunit disulfide bond, an obligatory step in the intoxication process, have been measured and compared under a variety of conditions. Intersubunit disulfide reduction of abrin I is fourfold faster than that of abrin II at pH 7.2. The rate of disulfide reduction in abrin I could be decreased 1 I-fold by adding lactose, compared to that without lactose. The observed differences in the efficiencies of A chains, the dose-dependent lag period, the modulating effect of lactose on the rates of disulfide reduction and similarity in binding properties make the variants a valuable tool to probe the processing events in toxin transport in detail.

Efficiency of nonhomologous DNA end joining varies among somatic tissues, despite similarity in mechanism

Failure to repair DNA double-strand breaks (DSBs) can lead to cell death or cancer. Although nonhomologous end joining (NHEJ) has been studied extensively in mammals, little is known about it in primary tissues. Using oligomeric DNA mimicking endogenous DSBs, NHEJ in cell-free extracts of rat tissues were studied. Results show that efficiency of NHEJ is highest in lungs compared to other somatic tissues. DSBs with compatible and blunt ends joined without modifications, while noncompatible ends joined with minimal alterations in lungs and testes. Thymus exhibited elevated joining, followed by brain and spleen, which could be correlated with NHEJ gene expression. However, NHEJ efficiency was poor in terminally differentiated organs like heart, kidney and liver. Strikingly, NHEJ junctions from these tissues also showed extensive deletions and insertions. Hence, for the first time, we show that despite mode of joining being generally comparable, efficiency of NHEJ varies among primary tissues of mammals.

Effect of Substrate Binding Loop Mutations on the Structure, Kinetics, and Inhibition of Enoyl Acyl Carrier Protein Reductase from Plasmodium falciparum

Enoyl acyl carrier protein reductase (ENR), which catalyzes the final and rate limiting step of fatty acid elongation, has been validated as a potential drug target. Triclosan is known to be an effective inhibitor for this enzyme. We mutated the substrate binding site residue Ala372 of the ENR of Plasmodium falciparum (PfENR) to Methionine and Valine which increased the affinity of the enzyme towards triclosan to almost double, close to that of Escherichia coli ENR (EcENR) which has a Methionine at the structurally similar position of Ala372 of PfENR. Kinetic studies of the mutants of PfENR and the crystal structure analysis of the A372M mutant revealed that a more hydrophobic environment enhances the affinity of the enzyme for the inhibitor. A triclosan derivative showed a threefold increase in the affinity towards the mutants compared to the wild type, due to additional interactions with the A372M mutant as revealed by the crystal structure. The enzyme has a conserved salt bridge which stabilizes the substrate binding loop and appears to be important for the active conformation of the enzyme. We generated a second set of mutants to check this hypothesis. These mutants showed loss of function, except in one case, where the crystal structure showed that the substrate binding loop is stabilized by a water bridge network. (C) 2011 IUBMB mum Life, 63(1): 30-41,2011

One protein�many functions

The concept of one enzyme-one activity had influenced biochemistry for over half a century. Over 1000 enzymes are now described. Many of them are highly 'specific'. Some of them are crystallized and their three-dimensional structures determined. They range from 12 to 1000 kDa in molecular weight and possess 124 to several hundreds of amino acids. They occur as single polypeptides or multiple-subunit proteins. The active sites are assembled on these by appropriate tertiary folding of the polypeptide chain, or by interaction of the constituent subunits. The substrate is held by the side-chains of a few amino acids at the active site on the surface, occupying a tiny fraction of the total area. What is the bulk of the protein behind the active site doing? Do all proteins have only one function each? Why not a protein have more than one active site on its large surface? Will we discover more than one activity for some proteins? These newer possibilities are emerging and are finding experimental support. Some proteins purified to homogeneity using assay methods for different activities are now recognized to have the same molecular weight and a high degree of homology of amino acid sequence. Obviously they are identical. They represent the phenomenon of one protein-many functions.

Primary structure of sesbania mosaic virus coat protein: its implications to the assembly and architecture of the virus

Sesbania mosaic virus (SMV) is a plant virus that infects Sesbania grandiflora plants in Andhra Pradesh, India. The amino acid sequence of the coat protein of SMV was determined using purified peptides generated by cleavage with trypsin, chymotrypsin, V8 protease and clostripain. The 230 residues so far determined were compared to the corresponding residues of southern bean mosaic virus (SBMV), the type member of sobemoviruses. The overall identity between the sequences is 61.7%. The amino terminal 64 residues, which constitute an independent domain (R-domain) known to interact with RNA, are conserved to a lower extent (52.5%). Comparison of the positively charged residues in this domain suggests that the RNA-protein interactions are considerably weaker in SMV. The residues that constitute the major domain of the coat protein, the surface domain (S-domain, residues 65-260), are better conserved (66.5%). The positively charged residues of this domain that face the nucleic acid are well conserved. The longest conserved stretch of residues (131-142) corresponds to the loop involved in intersubunit interactions between subunits related by the quasi 3-fold symmetry. A unique cation binding site located on the quasi 3-fold axis contributes to the stability of SMV. These differences are reflected in the increased stability of the SMV coat protein and its ability to be reconstituted with RNA at pH 7.5. A major epitope was identified using monoclonal antibodies to SMV in the segment 201-223 which contains an exposed helix in the capsid structure. This region is highly conserved between SMV and SBMV (70%) suggesting that it could represent the site of an important function such as vector recognition.

Using Correlated Parameters for Improved Ranking of Protein-Protein Docking Decoys

A successful protein-protein docking study culminates in identification of decoys at top ranks with near-native quaternary structures. However, this task remains enigmatic because no generalized scoring functions exist that effectively infer decoys according to the similarity to near-native quaternary structures. Difficulties arise because of the highly irregular nature of the protein surface and the significant variation of the nonbonding and solvation energies based on the chemical composition of the protein-protein interface. In this work, we describe a novel method combining an interface-size filter, a regression model for geometric compatibility (based on two correlated surface and packing parameters), and normalized interaction energy (calculated from correlated nonbonded and solvation energies), to effectively rank decoys from a set of 10,000 decoys. Tests on 30 unbound binary protein-protein complexes show that in 16 cases we can identify at least one decoy in top three ranks having <= 10 angstrom backbone root mean square deviation from true binding geometry. Comparisons with other state-of-art methods confirm the improved ranking power of our method without the use of any experiment-guided restraints, evolutionary information, statistical propensities, or modified interaction energy equations. Tests on 118 less-difficult bound binary protein-protein complexes with <= 35% sequence redundancy at the interface showed that in 77% cases, at least 1 in 10,000 decoys were identified with <= 5 angstrom backbone root mean square deviation from true geometry at first rank. The work will promote the use of new concepts where correlations among parameters provide more robust scoring models. It will facilitate studies involving molecular interactions, including modeling of large macromolecular assemblies and protein structure prediction. (C) 2010 Wiley Periodicals, Inc. J Comput Chem 32: 787-796, 2011.

Characterization of the DNA-binding domain of beta protein, a component of phage lambda Red-pathway by UV catalyzed cross-linking

beta protein, a key component of Red-pathway of phage lambda is necessary for its growth and general genetic recombination in recombination-deficient mutants of Escherichia coli. To facilitate studies on structure-function relationships, we overexpressed beta protein and purified it to homogeneity. A chemical cross-linking reagent, glutaraldehyde, was used to stabilize the physical association of beta protein in solution. A 67-kDa band, corresponding to homodimer, was identified after separation by SDS-polyacrylamide gel electrophoresis. Stoichiometric measurements indicated a site-size of 1 monomer of beta protein/5 nucleotide residues. Electrophoretic gel mobility shift assays suggested that beta protein formed stable nucleoprotein complexes with 36-mer, but not with 27- or 17-mer DNA. Interestingly, the interaction of beta protein with DNA and the stability of nucleoprotein complexes was dependent on the presence of MgCl2, and the binding was abolished by 250 mM NaCl. The K-d of beta protein binding to 36-mer DNA was on the order of 1.8 x 10(-6) M. Photochemical cross-linking of native beta protein or its fragments, generated by chymotrypsin, to 36-mer DNA was performed to identify its DNA-binding domain. Characterization of the cross-linked peptide disclosed that amino acids required for DNA-binding specificity resided within a 20-kDa peptide at the N-terminal end. These findings provide a basis for further understanding oi the structure and function of beta protein.

Engineered dimer interface mutants of triosephosphate isomerase: the role of inter-subunit interactions in enzyme function and stability

The role of inter-subunit interactions in maintaining optimal catalytic activity in triosephosphate isomerase (TIM) has been probed, using the Plasmodium falciparum enzyme as a model. Examination of subunit interface contacts in the crystal structures suggests that residue 75 (Thr, conserved) and residue 13 (Cys, variable) make the largest number of inter-subunit contacts. The mutants Cys13Asp (C13D) and Cys13Glu (C13E) have been constructed and display significant reduction in catalytic activity when compared with wild-type (WT) enzyme (similar to 7.4-fold decrease in k(cat) for the C13D and similar to 3.3-fold for the C13E mutants). Analytical gel filtration demonstrates that the C13D mutant dissociates at concentrations < 1.25 mu M, whereas the WT and the C13E enzymes retain the dimeric structure. The order of stability of the mutants in the presence of chemical denaturants, like urea and guanidium chloride, is WT > Cys13Glu > Cys13Asp. Irreversible thermal precipitation temperatures follow the same order as well. Modeling studies establish that the Cys13Asp mutation is likely to cause a significantly greater structural perturbation than Cys13Glu. Analysis of sequence and structural data for TIMs from diverse sources suggests that residues 13 and 82 form a pair of proximal sites, in which a limited number of residue pairs may be accommodated.

Prediction of protein-protein interactions between human host and a pathogen and its application to three pathogenic bacteria

Molecular understanding of disease processes can be accelerated if all interactions between the host and pathogen are known. The unavailability of experimental methods for large-scale detection of interactions across host and pathogen organisms hinders this process. Here we apply a simple method to predict protein-protein interactions across a host and pathogen organisms. We use homology detection approaches against the protein-protein interaction databases. DIP and iPfam in order to predict interacting proteins in a host-pathogen pair. In the present work, we first applied this approach to the test cases involving the pairs phage T4 - Escherichia coli and phage lambda - E. coli and show that previously known interactions could be recognized using our approach. We further apply this approach to predict interactions between human and three pathogens E. coli, Salmonella enterica typhimurium and Yersinia pestis. We identified several novel interactions involving proteins of host or pathogen that could be thought of as highly relevant to the disease process. Serendipitously, many interactions involve hypothetical proteins of yet unknown function. Hypothetical proteins are predicted from computational analysis of genome sequences with no laboratory analysis on their functions yet available. The predicted interactions involving such proteins could provide hints to their functions. (C) 2011 Elsevier B.V. All rights reserved.

Angelman Syndrome Protein UBE3A Interacts with Primary Microcephaly Protein ASPM, Localizes to Centrosomes and Regulates Chromosome Segregation

Many proteins associated with the phenotype microcephaly have been localized to the centrosome or linked to it functionally. All the seven autosomal recessive primary microcephaly (MCPH) proteins localize at the centrosome. Microcephalic osteodysplastic primordial dwarfism type II protein PCNT and Seckel syndrome (also characterized by severe microcephaly) protein ATR are also centrosomal proteins. All of the above findings show the importance of centrosomal proteins as the key players in neurogenesis and brain development. However, the exact mechanism as to how the loss-of-function of these proteins leads to microcephaly remains to be elucidated. To gain insight into the function of the most commonly mutated MCPH gene ASPM, we used the yeast two-hybrid technique to screen a human fetal brain cDNA library with an ASPM bait. The analysis identified Angelman syndrome gene product UBE3A as an ASPM interactor. Like ASPM, UBE3A also localizes to the centrosome. The identification of UBE3A as an ASPM interactor is not surprising as more than 80% of Angelman syndrome patients have microcephaly. However, unlike in MCPH, microcephaly is postnatal in Angelman syndrome patients. Our results show that UBE3A is a cell cycle regulated protein and its level peaks in mitosis. The shRNA knockdown of UBE3A in HEK293 cells led to many mitotic abnormalities including chromosome missegregation, abnormal cytokinesis and apoptosis. Thus our study links Angelman syndrome protein UBE3A to ASPM, centrosome and mitosis for the first time. We suggest that a defective chromosome segregation mechanism is responsible for the development of microcephaly in Angelman syndrome.

Combinatorial logic synthesis using technology directed decomposition

The aim of logic synthesis is to produce circuits which satisfy the given boolean function while meeting timing constraints and requiring the minimum silicon area. Logic synthesis involves two steps namely logic decomposition and technology mapping. Existing methods treat the two as separate operation. The traditional approach is to minimize the number of literals without considering the target technology during the decomposition phase. The decomposed expressions are then mapped on to the target technology to optimize the area, Timing optimization is carried out subsequently, A new approach which treats logic decomposition and technology maping as a single operation is presented. The logic decomposition is based on the parameters of the target technology. The area and timing optimization is carried out during logic decomposition phase itself. Results using MCNC circuits are presented to show that this method produces circuits which are 38% faster while requiring 14% increase in area.

Water-mediated ionic interactions in protein structures

It is well known that water molecules play an indispensable role in the structure and function of biological macromolecules. The water-mediated ionic interactions between the charged residues provide stability and plasticity and in turn address the function of the protein structures. Thus, this study specifically addresses the number of possible water-mediated ionic interactions, their occurrence, distribution and nature found in 90% non-redundant protein chains. Further, it provides a statistical report of different charged residue pairs that are mediated by surface or buried water molecules to form the interactions. Also, it discusses its contributions in stabilizing various secondary structural elements of the protein. Thus, the present study shows the ubiquitous nature of the interactions that imparts plasticity and flexibility to a protein molecule.

Absolute Entropy and Energy of Carbon Dioxide Using the Two-Phase Thermodynamic Model

The two-phase thermodynamic (2PT) model is used to determine the absolute entropy and energy of carbon dioxide over a wide range of conditions from molecular dynamics trajectories. The 2PT method determines the thermodynamic properties by applying the proper statistical mechanical partition function to the normal modes of a fluid. The vibrational density of state (DoS), obtained from the Fourier transform of the velocity autocorrelation function, converges quickly, allowing the free energy, entropy, and other thermodynamic properties to be determined from short 20-ps MD trajectories. The anharmonic effects in the vibrations are accounted for by the broadening of the normal modes into bands from sampling the velocities over the trajectory. The low frequency diffusive modes, which lead to finite DoS at zero frequency, are accounted for by considering the DoS as a superposition of gas-phase and solid-phase components (two phases). The analytical decomposition of the DoS allows for an evaluation of properties contributed by different types of molecular motions. We show that this 2PT analysis leads to accurate predictions of entropy and energy of CO2 over a wide range of conditions (from the triple point to the critical point of both the vapor and the liquid phases along the saturation line). This allows the equation of state of CO2 to be determined, which is limited only by the accuracy of the force field. We also validated that the 2PT entropy agrees with that determined from thermodynamic integration, but 2PT requires only a fraction of the time. A complication for CO2 is that its equilibrium configuration is linear, which would have only two rotational modes, but during the dynamics it is never exactly linear, so that there is a third mode from rotational about the axis. In this work, we show how to treat such linear molecules in the 2PT framework.