A 3-D computational model predicts that cell deformation affects selectin-mediated leukocyte rolling

Leukocyte recruitment to sites of inflammation is initiated by their tethering and rolling on the activated endothelium under flow. Even though the fast kinetics and high tensile strength of selectin-ligand bonds are primarily responsible for leukocyte rolling, experimental evidence suggests that cellular properties such as cell deformability and microvillus elasticity actively modulate leukocyte rolling behavior. Previous theoretical models either assumed cells as rigid spheres or were limited to two-dimensional representations of deformable cells with deterministic receptor-ligand kinetics, thereby failing to accurately predict leukocyte rolling. We therefore developed a three-dimensional computational model based on the immersed boundary method to predict receptor-mediated rolling of deformable cells in shear flow coupled to a Monte Carlo method simulating the stochastic receptor-ligand interactions. Our model predicts for the first time that the rolling of more compliant cells is relatively smoother and slower compared to cells with stiffer membranes, due to increased cell-substrate contact area. At the molecular level, we show that the average number of bonds per cell as well as per single microvillus decreases with increasing membrane stiffness. Moreover, the average bond lifetime decreases with increasing shear rate and with increasing membrane stiffness, due to higher hydrodynamic force experienced by the cell. Taken together, our model captures the effect of cellular properties on the coupling between hydrodynamic and receptor-ligand bond forces, and successfully explains the stable leukocyte rolling at a wide range of shear rates over that of rigid microspheres.     

In vitro evolution of a hyperstable Gbeta1 variant

An in-vitro selection strategy was used to obtain strongly stabilized variants of the beta1 domain of protein G (Gbeta1). In a two-step approach, first candidate positions with a high potential for stabilization were identified in Gbeta1 libraries that were created by error-prone PCR, and then, after randomization of these positions by saturation mutagenesis, strongly stabilized variants were selected. For both steps the in-vitro selection method Proside was employed. Proside links the stability of a protein with the infectivity of a filamentous phage. Ultimately, residues from the two best selected variants were combined in a single Gbeta1 molecule. This variant with the four mutations E15V, T16L, T18I, and N37L showed an increase of 35.1 degrees C in the transition midpoint and of 28.5 kJ mol(-1) (at 70 degrees C) in the Gibbs free energy of stabilization. It was considerably more stable than the best variant from a previous Proside selection, in which positions were randomized that had originally been identified by computational design. Only a single substitution (T18I) was found in both selections. The best variants from the present selection showed a higher cooperativity of thermal unfolding, as indicated by an increase in the enthalpy of unfolding by about 60 kJ mol(-1). This increase is apparently correlated with the presence of Leu residues that were selected at the positions 16 and 37.     

Key role of coulombic interactions for the folding transition state of the cold shock protein

The cold shock protein CspB shows a five-stranded beta-sheet structure, and it folds rapidly via a native-like transition state. A previous Phi value analysis showed that most of the residues with Phi values close to one reside in strand beta1, and two of them, Lys5 and Lys7 are partially exposed charged residues. To elucidate how coulombic interactions of these two residues contribute to the energetic organisation of the folding transition state we performed comparative folding experiments in the presence of an ionic denaturant (guanidinium chloride) and a non-ionic denaturant (urea) and a double-mutant analysis. Lys5 contributes 6.6 kJ mol(-1) to the stability of the transition state, and half of it originates from screenable coulombic interactions. Lys7 contributes 5.3 kJ mol(-1), and 3.4 kJ mol(-1) of it are screened by salt. In the folded protein Lys7 interacts with Asp25, and the screenable coulombic interaction between these two residues is fully formed in the transition state. This suggests that long-range coulombic interactions such as those originating from Lys5 and Lys7 of CspB can be important for organizing and stabilizing native-like structure early in protein folding.     

The role of organ of Corti mass in passive cochlear tuning

