Double-exponential kinetics of binding and redistribution of the fluorescent dyes in cell membranes witness for the existence of lipid microdomains.

New technique of detecting lateral heterogeneity of the plasma membrane of living cells by means of membrane-binding fluorescent dyes is proposed. The kinetics of dye incorporation into the membrane or its lateral diffusion inside the membrane is measured and decomposed into exponential components by means of the Maximum Entropy Method. Two distinct exponential components are obtained consistently in all cases for several fluorescent dyes, two different cell lines and in different types of experiments including spectroscopy, flow cytometry and fluorescence recovery after photobleaching. These components are attributed to the liquid-ordered and disordered phases in the plasma membrane of studied cells in their dynamic equilibrium.

Low-Density Lipoproteins and Human Serum Albumin as Carriers of Squalenoylated Drugs: Insights from Molecular Simulations.

We have studied the interaction of three clinically promising squalenoylated drugs (gemcitabine-squalene, adenine-squalene, and doxorubicin-squalene) with low-density lipoproteins (LDL) by means of atomistic molecular dynamics simulations. It is shown that all studied squalenoylated drugs accumulate inside the LDL particles. This effect is promoted by the squalene moiety, which acts as an anchor and drives the hydrophilic drugs into the hydrophobic core of the LDL lipid droplet. Our data suggest that LDL particles could be a universal carriers of squalenoylated drugs in the bloodstream. Interaction of gemcitabine-squalene with human serum albumin (HSA) was also studied by ensemble of docking simulations. It is shown that HSA could also act as a passive carrier of this bioconjugate. It should be noted that the binding of squalene moiety to HSA was unspecific and did not occur in the binding pockets devoted to fatty acids.

Caspase-3 activation decreases lipid order in the outer plasma membrane leaflet during apoptosis: A fluorescent probe study.

In this research we investigate the connection between the cytoplasmic machinery of apoptosis and the plasma membrane organization by studying the coupling of caspase-3 activation and inhibition with PS exposure and the change of lipid order in plasma membrane sensed by a fluorescent membrane probe NR12S. First, we performed in silico molecular dynamics simulations, which suggest that the mechanism of response of NR12S to lipid order may combine both sensitivity to membrane polarity/hydration and change in the fluorophore orientation. Second, cellular studies revealed that upon triggering apoptosis with IPA-3 and camptothecin the NR12S response is similar to that observed after decrease of lipid order induced by cholesterol depletion, 7-ketocholesterol enrichment or sphingomyelin hydrolysis. NR12S response can be influenced by a caspase-3 inhibitor Z-DEVD-FMK. Flow cytometry data further indicate that the NR12S response correlates with the response of FITC-labeled DEVD-FMK peptide and GFP-labeled Annexin V on the whole time scale (0-24h) of apoptosis induction by camptothecin. We conclude that fine changes in lipid order observed by NR12S are coupled with early steps of cellular events in apoptosis.

Selective Inhibition of STAT3 with Respect to STAT1: Insights from Molecular Dynamics and Ensemble Docking Simulations.

STAT3 protein, which is known to be involved in cancer development, is a promising target for anticancer therapy. Successful inhibitors of STAT3 should not affect an activity of closely related protein STAT1, which makes their development challenging. The mechanisms of selectivity of several existing STAT3 inhibitors are not clear. In this work, we studied molecular mechanisms of selectivity of 13 experimentally tested STAT3 inhibitors by means of extensive molecular dynamics and ensemble docking simulations. It is shown that all studied inhibitors bind to the large part of the protein surface in an unspecific statistical manner. The binding to the dimerization interface of the SH2 domain, which is usually considered as the main target region, is not energetically preferable. Binding in this region is remarkably similar for STAT1 and STAT3 proteins and cannot explain experimentally observed selectivity toward STAT3. We propose a new mechanism of selectivity called "selectivity by distraction" for existing STAT3 inhibitors. This mechanism is based on equilibrium statistical partitioning of inhibitor molecules between protein domains. The unspecific binding of inhibitors to the DNA-binding and the coil-coil domains is stronger in STAT1 in comparison to STAT3 while the energies of their binding to SH2 domains are comparable. This "distracts" inhibitor molecules from the SH2 domain of STAT1 and leads to higher effective concentration of inhibitors in the vicinity of the SH2 domain of STAT3.

