Showing papers in "arXiv: Biomolecules in 2011"

Crowded Charges in Ion Channels

Abstract: . Ions in water are the liquid of life. Life occurs almost entirely in ‘salt water’. Life began in salty oceans. Animals kept that salt water within them when they moved out of the ocean to drier surroundings. The plasma and blood that surrounds all cells are electrolytes more or less resembling sea water. The plasma inside cells is an electrolyte solution that more or less resembles the sea water in which life began. Water itself (without ions) is lethal to animal cells and damaging for most proteins. Water must contain the right ions in the right amounts if it is to sustain life. Physical chemistry is the language of electrolyte solutions and so physical chemistry, and biology, particularly physiology, have been intertwined since physical chemistry was developed some one hundred fifty years ago. Physiology, of course, was studied by the Greeks some millennia earlier, but the biological role of electrolyte solutionscould not be understood until ions were discovered by chemists some 2,000 years later.

Classes of fast and specific search mechanisms for proteins on DNA

Abstract: Problems of search and recognition appear over different scales in biological systems. In this review we focus on the challenges posed by interactions between proteins, in particular transcription factors, and DNA and possible mechanisms which allow for a fast and selective target location. Initially we argue that DNA-binding proteins can be classified, broadly, into three distinct classes which we illustrate using experimental data. Each class calls for a different search process and we discuss the possible application of different search mechanisms proposed over the years to each class. The main thrust of this review is a new mechanism which is based on barrier discrimination. We introduce the model and analyze in detail its consequences. It is shown that this mechanism applies to all classes of transcription factors and can lead to a fast and specific search. Moreover, it is shown that the mechanism has interesting transient features which allow for stability at the target despite rapid binding and unbinding of the transcription factor from the target.

Mechanisms of kinetic trapping in self-assembly and phase transformation

Abstract: In self-assembly processes, kinetic trapping effects often hinder the formation of thermodynamically stable ordered states. In a model of viral capsid assembly and in the phase transformation of a lattice gas, we show how simulations in a self-assembling steady state can be used to identify two distinct mechanisms of kinetic trapping. We argue that one of these mechanisms can be adequately captured by kinetic rate equations, while the other involves a breakdown of theories that rely on cluster size as a reaction coordinate. We discuss how these observations might be useful in designing and optimising self-assembly reactions.

Crowding Promotes the Switch from Hairpin to Pseudoknot Conformation in Human Telomerase RNA

Abstract: Formation of a pseudoknot in the conserved RNA core domain in the ribonucleoprotein human telomerase is required for function. In vitro experiments show that the pseudoknot (PK) is in equilibrium with an extended hairpin (HP) structure. We use molecular simulations of a coarse-grained model, which reproduces most of the salient features of the experimental melting profiles of PK and HP, to show that crowding enhances the stability of PK relative to HP in the wild type and in a mutant associated with dyskeratosis congenita. In monodisperse suspensions, small crowding particles increase the stability of compact structures to a greater extent than larger crowders. If the sizes of crowders in a binary mixture are smaller than the unfolded RNA, the increase in melting temperature due to the two components is additive. In a ternary mixture of crowders that are larger than the unfolded RNA, which mimics the composition of ribosome, large enzyme complexes and proteins in E. coli, the marginal increase in stability is entirely determined by the smallest component. We predict that crowding can restore partially telomerase activity in mutants, which dramatically decrease the PK stability.

Dominant Folding Pathways of a WW Domain

Abstract: We investigate the folding mechanism of the WW domain Fip35 using a realistic atomistic force field by applying the Dominant Reaction Pathways (DRP) approach. We find evidence for the existence of two folding pathways, which differ by the order of formation of the two hairpins. This result is consistent with the analysis of the experimental data on the folding kinetics of WW domains and with the results obtained from large-scale molecular dynamics (MD) simulations of this system. Free-energy calculations performed in two coarse-grained models support the robustness of our results and suggest that the qualitative structure of the dominant paths are mostly shaped by the native interactions. Computing a folding trajectory in atomistic detail only required about one hour on 48 CPU's. The gain in computational efficiency opens the door to a systematic investigation of the folding pathways of a large number of globular proteins.

Mechanical Characterization of Amyloid Fibrils Using Coarse-Grained Normal Mode Analysis