The mechanism for passive cochlear tuning remains unsettled. Early models considered the organ of Corti complex (OCC) as a succession of spring-mass resonators. Later, traveling wave models showed that passive tuning could arise through the interaction of cochlear fluid mass and OCC stiffness without local resonators. However, including enough OCC mass to produce local resonance enhanced the tuning by slowing and thereby growing the traveling wave as it approached its resonant segment. To decide whether the OCC mass plays a role in tuning, the frequency variation of the wavenumber of the cochlear traveling wave was measured (in vivo, passive cochleae) and compared to theoretical predictions. The experimental wavenumber was found by taking the phase difference of basilar membrane motion between two longitudinally spaced locations and dividing by the distance between them. The theoretical wavenumber was a solution of the dispersion relation of a three-dimensional cochlear model with OCC mass and stiffness as the free parameters. The experimental data were only well fit by a model that included OCC mass. However, as the measurement position moved from a best-frequency place of 40 to 12 kHz, the role of mass was diminished. The notion of local resonance seems to only apply in the very high-frequency region of the cochlea.     

The single chlorophyll a molecule in the cytochrome b6f complex: unusual optical properties protect the complex against singlet oxygen

The cytochrome b(6)f complex of oxygenic photosynthesis mediates electron transfer between the reaction centers of photosystems I and II and facilitates coupled proton translocation across the membrane. High-resolution x-ray crystallographic structures (Kurisu et al., 2003; Stroebel et al., 2003) of the cytochrome b(6)f complex unambiguously show that a Chl a molecule is an intrinsic component of the cytochrome b(6)f complex. Although the functional role of this Chl a is presently unclear (Kuhlbrandt, 2003), an excited Chl a molecule is known to produce toxic singlet oxygen as the result of energy transfer from the excited triplet state of the Chl a to oxygen molecules. To prevent singlet oxygen formation in light-harvesting complexes, a carotenoid is typically positioned within approximately 4 A of the Chl a molecule, effectively quenching the triplet excited state of the Chl a. However, in the cytochrome b(6)f complex, the beta-carotene is too far (> or =14 Angstroms) from the Chl a for effective quenching of the Chl a triplet excited state. In this study, we propose that in this complex, the protection is at least partly realized through special arrangement of the local protein structure, which shortens the singlet excited state lifetime of the Chl a by a factor of 20-25 and thus significantly reduces the formation of the Chl a triplet state. Based on optical ultrafast absorption difference experiments and structure-based calculations, it is proposed that the Chl a singlet excited state lifetime is shortened due to electron exchange transfer with the nearby tyrosine residue. To our knowledge, this kind of protection mechanism against singlet oxygen has not yet been reported for any other chlorophyll-containing protein complex. It is also reported that the Chl a molecule in the cytochrome b(6)f complex does not change orientation in its excited state.     

Characterization of active site residues of nitroalkane oxidase

The flavoenzyme nitroalkane oxidase catalyzes the oxidation of primary and secondary nitroalkanes to the corresponding aldehydes and ketones plus nitrite. The structure of the enzyme shows that Ser171 forms a hydrogen bond to the flavin N5, suggesting that it plays a role in catalysis. Cys397 and Tyr398 were previously identified by chemical modification as potential active site residues. To more directly probe the roles of these residues, the S171A, S171V, S171T, C397S, and Y398F enzymes have been characterized with nitroethane as substrate. The C397S and Y398 enzymes were less stable than the wild-type enzyme, and the C397S enzyme routinely contained a substoichiometric amount of FAD. Analysis of the steady-state kinetic parameters for the mutant enzymes, including deuterium isotope effects, establishes that all of the mutations result in decreases in the rate constants for removal of the substrate proton by ∼5-fold and decreases in the rate constant for product release of ∼2-fold. Only the S171V and S171T mutations alter the rate constant for flavin oxidation. These results establish that these residues are not involved in catalysis, but rather are required for maintaining the protein structure.     

Cholesterol-induced protein sorting: an analysis of energetic feasibility