Interaction of C60 fullerenes with asymmetric and curved lipid membranes: a molecular dynamics study.

Interaction of fullerenes with asymmetric and curved DOPC/DOPS bicelles is studied by means of coarse-grained molecular dynamics simulations. The effects caused by asymmetric lipid composition of the membrane leaflets and the curvature of the membrane are analyzed. It is shown that the aggregates of fullerenes prefer to penetrate into the membrane in the regions of the moderately positive mean curvature. Upon penetration into the hydrophobic core of the membrane fullerenes avoid the regions of the extreme positive or the negative curvature. Fullerenes increase the ordering of lipid tails, which are in direct contact with them, but do not influence other lipids significantly. Our data suggest that the effects of the membrane curvature should be taken into account in the studies concerning permeability of the membranes to fullerenes and fullerene-based drug delivery systems.

Pteros 2.0: Evolution of the fast parallel molecular analysis library for C++ and python.

Pteros is the high-performance open-source library for molecular modeling and analysis of molecular dynamics trajectories. Starting from version 2.0 Pteros is available for C++ and Python programming languages with very similar interfaces. This makes it suitable for writing complex reusable programs in C++ and simple interactive scripts in Python alike. New version improves the facilities for asynchronous trajectory reading and parallel execution of analysis tasks by introducing analysis plugins which could be written in either C++ or Python in completely uniform way. The high level of abstraction provided by analysis plugins greatly simplifies prototyping and implementation of complex analysis algorithms. Pteros is available for free under Artistic License from http://sourceforge.net/projects/pteros/.

Shape-independent model (SHIM) approach for studying aggregation by NMR diffusometry.

NMR diffusometry has been gaining wide popularity in various areas of applied chemistry for investigating diffusion and complexation processes in solid and aqueous phases. To date, the application of this method to study aggregation phenomena proceeding beyond the dimer stage of assembly has been restricted by the need for a priori knowledge of the aggregates' shape, commonly difficult to know in practice. We describe here a comprehensive analysis of aggregation parameter-dependency on the type and shape selected for modeling assembly processes, and report for the first time a shape-independent model (designated the SHIM approach), which may be used as an alternative in cases when information on aggregates' shapes is unavailable. The model can be used for determining equilibrium aggregation parameters from self-diffusion NMR data including equilibrium self-association constant and changes in enthalpy, ΔH, and entropy, ΔS.

Efficiency of the monofunctionalized C60 fullerenes as membrane targeting agents studied by all-atom molecular dynamics simulations.

Transmembrane translocation of C60 fullerenes functionalized by the single amino-derivative in neutral and charged forms was studies by extensive all-atom molecular dynamics simulations. It is shown that these complexes exhibit very strong affinity to the membrane core, but their spontaneous translocation through the membrane is not possible at practical time scale. In contrast, free amino derivatives translocate through the membrane much easier than their complexes with fullerenes, but do not have pronounced affinity to the membrane interior. Our results suggest that monofunctionalized C60 could be extremely efficient membrane targeting agents, which facilitate accumulation of the water-soluble compounds in the hydrophobic core of lipid bilayer.

Asymmetric structure and domain binding interfaces of human tyrosyl-tRNA synthetase studied by molecular dynamics simulations.

Human tyrosyl-tRNA synthetase (HsTyrRS) is composed of two structural modules: N-terminal catalytic core and an EMAP II-like C-terminal domain. The structures of these modules are known, but no crystal structure of the full-length HsTyrRS is currently available. An all-atom model of the full-length HsTyrRS was developed in this work. The structure, dynamics, and domain binding interfaces of HsTyrRS were investigated by extensive molecular dynamics (MD) simulations. Our data suggest that HsTyrRS in solution consists of a number of compact asymmetric conformations, which differ significantly by their rigidity, internal mobility, orientation of C-terminal modules, and the strength of interdomain binding. Interfaces of domain binding obtained in MD simulations are in perfect agreement with our previous coarse-grained hierarchical rotations technique simulations. Formation of the hydrogen bonds between R93 residue of the ELR cytokine motif and the residues A340 and E479 in the C-module was observed. This observation supports the idea that the lack of cytokine activity in the full-length HsTyrRS is explained by interactions between N-modules and C-modules, which block the ELR motif.

Pteros: fast and easy to use open-source C++ library for molecular analysis.

