Observe without disturbing: tracer particles sense local stresses in cell collectives without affecting the cancer cell dynamics.

Measurements of local stresses on the cancer cells (CCs), inferred by embedding inert compressible tracer particles (TPs) in a growing multicellular spheroid (MCS), show that pressure decreases monotonically as the distance from the core of the MCS increases. How faithfully do the TPs report the local stresses in the CCs is an important question because pressure buildup in the MCS is dynamically generated due to CC division, which implies that the CC dynamics should be minimally altered by the TPs. Here using theory and simulations, we show that although the TP dynamics is unusual, exhibiting sub-diffusive behavior on times less than the CC division times and hyper-diffusive dynamics in the long-time limit, they do not affect the long-time CC dynamics. The CC pressure profile within the MCS, which decays from a high value at the core to the periphery, is almost identical with and without the TPs. That the TPs have a small effect on the local stresses in the MCS implies that they are reasonale reporters of the CC microenvironment.

Effect of pectin and attapulgite filler on swelling, network parameters and controlled release of diltiazem hydrochloride from polyacrylic copolymer gel.

Pectin grafted polyacrylic copolymer hydrogels were made by free radical crosslink copolymerization of acrylic acid (AA) and acrylamide (AM) in an aqueous solution of pectin. N'N-methylene bis acrylamide (MBA) was used as a crosslinker. During the polymerization reaction the attapulgite (APG) filler was also incorporated in situ into the network of the copolymer gel. Several filled hydrogels were prepared by varying the amount of pectin and APG filler. These hydrogels were characterized by FTIR, 13C NMR, XRD, TGA, SEM, mechanical properties, DMA, swelling, diffusion characteristics and network parameters. The release kinetics of a model drug diltiazem hydrochloride (DT) was studied with these hydrogels. The wt% of pectin, APG and MBA was optimized with a central composite design (CCD) model of response surface methodology (RSM) with equilibrium swelling ratio (ESR), drug adsorption (mg/100 mg gel) and drug release% in 16 h as response. Accordingly, the hydrogel prepared with 5:1 AA:AM molar ratio, 25 wt% monomer concentration, 1% each of initiator and MBA concentration, 18 wt% pectin and 2 wt% APG showed an optimized ESR of 17.75, drug loading of 27.58 and a drug release % of 92.5 in 16 h at a solution pH of 7.4.

Correction: Spatially heterogeneous dynamics of cells in a growing tumor spheroid: comparison between theory and experiments.

Correction for 'Spatially heterogeneous dynamics of cells in a growing tumor spheroid: comparison between theory and experiments' by Sumit Sinha et al., Soft Matter, 2020, 16, 5294-5304, DOI: .

Spatially heterogeneous dynamics of cells in a growing tumor spheroid: comparison between theory and experiments.

Collective cell movement, characterized by multiple cells that are in contact for substantial periods of time and undergo correlated motion, plays a central role in cancer and embryogenesis. Recent imaging experiments have provided time-dependent traces of individual cells, thus providing an unprecedented picture of tumor spheroid growth. By using simulations of a minimal cell model, we analyze the experimental data that map the movement of cells in a fibrosarcoma tumor spheroid embedded in a collagen matrix. Both simulations and experiments show that cells in the core of the spheroid exhibit subdiffusive glassy dynamics (mean square displacement, Δ(t) ≈ tα with α < 1), whereas cells in the periphery exhibit superdiffusive motion, Δ(t) ≈ tα with α > 1. The motion of most of the cells near the periphery is highly persistent and correlated directional motion due to cell doubling and apoptosis rates, thus explaining the observed superdiffusive behavior. The α values for cells in the core and periphery, extracted from simulations and experiments, are in near quantitative agreement with each other, which is surprising given that no parameter in the model was used to fit the measurements. The qualitatively different dynamics of cells in the core and periphery is captured by the fourth order susceptibility, introduced to characterize metastable states in glass forming systems. Analyses of the velocity autocorrelation of individual cells show remarkable spatial heterogeneity with no two cells exhibiting similar behavior. The prediction that α should depend on the location of the cells in the tumor is amenable to experimental testing. The highly heterogeneous dynamics of cells in the tumor spheroid provides a plausible mechanism for the origin of intratumor heterogeneity.

Heterocycle amide isosteres: An approach to overcoming resistance for HIV-1 integrase strand transfer inhibitors.