The mechanism(s) underlying the sorting of integral membrane proteins between the Golgi complex and the plasma membrane remain uncertain because no specific Golgi retention signal has been found. Moreover one can alter a protein's eventual localization simply by altering the length of its transmembrane domain (TMD). M. S. Bretscher and S. Munro (SCIENCE: 261:1280-1281, 1993) therefore proposed a physical sorting mechanism based on the hydrophobic match between the proteins' TMD and the bilayer thickness, in which cholesterol would regulate protein sorting by increasing the lipid bilayer thickness. In this model, Golgi proteins with short TMDs would be excluded from cholesterol-enriched domains (lipid rafts) that are incorporated into transport vesicles destined for the plasma membrane. Although attractive, this model remains unproven. We therefore evaluated the energetic feasibility of a cholesterol-dependent sorting process using the theory of elastic liquid crystal deformations. We show that the distribution of proteins between cholesterol-enriched and cholesterol-poor bilayer domains can be regulated by cholesterol-induced changes in the bilayer physical properties. Changes in bilayer thickness per se, however, have only a modest effect on sorting; the major effect arises because cholesterol changes also the bilayer material properties, which augments the energetic penalty for incorporating short TMDs into cholesterol-enriched domains. We conclude that cholesterol-induced changes in the bilayer physical properties allow for effective and accurate sorting which will be important generally for protein partitioning between different membrane domains.     

A simple model of backbone flexibility improves modeling of side-chain conformational variability

The considerable flexibility of side-chains in folded proteins is important for protein stability and function, and may have a role in mediating allosteric interactions. While sampling side-chain degrees of freedom has been an integral part of several successful computational protein design methods, the predictions of these approaches have not been directly compared to experimental measurements of side-chain motional amplitudes. In addition, protein design methods frequently keep the backbone fixed, an approximation that may substantially limit the ability to accurately model side-chain flexibility. Here, we describe a Monte Carlo approach to modeling side-chain conformational variability and validate our method against a large dataset of methyl relaxation order parameters derived from nuclear magnetic resonance (NMR) experiments (17 proteins and a total of 530 data points). We also evaluate a model of backbone flexibility based on Backrub motions, a type of conformational change frequently observed in ultra-high-resolution X-ray structures that accounts for correlated side-chain backbone movements. The fixed-backbone model performs reasonably well with an overall rmsd between computed and predicted side-chain order parameters of 0.26. Notably, including backbone flexibility leads to significant improvements in modeling side-chain order parameters for ten of the 17 proteins in the set. Greater accuracy of the flexible backbone model results from both increases and decreases in side-chain flexibility relative to the fixed-backbone model. This simple flexible-backbone model should be useful for a variety of protein design applications, including improved modeling of protein-protein interactions, design of proteins with desired flexibility or rigidity, and prediction of correlated motions within proteins.     

A mesoscale mechanical model of cellular interactions

Computational models of cell mechanics allow the precise interrogation of cell shape change. These morphological changes are required for cells to survive in diverse tissue environments. Here, we present a mesoscale mechanical model of cell-substrate interactions using the level set method based on experimentally measured parameters. By implementing a viscoelastic mechanical equivalent circuit, we accurately model whole-cell deformations that are important for a variety of cellular processes. To effectively model shape changes as a cell interacts with a substrate, we have included receptor-mediated adhesion, which is governed by catch-slip bond behavior. The effect of adhesion was explored by subjecting cells to a variety of different substrates including flat, curved, and deformable surfaces. Finally, we increased the accuracy of our simulations by including a deformable nucleus in our cells. This model sets the foundation for further exploration into computational analyses of multicellular interactions.     

Mutations in cytochrome b that affect kinetics of the electron transfer reactions at center N in the yeast cytochrome bc1 complex

