Intracellular transport is crucial for many cellular processes where a large fraction of the cargo is transferred by motor-proteins over a network of microtubules. Malfunctions in the transport mechanism underlie a number of medical maladies.
Existing methods for studying how motor-proteins coordinate the transfer of a shared cargo over a microtubule are either analytical or are based on Monte-Carlo simulations. Approaches that yield analytical results, while providing unique insights into transport mechanism, make simplifying assumptions, where a detailed characterization of important transport modalities is difficult to reach. On the other hand, Monte-Carlo based simulations can incorporate detailed characteristics of the transport mechanism; however, the quality of the results depend on the number and quality of simulation runs used in arriving at results. Here, for example, it is difficult to simulate and study rare-events that can trigger abnormalities in transport.
In this article, a semi-analytical methodology that determines the probability distribution function of motor-protein behavior in an exact manner is developed. The method utilizes a finite-dimensional projection of the underlying infinite-dimensional Markov model, which retains the Markov property, and enables the detailed and exact determination of motor configurations, from which meaningful inferences on transport characteristics of the original model can be derived.
Under this novel probabilistic approach new insights about the mechanisms of action of these proteins are found, suggesting hypothesis about their behavior and driving the design and realization of new experiments.
The advantages provided in accuracy and efficiency make it possible to detect rare events in the motor protein dynamics, that could otherwise pass undetected using standard simulation methods. In this respect, the model has allowed to provide a possible explanation for possible mechanisms under which motor proteins could coordinate their motion.
Molecular motors; Rare event detection; Markov models
Keratins are important structural proteins found in skin, hair and nails. Keratin Intermediate Filaments are major components of corneocytes, nonviable horny cells of the Stratum Corneum, the outermost layer of skin. It is considered that interactions between unstructured domains of Keratin Intermediate Filaments are the key factor in maintaining the elasticity of the skin.
We have developed a model for the interactions between keratin intermediate filaments based on self-consistent field theory. The intermediate filaments are represented by charged surfaces, and the disordered terminal domains of the keratins are represented by charged heteropolymers grafted to these surfaces. We estimate the system is close to a charge compensation point where the heteropolymer grafting density is matched to the surface charge density. Using a protein model with amino acid resolution for the terminal domains, we find that the terminal chains can mediate a weak attraction between the keratin surfaces. The origin of the attraction is a combination of bridging and electrostatics. The attraction disappears when the system moves away from the charge compensation point, or when excess small ions and/or NMF-representing free amino acids are added.
These results are in concordance with experimental observations, and support the idea that the interaction between keratin filaments, and ultimately in part the elastic properties of the keratin-containing tissue, is controlled by a combination of the physico-chemical properties of the disordered terminal domains and the composition of the medium in the inter-filament region.
Stratum corneum; Skin keratins; Intermediate filaments; Unstructured terminal domains; Bridging attraction
Increasing applications of titanium dioxide (TiO2) fine particles (FPs) and nanoparticles (NPs) require coupled knowledge improvement concerning their biokinetic effects. Neutrophils are quickly recruited to titanium implantation areas. Neutrophils mechanical properties display a crucial role on cell physiology and immune responsive functions. Then, micro and nanomechanical characterization assessed by force spectroscopy (FS) technique has been largely applied in this field.
Scanning electron microscopy (SEM) images highlighted neutrophils morphological changes along TiO2 FPs and NPs aggregates exposure time (1, 5, and 30 min) compared to controls. FS approaches showed an increasing on attraction forces to TiO2 FPs and NPs treated neutrophils. This group depicted stronger stiffness features than controls just at 1 min of exposure. Treated neutrophils showed a tendency to increase adhesive properties after 1 and 5 min of exposure. These cells maintained comparatively higher elasticity behavior for a longer time possibly due to intense phagocytosis and cell stiffness opposing to the tip indentation. Neutrophils activation caused by FPs and NPs uptake could be related to increasing dissipated energy results.
Mechanical modifications resulted from TiO2 FPs and NPs aggregates interaction with neutrophils showed increasing stiffness and also cell morphology alteration. Cells treatment by this metal FPs and NPs caused an increase in attractive forces. This event was mainly observed on the initial exposure times probably regarding to the interaction of neutrophils membrane and phagocytosis. Similar results were found to adhesion forces and dissipated energy outcomes. Treated cells presented comparatively higher elasticity behavior for a longer time. SEM images clearly suggested cell morphology alteration along time course probably related to activation, cytoskeleton rearrangement and phagocytosis. This scenario with increase in stiffness strongly suggests a direct relationship over neutrophil rolling, arrest, and transmigration. Scrutinizing these interactions represents an essential step to clarify the mechanisms involved on treatments containing micro and nanomaterials and their fates on the organisms.
Force spectroscopy; Neutrophil nanomechanics; Titanium dioxide microparticles; Titanium dioxide nanoparticles
PQS (PseudomonasQuinolone Signal) and its precursor HHQ are signal molecules of the P. aeruginosa quorum sensing system. They explicate their role in mammalian pathogenicity by binding to the receptor PqsR that induces virulence factor production and biofilm formation. The enzyme PqsD catalyses the biosynthesis of HHQ.
Enzyme kinetic analysis and surface plasmon resonance (SPR) biosensor experiments were used to determine mechanism and substrate order of the biosynthesis. Comparative analysis led to the identification of domains involved in functionality of PqsD. A kinetic cycle was set up and molecular dynamics (MD) simulations were used to study the molecular bases of the kinetics of PqsD. Trajectory analysis, pocket volume measurements, binding energy estimations and decompositions ensured insights into the binding mode of the substrates anthraniloyl-CoA and β-ketodecanoic acid.
Enzyme kinetics and SPR experiments hint at a ping-pong mechanism for PqsD with ACoA as first substrate. Trajectory analysis of different PqsD complexes evidenced ligand-dependent induced-fit motions affecting the modified ACoA funnel access to the exposure of a secondary channel. A tunnel-network is formed in which Ser317 plays an important role by binding to both substrates. Mutagenesis experiments resulting in the inactive S317F mutant confirmed the importance of this residue. Two binding modes for β-ketodecanoic acid were identified with distinct catalytic mechanism preferences.
Water is essential for life, but some organisms can survive complete desiccation, while many more survive partial dehydration during drying or freezing. The function of some protective molecules, such as sugars, has been extensively studied, but much less is known about the effects of amphiphiles such as flavonoids and other aromatic compounds. Amphiphiles may be largely soluble under fully hydrated conditions, but will partition into membranes upon removal of water. Little is known about the effects of amphiphiles on membrane stability and how amphiphile structure and function are related. Here, we have used two of the most intensively studied amphiphiles, tryptophan (Trp) and arbutin (Arb), along with their isolated hydrophilic moieties glycine (Gly) and glucose (Glc) to better understand structure-function relationships in amphiphile-membrane interactions in the dry state.
Fourier-transform infrared (FTIR) spectroscopy was used to measure gel-to-liquid crystalline phase transition temperatures (Tm) of liposomes formed from phosphatidylcholine and phosphatidylethanolamine in the presence of the different additives. In anhydrous samples, both Glc and Arb strongly depressed Tm, independent of lipid composition, while Gly had no measurable effect. Trp, on the other hand, either depressed or increased Tm, depending on lipid composition. We found no evidence for strong interactions of any of the compounds with the lipid carbonyl or choline groups, while all additives except Gly seemed to interact with the phosphate groups. In the case of Arb and Glc, this also had a strong effect on the sugar OH vibrations in the FTIR spectra. In addition, vibrations from the hydrophobic indole and phenol moieties of Trp and Arb, respectively, provided evidence for interactions with the lipid bilayers.
The two amphiphiles Arb and Trp interact differently with dry bilayers. The interactions of Arb are dominated by contributions of the Glc moiety, while the indole governs the effects of Trp. In addition, only Trp-membrane interactions showed a strong influence of lipid composition. Further investigations, using the large structural diversity of plant amphiphiles will help to understand how their structure determines the interaction with membranes and how that influences their biological functions, for example under freezing or dehydration conditions.
Amphiphiles; Arbutin; Desiccation; Fourier-transform infrared spectroscopy; Lipid phase transition; Model membranes; Tryptophan
Drop drying is a key factor in a wide range of technical applications, including spotted microarrays. The applied nL liquid volume provides specific reaction conditions for the immobilization of probe molecules to a chemically modified surface.
We investigated the influence of nL and μL liquid drop volumes on the process of probe immobilization and compare the results obtained to the situation in liquid solution. In our data, we observe a strong relationship between drop drying effects on immobilization and surface chemistry. In this work, we present results on the immobilization of dye labeled 20mer oligonucleotides with and without an activating 5′-aminoheptyl linker onto a 2D epoxysilane and a 3D NHS activated hydrogel surface.
Our experiments identified two basic processes determining immobilization. First, the rate of drop drying that depends on the drop volume and the ambient relative humidity. Oligonucleotides in a dried spot react unspecifically with the surface and long reaction times are needed. 3D hydrogel surfaces allow for immobilization in a liquid environment under diffusive conditions. Here, oligonucleotide immobilization is much faster and a specific reaction with the reactive linker group is observed. Second, the effect of increasing probe concentration as a result of drop drying. On a 3D hydrogel, the increasing concentration of probe molecules in nL spotting volumes accelerates immobilization dramatically. In case of μL volumes, immobilization depends on whether the drop is allowed to dry completely. At non-drying conditions, very limited immobilization is observed due to the low oligonucleotide concentration used in microarray spotting solutions. The results of our study provide a general guideline for microarray assay development. They allow for the initial definition and further optimization of reaction conditions for the immobilization of oligonucleotides and other probe molecule classes to different surfaces in dependence of the applied spotting and reaction volume.
A recent study by Dietz et al. using single-molecule fluorescence microscopy techniques demonstrates that, in the absence of the ligand InlB, the MET receptor exists as both a monomer and a dimer on the cell membrane, and addition of the ligand leads to increased MET dimerization. Under the crowded conditions of the cell membrane, dimer formation may be a common phenomenon for cell surface receptors. Ligand binding to both monomeric and dimeric receptors may provide parallel routes to receptor activation.
The human receptor tyrosine kinase MET and its ligand hepatocyte growth factor/scatter factor are essential during embryonic development and play an important role during cancer metastasis and tissue regeneration. In addition, it was found that MET is also relevant for infectious diseases and is the target of different bacteria, amongst them Listeria monocytogenes that induces bacterial uptake through the surface protein internalin B. Binding of ligand to the MET receptor is proposed to lead to receptor dimerization. However, it is also discussed whether preformed MET dimers exist on the cell membrane.
To address these issues we used single-molecule fluorescence microscopy techniques. Our photobleaching experiments show that MET exists in dimers on the membrane of cells in the absence of ligand and that the proportion of MET dimers increases significantly upon ligand binding.
Our results indicate that partially preformed MET dimers may play a role in ligand binding or MET signaling. The addition of the bacterial ligand internalin B leads to an increase of MET dimers which is in agreement with the model of ligand-induced dimerization of receptor tyrosine kinases.
MET receptor; Dimerization; Single-molecule photobleaching; Fluorescence correlation spectroscopy; Fluorescence; Signal transduction
Biomolecules or other complex macromolecules undergo conformational transitions upon exposure to an external perturbation such as ligand binding or mechanical force. To follow fluctuations in pairwise forces between atoms or residues during such conformational changes as observed in Molecular Dynamics (MD) simulations, we developed Time-Resolved Force Distribution Analysis (TRFDA).
The implementation focuses on computational efficiency and low-memory usage and, along with the wide range of output options, makes possible time series analysis of pairwise forces variation in long MD simulations and for large molecular systems. It also provides an exact decomposition of pairwise forces resulting from 3- and 4-body potentials and a unified treatment of pairwise forces between atoms or residues. As a proof of concept, we present a stress analysis during unfolding of ubiquitin in a force-clamp MD simulation.
TRFDA can be used, among others, in tracking signal propagation at atomic level, for characterizing dynamical intermolecular interactions (e.g. protein-ligand during flexible docking), in development of force fields and for following stress distribution during conformational changes.
UV–vis spectrophotometric optical density (OD) is the most commonly-used technique for estimating chromophore formation and cell concentration in liquid culture. OD wavelength is often chosen with little thought given to its effect on the quality of the measurement. Analysis of the contributions of absorption and scattering to the measured optical density provides a basis for understanding variability among spectrophotometers and enables a quantitative evaluation of the applicability of the Beer-Lambert law. This provides a rational approach for improving the accuracy of OD measurements used as a proxy for direct dry weight (DW), cell count, and pigment levels.
For pigmented organisms, the choice of OD wavelength presents a tradeoff between the robustness and the sensitivity of the measurement. The OD at a robust wavelength is primarily the result of light scattering and does not vary with culture conditions; whereas, the OD at a sensitive wavelength is additionally dependent on light absorption by the organism’s pigments. Suitably robust and sensitive wavelengths are identified for a wide range of organisms by comparing their spectra to the true absorption spectra of dyes. The relative scattering contribution can be reduced either by measurement at higher OD, or by the addition of bovine serum albumin. Reduction of scattering or correlation with off-peak light attenuation provides for more accurate assessment of chromophore levels within cells. Conversion factors between DW, OD, and colony-forming unit density are tabulated for 17 diverse organisms to illustrate the scope of variability of these correlations. Finally, an inexpensive short pathlength LED-based flow cell is demonstrated for the online monitoring of growth in a bioreactor at culture concentrations greater than 5 grams dry weight per liter which would otherwise require off-line dilutions to obtain non-saturated OD measurements.
OD is most accurate as a time-saving proxy measurement for biomass concentration when light attenuation is dominated by scattering. However, the applicability of OD-based correlations is highly dependent on the measurement specifications (spectrophotometer model and wavelength) and culture conditions (media type; growth stage; culture stress; cell/colony geometry; presence and concentration of secreted compounds). These variations highlight the importance of treating literature conversion factors as rough approximations as opposed to concrete constants. There is an opportunity to optimize measurements of cell pigment levels by considering scattering and absorption-dependent wavelengths of the OD spectrum.
Optical density; Light scattering; Algae; Conversion factors; Dry weight; Flow cell
Fluctuation-induced phenomena caused by both random and deterministic stimuli have been previously studied in a variety of contexts. They are based on the interplay between the spectro-temporal patterns of the signal and the kinetics of the system it is applied to. The aim of this study was to develop a method for designing fluctuating inputs into nonlinear system which would elicit the most desired system output and to implement the method to studies of ion channels.
We describe an algorithm based on constructing the input as a superposition of wavelets and optimizing it according to a selected cost functional. The algorithm is applied to ion channel electrophysiology where the input is the fluctuating voltage delivered through a patch-clamp experimental apparatus and the output is the whole-cell ionic current. The algorithm is optimized to aid selection of Markov models of the gating kinetics of the voltage-gated Shaker K+ channel and tested by comparison of numerically obtained ionic currents predicted by different models with experimental data obtained from the Shaker K+ channels. Other applications and optimization criteria are also suggested.
The method described in this paper can be useful in development and testing of models of ion channel gating kinetics, developing voltage inputs that optimize certain nonequilibrium phenomena in ion channels, such as the kinetic focusing, and potentially has applications to other fields.
Voltage-gated Shaker potassium channel; Markov model; Wavelet; Patch clamping
Transcription factories are nuclear domains where gene transcription takes place
although the molecular basis for their formation and maintenance are unknown. In this
study, we explored how the properties of chromatin as a polymer may contribute to the
structure of transcription factories. We found that transcriptional active chromatin
contains modifications like histone H4 acetylated at Lysine 16 (H4K16ac). Single
fibre analysis showed that this modification spans the entire body of the gene.
Furthermore, H4K16ac genes cluster in regions up to 500 Kb alternating active and
inactive chromatin. The introduction of H4K16ac in chromatin induces stiffness in the
chromatin fibre. The result of this change in flexibility is that chromatin could
behave like a multi-block copolymer with repetitions of stiff-flexible
(active-inactive chromatin) components. Copolymers with such structure self-organize
through spontaneous phase separation into microdomains. Consistent with such model
H4K16ac chromatin form foci that associates with nascent transcripts. We propose that
transcription factories are the result of the spontaneous concentration of H4K16ac
chromatin that are in proximity, mainly in cis.
Epigenetics; Biophysics; H4K16Ac; BrUTP; Transcription Factories; RNA pol II; Nuclear organization
Background and objective
Regulating protein function in the cell by small molecules, provide a rapid, reversible and tunable tool of metabolic control. However, due to its complexity the issue is poorly studied so far. The effects of small solutes on protein behavior can be studied by examining changes of protein secondary structure, in its hydrodynamic radius as well as its thermal aggregation. The study aim was to investigate effects of adenosine-5’-triphosphate (ATP), spermine NONOate (NO donor) as well as sodium/potassium ions on thermal aggregation of albumin and hemoglobin. To follow aggregation of the proteins, their diffusion coefficients were measured by quasi-elastic light scattering (QELS) at constant pH (7.4) in the presence of solutes over a temperature range from 25°C to 80°C.
Results and discussion
1) Spermine NONOate persistently decreased the hemoglobin aggregation temperature Tairrespectively of the Na+/K+ environment, 2) ATP alone had no effect on the protein’s thermal stability but it facilitated protein’s destabilization in the presence of spermine NONOate and 3) mutual effects of ATP and NO were strongly influenced by particular buffer ionic compositions.
The ATP effect on protein aggregation was ambiguous: ATP alone had no effect on the protein’s thermal stability but it facilitated protein’s destabilization in the presence of nitric oxide. The magnitude and direction of the observed effects strongly depended on concentrations of K+ and Na+ in the solution.
In order to replicate within their cellular host, many viruses have developed self-assembly strategies for their capsids which are sufficiently robust as to be reconstituted in vitro. Mathematical models for virus self-assembly usually assume that the bonds leading to cluster formation have constant reactivity over the time course of assembly (direct assembly). In some cases, however, binding sites between the capsomers have been reported to be activated during the self-assembly process (hierarchical assembly).
In order to study possible advantages of such hierarchical schemes for icosahedral virus capsid assembly, we use Brownian dynamics simulations of a patchy particle model that allows us to switch binding sites on and off during assembly. For T1 viruses, we implement a hierarchical assembly scheme where inter-capsomer bonds become active only if a complete pentamer has been assembled. We find direct assembly to be favorable for reversible bonds allowing for repeated structural reorganizations, while hierarchical assembly is favorable for strong bonds with small dissociation rate, as this situation is less prone to kinetic trapping. However, at the same time it is more vulnerable to monomer starvation during the final phase. Increasing the number of initial monomers does have only a weak effect on these general features. The differences between the two assembly schemes become more pronounced for more complex virus geometries, as shown here for T3 viruses, which assemble through homogeneous pentamers and heterogeneous hexamers in the hierarchical scheme. In order to complement the simulations for this more complicated case, we introduce a master equation approach that agrees well with the simulation results.
Our analysis shows for which molecular parameters hierarchical assembly schemes can outperform direct ones and suggests that viruses with high bond stability might prefer hierarchical assembly schemes. These insights increase our physical understanding of an essential biological process, with many interesting potential applications in medicine and materials science.
Understanding the mechanical properties of chromatin is an essential step towards deciphering the physical rules of gene regulation. In the past ten years, many single molecule experiments have been carried out, and high resolution measurements of the chromatin fiber stiffness are now available. Simulations have been used in order to link those measurements with structural cues, but so far no clear agreement among different groups has been reached.
We revisit here some of the most precise experimental results obtained with carefully reconstituted fibers.
We show that the mechanical properties of the chromatin fiber can be quantitatively accounted for by the stiffness of the DNA molecule and the 3D structure of the chromatin fiber.
Cholesterol is an important membrane component, but our knowledge about its transport in cells is sparse. Previous imaging studies using dehydroergosterol (DHE), an intrinsically fluorescent sterol from yeast, have established that vesicular and non-vesicular transport modes contribute to sterol trafficking from the plasma membrane. Significant photobleaching, however, limits the possibilities for in-depth analysis of sterol dynamics using DHE. Co-trafficking studies with DHE and the recently introduced fluorescent cholesterol analog BODIPY-cholesterol (BChol) suggested that the latter probe has utility for prolonged live-cell imaging of sterol transport.
We found that BChol is very photostable under two-photon (2P)-excitation allowing the acquisition of several hundred frames without significant photobleaching. Therefore, long-term tracking and diffusion measurements are possible. Two-photon temporal image correlation spectroscopy (2P-TICS) provided evidence for spatially heterogeneous diffusion constants of BChol varying over two orders of magnitude from the cell interior towards the plasma membrane, where D ~ 1.3 μm2/s. Number and brightness (N&B) analysis together with stochastic simulations suggest that transient partitioning of BChol into convoluted membranes slows local sterol diffusion. We observed sterol endocytosis as well as fusion and fission of sterol-containing endocytic vesicles. The mobility of endocytic vesicles, as studied by particle tracking, is well described by a model for anomalous subdiffusion on short time scales with an anomalous exponent α ~ 0.63 and an anomalous diffusion constant of Dα = 1.95 x 10-3 μm2/sα. On a longer time scale (t > ~5 s), a transition to superdiffusion consistent with slow directed transport with an average velocity of v ~ 6 x 10-3 μm/s was observed. We present an analytical model that bridges the two regimes and fit this model to vesicle trajectories from control cells and cells with disrupted microtubule or actin filaments. Both treatments reduced the anomalous diffusion constant and the velocity by ~40-50%.
The mobility of sterol-containing vesicles on the short time scale could reflect dynamic rearrangements of the cytoskeleton, while directed transport of sterol vesicles occurs likely along both, microtubules and actin filaments. Spatially varying anomalous diffusion could contribute to fine-tuning and local regulation of intracellular sterol transport.
Cholesterol; Transport; Fluorescence microscopy; Endocytosis; Vesicle; Tracking; Cytoskeleton dynamics
Mollusc shells are commonly investigated using high-resolution imaging techniques based on
cryo-fixation. Less detailed information is available regarding the light-optical properties. Sea shells of Haliotis pulcherina were embedded for polishing in defined orientations in order to investigate the interface between prismatic calcite and nacreous aragonite by standard materialographic methods. A polished thin section of the interface was prepared with a defined thickness of 60 μm for quantitative birefringence analysis using polarized light and LC-PolScope microscopy. Scanning electron microscopy images were obtained for comparison. In order to study structural-mechanical relationships, nanoindentation experiments were performed.
Incident light microscopy revealed a super-structure in semi-transparent regions of the polished cross-section under a defined angle. This super-structure is not visible in transmitted birefringence analysis due to the blurred polarization of small nacre platelets and numerous organic interfaces. The relative orientation and homogeneity of calcite prisms was directly identified, some of them with their optical axes exactly normal to the imaging plane. Co-oriented "prism colonies" were identified by polarized light analyses. The nacreous super-structure was also visualized by secondary electron imaging under defined angles. The domains of the super-structure were interpreted to consist of crystallographically aligned platelet stacks. Nanoindentation experiments showed that mechanical properties changed with the same periodicity as the domain size.
In this study, we have demonstrated that insights into the growth mechanisms of nacre can be obtained by conventional light-optical methods. For example, we observed super-structures formed by co-oriented nacre platelets as previously identified using X-ray Photo-electron Emission Microscopy (X-PEEM) [Gilbert et al., Journal of the American Chemical Society 2008, 130:17519–17527]. Polarized optical microscopy revealed unprecedented super-structures in the calcitic shell part. This bears, in principle, the potential for in vivo studies, which might be useful for investigating the growth modes of nacre and other shell types.
Advances in computational biology allow systematic investigations to ascertain whether internal chemical gradients can be maintained in bacteria – an open question at the resolution limit of fluorescence microscopy. While it was previously believed that the small bacterial cell size and fast diffusion in the cytoplasm effectively remove any such gradient, a new computational study published in BMC Biophysics supports the emerging view that gradients can exist. The study arose from the recent observation that phosphorylated CtrA forms a gradient prior to cell division in Caulobacter crescentus, a bacterium known for its complicated cell cycle. Tropini et al. (2012) postulate that such gradients can provide an internal chemical compass, directing protein localization, cell division and cell development. More specifically, they describe biochemical and physical constraints on the formation of such gradients and explore a number of existing bacterial cell morphologies. These chemical gradients may limit in vitro analyses, and may ensure timing control and robustness to fluctuations during critical stages in cell development.
Bacteria dynamically regulate their intricate intracellular organization involving proteins that facilitate cell division, motility, and numerous other processes. Consistent with this sophisticated organization, bacteria are able to create asymmetries and spatial gradients of proteins by localizing signaling pathway components. We use mathematical modeling to investigate the biochemical and physical constraints on the generation of intracellular gradients by the asymmetric localization of a source and a sink.
We present a systematic computational analysis of the effects of other regulatory mechanisms, such as synthesis, degradation, saturation, and cell growth. We also demonstrate that gradients can be established in a variety of bacterial morphologies such as rods, crescents, spheres, branched and constricted cells.
Taken together, these results suggest that gradients are a robust and potentially common mechanism for providing intracellular spatial cues.
Minor changes in protein structure induced by small organic and inorganic molecules can result in significant metabolic effects. The effects can be even more profound if the molecular players are chemically active and present in the cell in considerable amounts. The aim of our study was to investigate effects of a nitric oxide donor (spermine NONOate), ATP and sodium/potassium environment on the dynamics of thermal unfolding of human hemoglobin (Hb). The effect of these molecules was examined by means of circular dichroism spectrometry (CD) in the temperature range between 25°C and 70°C. The alpha-helical content of buffered hemoglobin samples (0.1 mg/ml) was estimated via ellipticity change measurements at a heating rate of 1°C/min.
Major results were: 1) spermine NONOate persistently decreased the hemoglobin unfolding temperature Tuirrespectively of the Na + /K + environment, 2) ATP instead increased the unfolding temperature by 3°C in both sodium-based and potassium-based buffers and 3) mutual effects of ATP and NO were strongly influenced by particular buffer ionic compositions. Moreover, the presence of potassium facilitated a partial unfolding of alpha-helical structures even at room temperature.
The obtained data might shed more light on molecular mechanisms and biophysics involved in the regulation of protein activity by small solutes in the cell.
EF-hand proteins can be activated by the binding of various heavy metals other than calcium, and such complexes can disturb the calcium-signaling pathway and cause toxicity and disease causing state. So far, no comprehensive study has been done to understand different heavy metals binding to calcium signaling proteins.
In this work, the flexibility of the EF-hand motifs are examined by crystallographic and thermodynamic studies of binding of Pb2+, Ba2+ and Sr2+ to Calcium Binding Protein-1 from Entamoeba histolytica (EhCaBP1). The structures of the EhCaBP1- heavy metal complexes are found to be overall similar, nevertheless specific differences in metal coordination, and small differences in the coordination distances between the metal and the ligands in the metal binding loop. The largest such distances occur for the Ba2+- EhCaBP1 complex, where two bariums are bound with partial occupancy at the EF2 motif. Thermodynamic studies confirm that EhCaBP1 has five binding sites for Ba2+ compared to four binding sites for the other metals. These structures and thermodynamic studies reveal that the EF-hand motifs can accommodate several heavy atoms with similar binding affinities. The binding of Ca2+ to the 1st, 2nd and 4th sites and the binding of Ba2+ to the 1st, 2nd, 4th and 5th sites are both enthalpically and entropically driven, whereas the binding of Sr2+ to the 1st, 2nd and 4th sites are simply enthalpy driven, interestingly in agreement with ITC data, Sr2+ do not coordinate with water in this structure. For all the metals, binding to the 3rd site is only entropy driven.
Energetically, Ca2+ is preferred in three sites, while in one site Ba2+ has better binding energy. The Sr2+-coordination in the EF hand motifs is similar to that of the native Ca2+ bound structure, except for the lack of water coordination. Sr2+ coordination seems to be a pre-formed in nature since all seven coordinating atoms are from the protein itself, which also correlates with entropy contributions in Sr2+ binding. These findings improve our understanding of metal association with calcium binding proteins and of metal induced conformational changes.
Calcium sensor; Calcium binding protein; Coordination geometry; EF-hand motifs; Anthropogenic toxicant; Domain swapped manner; Anomalous signal
Molecular Dynamics (MD) simulations are a promising tool to generate molecular understanding of processes related to the purification of proteins. Polyethylene glycols (PEG) of various length are commonly used in the production and purification of proteins. The molecular mechanisms behind PEG driven precipitation, aqueous two-phase formation or the effects of PEGylation are however still poorly understood.
In this paper, we ran MD simulations of single PEG molecules of variable length in explicitly simulated water. The resulting structures are in good agreement with experimentally determined 3D structures of PEG. The increase in surface hydrophobicity of PEG of longer chain length could be explained on an atomic scale. PEG-water interactions as well as aqueous two-phase formation in the presence of PO4 were found to be correlated to PEG surface hydrophobicity.
We were able to show that the taken MD simulation approach is capable of generating both structural data as well as molecule descriptors in agreement with experimental data. Thus, we are confident of having a good in silico representation of PEG.
G protein coupled receptors (GPCRs) are seven helical transmembrane proteins that function as signal transducers. They bind ligands in their extracellular and transmembrane regions and activate cognate G proteins at their intracellular surface at the other side of the membrane. The relay of allosteric communication between the ligand binding site and the distant G protein binding site is poorly understood. In this study, GREMLIN
, a recently developed method that identifies networks of co-evolving residues from multiple sequence alignments, was used to identify those that may be involved in communicating the activation signal across the membrane. The GREMLIN-predicted long-range interactions between amino acids were analyzed with respect to the seven GPCR structures that have been crystallized at the time this study was undertaken.
GREMLIN significantly enriches the edges containing residues that are part of the ligand binding pocket, when compared to a control distribution of edges drawn from a random graph. An analysis of these edges reveals a minimal GPCR binding pocket containing four residues (T1183.33, M2075.42, Y2686.51 and A2927.39). Additionally, of the ten residues predicted to have the most long-range interactions (A1173.32, A2726.55, E1133.28, H2115.46, S186EC2, A2927.39, E1223.37, G902.57, G1143.29 and M2075.42), nine are part of the ligand binding pocket.
We demonstrate the use of GREMLIN to reveal a network of statistically correlated and functionally important residues in class A GPCRs. GREMLIN identified that ligand binding pocket residues are extensively correlated with distal residues. An analysis of the GREMLIN edges across multiple structures suggests that there may be a minimal binding pocket common to the seven known GPCRs. Further, the activation of rhodopsin involves these long-range interactions between extracellular and intracellular domain residues mediated by the retinal domain.
GPCR; GREMLIN; Long-range interactions; Ligand binding pocket; Graphical model
Human carbonic anhydrases (CAs) play crucial role in various physiological processes including carbon dioxide and hydrocarbon transport, acid homeostasis, biosynthetic reactions, and various pathological processes, especially tumor progression. Therefore, CAs are interesting targets for pharmaceutical research. The structure-activity relationships (SAR) of designed inhibitors require detailed thermodynamic and structural characterization of the binding reaction. Unfortunately, most publications list only the observed thermodynamic parameters that are significantly different from the intrinsic parameters. However, only intrinsic parameters could be used in the rational design and SAR of the novel compounds.
Intrinsic binding parameters for several inhibitors, including ethoxzolamide, trifluoromethanesulfonamide, and acetazolamide, binding to recombinant human CA XIII isozyme were determined. The parameters were the intrinsic Gibbs free energy, enthalpy, entropy, and the heat capacity. They were determined by titration calorimetry and thermal shift assay in a wide pH and temperature range to dissect all linked protonation reaction contributions.
Precise determination of the inhibitor binding thermodynamics enabled correct intrinsic affinity and enthalpy ranking of the compounds and provided the means for SAR analysis of other rationally designed CA inhibitors.
During elongation, multi-subunit RNA polymerases (RNAPs) cycle between phosphodiester bond formation and nucleic acid translocation. In the conformation associated with catalysis, the mobile “trigger loop” of the catalytic subunit closes on the nucleoside triphosphate (NTP) substrate. Closing of the trigger loop is expected to exclude water from the active site, and dehydration may contribute to catalysis and fidelity. In the absence of a NTP substrate in the active site, the trigger loop opens, which may enable translocation. Another notable structural element of the RNAP catalytic center is the “bridge helix” that separates the active site from downstream DNA. The bridge helix may participate in translocation by bending against the RNA/DNA hybrid to induce RNAP forward movement and to vacate the active site for the next NTP loading. The transition between catalytic and translocation conformations of RNAP is not evident from static crystallographic snapshots in which macromolecular motions may be restrained by crystal packing.
All atom molecular dynamics simulations of Thermus thermophilus (Tt) RNAP reveal flexible hinges, located within the two helices at the base of the trigger loop, and two glycine hinges clustered near the N-terminal end of the bridge helix. As simulation progresses, these hinges adopt distinct conformations in the closed and open trigger loop structures. A number of residues (described as “switch” residues) trade atomic contacts (ion pairs or hydrogen bonds) in response to changes in hinge orientation. In vivo phenotypes and in vitro activities rendered by mutations in the hinge and switch residues in Saccharomyces cerevisiae (Sc) RNAP II support the importance of conformational changes predicted from simulations in catalysis and translocation. During simulation, the elongation complex with an open trigger loop spontaneously translocates forward relative to the elongation complex with a closed trigger loop.
Switching between catalytic and translocating RNAP forms involves closing and opening of the trigger loop and long-range conformational changes in the atomic contacts of amino acid side chains, some located at a considerable distance from the trigger loop and active site. Trigger loop closing appears to support chemistry and the fidelity of RNA synthesis. Trigger loop opening and limited bridge helix bending appears to promote forward nucleic acid translocation.