Abstract: Recent experimental studies have shown that amyloid fibril formed by aggregation of {\beta} peptide exhibits excellent mechanical properties comparable to other protein materials such as actin filaments and microtubules. These excellent mechanical properties of amyloid fibrils are related to their functional role in disease expression. This indicates the necessity to understand how an amyloid fibril achieves the remarkable mechanical properties through self-aggregation with structural hierarchy, highlighting the structure-property-function relationship for amyloids, whereas such relationship still remains elusive. In this work, we have studied the mechanical properties of human islet amyloid polypeptide (hIAPP) with respect to its structural hierarchies and structural shapes by coarse-grained normal mode analysis. Our simulation shows that hIAPP fibril can achieve the excellent bending rigidity via specific aggregation pattern such as antiparallel stacking of {\beta} peptides. Moreover, we have found the length-dependent mechanical properties of amyloids. This length-dependent property has been elucidated from Timoshenko beam model that takes into account the shear effect on the bending of amyloids. In summary, our study sheds light on the importance of not only the molecular architecture, which encodes the mechanical properties of the fibril, but also the shear effect on the mechanical (bending) behavior of the fibril.

Reconstructing an Icosahedral Virus from Single-Particle Diffraction Experiments

Abstract: The first experimental data from single-particle scattering experiments from free electron lasers (FELs) are now becoming available. The first such experiments are being performed on relatively large objects such as viruses, which produce relatively low-resolution, low-noise diffraction patterns in so-called "diffract-and-destroy" experiments. We describe a very simple test on the angular correlations of measured diffraction data to determine if the scattering is from an icosahedral particle. If this is confirmed, the efficient algorithm proposed can then combine diffraction data from multiple shots of particles in random unknown orientations to generate a full 3D image of the icosahedral particle. We demonstrate this with a simulation for the satellite tobacco necrosis virus (STNV), the atomic coordinates of whose asymmetric unit is given in Protein Data Bank entry 2BUK.

Medicinal Plants Database and Three Dimensional Structure of the Chemical Compounds from Medicinal Plants in Indonesia

Abstract: During this era of new drug designing, medicinal plants had become a very interesting object of further research. Pharmacology screening of active compound of medicinal plants would be time consuming and costly. Molecular docking is one of the in silico method which is more efficient compare to in vitro or in vivo method for its capability of finding the active compound in medicinal plants. In this method, three-dimensional structure becomes very important in the molecular docking methods, so we need a database that provides information on three-dimensional structures of chemical compounds from medicinal plants in Indonesia. Therefore, this study will prepare a database which provides information of the three dimensional structures of chemical compounds of medicinal plants. The database will be prepared by using MySQL format and is designed to be placed in http://herbaldb.farmasi.ui.ac.id website so that eventually this database can be accessed quickly and easily by users via the Internet.

Thermodynamics of Twisted DNA with Solvent Interaction

Abstract: The imaginary time path integral formalism is applied to a nonlinear Hamiltonian for a short fragment of heterogeneous DNA with a stabilizing solvent interaction term. Torsional effects are modeled by a twist angle between neighboring base pairs stacked along the molecule backbone. The base pair displacements are described by an ensemble of temperature dependent paths thus incorporating those fluctuational effects which shape the multisteps thermal denaturation. By summing over $\sim 10^7 - 10^8$ base pair paths, a large number of double helix configurations is taken into account consistently with the physical requirements of the model potential. The partition function is computed as a function of the twist. It is found that the equilibrium twist angle, peculiar of B-DNA at room temperature, yields the stablest helicoidal geometry against thermal disruption of the base pair hydrogen bonds. This result is corroborated by the computation of thermodynamical properties such as fractions of open base pairs and specific heat.

The nature of the low energy band of the Fenna-Matthews-Olson complex: vibronic signatures

