Solvent effects on the energy landscapes and folding kinetics of polyalanine

The effect of a solvation on the thermodynamics and kinetics of polyalanine (Ala12) is explored on the basis of its energy landscapes in vacuum and in an aqueous solution. Both energy landscapes are characterized by two basins, one associated with α-helical structures and the other with coil and β-structures of the peptide. In both environments, the basin that corresponds to the α-helical structure is considerably narrower than the basin corresponding to the β-state, reflecting their different contributions to the entropy of the peptide. In vacuum, the α-helical state of Ala12 constitutes the native state, in agreement with common helical propensity scales, whereas in the aqueous medium, the α-helical state is destabilized, and the β-state becomes the native state. Thus solvation has a dramatic effect on the energy landscape of this peptide, resulting in an inverted stability of the two states. Different folding and unfolding time scales for Ala12 in hydrophilic and hydrophobic chemical environments are caused by the higher entropy of the native state in water relative to vacuum. The concept of a helical propensity has to be extended to incorporate environmental solvent effects.

Recent improvements in prediction of protein structure by global optimization of a potential energy function

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsd's are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.

Point mutations and sequence variability in proteins: Redistributions of preexisting populations

Here we study the effect of point mutations in proteins on the redistributions of the conformational substates. We show that regardless of the location of a mutation in the protein structure and of its type, the observed movements of the backbone recur largely at the same positions in the structures. Despite the different interactions that are disrupted and formed by the residue substitution, not only are the conformations very similar, but the regions that move are also the same, regardless of their sequential or spatial distance from the mutation. This observation leads us to conclude that, apart from some extreme cases, the details of the interactions are not critically important in determining the protein conformation or in specifying which parts of the protein would be more prone to take on different local conformations in response to changes in the sequence. This finding further illustrates why proteins manifest a robustness toward many mutational events. This nonuniform distribution of the conformer population is consistently observed in a variety of protein structural types. Topology is critically important in determining folding pathways, kinetics, building block cutting, and anatomy trees. Here we show that topology is also very important in determining which regions of the protein structure will respond to sequence changes, regardless of the sequential or spatial location of the mutation.

Endogenous Msx1 antisense transcript: In vivo and in vitro evidences, structure, and potential involvement in skeleton development in mammals

Msx1 is a key factor for the development of tooth and craniofacial skeleton and has been proposed to play a pivotal role in terminal cell differentiation. In this paper, we demonstrated the presence of an endogenous Msx1 antisense RNA (Msx1-AS RNA) in mice, rats, and humans. In situ analysis revealed that this RNA is expressed only in differentiated dental and bone cells with an inverse correlation with Msx1 protein. These in vivo data and overexpression of Msx1 sense and AS RNA in an odontoblastic cell line (MO6-G3) showed that the balance between the levels of the two Msx1 RNAs is related to the expression of Msx1 protein. To analyze the impact of this balance in the Msx-Dlx homeoprotein pathway, we analyzed the effect of Msx1, Msx2, and Dlx5 overexpression on proteins involved in skeletal differentiation. We showed that the Msx1-AS RNA is involved in crosstalk between the Msx-Dlx pathways because its expression was abolished by Dlx5. Msx1 was shown to down-regulate a master gene of skeletal cells differentiation, Cbfa1. All these data strongly suggest that the ratio between Msx1 sense and antisense RNAs is a very important factor in the control of skeletal terminal differentiation. Finally, the initiation site for Msx1-AS RNA transcription was located by primer extension in both mouse and human in an identical region, including a consensus TATA box, suggesting an evolutionary conservation of the AS RNA-mediated regulation of Msx1 gene expression.

Structure and Expression of a Dhurrinase (β-Glucosidase) from Sorghum1

Sorghum (Sorghum bicolor L. Moench) has two isozymes of the cyanogenic β-glucosidase dhurrinase: dhurrinase-1 (Dhr1) and dhurrinase-2 (Dhr2). A nearly full-length cDNA encoding dhurrinase was isolated from 4-d-old etiolated seedlings and sequenced. The cDNA has a 1695-nucleotide-long open reading frame, which codes for a 565-amino acid-long precursor and a 514-amino acid-long mature protein, respectively. Deduced amino acid sequence of the sorghum Dhr showed 70% identity with two maize (Zea mays) β-glucosidase isozymes. Southern-blot data suggested that β-glu-cosidase is encoded by a small multigene family in sorghum. Northern-blot data indicated that the mRNA corresponding to the cloned Dhr cDNA is present at high levels in the node and upper half of the mesocotyl in etiolated seedlings but at low levels in the root—only in the zone of elongation and the tip region. Light-grown seedling parts had lower levels of Dhr mRNA than those of etiolated seedlings. Immunoblot analysis performed using maize-anti-β-glucosidase sera detected two distinct dhurrinases (57 and 62 kD) in sorghum. The distribution of Dhr activity in different plant parts supports the mRNA and immunoreactive protein data, suggesting that the cloned cDNA corresponds to the Dhr1 (57 kD) isozyme and that the dhr1 gene shows organ-specific expression.