An open-source Pteros library for molecular modeling and analysis of molecular dynamics trajectories for C++ programming language is introduced. Pteros provides a number of routine analysis operations ranging from reading and writing trajectory files and geometry transformations to structural alignment and computation of nonbonded interaction energies. The library features asynchronous trajectory reading and parallel execution of several analysis routines, which greatly simplifies development of computationally intensive trajectory analysis algorithms. Pteros programming interface is very simple and intuitive while the source code is well documented and easily extendible. Pteros is available for free under open-source Artistic License from http://sourceforge.net/projects/pteros/.

How cholesterol is distributed between monolayers in asymmetric lipid membranes.

The distribution of cholesterol in asymmetric lipid bilayers was studied by extensive coarse-grained molecular dynamics simulations. The effects of the lipid head group charge, acyl chain saturation, spontaneous membrane curvature and surface tension of the membrane were investigated. Four asymmetric bilayers containing DOPC, DOPS, DSPC or DSPS lipids were simulated on a time scale extended to tens of microseconds. We show that cholesterol strongly prefers anionic lipids to neutral and saturated lipid tails to unsaturated with a distribution ratio of ~0.7 in neutral/anionic bilayers and of ~0.4 in unsaturated/saturated bilayers. Multiple flip-flop transitions of cholesterol were observed directly, and their mean times ranged from 80 to 250 ns. It was shown that the distribution of cholesterol in the asymmetric membrane depends not only on the type of lipid, but also on the local membrane curvature and the surface tension. The membrane curvature enhances the influence of the lipid head groups on cholesterol distribution, while non-optimal surface tension caused by different areas per lipid in different monolayers increases the effect of the lipid tail saturation. It was clearly seen that the monolayers of asymmetric bilayers are interdependent. Mean distances from the bilayer center to cholesterol molecules depend not only on the type of the lipid in the considered monolayer but also on the composition of the opposite monolayer.

Polarizable water model for the coarse-grained MARTINI force field.

Coarse-grained (CG) simulations have become an essential tool to study a large variety of biomolecular processes, exploring temporal and spatial scales inaccessible to traditional models of atomistic resolution. One of the major simplifications of CG models is the representation of the solvent, which is either implicit or modeled explicitly as a van der Waals particle. The effect of polarization, and thus a proper screening of interactions depending on the local environment, is absent. Given the important role of water as a ubiquitous solvent in biological systems, its treatment is crucial to the properties derived from simulation studies. Here, we parameterize a polarizable coarse-grained water model to be used in combination with the CG MARTINI force field. Using a three-bead model to represent four water molecules, we show that the orientational polarizability of real water can be effectively accounted for. This has the consequence that the dielectric screening of bulk water is reproduced. At the same time, we parameterized our new water model such that bulk water density and oil/water partitioning data remain at the same level of accuracy as for the standard MARTINI force field. We apply the new model to two cases for which current CG force fields are inadequate. First, we address the transport of ions across a lipid membrane. The computed potential of mean force shows that the ions now naturally feel the change in dielectric medium when moving from the high dielectric aqueous phase toward the low dielectric membrane interior. In the second application we consider the electroporation process of both an oil slab and a lipid bilayer. The electrostatic field drives the formation of water filled pores in both cases, following a similar mechanism as seen with atomistically detailed models.

Dietary polyphenols inhibit plasma protein arabinosylation: Biomolecular interaction of genistein and ellagic acid with serum albumins.

The mode of interaction of polyphenolic compounds like genistein (GTN) and ellagic acid (EGA) with human and bovine serum albumin (HSA and BSA, respectively) was found to differ significantly. Stern-Volmer (SV) analysis of the fluorescence quenching data revealed that the binding strength of EGA (1.9 ± 0.09 × 105 M-1) to HSA is about one order of magnitude higher than GTN (2.24 ± 0.06 × 104 M-1). While the static quenching of HSA fluorescence was found to proceed through simple Stern-Volmer (SV) mechanism, a quenching sphere-of-action model was indispensable for BSA. Temperature dependent fluorescence along with a series of other biophysical experiments and ensemble docking calculation revealed that EGA and GTN bind to the serum proteins primarily through the entropy driven process. The α-helical content and the microenvironment near Trp residue of HSA and BSA did not show any appreciable change due to the binding of either GTN or EGA. Interestingly, both GTN and EGA were found to inhibit the formation of advanced glycated end (AGE) product of serum proteins up to the extent of 70-90% within 12-24 h. Relatively moderate binding propensity along with the anti-glycation ability of the polyphenols confirmed that GTN and EGA can be used either as an alternative or towards development of suitable drugs in the prevention of many diabetic-related complications.

