Search PubMed⌕ Search

Biomedical subjects

Juan J de Pablo

Publications and source records attributed to Juan J de Pablo.

At least 19 recordsLinked to original sources

Chromatin structures from integrated AI and polymer physics model.

The physical organization of the genome in three-dimensional space regulates many biological processes, including gene expression and cell differentiation. Three-dimensional characterization of genome structure is critical to understanding these biological processes. Direct experimental measurements of genome structure are challenging; computational models of chromatin structure are therefore necessary. We develop an approach that combines a particle-based chromatin polymer model, molecular simulation, and machine learning to efficiently and accurately estimate chromatin structure from indirect measures of genome structure. More specifically, we introduce a new approach where the interaction parameters of the polymer model are extracted from experimental Hi-C data using a graph neural network (GNN). We train the GNN on simulated data from the underlying polymer model, avoiding the need for large quantities of experimental data. The resulting approach accurately estimates chromatin structures across all chromosomes and across several experimental cell lines despite being trained almost exclusively on simulated data. The proposed approach can be viewed as a general framework for combining physical modeling with machine learning, and it could be extended to integrate additional biological data modalities. Ultimately, we achieve accurate and high-throughput estimations of chromatin structure from Hi-C data, which will be necessary as experimental methodologies, such as single-cell Hi-C, improve.

Chromatin↗

N log N method for hydrodynamic interactions of confined polymer systems: Brownian dynamics.

A Brownian dynamics simulation technique is presented where a Fourier-based N log N approach is used to calculate hydrodynamic interactions in confined flowing polymer systems between two parallel walls. A self-consistent coarse-grained Langevin description of the polymer dynamics is adopted in which the polymer beads are treated as point forces. Hydrodynamic interactions are therefore included in the diffusion tensor through a Green's function formalism. The calculation of Green's function is based on a generalization of a method developed for sedimenting particles by Mucha et al. [J. Fluid Mech. 501, 71 (2004)]. A Fourier series representation of the Stokeslet that satisfies no-slip boundary conditions at the walls is adopted; this representation is arranged in such a way that the total O(N2) contribution of bead-bead interactions is calculated in an O(N log N) algorithm. Brownian terms are calculated using the Chebyshev polynomial approximation proposed by Fixman [Macromolecules 19, 1195 (1986); 19, 1204 (1986)] for the square root of the diffusion tensor. The proposed Brownian dynamics simulation methodology scales as O(N1.25 log N). Results for infinitely dilute systems of dumbbells are presented to verify past predictions and to examine the performance and numerical consistency of the proposed method.

Algorithms↗

Thermodynamic stability of beta-peptide helices and the role of cyclic residues.

Beta-peptides are emerging as an attractive class of peptidomimetic molecules. In contrast to naturally occurring alpha-peptides, short oligomers of beta-amino acids (comprising just 4-6 monomers) exhibit stable secondary structures that make them amenable for quantitative, concerted experimental and theoretical studies of the effects of particular chemical interactions on structure. In this work, molecular simulations are used to study the thermodynamic stability of helical conformations formed by beta-peptides containing varying proportions of acyclic (beta(3)) and cyclic (ACH) residues. More specifically, several beta-peptides differing only in their content of cyclic residues are considered in this work. Previous computational studies of beta-peptides have relied mostly on energy minimization of molecular dynamics simulations. In contrast, our study relies on density-of-states based Monte Carlo simulations to calculate the free energy and examine the stability of various folded structures of these molecules along a well-defined order parameter. By resorting to an expanded-ensemble formalism, we are able to determine the free energy required to unfold specific molecules, a quantity that could be measured directly through single-molecule force spectroscopy. Simulations in both implicit and explicit solvents have permitted a systematic study of the role of cyclic residues and electrostatics on the stability of secondary structures. The molecules considered in this work are shown to exhibit stable H-14 helical conformations and, in some cases, relatively stable H-12 conformations, thereby suggesting that solvent quality may be used to manipulate the hydrogen-bonding patterns and structure of these peptides.

Amino Acids↗

3-D microwell culture of human embryonic stem cells.