A series of heterocyclic pyrimidinedione-based HIV-1 integrase inhibitors was prepared and screened for activity against purified integrase enzyme and/or viruses modified with the following mutations within integrase: Q148R, Q148H/G140S and N155H. These are mutations that result in resistance to the first generation integrase inhibitors raltegravir and elvitegravir. Based on consideration of drug-target interactions, an approach was undertaken to replace the amide moiety of the first generation pyrimidinedione inhibitor with azole heterocycles that could retain potency against these key resistance mutations. An imidazole moiety was found to be the optimal amide substitute and the observed activity was rationalized with the use of calculated properties and modeling. Rat pharmacokinetic (PK) studies of the lead imidazole compounds demonstrated moderate clearance and moderate exposure.

Universal Nature of Collapsibility in the Context of Protein Folding and Evolution.

Theory and simulations predicted that the sizes of the unfolded states of globular proteins should decrease as the denaturant concentration is reduced from a high to a low value. However, small angle X-ray scattering (SAXS) data were used to assert the opposite, while interpretation of single molecule Förster resonance energy transfer experiments (FRET) supported the theoretical predictions. The disagreement between the two experiments is the SAXS-FRET controversy. By harnessing recent advances in SAXS and FRET experiments and setting these findings in the context of a general theory and simulations, which do not rely on experimental data, we establish that compaction of unfolded states under native conditions is universal. The theory also predicts that proteins rich in ß-sheets are more collapsible than α-helical proteins. Because the extent of compaction is small, experiments have to be accurate and their interpretations should be as model-free as possible. Theory also suggests that collapsibility itself could be a physical restriction on the evolution of foldable sequences, and also provides a physical basis for the origin of multidomain proteins.

Giant Casimir Nonequilibrium Forces Drive Coil to Globule Transition in Polymers.

We develop a theory to probe the effect of nonequilibrium fluctuation-induced forces on the size of a polymer confined between two horizontal, thermally conductive plates subject to a constant temperature gradient, ∇ T. We assume that (a) the solvent is good and (b) the distance between the plates is large so that in the absence of a thermal gradient the polymer is a coil, whose size scales with the number of monomers as Nν, with ν ≈ 0.6. We find that above a critical temperature gradient, ∇ Tc ≈ N-5/4, a favorable attractive monomer-monomer interaction due to the giant Casimir force (GCF) overcomes the chain conformational entropy, resulting in a coil-globule transition. Our predictions can be verified using light-scattering experiments with polymers, such as polystyrene or polyisoprene in organic solvents in which the GCF is attractive.

Origin of superdiffusive behavior in a class of nonequilibrium systems.

Experiments and simulations have established that dynamics in a class of living and abiotic systems that are far from equilibrium exhibit superdiffusive behavior at long times, which in some cases (for example, an evolving tumor) is preceded by slow glass-like dynamics. By using the evolution of a collection of tumor cells, driven by mechanical forces and subject to cell birth and apoptosis, as a case study we show theoretically that on short timescales the mean-square displacement is subdiffusive due to jamming, whereas at long times it is superdiffusive. The results obtained by using a stochastic quantization method, which is needed because of the absence of the fluctuation-dissipation theorem, show that the superdiffusive behavior is universal and impervious to the nature of cell-cell interactions. Surprisingly, the theory also quantitatively accounts for the nontrivial dynamics observed in simulations of a model soap foam characterized by creation and destruction of spherical bubbles, which suggests that the two nonequilibrium systems belong to the same universality class. The theoretical prediction for the superdiffusion exponent is in excellent agreement with simulations for collective motion of tumor cells and dynamics associated with soap bubbles.

Charge fluctuation effects on the shape of flexible polyampholytes with applications to intrinsically disordered proteins.

Random polyampholytes (PAs) contain positively and negatively charged monomers that are distributed randomly along the polymer chain. The interaction between charges is assumed to be given by the Debye-Huckel potential. We show that the size of the PA is determined by an interplay between electrostatic interactions, giving rise to the polyelectrolyte effect due to net charge per monomer (σ) and an effective attractive PA interaction due to charge fluctuations, δσ. The interplay between these terms gives rise to non-monotonic dependence of the radius of gyration, R g , on the inverse Debye length, κ, when PA effects are important ( δ σ σ > 1 ). In the opposite limit, R g decreases monotonically with increasing κ. Simulations of PA chains, using a charged bead-spring model, further corroborate our theoretical predictions. The simulations unambiguously show that conformational heterogeneity manifests itself among sequences that have identical PA parameters. A clear implication is that the phases of PA sequences, and by inference intrinsically disordered proteins (IDPs), cannot be determined using only the bare PA parameters (σ and δσ). The theory is used to calculate the changes in R g on N, the number of residues for a set of IDPs. For a certain class of IDPs, with N between 24 and 441, the size grows as R g ∼ N 0.6, which agrees with data from small angle X-ray scattering experiments.