Fourier Transform Infrared Microspectroscopy Detects Changes in Protein Secondary Structure Associated with Desiccation Tolerance in Developing Maize Embryos1

Isolated immature maize (Zea mays L.) embryos have been shown to acquire tolerance to rapid drying between 22 and 25 d after pollination (DAP) and to slow drying from 18 DAP onward. To investigate adaptations in protein profile in association with the acquisition of desiccation tolerance in isolated, immature maize embryos, we applied in situ Fourier transform infrared microspectroscopy. In fresh, viable, 20- and 25-DAP embryo axes, the shapes of the different amide-I bands were identical, and this was maintained after flash drying. On rapid drying, the 20-DAP axes had a reduced relative proportion of α-helical protein structure and lost viability. Rapidly dried 25-DAP embryos germinated (74%) and had a protein profile similar to the fresh control axes. On slow drying, the α-helical contribution in both the 20- and 25-DAP embryo axes increased compared with that in the fresh control axes, and survival of desiccation was high. The protein profile in dry, mature axes resembled that after slow drying of the immature axes. Rapid drying resulted in an almost complete loss of membrane integrity in the 20-DAP embryo axes and much less so in the 25-DAP axes. After slow drying, low plasma membrane permeability ensued in both the 20- and 25-DAP axes. We conclude that slow drying of excised, immature embryos leads to an increased proportion of α-helical protein structures in their axes, which coincides with additional tolerance of desiccation stress.

Domain structure and dynamics in the helical filaments formed by RecA and Rad51 on DNA

Both the bacterial RecA protein and the eukaryotic Rad51 protein form helical nucleoprotein filaments on DNA that catalyze strand transfer between two homologous DNA molecules. However, only the ATP-binding cores of these proteins have been conserved, and this same core is also found within helicases and the F1-ATPase. The C-terminal domain of the RecA protein forms lobes within the helical RecA filament. However, the Rad51 proteins do not have the C-terminal domain found in RecA, but have an N-terminal extension that is absent in the RecA protein. Both the RecA C-terminal domain and the Rad51 N-terminal domain bind DNA. We have used electron microscopy to show that the lobes of the yeast and human Rad51 filaments appear to be formed by N-terminal domains. These lobes are conformationally flexible in both RecA and Rad51. Within RecA filaments, the change between the “active” and “inactive” states appears to mainly involve a large movement of the C-terminal lobe. The N-terminal domain of Rad51 and the C-terminal domain of RecA may have arisen from convergent evolution to play similar roles in the filaments.

Structure and inhibition of plasmepsin II, a hemoglobin-degrading enzyme from Plasmodium falciparum.

Plasmodium falciparum is the major causative agent of malaria, a disease of worldwide importance. Resistance to current drugs such as chloroquine and mefloquine is spreading at an alarming rate, and our antimalarial armamentarium is almost depleted. The malarial parasite encodes two homologous aspartic proteases, plasmepsins I and II, which are essential components of its hemoglobin-degradation pathway and are novel targets for antimalarial drug development. We have determined the crystal structure of recombinant plasmepsin II complexed with pepstatin A. This represents the first reported crystal structure of a protein from P. falciparum. The crystals contain molecules in two different conformations, revealing a remarkable degree of interdomain flexibility of the enzyme. The structure was used to design a series of selective low molecular weight compounds that inhibit both plasmepsin II and the growth of P. falciparum in culture.

Primary structure and tissue distribution of the orphanin FQ precursor.

The heptadecapeptide orphanin FQ (OFQ) is a recently discovered neuropeptide that exhibits structural features reminiscent of the opioid peptides and that is an endogenous ligand to a G protein-coupled receptor sequentially related to the opioid receptors. We have cloned both the human and rat cDNAs encoding the OFQ precursor proteins, to investigate whether the sequence relationships existing between the opioid and OFQ systems are also found at the polypeptide precursor level, in particular whether the OFQ precursor would encode several bioactive peptides as do the opioid precursors, and to study the regional distribution of OFQ sites of synthesis. The entire precursor protein displays structural homology to the opioid peptide precursors, especially preprodynorphin and preproenkephalin. The predicted amino acid sequence of the OFQ precursor contains a putative signal peptide and one copy of the OFQ sequence flanked by pairs of basic amino acid residues. Carboxyl-terminal to the OFQ sequence, the human and rat precursors contain a stretch of 28 amino acids that is 100% conserved and thus may encode novel bioactive peptides. Two peptides derived from this stretch were synthesized but were found to be unable to activate the OFQ receptor, suggesting that if they are produced in vivo, these peptides would likely recognize receptors different from the OFQ receptor. To begin analyzing the sites of OFQ mRNA synthesis, Northern analysis of human and rat tissues were carried out and showed that the OFQ precursor mRNA is mainly expressed in the brain. In situ hybridization of rat brain slices demonstrated a regional distribution pattern of the OFQ precursor mRNA, which is distinct from that of the opioid peptide precursors. These data confirm that the OFQ system differs from the opioid system at the molecular level, although the OFQ and opioid precursors may have arisen from a common ancestral gene.