Abstract: Based entirely upon actual experimental observations on electron-phonon coupling, we develop a theoretical framework to show that the lowest energy band of the Fenna- Matthews-Olson (FMO) complex exhibits observable features due to the quantum nature of the vibrational manifolds present in its chromophores. The study of linear spectra provides us with the basis to understand the dynamical features arising from the vibronic structure in non-linear spectra in a progressive fashion, starting from a microscopic model to finally performing an inhomogenous average. We show that the discreteness of the vibronic structure can be witnessed by probing the diagonal peaks of the non-linear spectra by means of a relative phase shift in the waiting time resolved signal. Moreover, we demonstrate the photon-echo and non-rephasing paths are sensitive to different harmonics in the vibrational manifold when static disorder is taken into account. Supported by analytical and numerical calculations, we show that nondiagonal resonances in the 2D spectra in the waiting time, further capture the discreteness of vibrations through a modulation of the amplitude without any effect in the signal intrinsic frequency. This fact generates a signal that is highly sensitive to correlations in the static disorder of the excitonic energy albeit protected against dephasing due to inhomogeneities of the vibrational ensemble.

Elastic Network Models: Theoretical and Empirical Foundations

Abstract: Fifteen years ago Monique Tirion showed that the low-frequency normal modes of a protein are not significantly altered when non-bonded interactions are replaced by Hookean springs, for all atom pairs whose distance is smaller than a given cutoff value. Since then, it has been shown that coarse-grained versions of Tirion's model are able to provide fair insights on many dynamical properties of biological macromolecules. In this text, theoretical tools required for studying these so-called Elastic Network Models are described, focusing on practical issues and, in particular, on possible artifacts. Then, an overview of some typical results that have been obtained by studying such models is given.

Monomer abundance distribution patterns as a universal biosignature: Examples from terrestrial and digital life

Abstract: Organisms leave a distinctive chemical signature in their environment because they synthesize those molecules that maximize their fitness. As a result, the relative concentrations of related chemical monomers in life-bearing environmental samples reflect, in part, those compounds' adaptive utility. In contrast, rates of molecular synthesis in a lifeless environment are dictated by reaction kinetics and thermodynamics, so concentrations of related monomers in abiotic samples tend to exhibit specific patterns dominated by small, easily formed, low-formation-energy molecules. We contend that this distinction can serve as a universal biosignature: the measurement of chemical concentration ratios that belie formation kinetics or equilibrium thermodynamics indicates the likely presence of life. We explore the features of this biosignature as observed in amino acids and carboxylic acids, using published data from numerous studies of terrestrial sediments, abiotic (spark, UV, and high-energy proton) synthesis experments, and meteorite bodies. We then compare these data to the results of experimental studies of an evolving digital life system. We observe the robust and repeatable evolution of an analogous biosignature in a digital lifeform, suggesting that evolutionary selection necessarily constrains organism composition and that the monomer abundance biosignature phenomenon is universal to evolved biosystems.

Subunit rotation in a single FoF1-ATP synthase in a living bacterium monitored by FRET

Abstract: FoF1-ATP synthase is the ubiquitous membrane-bound enzyme in mitochondria, chloroplasts and bacteria which provides the 'chemical energy currency' adenosine triphosphate (ATP) for cellular processes. In Escherichia coli ATP synthesis is driven by a proton motive force (PMF) comprising a proton concentration difference {\Delta}pH plus an electric potential {\Delta}{\Psi} across the lipid membrane. Single-molecule in vitro experiments have confirmed that proton-driven subunit rotation within FoF1-ATP synthase is associated with ATP synthesis. Based on intramolecular distance measurements by single-molecule fluorescence resonance energy transfer (FRET) the kinetics of subunit rotation and the step sizes of the different rotor parts have been unraveled. However, these experiments were accomplished in the presence of a PMF consisting of a maximum {\Delta}pH ~ 4 and an unknown {\Delta}{\Psi}. In contrast, in living bacteria the maximum {\Delta}pH across the plasma membrane is likely 0.75, and {\Delta}{\Psi} has been measured between -80 and -140 mV. Thus the problem of in vivo catalytic turnover rates, or the in vivo rotational speed in single FoF1-ATP synthases, respectively, has to be solved. In addition, the absolute number of functional enzymes in a single bacterium required to maintain the high ATP levels has to be determined. We report our progress of measuring subunit rotation in single FoF1-ATP synthases in vitro and in vivo, which was enabled by a new labeling approach for single-molecule FRET measurements.

Determining the DNA stability parameters for the breathing dynamics of heterogeneous DNA by stochastic optimization