We have examined the pre-steady-state kinetics and thermodynamic properties of the b hemes in variants of the yeast cytochrome bc₁ complex that have mutations in the quinone reductase site (center N). Trp-30 is a highly conserved residue, forming a hydrogen bond with the propionate on the high potential b heme (b_H heme). The substitution by a cysteine (W30C) lowers the redox potential of the heme and an apparent consequence is a lower rate of electron transfer between quinol and heme at center N. Leu-198 is also in close proximity to the b_H heme and a L198F mutation alters the spectral properties of the heme but has only minor effects on its redox properties or the electron transfer kinetics at center N. Substitution of Met-221 by glutamine or glutamate results in the loss of a hydrophobic interaction that stabilizes the quinone ligands. Ser-20 and Gln-22 form a hydrogen-bonding network that includes His-202, one of the carbonyl groups of the ubiquinone ring, and an active-site water. A S20T mutation has long-range structural effects on center P and thermodynamic effects on both b hemes. The other mutations (M221E, M221Q, Q22E and Q22T) do not affect the ubiquinol oxidation kinetics at center P, but do modify the electron transfer reactions at center N to various extents. The pre-steady reduction kinetics suggest that these mutations alter the binding of quinone ligands at center N, possibly by widening the binding pocket and thus increasing the distance between the substrate and the b_H heme. These results show that one can distinguish between the contribution of structural and thermodynamic factors to center N function.     

对氢键网络的干扰降低细胞色素b_6f复合体中Chla的光稳定性

细胞色素b6f蛋白复合体(Cyt b6f)参与光合膜上电子传递和跨膜质子转移,在体内以二体形式存在,每个单体只结合1分子叶绿素a(Chla).该Chla性质独特,光照条件下十分稳定,是甲醇中游离Chla的120~130倍,然而其光稳定性的机制仍未彻底阐明.Cyt b6f 2.7的晶体结构显示,Chla中心的Mg离子和H2O分子配位,并且该H2O分子通过氢键与复合体亚基Ⅳ的氨基酸G136和T137相互作用.本研究基于这一结构特点,对上述2个氨基酸进行了定点突变,以干扰、破坏氢键网络.结果发现,突变不仅导致蛋白与Chla结合能力下降,而且显著降低了Chla的光稳定性,这一结果表明,Cyt b6f中Chla相关的氢键网络对其稳定性具有重要的作用.     

Molecular Determinants of a Native-State Prolyl Isomerization

Prolyl cis/trans isomerizations determine the rates of many protein-folding reactions, and they can serve as molecular switches and timers. The energy required to shift the prolyl cis/trans equilibrium during these processes originates from conformational reactions that are linked structurally and energetically with prolyl isomerization. We used the N2 domain of the gene-3-protein of phage fd to elucidate how such an energetic linkage develops in the course of folding. The Asp160-Pro161 bond at the tip of a β hairpin of N2 is cis in the crystal structure, but in fact, it exists as a mixture of conformers in folded N2. During refolding, about 10 kJ mol^− 1 of conformational energy becomes available for a 75-fold shift of the cis/trans equilibrium constant at Pro161, from 7/93 in the unfolded to 90/10 in the folded form. We combined single- and double-mixing kinetic experiments with a mutational analysis to identify the structural origin of this proline shift energy and to elucidate the molecular path for the transfer of this energy to Pro161. It originates largely, if not entirely, from the two-stranded β sheet at the base of the Pro161 hairpin. The two strands improve their stabilizing interactions when Pro161 is cis, and this stabilization is propagated to Pro161, because the connector peptides between the β strands and Pro161 are native-like folded when Pro161 is cis. In the presence of a trans-Pro161, the connector peptides are locally unfolded, and thus, Pro161 is structurally and energetically uncoupled from the β sheet. Such interrelations between local folding and prolyl isomerization and the potential modulation by prolyl isomerases might also be used to break and reestablish slow communication pathways in proteins.     

On the evolution of reproductive restraint in malaria