Human embryonic stem cells (hESCs) have the ability to proliferate indefinitely and differentiate into each of the embryonic cell lineages. Great care is required to maintain undifferentiated hESC cultures since spontaneous differentiation often occurs in culture, presumably resulting from soluble factors, cell-cell contact, and/or cell-matrix signaling. hESC differentiation is typically stimulated via generation of embryoid bodies (EBs) and lineage commitment of individual cells depends upon numerous cues throughout the EB environment, including EB shape and size. Common EB formation protocols, however, produce a very heterogeneous size distribution, perhaps reducing efficiency of directed differentiation. We have developed a 3-D microwell-based method to maintain undifferentiated hESC cultures for weeks without passaging using physical and extracellular matrix patterning constraints to limit colony growth. Over 90% of hESCs cultured in microwells for 2-3 weeks were viable and expressed the hESC transcription marker Oct-4. Upon passaging to Matrigel-coated tissue culture-treated polystyrene dishes (TCPS), microwell cultured hESCs maintained undifferentiated proliferation. Microwell culture also permits formation of hESC colonies with a defined size, which can then be used to form monodisperse EBs. When cultured in this system, hESCs retained pluripotency and self-renewal, and were able to be passaged to standard unconstrained culture conditions.

Biocompatible Materials↗

Anisotropic nanoparticles immersed in a nematic liquid crystal: defect structures and potentials of mean force.

We report results for the potential of mean force (PMF) and the defect structures that arise when spherocylindrical nanoparticles are immersed in a nematic liquid crystal. Using a dynamic field theory for the tensor order parameter Q of the liquid crystal, we analyzed configurations, including one, two, and three elongated particles, with strong homeotropic anchoring at their surfaces. For systems with one nanoparticle, the most stable configuration is achieved when the spherocylinder is placed with its long axis perpendicular to the far-field director, for which the defect structure consists of an elongated Saturn ring. For systems with two or three nanoparticles with their long axes placed perpendicular to the far-field director, at small separations the defect structures consist of incomplete Saturn rings fused with new disclination rings orthogonal to the original ones, in analogy to results previously observed for spherical nanoparticles. The shape of these orthogonal rings depends on the nanoparticles' configuration, i.e., triangular, linear, or parallel with respect to their long axis. A comparison of the PMFs indicates that the latter configuration is the most stable. The stability of the different arrays depends on whether orthogonal disclination rings form or not, their size, and the curvature effects in the interparticle regions. Our results suggest that the one-elastic-constant approximation is valid for the considered systems; similar results were obtained when a three-constant expression is used to represent the elastic free energy. The attractive interactions between the elongated particles were compared to those observed for spheres of similar diameters. Similar interparticle energies were observed for linear arrays; in contrast, parallel and triangular arrays of spherocylinders yielded interactions that were up to 3.4 times stronger than those observed for spherical particles.

Journal Article↗

Influence of confinement on the fragility of antiplasticized and pure polymer films.

We investigate by molecular dynamics simulation how thin film confinement modifies the fragility of a model glass-forming liquid characterized previously in the bulk. Film confinement is found to reduce the relative fragility of the polymer fluid, leading to effects similar to simulations of the addition of an antiplasticizer additive. A reduction of fragility is not observed for the antiplasticized polymer film. These effects are interpreted in terms of variation in the string (cooperatively moving segments) concentration with film confinement and the addition of antiplasticizing additives.

Journal Article↗

Cross-stream migration of flexible molecules in a nanochannel.

Molecular dynamics simulations are used to examine the cross-stream chain migration phenomenon in dilute polymer solutions that are flowing in nanochannels. In particular, both uniform planar shear (Couette) and pressure driven (Poiseuille) flows of a dilute polymer solution are studied using a bead-spring representation of polymer chains and a coarse grained model for the solvent. Our results show that three mechanisms govern the migration of deformable molecules in a nanochannel: (1) chain-wall hydrodynamic interactions, (2) thermal diffusion, and (3) gradient in chain mobility. These results are discussed in the context of recent experimental, numerical and theoretical work.

Bacteriophage lambda↗

Measurement of the azimuthal anchoring energy of liquid crystals in contact with oligo(ethylene glycol)-terminated self-assembled monolayers supported on obliquely deposited gold films.

