Structure of molten TbCl3 measured by neutron diffraction

The total structure factor of molten TbCl3 at 617ºC was measured by using neutron diffraction. The data are in agreement with results from previous experimental work but the use of a diffractometer having an extended reciprocal-space measurement window leads to improved resolution in real space. Significant discrepancies with the results obtained from recent molecular dynamics simulations carried out using a polarizable ion model, in which the interaction potentials were optimized to enhance agreement with previous diffraction data, are thereby highlighted. It is hence shown that there is considerable scope for the development of this model for TbCl3 and for other trivalent metal halide systems spanning a wide range of ion size ratios.

Visualising and controlling the flow in biomolecular systems at and between multiple scales:from atoms to hydrodynamics at different locations in time and space

A novel framework for modelling biomolecular systems at multiple scales in space and time simultaneously is described. The atomistic molecular dynamics representation is smoothly connected with a statistical continuum hydrodynamics description. The system behaves correctly at the limits of pure molecular dynamics (hydrodynamics) and at the intermediate regimes when the atoms move partly as atomistic particles, and at the same time follow the hydrodynamic flows. The corresponding contributions are controlled by a parameter, which is defined as an arbitrary function of space and time, thus, allowing an effective separation of the atomistic 'core' and continuum 'environment'. To fill the scale gap between the atomistic and the continuum representations our special purpose computer for molecular dynamics, MDGRAPE-4, as well as GPU-based computing were used for developing the framework. These hardware developments also include interactive molecular dynamics simulations that allow intervention of the modelling through force-feedback devices.

Development of potent and selective tissue transglutaminase inhibitors:their effect on TG2 function and application in pathological conditions

Potent-selective peptidomimetic inhibitors of tissue transglutaminase (TG2) were developed through a combination of protein-ligand docking and molecular dynamic techniques. Derivatives of these inhibitors were made with the aim of specific TG2 targeting to the intra- and extracellular space. A cell-permeable fluorescently labeled derivative enabled detection of in situ cellular TG2 activity in human umbilical cord endothelial cells and TG2-transduced NIH3T3 cells, which could be enhanced by treatment of cells with ionomycin. Reaction of TG2 with this fluorescent inhibitor in NIH3T3 cells resulted in loss of binding of TG2 to cell surface syndecan-4 and inhibition of translocation of the enzyme into the extracellular matrix, with a parallel reduction in fibronectin deposition. In human umbilical cord endothelial cells, this same fluorescent inhibitor also demonstrated a reduction in fibronectin deposition, cell motility, and cord formation in Matrigel. Use of the same inhibitor in a mouse model of hypertensive nephrosclerosis showed over a 40% reduction in collagen deposition.

Structural determinants of oligomerization of the aquaporin-4 channel

The aquaporins (AQP) family of integral membrane protein channels mediate cellular water and solute flow. Although qualitative and quantitative differences in channel permeability, selectivity, subcellular localization and trafficking responses have been observed for different members of the AQP family, the signature homotetrameric quaternary structure is conserved. Using a variety of biophysical techniques, we show that mutations to an intracellular loop (loop D) of human AQP4 reduce oligomerization. Non-tetrameric AQP4 mutants are unable to relocalize to the plasma membrane in response to changes in extracellular tonicity, despite equivalent constitutive surface expression levels and water permeability to wild-type AQP4. A network of AQP4 loop D hydrogen bonding interactions, identified using molecular dynamics simulations and based on a comparative mutagenic analysis of AQPs 1, 3 and 4, suggest that loop D interactions may provide a general structural framework for tetrameric assembly within the AQP family.

Movement Effects on the Flow Physics and Nutrient Delivery in Engineered Valvular Tissues