Malaria is one of the leading causes of death among infectious diseases in the world, claiming over one million lives every year. By these standards, this highly complex parasite is extremely successful at generating new infections. Somewhat surprisingly, however, many malaria species seem to invest relatively little in gametocytes, converting only a small percentage of circulating asexual parasite forms into this transmissible form. In this article, we use mathematical models to explore three of the hypotheses that have been proposed to explain this apparent 'reproductive restraint' and develop a novel, fourth hypothesis. We find that only one of the previous three hypotheses we explore can explain such low gametocyte conversion rates, and this hypothesis involves a very specific form of density-dependent transmission-blocking immunity. Our fourth hypothesis also provides a potential explanation and is based on the occurrence of multiple infections and the resultant within-host competition between malaria strains that this entails. Further experimental work is needed to determine which of these two hypotheses provides the most likely explanation.     

Inhibition of autophagy attenuates pancreatic cancer growth independent of TP53/TRP53 status

Basal levels of autophagy are elevated in most pancreatic ductal adenocarcinomas (PDAC). Suppressing autophagy pharmacologically using chloroquine (CQ) or genetically with RNAi to essential autophagy genes inhibits human pancreatic cancer growth in vitro and in vivo, which presents possible treatment opportunities for PDAC patients using the CQ-derivative hydroxychloroquine (HCQ). Indeed, such clinical trials are ongoing. However, autophagy is a complex cellular mechanism to maintain cell homeostasis under stress. Based on its biological role, a dual role of autophagy in tumorigenesis has been proposed: at tumor initiation, autophagy helps maintain genomic stability and prevent tumor initiation; while in advanced disease, autophagy degrades and recycles cellular components to meet the metabolic needs for rapid growth. This model was proven to be the case in mouse lung tumor models. However, in contrast to prior work in various PDAC model systems, loss of autophagy in PDAC mouse models with embryonic homozygous Trp53 deletion does not inhibit tumor growth and paradoxically increases progression. This raised concerns whether there may be a genotype-dependent reliance of PDAC on autophagy. In a recent study, our group used a Trp53 heterozygous mouse PDAC model and human PDX xenografts to address the question. Our results demonstrate that autophagy inhibition was effective against PDAC tumors irrespective of TP53/TRP53 status.     

Biological insights from hydrogen exchange mass spectrometry

Over the past two decades, hydrogen exchange mass spectrometry (HXMS) has achieved the status of a widespread and routine approach in the structural biology toolbox. The ability of hydrogen exchange to detect a range of protein dynamics coupled with the accessibility of mass spectrometry to mixtures and large complexes at low concentrations result in an unmatched tool for investigating proteins challenging to many other structural techniques. Recent advances in methodology and data analysis are helping HXMS deliver on its potential to uncover the connection between conformation, dynamics and the biological function of proteins and complexes. This review provides a brief overview of the HXMS method and focuses on four recent reports to highlight applications that monitor structure and dynamics of proteins and complexes, track protein folding, and map the thermodynamics and kinetics of protein unfolding at equilibrium. These case studies illustrate typical data, analysis and results for each application and demonstrate a range of biological systems for which the interpretation of HXMS in terms of structure and conformational parameters provides unique insights into function. This article is part of a Special Issue entitled: Mass spectrometry in structural biology.     

Intramolecular hydrogen bonding in articaine can be related to superior bone tissue penetration: a molecular dynamics study

Local anesthetics (LAs) are drugs that cause reversible loss of nociception during surgical procedures. Articaine is a commonly used LA in dentistry that has proven to be exceptionally effective in penetrating bone tissue and induce anesthesia on posterior teeth in maxilla and mandibula. In the present study, our aim was to gain a deeper understanding of the penetration of articaine through biological membranes by studying the interactions of articaine with a phospholipid membrane. Our approach involves Langmuir monolayer experiments combined with molecular dynamics simulations. Membrane permeability of LAs can be modulated by pH due to a titratable amine group with a pKa value close to physiological pH. A change in protonation state is thus known to act as a lipophilicity switch in LAs. Our study shows that articaine has an additional unique lipophilicity switch in its ability to form an intramolecular hydrogen bond. We suggest this intramolecular hydrogen bond as a novel and additional solvent-dependent mechanism for modulation of lipophilicity of articaine which may enhance its diffusion through membranes and connective tissue.     