Design, Synthesis, and SAR of C-3 Benzoic Acid, C-17 Triterpenoid Derivatives. Identification of the HIV-1 Maturation Inhibitor 4-((1 R,3a S,5a R,5b R,7a R,11a S,11b R,13a R,13b R)-3a-((2-(1,1-Dioxidothiomorpholino)ethyl)amino)-5a,5b,8,8,11a-pentamethyl-1-(prop-1-en-2-yl)-2,3,3a,4,5,5a,5b,6,7,7a,8,11,11a,11b,12,13,13a,13b-octadecahydro-1 H-cyclopenta[ a]chrysen-9-yl)benzoic Acid (GSK3532795, BMS-955176).

GSK3532795, formerly known as BMS-955176 (1), is a potent, orally active, second-generation HIV-1 maturation inhibitor (MI) that advanced through phase IIb clinical trials. The careful design, selection, and evaluation of substituents appended to the C-3 and C-17 positions of the natural product betulinic acid (3) was critical in attaining a molecule with the desired virological and pharmacokinetic profile. Herein, we highlight the key insights made in the discovery program and detail the evolution of the structure-activity relationships (SARs) that led to the design of the specific C-17 amine moiety in 1. These modifications ultimately enabled the discovery of 1 as a second-generation MI that combines broad coverage of polymorphic viruses (EC50 <15 nM toward a panel of common polymorphisms representative of 96.5% HIV-1 subtype B virus) with a favorable pharmacokinetic profile in preclinical species.

Protein collapse is encoded in the folded state architecture.

Folded states of single domain globular proteins are compact with high packing density. The radius of gyration, Rg, of both the folded and unfolded states increase as Nν where N is the number of amino acids in the protein. The values of the Flory exponent ν are, respectively, ≈⅓ and ≈0.6 in the folded and unfolded states, coinciding with those for homopolymers. However, the extent of compaction of the unfolded state of a protein under low denaturant concentration (collapsibility), conditions favoring the formation of the folded state, is unknown. We develop a theory that uses the contact map of proteins as input to quantitatively assess collapsibility of proteins. Although collapsibility is universal, the propensity to be compact depends on the protein architecture. Application of the theory to over two thousand proteins shows that collapsibility depends not only on N but also on the contact map reflecting the native structure. A major prediction of the theory is that ß-sheet proteins are far more collapsible than structures dominated by α-helices. The theory and the accompanying simulations, validating the theoretical predictions, provide insights into the differing conclusions reached using different experimental probes assessing the extent of compaction of proteins. By calculating the criterion for collapsibility as a function of protein length we provide quantitative insights into the reasons why single domain proteins are small and the physical reasons for the origin of multi-domain proteins. Collapsibility of non-coding RNA molecules is similar ß-sheet proteins structures adding support to "Compactness Selection Hypothesis".

Optimal information transfer in enzymatic networks: A field theoretic formulation.

Signaling in enzymatic networks is typically triggered by environmental fluctuations, resulting in a series of stochastic chemical reactions, leading to corruption of the signal by noise. For example, information flow is initiated by binding of extracellular ligands to receptors, which is transmitted through a cascade involving kinase-phosphatase stochastic chemical reactions. For a class of such networks, we develop a general field-theoretic approach to calculate the error in signal transmission as a function of an appropriate control variable. Application of the theory to a simple push-pull network, a module in the kinase-phosphatase cascade, recovers the exact results for error in signal transmission previously obtained using umbral calculus [Hinczewski and Thirumalai, Phys. Rev. X 4, 041017 (2014)2160-330810.1103/PhysRevX.4.041017]. We illustrate the generality of the theory by studying the minimal errors in noise reduction in a reaction cascade with two connected push-pull modules. Such a cascade behaves as an effective three-species network with a pseudointermediate. In this case, optimal information transfer, resulting in the smallest square of the error between the input and output, occurs with a time delay, which is given by the inverse of the decay rate of the pseudointermediate. Surprisingly, in these examples the minimum error computed using simulations that take nonlinearities and discrete nature of molecules into account coincides with the predictions of a linear theory. In contrast, there are substantial deviations between simulations and predictions of the linear theory in error in signal propagation in an enzymatic push-pull network for a certain range of parameters. Inclusion of second-order perturbative corrections shows that differences between simulations and theoretical predictions are minimized. Our study establishes that a field theoretic formulation of stochastic biological signaling offers a systematic way to understand error propagation in networks of arbitrary complexity.