Primary structure and functional expression of a cGMP-gated potassium channel.

Cyclic nucleotides modulate potassium (K) channel activity in many cells and are thought to act indirectly by inducing channel protein phosphorylation. Herein we report the isolation from rabbit of a gene encoding a K channel (Kcn1) that is specifically activated by cGMP and not by cAMP. Analysis of the deduced amino acid sequence (725 amino acids) indicates that, in addition to a core region that is highly homologous to Shaker K channels, Kcn1 also contains a cysteine-rich region similar to that of ligand-gated ion channels and a cyclic nucleotide-binding region. Northern blot analysis detects gene expression in kidney, aorta, and brain. Kcn1 represents a class of K channels that may be specifically regulated by cGMP and could play an important role in mediating the effects of substances, such as nitric oxide, that increase intracellular cGMP.

Protein structure-based design of potent orally bioavailable, nonpeptide inhibitors of human immunodeficiency virus protease.

A class of potent nonpeptidic inhibitors of human immunodeficiency virus protease has been designed by using the three-dimensional structure of the enzyme as a guide. By employing iterative protein cocrystal structure analysis, design, and synthesis the binding affinity of the lead compound was incrementally improved by over four orders of magnitude. An inversion in inhibitor binding mode was observed crystallographically, providing information critical for subsequent design and highlighting the utility of structural feedback in inhibitor optimization. These inhibitors are selective for the viral protease enzyme, possess good antiviral activity, and are orally available in three species.

New multivalent calixarene-based ligands for protein recognition and inhibition

One of the challenges that concerns chemistry is the design of molecules able to modulate protein-protein and protein-ligand interactions, since these are involved in many physiological and pathological processes. The interactions occurring between proteins and their natural counterparts can take place through reciprocal recognition of rather large surface areas, through recognition of single contact points and single residues, through inclusion of the substrates in specific, more or less deep binding sites. In many cases, the design of synthetic molecules able to interfere with the processes involving proteins can benefit from the possibility of exploiting the multivalent effect. Multivalency, widely spread in Nature, consists in the simultaneous formation between two entities (cell-cell, cell-protein, protein-protein) of multiple equivalent ligand-recognition site complexes. In this way the whole interaction results particularly strong and specific. Calixarenes furnish a very interesting scaffold for the preparation of multivalent ligands and in the last years calixarene-based ligands demonstrated their remarkable capability to recognize and inhibit or restore the activity of different proteins, with a high efficiency and selectivity in several recognition phenomena. The relevance and versatility of these ligands is due to the different exposition geometries of the binding units that can be explored exploiting the conformational properties of these macrocycles, the wide variety of functionalities that can be linked to their structure at different distances from the aromatic units and to their intrinsic multivalent nature. With the aim of creating new multivalent systems for protein targeting, the work reported in this thesis regards the synthesis and properties of glycocalix[n]arenes and guanidino calix[4]arenes for different purposes. Firstly, a new bolaamphiphile glycocalix[4]arene in 1,3-alternate geometry, bearing cellobiose, was synthesized for the preparation of targeted drug delivery systems based on liposomes. The formed stable mixed liposomes obtained by mixing the macrocycle with DOPC were shown to be able of exploiting the sugar units emerging from the lipid bilayer to agglutinate Concanavalin A, a lectin specific for glucose. Moreover, always thanks to the presence of the glycocalixarene in the layer, the same liposomes demonstrated through preliminary experiments to be uptaken by cancer cells overexpressing glucose receptors on their exterior surface more efficiently respect to simple DOPC liposomes lacking glucose units in their structure. Then a small library of glycocalix[n]arenes having different valency and geometry was prepared, for the creation of potentially active immunostimulants against Streptococcus pneumoniae, particularly the 19F serotype, one of the most virulent. These synthesized glycocalixarenes bearing β-N-acetylmannosamine as antigenic unit were compared with the natural polysaccharide on the binding to the specific anti-19F human polyclonal antibody, to verify their inhibition potency. Among all, the glycocalixarene based on the conformationally mobile calix[4]arene resulted the more efficient ligand, probably due its major possibility to explore the antibody surface and dispose the antigenic units in a proper arrangement for the interaction process. These results pointed out the importance of how the different multivalent presentation in space of the glycosyl units can influence the recognition phenomena. At last, NMR studies, using particularly 1H-15N HSQC experiments, were performed on selected glycocalix[6]arenes and guanidino calix[4]arenes blocked in the cone geometry, in order to better understand protein-ligand interactions. The glycosylated compounds were studied with Ralstonia solanacearum lectin, in order to better understand the nature of the carbohydrate‐lectin interactions in solution. The series of cationic calixarene was employed with three different acidic proteins: GB1, Fld and alpha synuclein. Particularly GB1 and Fld were observed to interact with all five cationic calix[4]arenes but showing different behaviours and affinities.