Intrinsic curvature in the VP1 gene of SV40: comparison of solution and gel results

DNA restriction fragments that are stably curved are usually identified by polyacrylamide gel electrophoresis because curved fragments migrate more slowly than normal fragments containing the same number of basepairs. In free solution, curved DNA molecules can be identified by transient electric birefringence (TEB) because they exhibit rotational relaxation times that are faster than those of normal fragments of the same size. In this article, the results observed in free solution and in polyacrylamide gels are compared for a highly curved 199-basepair (bp) restriction fragment taken from the VP1 gene in Simian Virus 40 (SV40) and various sequence mutants and insertion derivatives. The TEB method of overlapping fragments was used to show that the 199-bp fragment has an apparent bend angle of 46 +/- 2 degrees centered at sequence position 1922 +/- 2 bp. Four unphased A- and T-tracts and a mixed A3T4-tract occur within a span of approximately 60 bp surrounding the apparent bend center; for brevity, this 60-bp sequence element is called a curvature module. Modifying any of the A- or T-tracts in the curvature module by site-directed mutagenesis decreases the curvature of the fragment; replacing all five A- and T-tracts by random-sequence DNA causes the 199-bp mutant to adopt a normal conformation, with normal electrophoretic mobilities and birefringence relaxation times. Hence, stable curvature in this region of the VP1 gene is due to the five unphased A- and T- tracts surrounding the apparent bend center. Discordant solution and gel results are observed when long inverted repeats are inserted within the curvature module. These insertion derivatives migrate anomalously slowly in polyacrylamide gels but have normal, highly flexible conformations in free solution. Discordant solution and gel results are not observed if the insert does not contain a long inverted repeat or if the long inverted repeat is added to the 199-bp fragment outside the curvature module. The results suggest that long inverted repeats can form hairpins or cruciforms when they are located within a region of the helix backbone that is intrinsically curved, leading to large mobility anomalies in polyacrylamide gels. Hairpin/cruciform formation is not observed in free solution, presumably because of rapid conformational exchange. Hence, DNA restriction fragments that migrate anomalously slowly in polyacrylamide gels are not necessarily stably curved in free solution.     

Activation of PI3Kα by physiological effectors and by oncogenic mutations: structural and dynamic effects

PI3Kα, a heterodimeric lipid kinase, catalyzes the conversion of phosphoinositide-4,5-bisphosphate (PIP₂) to phosphoinositide-3,4,5-trisphosphate (PIP₃), a lipid that recruits to the plasma membrane proteins that regulate signaling cascades that control key cellular processes such as cell proliferation, carbohydrate metabolism, cell motility, and apoptosis. PI3Kα is composed of two subunits, p110α and p85, that are activated by binding to phosphorylated receptor tyrosine kinases (RTKs) or their substrates. The gene coding for p110α, PIK3CA, has been found to be mutated in a large number of tumors; these mutations result in increased PI3Kα kinase activity. The structure of the complex of p110α with a fragment of p85 containing the nSH2 and the iSH2 domains has provided valuable information about the mechanisms underlying the physiological activation of PI3Kα and its pathological activation by oncogenic mutations. This review discusses information derived from x-ray diffraction and theoretical calculations regarding the structural and dynamic effects of mutations in four highly mutated regions of PI3K p110α, as well as the proposed mechanisms by which these mutations increase kinase activity. During the physiological activation of PI3Kα, the phosphorylated tyrosine of RTKs binds to the nSH2 domain of p85, dislodging an inhibitory interaction between the p85 nSH2 and a loop of the helical domain of p110α. Several of the oncogenic mutations in p110α activate the enzyme by weakening this autoinhibitory interaction. These effects involve structural changes as well as changes in the dynamics of the enzyme. One of the most common p110α mutations, H1047R, activates PI3Kα by a different mechanism: it increases the interaction of the enzyme with the membrane, maximizing the access of the PI3Kα to its substrate PIP₂, a membrane lipid.