Mechanical conditioning has been shown to promote tissue formation in a wide variety of tissue engineering efforts. However the underlying mechanisms by which external mechanical stimuli regulate cells and tissues are not known. This is particularly relevant in the area of heart valve tissue engineering (HVTE) owing to the intense hemodynamic environments that surround native valves. Some studies suggest that oscillatory shear stress (OSS) caused by steady flow and scaffold flexure play a critical role in engineered tissue formation derived from bone marrow derived stem cells (BMSCs). In addition, scaffold flexure may enhance nutrient (e.g. oxygen, glucose) transport. In this study, we computationally quantified the i) magnitude of fluid-induced shear stresses; ii) the extent of temporal fluid oscillations in the flow field using the oscillatory shear index (OSI) parameter, and iii) glucose and oxygen mass transport profiles. Noting that sample cyclic flexure induces a high degree of oscillatory shear stress (OSS), we incorporated moving boundary computational fluid dynamic simulations of samples housed within a bioreactor to consider the effects of: 1) no flow, no flexure (control group), 2) steady flow-alone, 3) cyclic flexure-alone and 4) combined steady flow and cyclic flexure environments. We also coupled a diffusion and convention mass transport equation to the simulated system. We found that the coexistence of both OSS and appreciable shear stress magnitudes, described by the newly introduced parameter OSI-t , explained the high levels of engineered collagen previously observed from combining cyclic flexure and steady flow states. On the other hand, each of these metrics on its own showed no association. This finding suggests that cyclic flexure and steady flow synergistically promote engineered heart valve tissue production via OSS, so long as the oscillations are accompanied by a critical magnitude of shear stress. In addition, our simulations showed that mass transport of glucose and oxygen is enhanced by sample movement at low sample porosities, but did not play a role in highly porous scaffolds. Preliminary in-house in vitro experiments showed that cell proliferation and phenotype is enhanced in OSI-t environments.

Avaliação funcional e estrutural de um novo peptídeo antimicrobiano do escorpião Tityus stigmurus

In Brazil, there is a high incidence of venomous animals. Among them, scorpions are highlighted by their medical importance, and for being their venom a source of several molecules with biological and pharmacological activity not yet fully understood, including several bioactive peptides. Antimicrobial peptides (AMPs) are components of the immune system in prokaryotes and eukaryotes, used in the first line of defense against microorganisms. In the present study, we characterized the first PAM previously identified through transcriptome of the venom gland of the scorpion Tityus stigmurus, named Stigmurin. The characteristics of Stigmurin were investigated by computational modeling and construction of dendrogram. In vitro tests investigated the antibacterial, antifungal, haemolytic and cytotoxic effects of crude venom and Stigmurin. In addition, the structural characteristics of Stigmurin were investigated by circular dochroism in water, 2, 2 , 2- trifluoethanol (TFE) and sodium dodecyl sulfate (SDS) and the models were refined by molecular dynamics simulations. The results showed that the selected sequence encodes a mature protein of 17 amino acid residues and the dendrogram reveals a case of convergent evolution. The crude venom showed no antimicrobial activity, however, Stigmurin exhibited a broad spectrum of antibacterial activity, with minimal inhibitory concentrations (MIC) ranging from 31.25 and 250 µg/mL for different strains, while the hemolytic activity at these concentrations was low. In cytotoxicity studies, the crude venom was unable to reduce cell viability in VERO E6 cells; in contrast, its activity in SiHa cells was significantly higher, corresponding to IC50 of 3.6 µg/mL. For Stigmurin the concentration sable to decrease cell viability of Vero E6 and SiHa cells in 50% were 275.67 µg/mL and 212.54 µg/mL, respectively. The dichroism spectra revealed the conformational flexibility, with predominating extended and β–sheet structures, as well as a remark able renaturation ability. The results suggest that Stigmurin could be considered as a potential antiinfective drug

An allosteric role for receptor activity-modifying proteins in defining GPCR pharmacology

G protein-coupled receptors are allosteric proteins that control transmission of external signals to regulate cellular response. Although agonist binding promotes canonical G protein signalling transmitted through conformational changes, G protein-coupled receptors also interact with other proteins. These include other G protein-coupled receptors, other receptors and channels, regulatory proteins and receptor-modifying proteins, notably receptor activity-modifying proteins (RAMPs). RAMPs have at least 11 G protein-coupled receptor partners, including many class B G protein-coupled receptors. Prototypic is the calcitonin receptor, with altered ligand specificity when co-expressed with RAMPs. To gain molecular insight into the consequences of this protein–protein interaction, we combined molecular modelling with mutagenesis of the calcitonin receptor extracellular domain, assessed in ligand binding and functional assays. Although some calcitonin receptor residues are universally important for peptide interactions (calcitonin, amylin and calcitonin gene-related peptide) in calcitonin receptor alone or with receptor activity-modifying protein, others have RAMP-dependent effects, whereby mutations decreased amylin/calcitonin gene-related peptide potency substantially only when RAMP was present. Remarkably, the key residues were completely conserved between calcitonin receptor and AMY receptors, and between subtypes of AMY receptor that have different ligand preferences. Mutations at the interface between calcitonin receptor and RAMP affected ligand pharmacology in a RAMP-dependent manner, suggesting that RAMP may allosterically influence the calcitonin receptor conformation. Supporting this, molecular dynamics simulations suggested that the calcitonin receptor extracellular N-terminal domain is more flexible in the presence of receptor activity-modifying protein 1. Thus, RAMPs may act in an allosteric manner to generate a spectrum of unique calcitonin receptor conformational states, explaining the pharmacological preferences of calcitonin receptor-RAMP complexes. This provides novel insight into our understanding of G protein-coupled receptor-protein interaction that is likely broadly applicable for this receptor class.