Abstract: We suggest that the thermodynamic stability parameters (nearest neighbor stacking and hydrogen bonding free energies) of double-stranded DNA molecules can be inferred reliably from time series of the size fluctuations (breathing) of local denaturation zones (bubbles). On the basis of the reconstructed bubble size distribution, this is achieved through stochastic optimization of the free energies in terms of Simulated Annealing. In particular, it is shown that even noisy time series allow the identification of the stability parameters at remarkable accuracy. This method will be useful to obtain the DNA stacking and hydrogen bonding free energies from single bubble breathing assays rather than equilibrium data.

Temperature Dependence of Protein Folding Deduced from Quantum Transition

Abstract: A quantum theory on conformation-electron system is presented. Protein folding is regarded as the quantum transition between torsion states on polypeptide chain, and the folding rate is calculated by nonadiabatic operator method. The theory is used to study the temperature dependences of folding rate of 15 proteins and their non-Arrhenius behavior can all be deduced in a natural way. A general formula on the rate-temperature dependence has been deduced which is in good accordance with experimental data. These temperature dependences are further analyzed in terms of torsion potential parameters. Our results show it is necessary to move outside the realm of classical physics when the temperature dependence of protein folding is studied quantitatively.

Structure and Thermodynamics of Molecular Hydration via Grid Inhomogeneous Solvation Theory

Abstract: Changes in hydration are central to the phenomenon of biomolecular recognition, but it has been difficult to properly frame and answer questions about their precise thermodynamic role. We address this problem by introducing Grid Inhomogeneous Solvation Theory (GIST), which discretizes the equations of Inhomogeneous Solvation Theory on a 3D grid in a volume of interest. Here, the solvent volume is divided into small grid boxes and localized thermodynamic entropies, energies and free energies are defined for each grid box. Thermodynamic solvation quantities are defined in such a manner that summing the quantities over all the grid boxes yields the desired total quantity for the system. This approach smoothly accounts for the thermodynamics of not only highly occupied water sites but also partly occupied and water depleted regions of the solvent, without the need for ad hoc terms drawn from other theories. The GIST method has the further advantage of allowing a rigorous end-states analysis that, for example in the problem of molecular recognition, can account for not only the thermodynamics of displacing water from the surface but also for the thermodynamics of solvent reorganization around the bound complex. As a preliminary application, we present GIST calculations at the 1-body level for the host cucurbit[7]uril, a low molecular weight receptor molecule which represents a tractable model for biomolecular recognition. One of the most striking results is the observation of a toroidal region of water density, at the center of the host's nonpolar cavity, which is significantly disfavored entropically, and hence may contribute to the ability of this small receptor to bind guest molecules with unusually high affinities.

Molecular jamming - the cystine slipknot mechanical clamp in all-atom simulations

Abstract: A recent survey of 17 134 proteins has identified a new class of proteins which are expected to yield stretching induced force-peaks in the range of 1 nN. Such high force peaks should be due to forcing of a slip-loop through a cystine ring, i.e. by generating a cystine slipknot. The survey has been performed in a simple coarse grained model. Here, we perform all-atom steered molecular dynamics simulations on 15 cystine knot proteins and determine their resistance to stretching. In agreement with previous studies within a coarse grained structure based model, the level of resistance is found to be substantially higher than in proteins in which the mechanical clamp operates through shear. The large stretching forces arise through formation of the cystine slipknot mechanical clamp and the resulting steric jamming. We elucidate the workings of such a clamp in an atomic detail. We also study the behavior of five top strength proteins with the shear-based mechanostability in which no jamming is involved. We show that in the atomic model, the jamming state is relieved by moving one amino acid at a time and there is a choice in the selection of the amino acid that advances the first. In contrast, the coarse grained model also allows for a simultaneous passage of two amino acids.

Link between allosteric signal transduction and functional dynamics in a multi-subunit enzyme: S-adenosylhomocysteine hydrolase