We report measurements of the orientations and azimuthal anchoring energies of the nematic liquid crystal 4-cyano-4'-pentylbiphenyl (5CB) on polycrystalline gold films that are deposited from a vapor at an oblique angle of incidence and subsequently decorated with organized monolayers of oligomers of ethylene glycol. Whereas the gold films covered with monolayers presenting tetra(ethylene glycol) (EG4) lead to orientations of 5CB that are perpendicular to the plane of incidence of the gold, monolayers presenting tri(ethylene glycol) (EG3) direct 5CB to orient parallel to the plane of incidence of the gold during deposition of the gold film. We also measure the azimuthal anchoring energy of the 5CB to be smaller on the surfaces presenting EG3 (3.2 +/- 0.8 microJ/m2) as compared to EG4 (5.5 +/- 0.9 microJ/m2). These measurements, when combined with other results presented in this paper, are consistent with a physical model in which the orientation and anchoring energies of LCs on these surfaces are influenced by both (i) short-range interactions of 5CB with organized oligomers of ethylene glycol at these surfaces and (ii) long-range interactions of 5CB with the nanometer-scale topography of the obliquely deposited films. For surfaces presenting EG3, these short- and long-range interactions oppose each other, leading to small net values of anchoring energies that we predict are dependent on the level of order in the EG3 SAM. These measurements provide insights into the balance of interactions that control the orientational response of LCs to biological species (proteins, viruses, cells) on these surfaces.

Biphenyl Compounds↗

Effects of trehalose on the phase behavior of DPPC-cholesterol unilamellar vesicles.

A systematic study is presented of the effects of trehalose on the physical properties of extruded DPPC-cholesterol unilamellar vesicles. Particular emphasis is placed on examining how the interactions present in the hydrated state translate into those in the dehydrated state. Observations from HSDSC and DSC are used to examine the phase behavior of hydrated and dehydrated vesicles, respectively. The concentration of trehalose inside and outside the vesicles is manipulated, and is shown to affect the relative stability of the membranes. Our results show for the first time that a combination of high inner and low outer trehalose concentration is able to decrease the gel-to-liquid crystalline phase temperature (T(m)), while any other combination will not. Upon dehydration, the T(m) of all lipid mixtures increases. The extent of the increase depends on the trehalose distribution across the bilayer. The T(m) changes in the same direction with trehalose concentration for both freeze-dried and fully hydrated samples, suggesting that the trehalose distribution across the vesicle membrane, as well as the trehalose-phospholipid interaction, is maintained upon lyophilization. The results presented in this work may aid in the formulation of systems to be used in the lyophilization of liposomes for drug delivery applications.

1,2-Dipalmitoylphosphatidylcholine↗

Fabrication of complex three-dimensional nanostructures from self-assembling block copolymer materials on two-dimensional chemically patterned templates with mismatched symmetry.

A study is presented of the self-assembly of a lamella-forming blend of a diblock copolymer and its respective homopolymers on periodically patterned substrates consisting of square arrays of spots, that preferentially attract one component, as a function of pattern dimensions and film thickness. The blend morphology follows the pattern at the substrate and forms a single quadratically perforated lamella (QPL). At intermediate film thicknesses necks connect this QPL to the film surface, resulting in a bicontinuous morphology. The necks do not register with the underlying square lattice but exhibit a substantial amount of hexagonal short-range order. For thicker films we observe bicontinuous morphologies consisting of parallel lamellae with disordered perforations. These results demonstrate a promising strategy for the fabrication of complex interfacial nanostructures from two-dimensional chemically patterned templates.

Journal Article↗

Inhibition of human embryonic stem cell differentiation by mechanical strain.

Mechanical forces have been reported to induce proliferation and/or differentiation in many cell types, but the role of mechanotransduction during embryonic stem cell fate decisions is unknown. To ascertain the role of mechanical strain in human embryonic stem cell (hESC) differentiation, we measured the rate of hESC differentiation in the presence and absence of biaxial cyclic strain. Above a threshold of 10% cyclic strain, applied to a deformable elastic substratum upon which the hESC colonies were cultured, hESC differentiation was reduced and self-renewal was promoted without selecting against survival of differentiated or undifferentiated cells. Frequency of mechanical strain application had little effect on extent of differentiation. hESCs cultured under cyclic strain retained pluripotency, evidenced by their ability to differentiate to cell lineages in all three germ layers. Mechanical inhibition of hESC differentiation could not be traced to secretion of chemical factors into the media suggesting that mechanical forces may directly regulate hESC differentiation. Mechanical strain is not sufficient to inhibit differentiation, however, in unconditioned medium, hESCs grown under strain differentiated at the same rate as cells cultured in the absence of strain. Thus, while mechanical forces play a role in regulating hESC self-renewal and differentiation, they must act synergistically with chemical signals. These findings imply that application of mechanical forces may be useful, in combination with chemical and matrix-encoded signals, towards controlling differentiation of hESCs for therapeutic applications.