Fluorescence Probes Exhibit Photoinduced Structural Planarization: Sensing In Vitro and In Vivo Microscopic Dynamics of Viscosity Free from Polarity Interference.

We demonstrate the construction of wavelength λ-ratiometric images that allow visualizing the distribution of microscopic dynamics within living cells and tissues by using the newly developed principle of fluorescence response. The bent-to-planar motion in the excited state of incorporated fluorescence probes leads to elongation of the π-delocalization, resulting in microviscosity-dependent but polarity-insensitive interplay between well-separated blue and red bands in emission spectra. This allows constructing the exceptionally contrasted images of cellular dynamics. Moreover, the application of probes with increased affinity toward biological membranes allowed detecting the differences in dynamics between the plasma membrane and intracellular membrane structures. Such λ-ratiometric microviscosity imaging was extended for mapping the living tissues and observing their inflammation-dependent changes.

Fuzzy domains: new way of describing flexibility and interdependence of the protein domains.

We proposed the innovative method of domain identification based on the concept of the fuzzy domains. In this method each residue of the protein can belong to several domains simultaneously with certain weights, which reflect to what extent this residue shares the motion pattern of the given domain. Our method allows describing the fuzzy boundaries between the domains and the gradual changes of the motion pattern from one domain to the other. It provides the reasonable compromise between the continuous change of the protein dynamics from one residue to the other and the discrete description of the structure in terms of small number of domains. We suggested quantitative criterion, which shows the overall degree of domain flexibility in the protein. The concept of the fuzzy domains provides an innovative way of visualization of domain flexibility, which makes the gradual transitions between the domains clearly visible and comparable to available experimental and structural data. In the future, the concept of the fuzzy domains can be used in the coarse-grained simulations of the domain dynamics in order to account for internal protein flexibility.

Therapeutic potency of substituted chromones as Alzheimer's drug: Elucidation of acetylcholinesterase inhibitory activity through spectroscopic and molecular modelling investigation.

Introduction: Documentation on the potency of chromones as acetylcholinesterase (AChE) antagonists has paved the way for the design and usage of new chromone analogues as inhibitors of AChE modelled on the hypothesis based on cholinergic pathway of Alzheimer's disease (AD). Here, 2 minimally substituted chromones, namely 3-cyanochromone (CyC) and 7-amino-3- methylchromone (AMC), were checked for their AChE inhibition efficacies and plasma protein modulation. Methods: Colorimetric enzymatic assay as well as fluorescence measurements were performed for obtaining the experimental results, which were further corroborated by molecular docking and simulation studies. Results: The investigated systems exhibited strong inhibition activities against AChE, with CyC (IC50= 85.12 ± 6.70 nM) acting as better inhibitor than AMC (IC50 = 103.09 ± 11.90 nM) and both having IC50 values in the range of FDA approved cholinergic drug Donepezil (IC50 = 74.13 ± 8.30 nM). Non-competitive inhibition was observed in both the cases with the inhibitors binding near the peripheral anionic site (PAS) of the enzyme. Having one planar nitrile group in CyC as compared to sp3 hybridised substituents in AMC facilitated stacking interactions in the former, accounting for its higher inhibitory efficacy. A significant decrease in the inhibition potency of CyC (~32%) was noted in comparison with AMC (~5%) when the experiments were performed in presence of human serum albumin (HSA) instead of pure aqueous buffer. Conclusion: This comparative study affirms the importance of meticulous substitution in the chromone scaffold to promote maximum inhibition potency, while considering their usage as AD drugs.

Novel coumarin derivatives as potent acetylcholinesterase inhibitors: insight into efficacy, mode and site of inhibition.