Abstract: S-adenosylhomocysteine hydrolase (SAHH), a cellular enzyme that plays a key role in methylation reactions including those required for maturation of viral mRNA, is an important drug target in the discovery of antiviral agents. While targeting the active site is a straightforward strategy of enzyme inhibition, evidences of allosteric modulation of active site in many enzymes underscore the molecular origin of signal transduction. Information of co-evolving sequences in SAHH family and the key residues for functional dynamics that can be identified using native topology of the enzyme provide glimpses into how the allosteric signaling network, dispersed over the molecular structure, coordinates intra- and inter-subunit conformational dynamics. To study the link between the allosteric communication and functional dynamics of SAHHs, we performed Brownian dynamics simulations by building a coarse-grained model based on the holo and ligand-bound structures. The simulations of ligand-induced transition revealed that the signal of intra-subunit closure dynamics is transmitted to form inter-subunit contacts, which in turn invoke a precise alignment of active site, followed by the dimer-dimer rotation that compacts the whole tetrameric structure. Further analyses of SAHH dynamics associated with ligand binding provided evidence of both induced fit and population shift mechanisms, and also showed that the transition state ensemble is akin to the ligand-bound state. Besides the formation of enzyme-ligand contacts at the active site, the allosteric couplings from the residues distal to the active site is vital to the enzymatic function.

Comprehensive data-driven analysis of the impact of chemoinformatic structure on the genome-wide biological response profiles of cancer cells to 1159 drugs

Abstract: Detailed and systematic understanding of the biological effects of millions of available compounds on living cells is a significant challenge. As most compounds impact multiple targets and pathways, traditional methods for analyzing structure-function relationships are not comprehensive enough. Therefore more advanced integrative models are needed for predicting biological effects elicited by specific chemical features. As a step towards creating such computational links we developed a data-driven chemical systems biology approach to comprehensively study the relationship of 76 structural 3D-descriptors (VolSurf, chemical space) of 1159 drugs with the gene expression responses (biological space) they elicited in three cancer cell lines. The analysis covering 11350 genes was based on data from the Connectivity Map. We decomposed these biological response profiles into components, each linked to a characteristic chemical descriptor profile. The integrated quantitative analysis of the chemical and biological spaces was more informative about protein-target based drug similarity than either dataset separately. We identified ten major components that link distinct VolSurf features across multiple compounds to specific biological activity types. For example, component 2 (hydrophobic properties) strongly links to DNA damage response, while component 3 (hydrogen bonding) connects to metabolic stress. Individual structural and biological features were often linked to one cell line only, such as leukemia cells (HL-60) specifically responding to cardiac glycosides. In summary, our approach identified specific chemical structures shared across multiple drugs causing distinct biological responses. The decoding of such systematic chemical-biological relationships is necessary to build better models of drug effects, including unidentified types of molecular properties with strong biological effects.

Single-image diffusion coefficient measurements of proteins in free solution

Abstract: Diffusion coefficient measurements are important for many biological and material investigations, such as particle dynamics, kinetics, and size determinations. Amongst current measurement methods, single particle tracking (SPT) offers the unique capability of providing location and diffusion information of a molecule simultaneously while using only femptomoles of sample. However, the temporal resolution of SPT is limited to seconds for single-color labeled samples. By directly imaging three dimensional (3D) diffusing fluorescent proteins and studying the widths of their intensity profiles, we determine the proteins' diffusion coefficients using single protein images of sub-millisecond exposure times. This simple method improves the temporal resolution of diffusion coefficient measurements to sub-millisecond, and can be readily applied to a range of particle sizes in SPT investigations and applications where diffusion coefficient measurements are needed, such as reaction kinetics and particle size determinations.

Analysis of the Equilibrium and Kinetics of the Ankyrin Repeat Protein Myotrophin

Abstract: We apply the Wako-Saito-Munoz-Eaton model to the study of Myotrophin, a small ankyrin repeat protein, whose folding equilibrium and kinetics have been recently characterized experimentally. The model, which is a native-centric with binary variables, provides a finer microscopic detail than the Ising model, that has been recently applied to some different repeat proteins, while being still amenable for an exact solution. In partial agreement with the experiments, our results reveal a weakly three-state equilibrium and a two-state-like kinetics of the wild type protein despite the presence of a non-trivial free-energy profile. These features appear to be related to a careful "design" of the free-energy landscape, so that mutations can alter this picture, stabilizing some intermediates and changing the position of the rate-limiting step. Also the experimental findings of two alternative pathways, an N-terminal and a C-terminal one, are qualitatively confirmed, even if the variations in the rates upon the experimental mutations cannot be quantitatively reproduced. Interestingly, folding and unfolding pathway appear to be different, even if closely related: a property that is not generally considered in the phenomenological interpretation of the experimental data.

Life's Solutions are Not Ideal

