Search

Selection of pharmacological agents based on potency measurements performed at equilibrium fail to incorporate the kinetic aspects of the drug–target interaction. Here we describe a method for screening or characterization of enzyme inhibitors that allows the concomitant determination of the equilibrium inhibition constant in unison with rates of complex formation and dissociation. The assay is distinct from conventional enzymatic assays and is based on the analysis of inhibition curves recorded prior to full equilibration of the system. The methodology is illustrated using bicyclic peptide inhibitors of the serine protease plasma kallikrein.

Bio-catalysis is the outcome of a subtle interplay between internal motions in enzymes and chemical kinetics. Small-angle X-ray scattering (SAXS) investigation of an enzyme’s internal motions during catalysis offers an integral view of the protein’s structural plasticity, dynamics, and function, which is useful for understanding allosteric effects and developing novel medicines. Guanylate kinase (GMPK) is an essential enzyme involved in the guanine nucleotide metabolism of unicellular and multicellular organisms. It is also required for the intracellular activation of numerous antiviral and anticancer purine nucleoside analog prodrugs. Catalytically active recombinant human GMPK (hGMPK) was purified for the first time and changes in the size and shape of open/closed hGMPK were tracked by SAXS. The binding of substrates (GMP + AMPPNP or Ap5G or GMP + ADP) resulted in the compaction of size and shape of hGMPK. The structural changes between open and completely closed hGMPK conformation were confirmed by observing differences in the hGMPK secondary structures with circular dichroism spectroscopy. Graphical abstract Electronic supplementary material The online version of this article (doi:10.1007/s00249-015-1079-9) contains supplementary material, which is available to authorized users.

One of the goals of mechanosensitive channel (MSC) studies is to understand the underlying molecular and biophysical mechanisms of the mechano-gating process from force sensing to gate opening. We focus on the latter process and investigate the role of water in the bacterial MSC MscL, which is activated by membrane tension. We analyze the interplay between water and the gate-constituting amino acids, Leu19-Gly26, through molecular dynamics simulations. To highlight the role of water, specifically hydration of the gate, in MscL gating, we restrain lateral movements of the water molecules along the water-vapor interfaces at the top and bottom of the vapor bubble, plugging the closed gate. The gating behaviors in this model and the normal MscL model, in which water movements are unrestrained, are compared. In the normal model, increased membrane tension breaks the hydrogen bond between Leu19 and Val 23 of the inner helix, exposing the backbone carbonyl oxygen of Leu19 to the water-accessible lumen side of the gate. Associated with this activity, water comes to access the vapor region and stably interacts with the carbonyl oxygen to induce a dewetting to wetting transition that facilitates gate expansion toward channel opening. By contrast, in the water-restrained model, carbonyl oxygen is also exposed, but no further conformational changes occur at the gate. This suggests that gate opening relies on a conformational change initiated by wetting. The penetrated water weakens the hydrophobic interaction between neighboring transmembrane inner helices called the "hydrophobic lock" by wedging into the space between their interacting portions.

The addition of trace molecules into membranes can significantly alter the morphology of the co-existing liquid phases and lipid phase transition temperature. Membrane additives may affect lipid phase dynamics through preferentially partitioning to the boundary between lipid phases or preferentially mixing into one lipid phase. The characteristic differences between these mechanisms are demonstrated here in a minimalistic nearest neighbor model to provide a framework for how slight changes to membrane composition may affect lipid-phase-dependent processes, such as lipid-raft formation, immunological signaling, and molecular sorting preceding endocytosis with coexisting liquid phases. Within the low mole fractions explored here (≤3 mol%), increasing the additive concentration linearly changed the phase miscibility temperature. Rotationally asymmetric Janus particles reduced the miscibility transition temperature for all fractions and degree of phase polarization. Rotationally symmetric additives, however, either increased or decreased the phase miscibility temperature depending on the phase preference of the additive. While most experimental molecules may contain aspects of both of these idealized additives, this model provides a broad framework to quantify the effects of membrane additives in regard to lipid phase preference, lipid-raft association, and contribution to lipid phase-dependent molecular sorting. Electronic supplementary material The online version of this article (doi:10.1007/s00249-015-1017-x) contains supplementary material, which is available to authorized users.