Structure and mechanical properties of sodium and calcium caseinate edible active films with carvacrol

Edible active films based on sodium caseinate (SC) and calcium caseinate (CC) plasticized with glycerol (G) at three different concentrations and carvacrol (CRV) as active agent were prepared by solvent casting. Transparent films were obtained and their surfaces were analysed by optical microscopy and scanning electron microscopy (SEM). The influence of the addition of three different plasticizer concentrations was studied by determining tensile properties, while Fourier transformed infrared spectroscopy (FTIR) and thermogravimetric analysis (TGA) were used to evaluate the structural and thermal behavior of such films. The addition of glycerol resulted in a reduction in the elastic modulus and tensile strength, while some increase in the elongation at break was observed. In general terms, SC films showed flexibility higher than the corresponding CC counterparts. In addition, the presence of carvacrol caused further improvements in ductile properties suggesting the presence of stronger interactions between the protein matrix and glycerol, as it was also observed in thermal degradation studies. FTIR spectra of all films showed the characteristic bands and peaks corresponding to proteins as well as to primary and secondary alcohols. In summary, the best results regarding mechanical and structural properties for caseinates-based films containing carvacrol were found for the formulations with high glycerol concentrations.

On the initiation of boudinage and folding instabilities in layered elasto-visco-plastic solids - Insights from natural microstructures and numerical simulations

The present understanding of the initiation of boudinage and folding structures is based on viscosity contrasts and stress exponents, considering an intrinsically unstable state of the layer. The criterion of localization is believed to be prescribed by geometry-material interactions, which are often encountered in natural structures. An alternative localization phenomenon has been established for ductile materials, in which instability emerges for critical material parameters and loading rates from homogeneous conditions. In this thesis, conditions are sought under which this type of instability prevails and whether localization in geological materials necessarily requires a trigger by geometric imperfections. The relevance of critical deformation conditions, material parameters and the spatial configuration of instabilities are discussed in a geological context. In order to analyze boudinage geometries, a numerical eigenmode analysis is introduced. This method allows determining natural frequencies and wavelengths of a structure and inducing perturbations on these frequencies. In the subsequent coupled thermo-mechanical simulations, using a grain size evolution and end-member flow laws, localization emerges when material softening through grain size sensitive viscous creep sets in. Pinch-and-swell structures evolve along slip lines through a positive feedback between the matrix response and material bifurcations inside the layer, independent from the mesh-discretization length scale. Since boudinage and folding are considered to express the same general instability, both structures should arise independently of the sign of the loading conditions and for identical material parameters. To this end, the link between material to energy instabilities is approached by means of bifurcation analyses of the field equations and finite element simulations of the coupled system of equations. Boudinage and folding structures develop at the same critical energy threshold, where dissipative work by temperature-sensitive creep overcomes the diffusive capacity of the layer. This finding provides basis for a unified theory for strain localization in layered ductile materials. The numerical simulations are compared to natural pinch-and-swell microstructures, tracing the adaption of grain sizes, textures and creep mechanisms in calcite veins. The switch from dislocation to diffusion creep relates to strain-rate weakening, which is induced by dissipated heat from grain size reduction, and marks the onset of continuous necking. The time-dependent sequence uncovers multiple steady states at different time intervals. Microstructurally and mechanically stable conditions are finally expressed in the pinch-and-swell end members. The major outcome of this study is that boudinage and folding can be described as the same coupled energy-mechanical bifurcation, or as one critical energy attractor. This finding allows the derivation of critical deformation conditions and fundamental material parameters directly from localized structures in the field.

Investigating the Effect of Histone H4 Sumoylation on Chromatin Structure and Function

Thesis (Ph.D.)--University of Washington, 2016-06