The inhibitory efficacy of two substituted coumarin derivatives on the activity of neurodegenerative enzyme acetylcholinesterase (AChE) was assessed in aqueous buffer as well as in the presence of human serum albumin (HSA) and compared against standard cholinergic AD drug, Donepezil (DON). The experimental data revealed the inhibition to be of non-competitive type with both the systems showing substantial inhibitory activity on AChE. In fact, one of the tested compounds Chromenyl Coumarate (CC) was found to be better inhibitor (IC50 = 48.49 ± 5.6 nM) than the reference drug DON (IC50 = 74.13 ± 8.3 nM), unequivocally amplifying its importance. The structure of the compound was found to play a vital role in the inhibitory efficiency, validating previous Structure Activity Relationship (SAR) reviews for coumarin. The mechanism of inhibition remained impervious when the experimental medium was switched from aqueous buffer to HSA, albeit noticeable change in the inhibition potency of the compound 3, 3'- Methylene-bis (4-hydroxy coumarin) (MHC) (38%) and CC (35%). Both the coumarin derivatives were observed to bind to the peripheral anionic site (PAS) of AChE and also found to displace the fluorescence marker thioflavinT (ThT) from AChE binding pocket. All experimental observations were seconded by molecular docking and MD simulation results. The inferences drawn in this study form a foundation for further investigation on these compounds; magnifying the probability of their usage as AD drugs and re-emphasizes the significance of drug delivery media while considering the inhibition potencies of targeted drugs.

Stacking as a Key Property for Creating Nanoparticles with Tunable Shape: The Case of Squalenoyl-Doxorubicin.

The development of elongated nanoparticles for drug delivery is of growing interest in recent years, due to longer blood circulation and improved efficacy compared to spherical counterparts. Squalenoyl-doxorubicin (SQ-Dox) conjugate was previously shown to form elongated nanoparticles with improved therapeutic efficacy and decreased toxicity compared to free doxorubicin. By using experimental and computational techniques, we demonstrate here that the specific physical properties of SQ-Dox, which include stacking and electrostatic interactions of doxorubicin as well as hydrophobic interactions of squalene, are involved in the formation of nanoassemblies with diverse elongated structures. We show that SQ-Dox bioconjugate concentration, ionic strength, and anion nature can be used to modulate the shape and stiffness of SQ-Dox nanoparticles. As those parameters are involved in nanoparticle behavior in biological media, these findings could bring interesting opportunities for drug delivery and serve as an example for the design of original nanodrugs with stacking properties tuned for particular clinical purposes.

The blind search for the closed states of hinge-bending proteins.

The hinge-bending proteins provide the most pronounced example of the large-amplitude slow motions in a number of proteins, which are critical for their functioning. They are often used as the test ground for developing novel approaches to the simulation of slow protein dynamics. In the present study, we present the algorithm, which allows physically-consistent simulations of slow protein dynamics in globular proteins. Our algorithm is based on the hierarchical clustering of the correlation patterns (HCCP) technique of domain identification, which allows subdividing the protein into the hierarchy of the rigid-body-like clusters. The clusters are allowed to rotate relative to one another on the automatically identified hinges. The clusters are found in the course of automated, objective and well-tested procedure. In the present communication, our technique is applied to 10 hinge-bending proteins. For each of the proteins, we performed the blind search for the closed conformation, staring from the open one. Resulting closed conformations are compared with the closed states observed in crystallographic structures. It is shown that our technique produces realistic closed conformations for 8 out of 10 studied proteins. This demonstrates that HCCP technique can be used for finding alternative protein conformations and for sampling the slow motions in proteins.

The influence of curvature on the properties of the plasma membrane. Insights from atomistic molecular dynamics simulations.

In this work we conduct a systematic analysis of the influence of curvature on various properties of a realistic model of mammalian plasma membrane with asymmetric lipid content of monolayers and a realistic concentration of cholesterol. In order to do this we developed new technique for simulating membranes with the global membrane curvature restricted to any desirable value while keeping free lateral diffusion of the lipids and without introducing artifacts or perturbing the membrane structure. We show that the thickness of the hydrophobic core of the concave monolayer decreases by approximately 1.3 Å in comparison to that of the flat membrane, while the thickness of the convex monolayer does not change. The order parameter of the lipid tails decreases significantly in the certain layers of the curved membrane. The area per lipid increases in the convex monolayer and decreases in the concave monolayer. The cholesterol inclination angle decreases when the curvature of a particular monolayer changes from concave to convex. The amount of cholesterol in the minor fraction located between the membrane leaflets is zero in the membrane with positive curvature and increases to 1.7% in the flat membrane and to 2.5% in the membrane with negative curvature.