Discovery of BMS-955176, a Second Generation HIV-1 Maturation Inhibitor with Broad Spectrum Antiviral Activity.

HIV-1 maturation inhibition (MI) has been clinically validated as an approach to the control of HIV-1 infection. However, identifying an MI with both broad polymorphic spectrum coverage and good oral exposure has been challenging. Herein, we describe the design, synthesis, and preclinical characterization of a potent, orally active, second generation HIV-1 MI, BMS-955176 (2), which is currently in Phase IIb clinical trials as part of a combination antiretroviral regimen.

Identification and Characterization of BMS-955176, a Second-Generation HIV-1 Maturation Inhibitor with Improved Potency, Antiviral Spectrum, and Gag Polymorphic Coverage.

BMS-955176 is a second-generation human immunodeficiency virus type 1 (HIV-1) maturation inhibitor (MI). A first-generation MI, bevirimat, showed clinical efficacy in early-phase studies, but ∼50% of subjects had viruses with reduced susceptibility associated with naturally occurring polymorphisms in Gag near the site of MI action. MI potency was optimized using a panel of engineered reporter viruses containing site-directed polymorphic changes in Gag that reduce susceptibility to bevirimat (including V362I, V370A/M/Δ, and T371A/Δ), leading incrementally to the identification of BMS-955176. BMS-955176 exhibits potent activity (50% effective concentration [EC50], 3.9 ± 3.4 nM [mean ± standard deviation]) toward a library (n = 87) of gag/pr recombinant viruses representing 96.5% of subtype B polymorphic Gag diversity near the CA/SP1 cleavage site. BMS-955176 exhibited a median EC50 of 21 nM toward a library of subtype B clinical isolates assayed in peripheral blood mononuclear cells (PBMCs). Potent activity was maintained against a panel of reverse transcriptase, protease, and integrase inhibitor-resistant viruses, with EC50s similar to those for the wild-type virus. A 5.4-fold reduction in EC50 occurred in the presence of 40% human serum plus 27 mg/ml of human serum albumin (HSA), which corresponded well to an in vitro measurement of 86% human serum binding. Time-of-addition and pseudotype reporter virus studies confirm a mechanism of action for the compound that occurs late in the virus replication cycle. BMS-955176 inhibits HIV-1 protease cleavage at the CA/SP1 junction within Gag in virus-like particles (VLPs) and in HIV-1-infected cells, and it binds reversibly and with high affinity to assembled Gag in purified HIV-1 VLPs. Finally, in vitro combination studies showed no antagonistic interactions with representative antiretrovirals (ARVs) of other mechanistic classes. In conclusion, BMS-955176 is a second-generation MI with potent in vitro anti-HIV-1 activity and a greatly improved preclinical profile compared to that of bevirimat.

C-3 benzoic acid derivatives of C-3 deoxybetulinic acid and deoxybetulin as HIV-1 maturation inhibitors.

A series of C-3 phenyl- and heterocycle-substituted derivatives of C-3 deoxybetulinic acid and C-3 deoxybetulin was designed and synthesized as HIV-1 maturation inhibitors (MIs) and evaluated for their antiviral activity and cytotoxicity in cell culture. A 4-subsituted benzoic acid moiety was identified as an advantageous replacement for the 3'3'-dimethylsuccinate moiety present in previously disclosed MIs that illuminates new aspects of the topography of the pharmacophore. The new analogs exhibit excellent in vitro antiviral activity against wild-type (wt) virus and a lower serum shift when compared with the prototypical HIV-1 MI bevirimat (1, BVM), the first MI to be evaluated in clinical studies. Compound 9a exhibits comparable cell culture potency toward wt virus as 1 (WT EC50=16 nM for 9a compared to 10nM for 1). However, the potency of 9a is less affected by the presence of human serum, while the compound displays a similar pharmacokinetic profile in rats to 1. Hence 9a, the 4-benzoic acid derivative of deoxybetulinic acid, represents a new starting point from which to explore the design of a 2nd generation MI.