Formation of helical ion chains

We study the nonequilibrium dynamics of the linear to zigzag structural phase transition exhibited by an ion chain confined in a trap with periodic boundary conditions. The transition is driven by reducing the transverse confinement at a finite quench rate, which can be accurately controlled. This results in the formation of zigzag domains oriented along different transverse planes. The twists between different domains can be stabilized by the topology of the trap and under laser cooling the system has a chance to relax to a helical chain with nonzero winding number. Molecular dynamics simulations are used to obtain a large sample of possible trajectories for different quench rates. The scaling of the average winding number with different quench rates is compared to the prediction of the Kibble-Zurek theory, and a good quantitative agreement is found.

How curvature-generating proteins build scaffolds on membrane nanotubes

Bin/Amphiphysin/Rvs (BAR) domain proteins control the curvature of lipid membranes in endocytosis, trafficking, cell motility, the formation of complex sub-cellular structures, and many other cellular phenomena. They form three-dimensional assemblies, which act as molecular scaffolds to reshape the membrane and alter its mechanical properties. It is unknown, however, how a protein scaffold forms and how BAR domains interact in these assemblies at protein densities relevant for a cell. In this work, we employ various experimental, theoretical and simulation approaches to explore how BAR proteins organize to form a scaffold on a membrane nanotube. By combining quantitative microscopy with analytical modeling, we demonstrate that a highly curving BAR protein endophilin nucleates its scaffolds at the ends of a membrane tube, contrary to a weaker curving protein centaurin, which binds evenly along the tube’s length. Our work implies that the nature of local protein-membrane interactions can affect the specific localization of proteins on membrane-remodeling sites. Furthermore, we show that amphipathic helices are dispensable in forming protein scaffolds. Finally, we explore a possible molecular structure of a BAR-domain scaffold using coarse-grained molecular dynamics simulations. Together with fluorescence microscopy, the simulations show that proteins need only to cover 30–40% of a tube’s surface to form a rigid assembly. Our work provides mechanical and structural insights into the way BAR proteins may sculpt the membrane as a high-order cooperative assembly in important biological processes. 

Interactions of Cytochrome P450 3A4 with Phospholipid Bilayer Nanodiscs

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

Studies of Protein-Protein and Protein-Water Interactions by Small Angle X-Ray Scattering, Terahertz Spectroscopy, ASMOS, And Computer Simulation

The protein folding problem has been one of the most challenging subjects in biological physics due to its complexity. Energy landscape theory based on statistical mechanics provides a thermodynamic interpretation of the protein folding process. We have been working to answer fundamental questions about protein-protein and protein-water interactions, which are very important for describing the energy landscape surface of proteins correctly. At first, we present a new method for computing protein-protein interaction potentials of solvated proteins directly from SAXS data. An ensemble of proteins was modeled by Metropolis Monte Carlo and Molecular Dynamics simulations, and the global X-ray scattering of the whole model ensemble was computed at each snapshot of the simulation. The interaction potential model was optimized and iterated by a Levenberg-Marquardt algorithm. Secondly, we report that terahertz spectroscopy directly probes hydration dynamics around proteins and determines the size of the dynamical hydration shell. We also present the sequence and pH-dependence of the hydration shell and the effect of the hydrophobicity. On the other hand, kinetic terahertz absorption (KITA) spectroscopy is introduced to study the refolding kinetics of ubiquitin and its mutants. KITA results are compared to small angle X-ray scattering, tryptophan fluorescence, and circular dichroism results. We propose that KITA monitors the rearrangement of hydrogen bonding during secondary structure formation. Finally, we present development of the automated single molecule operating system (ASMOS) for a high throughput single molecule detector, which levitates a single protein molecule in a 10 µm diameter droplet by the laser guidance. I also have performed supporting calculations and simulations with my own program codes.