Membrane-active peptides include peptides that can cross cellular membranes and deliver macromolecular cargo as well as peptides that inhibit bacterial growth. Some of these peptides can act as both transporters and antibacterial agents. It is desirable to combine the knowledge from these two different fields of membrane-active peptides into design of new peptides with tailored actions, as transporters of cargo or as antibacterial substances, targeting specific membranes. We have previously shown that the position of the amino acid tryptophan in the peptide sequence of three arginine-tryptophan peptides affects their uptake and intracellular localization in live mammalian cells, as well as their ability to inhibit bacterial growth. Here, we use quartz crystal microbalance with dissipation monitoring to assess the induced changes caused by binding of the three peptides to supported model membranes composed of POPC, POPC/POPG, POPC/POPG/cholesterol or POPC/lactosyl PE. Our results indicate that the tryptophan position in the peptide sequence affects the way these peptides interact with the different model membranes and that the presence of cholesterol in particular seems to affect the membrane interaction of the peptide with an even distribution of tryptophans in the peptide sequence. These results give mechanistic insight into the function of these peptides and may aid in the design of membrane-active peptides with specified cellular targets and actions. Electronic supplementary material The online version of this article (doi:10.1007/s00249-014-0958-9) contains supplementary material, which is available to authorized users.

MultiSig is a newly developed mode of analysis of sedimentation equilibrium (SE) experiments in the analytical ultracentrifuge, having the capability of taking advantage of the remarkable precision (~0.1 % of signal) of the principal optical (fringe) system employed, thus supplanting existing methods of analysis through reducing the ‘noise’ level of certain important parameter estimates by up to orders of magnitude. Long-known limitations of the SE method, arising from lack of knowledge of the true fringe number in fringe optics and from the use of unstable numerical algorithms such as numerical differentiation, have been transcended. An approach to data analysis, akin to ‘spatial filtering’, has been developed, and shown by both simulation and practical application to be a powerful aid to the precision with which near-monodisperse systems can be analysed, potentially yielding information on protein-solvent interaction. For oligo- and poly-disperse systems the information returned includes precise average mass distributions over both cell radial and concentration ranges and mass-frequency histograms at fixed radial positions. The application of MultiSig analysis to various complex heterogenous systems and potentially multiply-interacting carbohydrate oligomers is described.

Formins are highly conserved proteins that are essential in the formation and regulation of the actin cytoskeleton. The formin homology 2 (FH2) domain is responsible for actin binding and acts as an important nucleating factor in eukaryotic cells. In this work EPR and DSC were used to investigate the properties of the mDia1-FH2 formin fragment and its interaction with actin. MDia1-FH2 was labeled with a maleimide spin probe (MSL). EPR results suggested that the MSL was attached to a single SH group in the FH2. In DSC and temperature-dependent EPR experiments we observed that mDia1-FH2 has a flexible structure and observed a major temperature-induced conformational change at 41 °C. The results also confirmed the previous observation obtained by fluorescence methods that formin binding can destabilize the structure of actin filaments. In the EPR experiments the intermolecular connection between the monomers of formin dimers proved to be flexible. Considering the complex molecular mechanisms underlying the cellular roles of formins this internal flexibility of the dimers is probably important for manifestation of their biological functions. Electronic supplementary material The online version of this article (doi:10.1007/s00249-013-0922-0) contains supplementary material, which is available to authorized users.