Abstract: Life occurs in ionic solutions, not pure water. The ionic mixtures of these solutions are very different from water and have dramatic effects on the cells and molecules of biological systems, yet theories and simulations cannot calculate their properties. I suggest the reason is that existing theories stem from the classical theory of ideal or simple gases in which (to a first approximation) atoms do not interact. Even the law of mass action describes reactants as if they were ideal. I propose that theories of ionic solutions should start with the theory of complex fluids because that theory is designed to deal with interactions from the beginning. The variational theory of complex fluids is particularly well suited to describe mixtures like the solutions in and outside biological cells. When a component or force is added to a solution, the theory derives - by mathematics alone - a set of partial differential equations that captures the resulting interactions self-consistently. Such a theory has been implemented and shown to be computable in biologically relevant systems but it has not yet been thoroughly tested in equilibrium or flow.

Channels as Enzymes: Oxymoron and Tautology

Abstract: Channels are Catalysts for Diffusion and in that sense are Enzymes. This idea is useful for the design and interpretation of experiments.

Nanodiamond as a vector for siRNA delivery to Ewing sarcoma cells

Abstract: We investigated the ability of diamond nanoparticles (nanodiamonds, NDs) to deliver small interfering RNA (siRNA) in Ewing sarcoma cells, in the perspective of in vivo anti-cancer nucleic acid drug delivery. siRNA was adsorbed onto NDs previously coated with cationic polymer. Cell uptake of NDs has been demonstrated by taking advantage of NDs intrinsic fluorescence coming from embedded color center defects. Cell toxicity of these coated NDs was shown to be low. Consistent with the internalization efficacy, we have shown a specific inhibition of EWS/Fli-1 gene expression at the mRNA and protein level by the ND vectorized siRNA in a serum containing medium.

Protein Photo-folding and Quantum Folding Theory