Controlled release of tinidazole and theophylline from chitosan based composite hydrogels.

Several composite hydrogels were synthesized by free radical crosslink copolymerization of acrylic acid (AA) and N' methylene bis-acrylamide (MBA) in the presence of chitosan (CS). During polymerization CS was incorporated in situ in the crosslinked polyacrylic acid gel to produce composite hydrogels. The structure and properties of the hydrogels were characterized by FTIR, (13)C NMR, DTA-TGA, XRD, swelling and diffusion characteristic and also network parameters. The loading and the in vitro release behaviours of theophylline and tinidazole model drugs were studied with these hydrogels. The wt% of CS and MBA and pH of the medium was found to strongly influence the drug release behaviour of the gels. Accordingly, the release rate of these two drugs was much faster at pH of 7.6 than at pH 1.5.

Synthesis, characterization, swelling and drug release behavior of semi-interpenetrating network hydrogels of sodium alginate and polyacrylamide.

Several semi interpenetrating network (SIPN) type hydrogels were synthesized by in-situ free radical crosslink copolymerization of acrylamide and crosslinker N,N'-methylene bisacrylamide (MBA) in aqueous solution of sodium alginate (SA).These SIPN hydrogels were characterized by FTIR, NMR SEM, DTA-TGA, XRD, PZC and also by swelling characteristics and network parameters. Adsorption (loading) and release of acetaminophen drug were studied with these hydrogels. Solution pH, crosslinker concentration and monomer to SA weight ratio of the hydrogels were found to have a strong effect on adsorption and in vitro release profile of the drug from the gel matrix.

On noise induced Poincaré-Andronov-Hopf bifurcation.

It has been numerically seen that noise introduces stable well-defined oscillatory state in a system with unstable limit cycles resulting from subcritical Poincaré-Andronov-Hopf (or simply Hopf) bifurcation. This phenomenon is analogous to the well known stochastic resonance in the sense that it effectively converts noise into useful energy. Herein, we clearly explain how noise induced imperfection in the bifurcation is a generic reason for such a phenomenon to occur and provide explicit analytical calculations in order to explain the typical square-root dependence of the oscillations' amplitude on the noise level below a certain threshold value. Also, we argue that the noise can bring forth oscillations in average sense even in the absence of a limit cycle. Thus, we bring forward the inherent general mechanism of the noise induced Hopf bifurcation naturally realisable across disciplines.

Exact solution of the Zwanzig-Lauritzen model of polymer crystallization under tension.

We solve a two-dimensional model for polymer chain folding in the presence of mechanical pulling force (f) exactly using equilibrium statistical mechanics. Using analytically derived expression for the partition function we determine the phase diagram for the model in the f-temperature (T) plane. A square root singularity in the susceptibility indicates a second order phase transition from a folded to an unfolded state at a critical force (fc) in the thermodynamic limit of infinitely long polymer chain. The temperature dependence of fc shows a reentrant phase transition, which is reflected in an increase in fc as T increases below a threshold value. As a result, for a range of f values, the unfolded state is stable at both low and high temperatures. The high temperature unfolded state is stabilized by entropy whereas the low temperature unfolded state is dominated by favorable energy. The exact calculation could serve as a benchmark for testing approximate theories that are used in analyzing single molecule pulling experiments.

The terminal (catalytic) adenosine of the HIV LTR controls the kinetics of binding and dissociation of HIV integrase strand transfer inhibitors.

Specific HIV integrase strand transfer inhibitors are thought to bind to the integrase active site, positioned to coordinate with two catalytic magnesium atoms in a pocket flanked by the end of the viral LTR. A structural role for the 3' terminus of the viral LTR in the inhibitor-bound state has not previously been examined. This study describes the kinetics of binding of a specific strand transfer inhibitor to integrase variants assembled with systematic changes to the terminal 3' adenosine. Kinetic experiments are consistent with a two-step binding model in which there are different functions for the terminal adenine base and the terminal deoxyribose sugar. Adenine seems to act as a "shield" which retards the rate of inhibitor association with the integrase active site, possibly by acting as an internal competitive inhibitor. The terminal deoxyribose is responsible for retarding the rate of inhibitor dissociation, either by sterically blocking inhibitor egress or by a direct interaction with the bound inhibitor. These findings further our understanding of the details of the inhibitor binding site of specific strand transfer inhibitors.