New self-assembled supramolecular hydrogels based on dehydropeptides

Supramolecular hydrogels rely on small molecules that self-assemble in water as a result of the cooperative effect of several relatively weak intermolecular interactions. Peptide-based low molecular weight hydrogelators have attracted enormous interest owing to the simplicity of small molecules combined with the versatility and biocompatibility of peptides. In this work, naproxen, a well known non-steroidal anti-inflammatory drug, was N-conjugated with various dehydrodipeptides to give aromatic peptide amphiphiles that resist proteolysis. Molecular dynamics simulations were used to obtain insight into the underlying molecular mechanism of self-assembly and to rationalize the design of this type of hydrogelators. The results obtained were in excellent agreement with the experimental observations. Only dehydrodipeptides having at least one aromatic amino acid gave hydrogels. The characterization of the hydrogels was carried out using transmission electron microscopy (TEM), circular dichroism (CD), fluorescence spectroscopy and also rheological assays.

Nanocomposites: the contribution of multiscale modelling to improve the dispersion levels

Tese de Doutoramento em Ciência e Engenharia de Polímeros e Compósitos

Synthesis of new peptide derivatives that self-assemble into nanostructured hydrogels for biomedical applications

Tese de Doutoramento em Ciências (área de especialização em Química)

Insights on the mechanical behavior of keratin fibrils

Supplementary data associated with this article can be found,in the online version, at http://dx.doi.org/10.1016/j.ijbiomac.2016.05.018.

Mutational analysis of the highly conserved arginine within the Glu/Asp-Arg-Tyr motif of the alpha(1b)-adrenergic receptor: effects on receptor isomerization and activation.

We have suggested previously that both the negatively and positively charged residues of the highly conserved Glu/Asp-Arg-Tyr (E/DRY) motif play an important role in the activation process of the alpha(1b)-adreneric receptor (AR). In this study, R143 of the E/DRY sequence in the alpha(1b)-AR was mutated into several amino acids (Lys, His, Glu, Asp, Ala, Asn, and Ile). The charge-conserving mutation of R143 into lysine not only preserved the maximal agonist-induced response of the alpha(1b)-AR, but it also conferred high degree of constitutive activity to the receptor. Both basal and agonist-induced phosphorylation levels were significantly increased for the R143K mutant compared with those of the wild-type receptor. Other substitutions of R143 resulted in receptor mutants with either a small increase in constitutive activity (R143H and R143D), impairment (R143H, R143D), or complete loss of receptor-mediated response (R143E, R143A, R143N, R143I). The R413E mutant displayed a small, but significant increase in basal phosphorylation despite being severely impaired in receptor-mediated response. Interestingly, all the arginine mutants displayed increased affinity for agonist binding compared with the wild-type alpha(1b)-AR. A correlation was found between the extent of the affinity shift and the intrinsic activity of the agonists. The analysis of the receptor mutants using the allosteric ternary complex model in conjunction with the results of molecular dynamics simulations on the receptor models support the hypothesis that mutations of R143 can drive the isomerization of the alpha(1b)-AR into different states, highlighting the crucial role of this residue in the activation process of the receptor.

Localization of breakage points in knotted strings

It is a common macroscopic observation that knotted ropes or fishing lines under tension easily break at the knot. However, a more precise localization of the breakage point in knotted macroscopic strings is a difficult task. In the present work, the tightening of knots was numerically simulated, a comparison of strength of different knots was experimentally performed and a high velocity camera was used to precisely localize the site where knotted macroscopic strings break. In the case of knotted spaghetti, the breakage occurs at the position with high curvature at the entry to the knot. This localization results from joint contributions of loading, bending and friction forces into the complex process of knot breakage. The present simulations and experiments are in agreement with recent molecular dynamics simulations of a knotted polymer chain and with experiments performed on actin and DNA filaments. The strength of the knotted string is greatly reduced (down to 50%) by the presence of a knot, therefore reducing the resistance to tension of all materials containing chains of any sort. The present work with macroscopic strings revels some important aspects, which are not accessible by experiments with microscopic chains.