The regulation of striated muscle contraction involves cooperative interactions between actin filaments, myosin-S1 (S1), tropomyosin (Tm), troponin (Tn), and calcium. These interactions are modeled by treating overlapping tropomyosins as a continuous flexible chain (CFC), weakly confined by electrostatic interactions with actin. The CFC is displaced locally in opposite directions on the actin surface by the binding of either S1 or Troponin I (TnI) to actin. The apparent rate constants for myosin and TnI binding to and detachment from actin are then intrinsically coupled via the CFC model to the presence of neighboring bound S1s and TnIs. Monte Carlo simulations at prescribed values of the CFC stiffness, the CFC’s degree of azimuthal confinement, and the angular displacements caused by the bound proteins were able to predict the stopped-flow transients of S1 binding to regulated F-actin. The transients collected over a large range of calcium concentrations could be well described by adjusting a single calcium-dependent parameter, the rate constant of TnI detachment from actin, k−I. The resulting equilibrium constant \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ K_{\text{B}} \equiv 1/K_{\text{I}} $$\end{document} varied sigmoidally with the free calcium, increasing from 0.12 at low calcium (pCa >7) to 12 at high calcium (pCa

Cytochrome P450 2D6 is a major drug-metabolising enzyme with a wide substrate range. A single-point mutation introduced in this enzyme induces stereoselective binding of R and S-propranolol whereas the wild type has no preference. The system has previously been studied both experimentally and computationally (de Graaf et al. in Eur Biophys J 36:589–599, 2007a). The in silico study reported hysteresis and significant deviations from closure of thermodynamic cycles, probably because of lack of sampling. Here, we focus on the effect of prolonged simulation time and enhanced sampling methods, such as Hamiltonian replica exchange, to reduce these problems and to improve the precision of free energy calculations. Finally we rationalize the results at a molecular level and compare data with experimental findings and previously estimated free energies. Graphical Abstract Propranolol (PPD) binds stereospecifically to cytochrome P450 2D6, if the mutation F483A is introduced, by reducing binding affinity to R-PPD. This observation has previously been studied by molecular dynamics, by use of several different thermodynamic cycles. Previous computational results suffered from internal inconsistencies, probably because of insufficient sampling. Here we use prolonged simulation time and Hamiltonian replica exchange to overcome the inconsistencies, and give a possible rationalization of the phenomenon at a molecular level.

Rat prostate cancer cells have been previously investigated using two cell lines: a highly metastatic one (Mat-Ly-Lu) and a nonmetastatic one (AT-2). It turns out that the highly metastatic Mat-Ly-Lu cells exhibit a phenomenon of cathodal galvanotaxis in an electric field which can be blocked by interrupting the voltage-gated sodium channel (VGSC) activity. The VGSC activity is postulated to be characteristic for metastatic cells and seems to be a reasonable driving force for motile behavior. However, the classical theory of cellular motion depends on calcium ions rather than sodium ions. The current research provides a theoretical connection between cellular sodium inflow and cathodal galvanotaxis of Mat-Ly-Lu cells. Electrical repulsion of intracellular calcium ions by entering sodium ions is proposed after depolarization starting from the cathodal side. The disturbance in the calcium distribution may then drive actin polymerization and myosin contraction. The presented modeling is done within a continuous one-dimensional Poisson-Nernst-Planck equation framework.

Short cationic antimicrobial peptides (AMPs) are believed to act either by inducing transmembrane pores or disrupting membranes in a detergent-like manner. For example, the antimicrobial peptides aurein 1.2, citropin 1.1, maculatin 1.1 and caerin 1.1, despite being closely related, appear to act by fundamentally different mechanisms depending on their length. Using molecular dynamics simulations, the structural properties of these four peptides have been examined in solution as well as in a variety of membrane environments. It is shown that each of the peptides has a strong preference for binding to regions of high membrane curvature and that the structure of the peptides is dependent on the degree of local curvature. This suggests that the shorter peptides aurein 1.2 and citropin 1.1 act via a detergent-like mechanism because they can induce high local, but not long-range curvature, whereas the longer peptides maculatin 1.1 and caerin 1.1 require longer range curvature to fold and thus bind to and stabilize transmembrane pores.

Parameterized models of biophysical and mechanical cell properties are important for predictive mathematical modeling of cellular processes. The concepts of turgor, cell wall elasticity, osmotically active volume, and intracellular osmolarity have been investigated for decades, but a consistent rigorous parameterization of these concepts is lacking. Here, we subjected several data sets of minimum volume measurements in yeast obtained after hyper-osmotic shock to a thermodynamic modeling framework. We estimated parameters for several relevant biophysical cell properties and tested alternative hypotheses about these concepts using a model discrimination approach. In accordance with previous reports, we estimated an average initial turgor of 0.6 ± 0.2 MPa and found that turgor becomes negligible at a relative volume of 93.3 ± 6.3% corresponding to an osmotic shock of 0.4 ± 0.2 Osm/l. At high stress levels (4 Osm/l), plasmolysis may occur. We found that the volumetric elastic modulus, a measure of cell wall elasticity, is 14.3 ± 10.4 MPa. Our model discrimination analysis suggests that other thermodynamic quantities affecting the intracellular water potential, for example the matrix potential, can be neglected under physiological conditions. The parameterized turgor models showed that activation of the osmosensing high osmolarity glycerol (HOG) signaling pathway correlates with turgor loss in a 1:1 relationship. This finding suggests that mechanical properties of the membrane trigger HOG pathway activation, which can be represented and quantitatively modeled by turgor. This work was supported via several projects funded by the European Commission: QUASI (Contract No. 503230 to SH, EK, FP and MP), CELLCOMPUT (Contract No. 043310 to SH, EK and FP), UNICELLSYS (Contract No. 201142 to SH, EK, FP, MP and MG), SYSTEMSBIOLOGY (Contract No. 514169 to SH and EK), and AMPKIN (Contract No. 518181 to SH and MG). In addition work was funded by grants from the Swedish Foundation for Strategic Research SSF (Bio-X to MG), the Swedish Research Council (project grants to SH and MG), the Carl Trygger Foundation (to MG), the Science Faculty, University of Gothenburg (to SH and MG), and the Swiss systemsX.ch (to MP)

The FtsZ protein is a self-polymerizing GTPase that plays a central role in bacterial cell division. Several C8-substituted GTP analogs are known to inhibit the polymerization of FtsZ by competing for the same binding site as its endogenous activating ligand GTP. Free energy calculations of the relative binding affinities to FtsZ for a set of five C8-substituted GTP analogs were performed. The calculated values agree well with the available experimental data, and the main contribution to the free energy differences is determined to be the conformational restriction of the ligands. The dihedral angle distributions around the glycosidic bond of these compounds in water are known to vary considerably depending on the physicochemical properties of the substituent at C8. However, within the FtsZ protein, this substitution has a negligible influence on the dihedral angle distributions, which fall within the narrow range of −140° to −90° for all investigated compounds. The corresponding ensemble average of the coupling constants 3 J(C4,H1′) is calculated to be 2.95 ± 0.1 Hz. The contribution of the conformational selection of the GTP analogs upon binding was quantified from the corresponding populations. The obtained restraining free energy values follow the same trend as the relative binding affinities to FtsZ, indicating their dominant contribution.

The p7 protein of hepatitis C virus (HCV) plays an important role in the viral lifecycle. Like other members of the viroporin family of small membrane proteins, the amino acid sequence of p7 is largely conserved over the entire range of genotypes, and it forms ion channels that can be blocked by a number of established channel-blocking compounds. Its characteristics as a membrane protein make it difficult to study by most structural techniques, since it requires the presence of lipids to fold and function properly. Purified p7 can be incorporated into phospholipid bilayers and micelles. Initial solid-state nuclear magnetic resonance (NMR) studies of p7 in 14-O-PC/6-O-PC bicelles indicate that the protein contains helical segments that are tilted approximately 10 degrees and 25 degrees relative to the bilayer normal. A truncated construct corresponding to the second transmembrane domain of p7 is shown to have properties similar to those of the full-length protein, and was used to determine that the helix segment tilted at 10 degrees is in the C-terminal portion of the protein. The addition of the channel blocker amantadine to the full-length protein resulted in selective chemical shift changes, demonstrating that NMR has a potential role in the development of drugs targeted to p7.

The conformation of a transmembrane peptide, sMTM7, encompassing the cytoplasmic hemi-channel domain of the seventh transmembrane section of subunit a from V-ATPase from Saccharomyces cerevisiae solubilized in SDS solutions was studied by circular dichroism (CD) spectroscopy and fluorescence spectroscopy of the single tryptophan residue of this peptide. The results show that the peptide adopts an alpha-helical conformation or aggregated beta-sheet depending on the peptide-to-SDS ratio used. The results are compared with published data about a longer version of the peptide (i.e., MTM7). It is concluded that the bulky, positively charged arginine residue located in the center of both peptides has a destabilizing effect on the helical conformation of the SDS-solubilized peptides, leading to beta-sheet formation and subsequent aggregation.

The sequential organization of genomes, i.e. the relations between distant base pairs and regions within sequences, and its connection to the three-dimensional organization of genomes is still a largely unresolved problem. Long-range power-law correlations were found using correlation analysis on almost the entire observable scale of 132 completely sequenced chromosomes of 0.5 x 10(6) to 3.0 x 10(7) bp from Archaea, Bacteria, Arabidopsis thaliana, Saccharomyces cerevisiae, Schizosaccharomyces pombe, Drosophila melanogaster, and Homo sapiens. The local correlation coefficients show a species-specific multi-scaling behaviour: close to random correlations on the scale of a few base pairs, a first maximum from 40 to 3,400 bp (for Arabidopsis thaliana and Drosophila melanogaster divided in two submaxima), and often a region of one or more second maxima from 10(5) to 3 x 10(5) bp. Within this multi-scaling behaviour, an additional fine-structure is present and attributable to codon usage in all except the human sequences, where it is related to nucleosomal binding. Computer-generated random sequences assuming a block organization of genomes, the codon usage, and nucleosomal binding explain these results. Mutation by sequence reshuffling destroyed all correlations. Thus, the stability of correlations seems to be evolutionarily tightly controlled and connected to the spatial genome organization, especially on large scales. In summary, genomes show a complex sequential organization related closely to their three-dimensional organization.

Localized induction of DNA damage is a valuable tool for studying cellular DNA damage responses. In recent decades, methods have been developed to generate DNA damage using radiation of various types, including photons and charged particles. Here we describe a simple ultrasoft X-ray multi-microbeam system for high dose-rate, localized induction of DNA strand breaks in cells at spatially and geometrically adjustable sites. Our system can be combined with fixed- and live-cell microscopy to study responses of cells to DNA damage.

Three well-known representatives of the cyclodextrin family were completely characterized by molecular hydrodynamics methods in three different solvents. For the first time the possibility of an estimation of velocity sedimentation coefficients s between 0.15 and 0.5 S by the numerical solution of the Lamm equation is shown. Comparison of the experimental hydrodynamic characteristics of the cyclodextrins with theoretical calculations for toroidal molecules allows an estimation of the thickness of the solvent layers on the surface of cyclodextrin molecules.

The function and dynamics of proteins depend on their direct environment, and much evidence has pointed to a strong coupling between water and protein motions. Recently however, neutron scattering measurements on deuterated and natural-abundance purple membrane (PM), hydrated in H2O and D2O, respectively, revealed that membrane and water motions on the ns–ps time scale are not directly coupled below 260 K (Wood et al. in Proc Natl Acad Sci USA 104:18049–18054, 2007). In the initial study, samples with a high level of hydration were measured. Here, we have measured the dynamics of PM and water separately, at a low-hydration level corresponding to the first layer of hydration water only. As in the case of the higher hydration samples previously studied, the dynamics of PM and water display different temperature dependencies, with a transition in the hydration water at 200 K not triggering a transition in the membrane at the same temperature. Furthermore, neutron diffraction experiments were carried out to monitor the lamellar spacing of a flash-cooled deuterated PM stack hydrated in H2O as a function of temperature. At 200 K, a sudden decrease in lamellar spacing indicated the onset of long-range translational water diffusion in the second hydration layer as has already been observed on flash-cooled natural-abundance PM stacks hydrated in D2O (Weik et al. in J Mol Biol 275:632–634, 2005), excluding thus a notable isotope effect. Our results reinforce the notion that membrane-protein dynamics may be less strongly coupled to hydration water motions than the dynamics of soluble proteins.

The molecular understanding of protein stabilization by the disaccharide trehalose in extreme temperature or hydration conditions is still debated. In the present study, we investigated the role of trehalose on the dynamics of the protein C-phycocyanin (C-PC) by neutron scattering. To single out the motions of C-PC hydrogen (H) atoms in various trehalose/water environments, measurements were performed in deuterated trehalose and heavy water (D2O). We report that trehalose decreases the internal C-PC dynamics, as shown by a reduced diffusion coefficient of protein H atoms. By fitting the Elastic Incoherent Structure Factor--which gives access to the "geometry" of the internal proton motions--with the model of diffusion inside a sphere, we found that the presence of trehalose induces a significantly higher proportion of immobile C-PC hydrogens. We investigated, by elastic neutron scattering, the mean square displacements (MSDs) of deuterated trehalose/D2O-embedded C-PC as a function of temperature in the range of 40-318 K. Between 40 and approximately 225 K, harmonic MSDs of C-PC are slightly smaller in samples containing trehalose. Above a transition temperature of approximately 225 K, we observed anharmonic motions in all trehalose/water-coated C-PC samples. In the hydrated samples, MSDs are not significantly changed by addition of 15% trehalose but are slightly reduced by 30% trehalose. In opposition, no dynamical transition was detected in dry trehalose-embedded C-PC, whose hydrogen motions remain harmonic up to 318 K. These results suggest that a role of trehalose would be to stabilize proteins by inhibiting some fluctuations at the origin of protein unfolding and denaturation.

Sugars play an important role in the desiccation tolerance of most anhydrobiotic organisms. It has been shown in previous studies that different structural families of oligosaccharides have different efficacies to interact with phospholipid headgroups and protect membranes from solute leakage during drying. Here, we have compared three families of linear oligosaccharides (fructans (inulins), malto-oligosaccharides, manno-oligosaccharides) for their chain-length dependent protection of egg phosphatidylcholine liposomes against membrane fusion. We found increased protection with chain length up to a degree of polymerization (DP) of 5 for malto-oligosaccharides, and a decrease for inulins and manno-oligosaccharides. Differential scanning calorimetry measurements showed that for all sugars the glass transition temperature (T(g)) increased with DP, although to different degrees for the different oligosaccharide families. Higher T(g) values resulted in reduced membrane fusion only for malto-oligosaccharides below DP5. Contrary to expectation, for inulins, manno-oligosaccharides and malto-oligosaccharides of a DP above five, fusion increased with increasing T(g), indicating that other physical parameters are more important in determining the ability of different sugars to protect membranes against fusion during drying. Further research will be necessary to experimentally define such parameters.

Several approaches to ion-channel gating modelling have been proposed. Although many models describe the dwell-time distributions correctly, they are incapable of predicting and explaining the long-term correlations between the lengths of adjacent openings and closings of a channel. In this paper we propose two simple random-walk models of the gating dynamics of voltage and Ca(2+)-activated potassium channels which qualitatively reproduce the dwell-time distributions, and describe the experimentally observed long-term memory quite well. Biological interpretation of both models is presented. In particular, the origin of the correlations is associated with fluctuations of channel mass density. The long-term memory effect, as measured by Hurst R/S analysis of experimental single-channel patch-clamp recordings, is close to the behaviour predicted by our models. The flexibility of the models enables their use as templates for other types of ion channel.

The paper presents experimental investigations of diffusion of antibiotics (ciprofloxacin or ampicillin) into the water phase from mixtures of neutral or negatively charged liposomes, and antibiotic–liposome interactions. Using the laser interferometry technique, the amounts and fluxes of released antibiotics, concentration field evolution, and the velocity of the concentration boundary layer’s “growth” were determined. To avoid the limitations of membranes, a measurement system without the artificial boundary of phases with a free water–solution interface has been proposed. It was found that the diffusion of anionic and neutral liposomes into the water phase was insignificant and mainly the diffusion of antibiotics was measured. Differences in the diffusion kinetics of ciprofloxacin and ampicillin from liposomal solutions to the water phase were observed. Ampicillin diffused more efficiently than ciprofloxacin regardless of the liposomal solution type. Moreover, the amount of ampicillin and ciprofloxacin released from the anionic liposomal phase was higher than that from the neutral one. Our results confirm that ciprofloxacin at neutral pH shows little tendency to bind neutral liposomes. Additionally, it was also observed that ciprofloxacin disrupts negatively charged liposomes as a final effect of antibiotic–lipid interactions.