Cell Culture Techniques↗

Generation and differentiation of human embryonic stem cell-derived keratinocyte precursors.

Human embryonic stem cells (hESC) hold tremendous potential in the future of tissue engineering, offering promise as a source of virtually unlimited quantities of desired cell and tissue types. We have identified soluble chemical and extracellular matrix factors that permit isolation of keratinocyte precursors from hESCs. Culturing embryoid bodies (EB) formed from hESCs in a defined serum-free keratinocyte growth medium on a gelatin matrix generated keratin 14 (K14) expressing cells with an epithelial morphology. These K14 expressing cells could be subcultured in medium supplemented with hydrocortisone and induced to stratify and terminally differentiate by addition of calcium. Optimum times for obtaining K14 expressing cells were found for EB formation and for differentiation and growth of cultures after EB plating. EB formation was not necessary to generate keratinocyte precursors; direct transfer of hESC colonies to keratinocyte growth medium permitted differentiation into the keratinocyte lineage. With further studies to optimize generation and purification of hESC-derived keratinocyte precursors, these cells could provide a source of epidermal cells for skin tissue engineering applications in vitro or in vivo.

Biomarkers↗

Modulating membrane properties: the effect of trehalose and cholesterol on a phospholipid bilayer.

The protective properties of trehalose on cholesterol-containing lipid dipalmitoylphosphatidylcholine (DPPC) bilayers are studied through molecular simulations. The ability of the disaccharide to interact with the phospholipid headgroups and stabilize the membrane persists even at high cholesterol concentrations and restricts some of the changes to the structure that would otherwise be imposed by cholesterol molecules. Predictions of bilayer properties such as area per lipid, tail ordering, and chain conformation support the notion that the disaccharide decreases the main melting transition in these multicomponent model membranes, which correspond more closely to common biological systems than pure bilayers. Molecular simulations indicate that the membrane dynamics are slowed considerably by the presence of trehalose, indicating that high sugar concentrations would serve to avert possible phase separations that could arise in mixed phospholipid systems. Various time correlation functions suggest that the character of the modifications in lipid dynamics induced by trehalose and cholesterol is different in the hydrophilic and hydrophobic regions of the membrane.

1,2-Dipalmitoylphosphatidylcholine↗

Confinement effects on the thermodynamics of protein folding: Monte Carlo simulations.

The effects of chaperonin-like cage-induced confinement on protein stability have been studied for molecules of varying sizes and topologies. Minimalist models based on Gō-like interactions are employed for the proteins, and density-of-states-based Monte Carlo simulations are performed to accurately characterize the thermodynamic transitions. This method permits efficient sampling of conformational space and yields precise estimates of free energy and entropic changes associated with protein folding. We find that confinement-driven stabilization is not only dependent on protein size and cage radius, but also on the specific topology. The choice of the confining potential is also shown to have an effect on the observed stabilization and the scaling behavior of the stabilization with respect to the cage size.

Computer Simulation↗

Stability of Na-, K-, and Ca-montmorillonite at high temperatures and pressures: a Monte Carlo simulation.

Monte Carlo grand canonical molecular simulations on the hydration of Na-, K-, and Ca-montmorillonite show that between 333 and 533 K and 300-1300 bar Na-montmorillonite forms stable one-layer hydrates of d(001) spacings 12.64-12.38 Angstroms, K-montmorillonite of 12.78-12.59 Angstroms, and Ca-montmorillonite of 12.48-12.32 Angstroms. A two-layer hydrate of 14.80 Angstroms occurs for Na-montmorillonite at 533 K and 1300 bar, for K-montmorillonite of 15.32 Angstroms at 533 K and 1300 bar and of 14.74 Angstroms at 533 K and 2000 bar, and for Ca-montmorillonite of 13.83 Angstroms at 473 K and 1000 bar. Three-layer hydrates may possibly form within these same ranges. Outside of them, one-layer hydrates simulate as the only stable hydrates. In sedimentary basins, the two-layer hydrate of Ca-montmorillonite will locate at 6.7 km depth and those of Na- and K-montmorillonite at 8.7 km depth; above and below these depths, the one-layer hydrates are the stable phases.

Journal Article↗

Structure and stability of a model three-helix-bundle protein on tailored surfaces.