Combined simulation and mutagenesis analyses reveal the involvement of key residues for peroxisome proliferator-activated receptor alpha helix 12 dynamic behavior.

The dynamic properties of helix 12 in the ligand binding domain of nuclear receptors are a major determinant of AF-2 domain activity. We investigated the molecular and structural basis of helix 12 mobility, as well as the involvement of individual residues with regard to peroxisome proliferator-activated receptor alpha (PPARalpha) constitutive and ligand-dependent transcriptional activity. Functional assays of the activity of PPARalpha helix 12 mutants were combined with free energy molecular dynamics simulations. The agreement between the results from these approaches allows us to make robust claims concerning the mechanisms that govern helix 12 functions. Our data support a model in which PPARalpha helix 12 transiently adopts a relatively stable active conformation even in the absence of a ligand. This conformation provides the interface for the recruitment of a coactivator and results in constitutive activity. The receptor agonists stabilize this conformation and increase PPARalpha transcription activation potential. Finally, we disclose important functions of residues in PPARalpha AF-2, which determine the positioning of helix 12 in the active conformation in the absence of a ligand. Substitution of these residues suppresses PPARalpha constitutive activity, without changing PPARalpha ligand-dependent activation potential.

Identification of potential small molecule binding pockets on Rho family GTPases

Rho GTPases are conformational switches that control a wide variety of signaling pathways critical for eukaryotic cell development and proliferation. They represent attractive targets for drug design as their aberrant function and deregulated activity is associated with many human diseases including cancer. Extensive high-resolution structures (.100) and recent mutagenesis studies have laid the foundation for the design of new structure-based chemotherapeutic strategies. Although the inhibition of Rho signaling with drug-like compounds is an active area of current research, very little attention has been devoted to directly inhibiting Rho by targeting potential allosteric non-nucleotide binding sites. By avoiding the nucleotide binding site, compounds may minimize the potential for undesirable off-target interactions with other ubiquitous GTP and ATP binding proteins. Here we describe the application of molecular dynamics simulations, principal component analysis, sequence conservation analysis, and ensemble small-molecule fragment mapping to provide an extensive mapping of potential small-molecule binding pockets on Rho family members. Characterized sites include novel pockets in the vicinity of the conformationaly responsive switch regions as well as distal sites that appear to be related to the conformations of the nucleotide binding region. Furthermore the use of accelerated molecular dynamics simulation, an advanced sampling method that extends the accessible time-scale of conventional simulations, is found to enhance the characterization of novel binding sites when conformational changes are important for the protein mechanism.

Understanding the two-component sensing System of virulent bacteria

Report for the scientific sojourn carried out at the l’ Institute for Computational Molecular Science of the Temple University, United States, from 2010 to 2012. Two-component systems (TCS) are used by pathogenic bacteria to sense the environment within a host and activate mechanisms related to virulence and antimicrobial resistance. A prototypical example is the PhoQ/PhoP system, which is the major regulator of virulence in Salmonella. Hence, PhoQ is an attractive target for the design of new antibiotics against foodborne diseases. Inhibition of the PhoQ-mediated bacterial virulence does not result in growth inhibition, presenting less selective pressure for the generation of antibiotic resistance. Moreover, PhoQ is a histidine kinase (HK) and it is absent in animals. Nevertheless, the design of satisfactory HK inhibitors has been proven to be a challenge. To compete with the intracellular ATP concentrations, the affinity of a HK inhibidor must be in the micromolar-nanomolar range, whereas the current lead compounds have at best millimolar affinities. Moreover, the drug selectivity depends on the conformation of a highly variable loop, referred to as the “ATP-lid, which is difficult to study by X-Ray crystallography due to its flexibility. I have investigated the binding of different HK inhibitors to PhoQ. In particular, all-atom molecular dynamics simulations have been combined with enhanced sampling techniques in order to provide structural and dynamic information of the conformation of the ATP-lid. Transient interactions between these drugs and the ATP-lid have been identified and the free energy of the different binding modes has been estimated. The results obtained pinpoint the importance of protein flexibility in the HK-inhibitor binding, and constitute a first step in developing more potent and selective drugs. The computational resources of the hosting institution as well as the experience of the members of the group in drug binding and free energy methods have been crucial to carry out this work.