Abstract: The rates of protein folding with photon absorption or emission and the cross section of photon -protein inelastic scattering are calculated from the quantum folding theory by use of standard field-theoretical method. All these protein photo-folding processes are compared with common protein folding without interaction of photons (nonradiative folding). It is demonstrated that there exists a common factor (thermo-averaged overlap integral of vibration wave function, TAOI) for protein folding and protein photo-folding. Based on this finding it is predicted that: 1) the stimulated photo-folding rates show the same temperature dependence as protein folding; 2) the spectral line of electronic transition is broadened to a band which includes abundant vibration spectrum without and with conformational transition and the width of the vibration spectral line is largely reduced; 3) the resonance fluorescence cross section changes with temperature obeying the same law (Luo-Lu's law). The particular form of the folding rate - temperature relation and the abundant spectral structure imply the existence of a set of quantum oscillators in the transition process and these oscillators are mainly of torsion type of low frequency, imply the quantum tunneling between protein conformations does exist in folding and photo-folding processes and the tunneling is rooted deeply in the coherent motion of the conformational-electronic system.

Quantum Proteomics

Abstract: We put forward the idea of establishing a novel interdisciplinary field of research at the interface between quantum mechanics and proteomics The new field, called quantum proteomics, is defined as the large-scale study of the electronic structure of the proteins that define an organism's proteome The electronic structure of proteins is unveiled with the aid of linear-scaling quantum mechanical calculations Such calculations provide information about the energy levels of the proteins, the charges of their amino acid side chains, their electrostatic potentials and permanent dipole moments ({\mu}) Since the magnitude of the electric dipole moment of any protein is not null ({\mu} eq0 Debye), the dipole moment can be employed to characterize the electronic structure of each protein that belongs to an organism's proteome As an example, we investigate six proteins from the thermophilic bacterium Methanobacterium thermoautotrophicum (Mth) whose atomic structures were characterized by solution NMR spectroscopy

The Boltzmann factor, DNA melting, and Brownian ratchets: Topics in an introductory physics sequence for biology and premedical students

Abstract: Three, interrelated biologically-relevant examples of biased random walks are presented: (1) A model for DNA melting, modelled as DNA unzipping, which provides a way to illustrate the role of the Boltzmann factor in a venue well-known to biology and pre-medical students; (2) the activity of helicase motor proteins in unzipping double-stranded DNA, for example, at the replication fork, which is an example of a Brownian ratchet; (3) force generation by actin polymerization, which is another Brownian ratchet, and for which the force and actin-concentration dependence of the velocity of actin polymerization is determined.

Nonequilibrium effects in DNA microarrays: a multiplatform study

Abstract: It has recently been shown that in some DNA microarrays the time needed to reach thermal equilibrium may largely exceed the typical experimental time, which is about 15h in standard protocols (Hooyberghs et al. Phys. Rev. E 81, 012901 (2010)). In this paper we discuss how this breakdown of thermodynamic equilibrium could be detected in microarray experiments without resorting to real time hybridization data, which are difficult to implement in standard experimental conditions. The method is based on the analysis of the distribution of fluorescence intensities I from different spots for probes carrying base mismatches. In thermal equilibrium and at sufficiently low concentrations, log I is expected to be linearly related to the hybridization free energy $\Delta G$ with a slope equal to $1/RT_{exp}$, where $T_{exp}$ is the experimental temperature and R is the gas constant. The breakdown of equilibrium results in the deviation from this law. A model for hybridization kinetics explaining the observed experimental behavior is discussed, the so-called 3-state model. It predicts that deviations from equilibrium yield a proportionality of $\log I$ to $\Delta G/RT_{eff}$. Here, $T_{eff}$ is an effective temperature, higher than the experimental one. This behavior is indeed observed in some experiments on Agilent arrays. We analyze experimental data from two other microarray platforms and discuss, on the basis of the results, the attainment of equilibrium in these cases. Interestingly, the same 3-state model predicts a (dynamical) saturation of the signal at values below the expected one at equilibrium.

Drug transport mechanism of P-glycoprotein monitored by single molecule fluorescence resonance energy transfer

Abstract: In this work we monitor the catalytic mechanism of P-glycoprotein (Pgp) using single-molecule fluorescence resonance energy transfer (FRET). Pgp, a member of the ATP binding cassette family of transport proteins, is found in the plasma membrane of animal cells where it is involved in the ATP hydrolysis driven export of hydrophobic molecules. When expressed in the plasma membrane of cancer cells, the transport activity of Pgp can lead to the failure of chemotherapy by excluding the mostly hydrophobic drugs from the interior of the cell. Despite ongoing effort, the catalytic mechanism by which Pgp couples MgATP binding and hydrolysis to translocation of drug molecules across the lipid bilayer is poorly understood. Using site directed mutagenesis, we have introduced cysteine residues for fluorescence labeling into different regions of the nucleotide binding domains (NBDs) of Pgp. Double-labeled single Pgp molecules showed fluctuating FRET efficiencies during drug stimulated ATP hydrolysis suggesting that the NBDs undergo significant movements during catalysis. Duty cycle-optimized alternating laser excitation (DCO-ALEX) is applied to minimize FRET artifacts and to select the appropriate molecules. The data show that Pgp is a highly dynamic enzyme that appears to fluctuate between at least two major conformations during steady state turnover.

3D Protein Structure Predicted from Sequence

Abstract: The evolutionary trajectory of a protein through sequence space is constrained by function and three-dimensional (3D) structure. Residues in spatial proximity tend to co-evolve, yet attempts to invert the evolutionary record to identify these constraints and use them to computationally fold proteins have so far been unsuccessful. Here, we show that co-variation of residue pairs, observed in a large protein family, provides sufficient information to determine 3D protein structure. Using a data-constrained maximum entropy model of the multiple sequence alignment, we identify pairs of statistically coupled residue positions which are expected to be close in the protein fold, termed contacts inferred from evolutionary information (EICs). To assess the amount of information about the protein fold contained in these coupled pairs, we evaluate the accuracy of predicted 3D structures for proteins of 50-260 residues, from 15 diverse protein families, including a G-protein coupled receptor. These structure predictions are de novo, i.e., they do not use homology modeling or sequence-similar fragments from known structures. The resulting low C{\alpha}-RMSD error range of 2.7-5.1A, over at least 75% of the protein, indicates the potential for predicting essentially correct 3D structures for the thousands of protein families that have no known structure, provided they include a sufficiently large number of divergent sample sequences. With the current enormous growth in sequence information based on new sequencing technology, this opens the door to a comprehensive survey of protein 3D structures, including many not currently accessible to the experimental methods of structural genomics. This advance has potential applications in many biological contexts, such as synthetic biology, identification of functional sites in proteins and interpretation of the functional impact of genetic variants.