Human Ezrin: Prediction of the 3D structure and interactions with mTOR

The protein Ezrin, is a member of the ERM family (Ezrin, Radixin and Moesin) that links the F-actin to the plasma membrane. The protein is made of three domains namely the FERM domain, a central α-helical domain and the CERMAD domain. The residues in Ezrin such as Ser66, Tyr145, Tyr353 and Tyr477 regulate the function of the protein through phosphorylation. The protein is found in two distinct conformations of active and dormant (inactive) state. The initial step during the conformation change is the breakage of intramolecular interaction in dormant Ezrin by phosphorylation of residue Thr567. The dormant structure of human Ezrin was predicted computationally since only partial active form structure was available. The validation analysis showed that 99.7% residues were positioned in favored, allowed and generously allowed regions of the Ramachandran plot. The Z-score of Ezrin was −7.36, G-factor was 0.1, and the QMEAN score of the model was 0.61 indicating a good model for human Ezrin. The comparison of the conformations of the activated and dormant Ezrin showed a major shift in the F2 lobe (residues 142-149 and 161-177) while changes in the conformation induced mobility shifts in lobe F3 (residues 261 to 267). The 3D positions of the phosphorylation sites Tyr145, Tyr353, Tyr477, Tyr482 and Thr567 were also located. Using targeted molecular dynamic simulation, the molecular movements during conformational change from active to dormant were visualized. The dormant Ezrin auto-inhibits itself by a head-to-tail interaction of the N-terminal and C-terminal residues. The trajectory shows the breakage of the interactions and mobility of the CERMAD domain away from the FERM domain. Protein docking and clustering analysis were used to predict the residues involved in the interaction between dormant Ezrin and mTOR. Residues Tyr477 and Tyr482 were found to be involved in interaction with mTOR.

Mechanism of pH-dependent activation of the sodium-proton antiporter NhaA

Data files to accompany the article in Nature Communications, in press.

Non-inclusion complexes between riboflavin and cyclodextrins

Objectives To investigate the molecular interaction between beta-cyclodextrin (beta CD) or hydroxypropyl-beta-cyclodextrin (HP beta CD) and riboflavin (RF), and to test the anticancer potential of these formulations. Methods The physicochemical characterization of the association between RF and CDs was performed by UV-vis absorption, fluorescence, differential scanning calorimetry and NMR techniques. Molecular dynamics simulation was used to shed light on the mechanism of interaction of RF and CDs. Additionally, in-vitro cell culture tests were performed to evaluate the cytotoxicity of the RFCD complexes against prostate cancer cells. Key findings Neither beta CD nor HP beta CD led to substantial changes in the physicochemical properties of RF (with the exception of solubility). Additionally, rotating frame Overhauser effect spectroscopy experiments detected no spatial correlations between hydrogens from the internal cavity of CDs and RF, while molecular dynamics simulations revealed out-of-ring RFCD interactions. Notwithstanding, both RF beta CD and RFHP beta CD complexes were cytotoxic to PC3 prostate cancer cells. Conclusions The interaction between RF and either beta CD or HP beta CD, at low concentrations, seems to be made through hydrogen bonding between the flavonoid and the external rim of both CDs. Regardless of the mechanism of complexation, our findings indicate that RFCD complexes significantly increase RF solubility and potentiate its antitumour effect.

Understanding the cation specific effects on the aqueous solubility of amino acids: from mono to polyvalent cations

The interactions established by mono and polyvalent cations in natural media have important implications on the structure formation, function and physico-chemical behavior of biomolecules, playing therefore a critical role in biochemical processes. In order to further elucidate the molecular phenomena behind the cation specific effects in biological environments, and clarify the influence of the charge of the ions, solubility measurements and molecular dynamics simulations were performed for aqueous solutions of three amino acids (alanine, valine and isoleucine), in the presence of a series of inorganic salts comprising mono-, di- and trivalent cations (LiCl, Li2SO4, K2SO4, CaCl2, AlCl3 and Al-2(SO4)(3)). The evidence gathered indicates that the mechanism by which (salting-in inducing) polyvalent cations affect the solubility of amino acids in aqueous solutions is different from that of monovalent cations. A consistent and refined molecular description of the effect of the cation on the solubility of amino acids based on specific interactions of the cations with the negatively charged moieties of the biomolecules is here proposed.