T-Cell Receptors Binding Orientation over Peptide/MHC Class I Is Driven by Long-Range Interactions.

Crystallographic data about T-Cell Receptor - peptide - major histocompatibility complex class I (TCRpMHC) interaction have revealed extremely diverse TCR binding modes triggering antigen recognition. Understanding the molecular basis that governs TCR orientation over pMHC is still a considerable challenge. We present a simplified rigid approach applied on all non-redundant TCRpMHC crystal structures available. The CHARMM force field in combination with the FACTS implicit solvation model is used to study the role of long-distance interactions between the TCR and pMHC. We demonstrate that the sum of the coulomb interactions and the electrostatic solvation energies is sufficient to identify two orientations corresponding to energetic minima at 0° and 180° from the native orientation. Interestingly, these results are shown to be robust upon small structural variations of the TCR such as changes induced by Molecular Dynamics simulations, suggesting that shape complementarity is not required to obtain a reliable signal. Accurate energy minima are also identified by confronting unbound TCR crystal structures to pMHC. Furthermore, we decompose the electrostatic energy into residue contributions to estimate their role in the overall orientation. Results show that most of the driving force leading to the formation of the complex is defined by CDR1,2/MHC interactions. This long-distance contribution appears to be independent from the binding process itself, since it is reliably identified without considering neither short-range energy terms nor CDR induced fit upon binding. Ultimately, we present an attempt to predict the TCR/pMHC binding mode for a TCR structure obtained by homology modeling. The simplicity of the approach and the absence of any fitted parameters make it also easily applicable to other types of macromolecular protein complexes.

Constitutively active mutants of the alpha 1B-adrenergic receptor: role of highly conserved polar amino acids in receptor activation.

Site-directed mutagenesis and molecular dynamics simulations of the alpha 1B-adrenergic receptor (AR) were combined to explore the potential molecular changes correlated with the transition from R (inactive state) to R (active state). Using molecular dynamics analysis we compared the structural/dynamic features of constitutively active mutants with those of the wild type and of an inactive alpha 1B-AR to build a theoretical model which defines the essential features of R and R. The results of site-directed mutagenesis were in striking agreement with the predictions of the model supporting the following hypothesis. (i) The equilibrium between R and R depends on the equilibrium between the deprotonated and protonated forms, respectively, of D142 of the DRY motif. In fact, replacement of D142 with alanine confers high constitutive activity to the alpha 1B-AR. (ii) The shift of R143 of the DRY sequence out of a conserved 'polar pocket' formed by N63, D91, N344 and Y348 is a feature common to all the active structures, suggesting that the role of R143 is fundamental for mediating receptor activation. Disruption of these intramolecular interactions by replacing N63 with alanine constitutively activates the alpha 1B-AR. Our findings might provide interesting generalities about the activation process of G protein-coupled receptors.

Enhanced EGFP Fluorescence Emission in Presence of PEG Aqueous Solutions and PIB1000-PEG6000-PIB1000 Copolymer Vesicles.