The interaction of protein molecules with surfaces is important in numerous applications. Theoretical work on protein adsorption has been limited. In particular, it is difficult to obtain quantitative predictions about the structure and stability of proteins on surfaces. In this study, density-of-states-based simulations were performed on a Gō-like model of a three-helix-bundle fragment from protein A (PDB ID: 1bdd). Both mechanical and thermal stability were investigated on neutral and attractive surfaces and compared to that in the absence of a surface. It was found that attaching the peptide to any type of surface decreases its melting temperature by as much as 9 K, depending upon orientation. Calorimetric cooperativity, as measured by van't Hoff to calorimetric enthalpy ratios, similarly decreased. It was also found that the mechanical strength of the peptide attached to surfaces is degraded to varying extents, depending upon the surface type and protein orientation. A comparison of mechanical and thermal stability showed that the two are not synonymous, but occur through different pathways, and that system configurations that are more thermally stable are not always so mechanically.

Adsorption↗

Molecular dynamics simulation of discontinuous volume phase transitions in highly-charged crosslinked polyelectrolyte networks with explicit counterions in good solvent.

The volumetric properties of highly-charged defect-free polyelectrolyte networks with tetrafunctional crosslinks are studied through molecular dynamics simulations in the canonical ensemble. The network backbone monomers, which are monovalent, and the counterions, which are mono-, di-, or trivalent, are modeled explicitly in the simulations, but the solvent is treated implicitly as a dielectric medium of good solvation quality. The osmotic pressure of the network-solvent system is found to depend greatly on the strength of electrostatic interactions. Discontinuous volume phase transitions are observed when the electrostatic interactions are strong, and the onset of these transitions shifts to higher solvent dielectricity as the counterion valency increases. The roles of the various virial contributions to the osmotic pressure are examined. The network elasticity entropy is found to behave nearly classically. As the network contracts and collapses with increasing strength of electrostatic interactions, the loss of counterion entropy leads to increased counterion osmotic pressure contributions via two mechanisms. The reduction in available configurational space increases the counterion translational entropy contribution to the ideal part of the osmotic pressure, and the greater number of counterion-monomer contacts formed due to counterion condensation and confinement increases the counterion excluded-volume entropy contribution to the excess part of the osmotic pressure. These observations contrast the decrease in the single ideal-gas-like counterion translational entropy contribution to the osmotic pressure predicted by the counterion condensation-charge renormalization theory. An accompanying decrease in the total electrostatic energy balances the loss of counterion excluded-volume entropy as the polyelectrolyte networks collapse in low-dielectric solvents. This interplay between the electrostatic energy and the counterion excluded-volume entropy appears to be responsible for the discontinuous volume phase transitions that are observed in polyelectrolyte networks. The structure of the polyelectrolyte network is also found to be affine in the swollen state, with constituent chains nearly fully extended, and nonaffine in the collapsed state, with the chains adopting a Gaussian conformation.

Journal Article↗

Platelet cryopreservation using a trehalose and phosphate formulation.

Long-term storage of platelets is infeasible due to platelet activation at low temperatures. In an effort to address this problem, we evaluated the effectiveness of a formulation combining trehalose and phosphate in protecting platelet structure and function following cryopreservation. An annexin V binding assay was used to quantify the efficacy of the trehalose and phosphate formulation in suppressing platelet activation during cryopreservation. Of the platelets cryopreserved with the trehalose plus phosphate formulation, 23% +/- 1.2% were nonactivated, compared with 9.8% +/- 0.26% nonactivated following cryopreservation with only trehalose. The presence of both trehalose and phosphate in the cryopreservation medium is critical for cell survival and preincubation in trehalose plus phosphate solutions further enhances viability. The effectiveness of trehalose plus phosphate in preserving platelets in a nonactivated state is comparable to 6% dimethyl sulfoxide (Me(2)SO). Measurements of platelet metabolic activity using an alamarBlue assay also established that trehalose plus phosphate is superior to trehalose alone. Finally, platelets protected by the trehalose plus phosphate formulation exhibit similar aggregation response upon thrombin addition as fresh platelets, but an increase of cytosolic calcium concentration upon thrombin addition was not observed in the cryopreserved platelets. These results suggest that trehalose and phosphate protect several aspects of platelet structure and function during cryopreservation, including an intact plasma membrane, metabolic activity, and aggregation in response to thrombin, but not intracellular calcium release in response to thrombin.

Annexin A5↗