An EGFP construct interacting with the PIB1000-PEG6000-PIB1000 vesicles surface reported a ~2-fold fluorescence emission enhancement. Because of the constructs nature with the amphiphilic peptide inserted into the PIB core, EGFP is expected to experience a "pure" PEG environment. To unravel this phenomenon PEG/water solutions at different molecular weights and concentrations were used. Already at ~1 : 10 protein/PEG molar ratio the increase in fluorescence emission is observed reaching a plateau correlating with the PEG molecular weight. Parallel experiments in presence of glycerol aqueous solutions did show a slight fluorescence enhancement however starting at much higher concentrations. Molecular dynamics simulations of EGFP in neat water, glycerol, and PEG aqueous solutions were performed showing that PEG molecules tend to "wrap" the protein creating a microenvironment where the local PEG concentration is higher compared to its bulk concentration. Because the fluorescent emission can be perturbed by the refractive index surrounding the protein, the clustering of PEG molecules induces an enhanced fluorescence emission already at extremely low concentrations. These findings can be important when related to the use of EGFP as reported in molecular biology experiments.

Heat transfer between nanoparticles: Thermal conductance for near-field interactions

We analyze the heat transfer between two nanoparticles separated by a distance lying in the near-field domain in which energy interchange is due to the Coulomb interactions. The thermal conductance is computed by assuming that the particles have charge distributions characterized by fluctuating multipole moments in equilibrium with heat baths at two different temperatures. This quantity follows from the fluctuation-dissipation theorem for the fluctuations of the multipolar moments. We compare the behavior of the conductance as a function of the distance between the particles with the result obtained by means of molecular dynamics simulations. The formalism proposed enables us to provide a comprehensive explanation of the marked growth of the conductance when decreasing the distance between the nanoparticles.

Metastable liquid-liquid phase transition in a single-component system with only one crystal phase and no density anomaly

We investigate the phase behavior of a single-component system in three dimensions with spherically-symmetric, pairwise-additive, soft-core interactions with an attractive well at a long distance, a repulsive soft-core shoulder at an intermediate distance, and a hard-core repulsion at a short distance, similar to potentials used to describe liquid systems such as colloids, protein solutions, or liquid metals. We showed [Nature (London) 409, 692 (2001)] that, even with no evidence of the density anomaly, the phase diagram has two first-order fluid-fluid phase transitions, one ending in a gas¿low-density-liquid (LDL) critical point, and the other in a gas¿high-density-liquid (HDL) critical point, with a LDL-HDL phase transition at low temperatures. Here we use integral equation calculations to explore the three-parameter space of the soft-core potential and perform molecular dynamics simulations in the interesting region of parameters. For the equilibrium phase diagram, we analyze the structure of the crystal phase and find that, within the considered range of densities, the structure is independent of the density. Then, we analyze in detail the fluid metastable phases and, by explicit thermodynamic calculation in the supercooled phase, we show the absence of the density anomaly. We suggest that this absence is related to the presence of only one stable crystal structure.

Liquid-liquid phase transition for soft-core attractive potential

Using event-driven molecular dynamics simulations, we study a three-dimensional one-component system of spherical particles interacting via a discontinuous potential combining a repulsive square soft core and an attractive square well. In the case of a narrow attractive well, it has been shown that this potential has two metastable gas-liquid critical points. Here we systematically investigate how the changes of the parameters of this potential affect the phase diagram of the system. We find a broad range of potential parameters for which the system has both a gas-liquid critical point C1 and a liquid-liquid critical point C2. For the liquid-gas critical point we find that the derivatives of the critical temperature and pressure, with respect to the parameters of the potential, have the same signs: they are positive for increasing width of the attractive well and negative for increasing width and repulsive energy of the soft core. This result resembles the behavior of the liquid-gas critical point for standard liquids. In contrast, for the liquid-liquid critical point the critical pressure decreases as the critical temperature increases. As a consequence, the liquid-liquid critical point exists at positive pressures only in a finite range of parameters. We present a modified van der Waals equation which qualitatively reproduces the behavior of both critical points within some range of parameters, and gives us insight on the mechanisms ruling the dependence of the two critical points on the potential¿s parameters. The soft-core potential studied here resembles model potentials used for colloids, proteins, and potentials that have been related to liquid metals, raising an interesting possibility that a liquid-liquid phase transition may be present in some systems where it has not yet been observed.