Showing papers in "arXiv: Biomolecules in 2005"

Quantifying Influenza Vaccine Efficacy and Antigenic Distance

Abstract: We introduce a new measure of antigenic distance between influenza A vaccine and circulating strains. The measure correlates well with efficacies of the H3N2 influenza A component of the annual vaccine between 1971 and 2004, as do results of a theory of the immune response to influenza following vaccination. This new measure of antigenic distance is correlated with vaccine efficacy to a greater degree than are current state-of-the-art phylogenetic sequence analyzes or ferret antisera inhibition assays. We suggest that this new measure of antigenic distance be used in the design of the annual influenza vaccine and in the interpretation of vaccine efficacy monitoring.

Functional Dynamics of PDZ Binding Domains: A Normal Mode Analysis

Abstract: PDZ (Post-synaptic density-95/discs large/zonula occludens-1) domains are relatively small (80 to 120 residues) protein binding modules central in the organization of receptor clusters and in the association of cellular proteins. Their main function is to bind C-terminals of selected proteins that are recognized through specific amino-acids in their carboxyl end. Binding is associated with a deformation of the PDZ native structure and is responsible for dynamical changes in regions not in direct contact with the target. We investigate how this deformation is related to the harmonic dynamics of the PDZ structure and show that one low-frequency collective normal mode, characterized by the concerted movements of different secondary structures, is involved in the binding process. Our results suggest that even minimal structural changes are responsible of communication between distant regions of the protein, in agreement with recent Nuclear Magnetic Resonance (NMR) experiments. Thus PDZ domains are a very clear example of how collective normal modes are able to characterize the relation between function and dynamics of proteins, and to provide indications on the precursors of binding/unbonding events.

Master equation approach to the assembly of viral capsids

Abstract: The distribution of inequivalent geometries occurring during self-assembly of the major capsid protein in thermodynamic equilibrium is determined based on a master equation approach. These results are implemented to characterize the assembly of SV40 virus and to obtain information on the putative pathways controlling the progressive build-up of the SV40 capsid. The experimental testability of the predictions is assessed and an analysis of the geometries of the assembly intermediates on the dominant pathways is used to identify targets for antiviral drug design.

Directed motion emerging from two coupled random processes: Translocation of a chain through a membrane nanopore driven by binding proteins

Abstract: We investigate the translocation of a stiff polymer consisting of M monomers through a nanopore in a membrane, in the presence of binding particles (chaperones) that bind onto the polymer, and partially prevent backsliding of the polymer through the pore. The process is characterized by the rates: k for the polymer to make a diffusive jump through the pore, q for unbinding of a chaperone, and the rate q kappa for binding (with a binding strength kappa); except for the case of no binding kappa=0 the presence of the chaperones give rise to an effective force that drives the translocation process. Based on a (2+1) variate master equation, we study in detail the coupled dynamics of diffusive translocation and (partial) rectification by the binding proteins. In particular, we calculate the mean translocation time as a function of the various physical parameters.

Activation and Protonation of Dinitrogen at the FeMo-Cofactor of Nitrogenase

Abstract: The protonation of N2 bound to the active center of nitrogenase has been investigated using state-of-the-art DFT calculations. Dinitrogen in the bridging mode is activated by forming two bonds to Fe sites, which results in a reduction of the energy for the first hydrogen transfer by 123 kJ/mol. The axial binding mode with open sulfur bridge is less reactive by 30 kJ/mol and the energetic ordering of the axial and bridged binding mode is reversed in favor of the bridging dinitrogen during the first protonation. Protonation of the central ligand is thermodynamically favorable but kinetically hindered. If the central ligand is protonated, the proton is transferred to dinitrogen following the second protonation. Protonation of dinitrogen at the Mo site does not lead to low-energy intermediates.

Equation of state of looped DNA

Abstract: We calculate the equation of state of DNA under tension for the case that the DNA features loops. Such loops occur transiently during DNA condensation in the presence of multivalent ions or sliding cationic protein linkers. The force-extension relation of such looped DNA modelled as a wormlike chain is calculated via path integration in the semiclassical limit. This allows us to determine rigorously the high stretching asymptotics. Notably the functional form of the force-extension curve resembles that of straight DNA, yet with a strongly renormalized apparent persistence length. That means that the experimentally extracted single molecule elasticity does not necessarily reflect the bare DNA stiffness only, but can also contain additional contributions that depend on the overall chain conformation and length.

Controlled DNA compaction within chromatin: the tail-bridging effect

Abstract: We study the mechanism underlying the attraction between nucleosomes, the fundamental packaging units of DNA inside the chromatin complex. We introduce a simple model of the nucleosome, the eight-tail colloid, consisting of a charged sphere with eight oppositely charged, flexible, grafted chains that represent the terminal histone tails. We demonstrate that our complexes are attracted via the formation of chain bridges and that this attraction can be tuned by changing the fraction of charged monomers on the tails. This suggests a physical mechanism of chromatin compaction where the degree of DNA condensation can be controlled via biochemical means, namely the acetylation and deacetylation of lysines in the histone tails.

Stretching of proteins in a force-clamp

Abstract: Mechanically induced protein unfolding in the force-clamp apparatus is shown, in a coarse-grained model of ubiquitin, to have lognormal statistics above a treshold force and exponential below it. Correspondingly, the mean unfolding time is slowly varying and exponentially decreases as a function of the force. The time dependencies of the end-to-end distances are also distinct. The time sequence of unfolding events weakly depends on force and much of it resembles that for stretching at constant speed. A more complicated time dependence arises for integrin.

Deconstructing the energy landscape: new algorithms for folding heteropolymers

Abstract: We apply the computational methodology of phase retrieval to the problem of folding heteropolymers. The ground state fold of the polymer is defined by the intersection of two sets in the configuration space of its constituent monomers: a geometrical chain constraint and a threshold constraint on the contact energy. A dynamical system is then defined in terms of the projections to these constraint sets, such that its fixed points solve the set intersection problem. We present results for two off-lattice HP models: one with only rotameric degrees of freedom, and one proposed by Stillinger et al. with flexible bond angles. Our phase retrieval inspired algorithm is competitive with more established algorithms and even finds lower energy folds for one of the longer polymer chains.

Predicting Secondary Structures, Contact Numbers, and Residue-wise Contact Orders of Native Protein Structure from Amino Acid Sequence by Critical Random Networks

Abstract: Prediction of one-dimensional protein structures such as secondary structures and contact numbers is useful for the three-dimensional structure prediction and important for the understanding of sequence-structure relationship. Here we present a new machine-learning method, critical random networks (CRNs), for predicting one-dimensional structures, and apply it, with position-specific scoring matrices, to the prediction of secondary structures (SS), contact numbers (CN), and residue-wise contact orders (RWCO). The present method achieves, on average, $Q_3$ accuracy of 77.8% for SS, correlation coefficients of 0.726 and 0.601 for CN and RWCO, respectively. The accuracy of the SS prediction is comparable to other state-of-the-art methods, and that of the CN prediction is a significant improvement over previous methods. We give a detailed formulation of critical random networks-based prediction scheme, and examine the context-dependence of prediction accuracies. In order to study the nonlinear and multi-body effects, we compare the CRNs-based method with a purely linear method based on position-specific scoring matrices. Although not superior to the CRNs-based method, the surprisingly good accuracy achieved by the linear method highlights the difficulty in extracting structural features of higher order from amino acid sequence beyond that provided by the position-specific scoring matrices.

Are structural biases at protein termini a signature of vectorial folding

Abstract: Experimental investigations of the biosynthesis of a number of proteins have pointed out that part of the native structure can be acquired already during translation. We carried out a comprehensive statistical analysis of some average structural properties of proteins that have been put forward as possible signatures of this progressive buildup process. Contrary to a widespread belief, it is found that there is no major propensity of the amino acids to form contacts with residues that are closer to the N terminus. Moreover, it is found that the C terminus is significantly more compact and locally-organized than the N one. Also this bias, though, is unlikely to be related to vectorial effects, since it correlates with subtle differences in the primary sequence. These findings indicate that even if proteins aquire their structure vectorially no signature of this seems to be detectable in their average structural properties.

A new series of polyhedra as blueprints for viral capsids in the family of Papovaviridae

Abstract: In a seminal paper Caspar and Klug established a theory that provides a family of polyhedra as blueprints for the structural organisation of viral capsids In particular, they encode the locations of the proteins in the shells that encapsulate, and hence provide protection for, the viral genome Despite of its huge success and numerous applications in virology experimental results have provided evidence for the fact that the theory is too restrictive to describe all known viruses Especially, the family of Papovaviridae, which contains cancer-causing viruses, falls out of the scope of this theory In a recent paper we have shown that certain members of the family of Papovaviridae can be described via tilings In this paper, we develop a comprehensive mathematical framework for the derivation of all surface structures of viral particles in this family We show that this formalism fixes the structure and relative sizes of all particles collectively so that there exists only one scaling factor that relates the sizes of all particles with their biological counterparts The series of polyhedra derived here complements the Caspar-Klug family of polyhedra It is the first mathematical result that provides a common organisational principle for different types of viral particles in the family of Papovaviridae and paves the way for an understanding of Papovaviridae polymorphism Moreover, it provides crucial input for the construction of assembly models

Vibrational energy relaxation (VER) of a CD stretching mode in cytochrome c

Abstract: We first review how to determine the rate of vibrational energy relaxation (VER) using perturbation theory. We then apply those theoretical results to the problem of VER of a CD stretching mode in the protein cytochrome c. We model cytochrome c in vacuum as a normal mode system with the lowest-order anharmonic coupling elements. We find that, for the “lifetime” width parameter = 3 ∼ 30 cm 1 , the VER time is 0.2 ∼ 0.3 ps, which agrees rather well with the previous classical calculation using the quantum correction factor method, and is consistent with spectroscopic experiments by Romesberg’s group. We decompose the VER rate into separate contributions from two modes, and find that the most significant contribution, which depends on the “lifetime” width parameter, comes from those modes most resonant with the CD vibrational mode.

Proteins and polymers

Abstract: Proteins, chain molecules of amino acids, behave in ways which are similar to each other yet quite distinct from standard compact polymers. We demonstrate that the Flory theorem, derived for polymer melts, holds for compact protein native state structures and is not incompatible with the existence of structured building blocks such as $\alpha$-helices and $\beta$-strands. We present a discussion on how the notion of the thickness of a polymer chain, besides being useful in describing a chain molecule in the continuum limit, plays a vital role in interpolating between conventional polymer physics and the phase of matter associated with protein structures.

Predicting Residue-wise Contact Orders of Native Protein Structure from Amino Acid Sequence

Abstract: Residue-wise contact order (RWCO) is a new kind of one-dimensional protein structures which represents the extent of long-range contacts. We have recently shown that a set of three types of one-dimensional structures (secondary structure, contact number, and RWCO) contains sufficient information for reconstructing the three-dimensional structure of proteins. Currently, there exist prediction methods for secondary structure and contact number from amino acid sequence, but none exists for RWCO. Also, the properties of amino acids that affect RWCO is not clearly understood. Here, we present a linear regression-based method to predict RWCO from amino acid sequence, and analyze the regression parameters to identify the properties that correlates with the RWCO. The present method achieves the significant correlation of 0.59 between the native and predicted RWCOs on average. An unusual feature of the RWCO prediction is the remarkably large optimal half window size of 26 residues. The regression parameters for the central and near-central residues of the local sequence segment highly correlate with those of the contact number prediction, and hence with hydrophobicity.

Correlations of Amino Acids with Secondary Structure Types: Connection with Amino Acid Structure

Abstract: The correlations of primary and secondary structures were analyzed using proteins with known structure from Protein Data Bank The correlation values of amino acid type and the eight secondary structure types at distant position were calculated for distances between -25 and 25 Shapes of the diagrams indicate that amino acids polarity and capability for hydrogen bonding have influence on the secondary structure at some distances Clear preference of most of the amino acids towards certain secondary structure type classifies amino acids into four groups: alpha-helix admirers, strand admirers, turn and bend admirers and the others Group four consists of His and Cis, the amino acids that do not show clear preference for any secondary structure Amino acids from a group have similar physicochemical properties, and the same structural characteristics The results suggest that amino acid preference for secondary structure type is based on the structural characteristics at Cb and Cg atoms of amino acid alpha-helix admirers do not have polar heteroatoms on Cb and Cg atoms, nor branching or aromatic group on Cb atom Amino acids that have aromatic groups or branching on Cb atom are strand admirers Turn and bend admirers have polar heteroatom on Cb or Cg atoms or do not have Cb atom at all Our results indicate that polarity and capability for hydrogen bonding have influence on the secondary structure at some distance, and that amino acid preference for secondary structure is caused by structural properties at Cb or Cg atoms

A Hydrophobic Gate in an Ion Channel: The Closed State of the Nicotinic Acetylcholine Receptor

Abstract: The nicotinic acetylcholine receptor (nAChR) is the prototypic member of the `Cys-loop' superfamily of ligand-gated ion channels which mediate synaptic neurotransmission, and whose other members include receptors for glycine, gamma-aminobutyric acid, and serotonin. Cryo-electron microscopy has yielded a three dimensional structure of the nAChR in its closed state. However, the exact nature and location of the channel gate remains uncertain. Although the transmembrane pore is constricted close to its center, it is not completely occluded. Rather, the pore has a central hydrophobic zone of radius about 3 A. Model calculations suggest that such a constriction may form a hydrophobic gate, preventing movement of ions through a channel. We present a detailed and quantitative simulation study of the hydrophobic gating model of the nicotinic receptor, in order to fully evaluate this hypothesis. We demonstrate that the hydrophobic constriction of the nAChR pore indeed forms a closed gate. Potential of mean force (PMF) calculations reveal that the constriction presents a barrier of height ca. 10 kT to the permeation of sodium ions, placing an upper bound on the closed channel conductance of 0.3 pS. Thus, a 3 A radius hydrophobic pore can form a functional barrier to the permeation of a 1 A radius Na+ ion. Using a united atom force field for the protein instead of an all atom one retains the qualitative features but results in differing conductances, showing that the PMF is sensitive to the detailed molecular interactions.

Chemical and Structural Disorder in Eumelanins - A Possible Explanation for Broad Band Absorbance

Abstract: We report the results of an experimental and theoretical study of the electronic and structural properties of a key eumelanin precursor - 5,6,-dihydroxyindole-2-carboxylic acid (DHICA) and its dimeric forms. We have used optical spectroscopy to follow the oxidative polymerization of DHICA to eumelanin, and observe red shifting and broadening of the absorption spectrum as the reaction proceeds. First principles density functional theory calculations indicate that DHICA oligomers (possible reaction products of oxidative polymerization) have red shifted HOMO-LUMO gaps with respect to the monomer. Furthermore, different bonding configurations (leading to oligomers with different structures) produce a range of gaps. These experimental and theoretical results lend support to the chemical disorder model where the broad band monotonic absorption characteristic of all melanins is a consequence of the superposition of a large number of inhomogeneously broadened Gaussian transitions associated with each of the components of a melanin ensemble. These results suggest that the traditional model of eumelanin as an amorphous organic semiconductor is not required to explain its optical properties, and should be thoroughly re-examined. These results have significant implications for our understanding of the physics, chemistry and biological function of these important biological macromolecules. Indeed, one may speculate that the robust functionality of melanins in vitro is a direct consequence of its heterogeneity, i.e. chemical disorder is a "low cost" natural resource in these systems.

Design of HIV-1-PR inhibitors which do not create resistance: blocking the folding of single monomers

Abstract: One of the main problems of drug design is that of optimizing the drug--target interaction. In the case in which the target is a viral protein displaying a high mutation rate, a second problem arises, namely the eventual development of resistance. We wish to suggest a scheme for the design of non--conventional drugs which do not face any of these problems and apply it to the case of HIV--1 protease. It is based on the knowledge that the folding of single--domain proteins, like e.g. each of the monomers forming the HIV--1--PR homodimer, is controlled by local elementary structures (LES), stabilized by local contacts among hydrophobic, strongly interacting and highly conserved amino acids which play a central role in the folding process. Because LES have evolved over myriads of generations to recognize and strongly interact with each other so as to make the protein fold fast as well as to avoid aggregation with other proteins, highly specific (and thus little toxic) as well as effective folding--inhibitor drugs suggest themselves: short peptides (or eventually their mimetic molecules), displaying the same amino acid sequence of that of LES (p--LES). Aside from being specific and efficient, these inhibitors are expected not to induce resistance: in fact, mutations which successfully avoid their action imply the destabilization of one or more LES and thus should lead to protein denaturation. Making use of Monte Carlo simulations within the framework of a simple although not oversimplified model, which is able to reproduce the main thermodynamic as well as dynamic properties of monoglobular proteins, we first identify the LES of the HIV--1--PR and then show that the corresponding p--LES peptides act as effective inhibitors of the folding of the protease which do not create resistance.

A toolkit for the construction of icosahedral particles with local symmetry axes

Abstract: A formalism is developed which allows to determine the locations of all local symmetry axes of three-dimensional particles with overall icosahedral symmetry. It relies on the fact that the root system of the non-crystallographic Coxeter group H_3 encodes the locations of the planes of reflection that generate the discrete rotational symmetries of the particles. Via an appropriate extension of the root system, new planes of reflection are introduced which determine local axes of rotational symmetry. An easy-to-implement formalism is derived that allows to compute the surface structure of any three-dimensional icosahedral particle with local symmetries. It can be used also for particles with overall octahedral and tetrahedral symmetry in conjunction with the root systems of the corresponding reflection groups. Applications to viruses are discussed explicitly. It is shown that the concept of quasi-equivalence in Caspar-Klug Theory corresponds to the special case of local six-fold symmetry axes contained in the theory developed here, and the corresponding geometries can hence be obtained with this formalism based on the root system of H_3. Moreover, as a by-product, the theory answers the long-standing open question why only certain types of capsomeres, i.e. clusters of protein subunits, are observed in the surface structures of viruses. Since the types of the capsomeres are determined by the orders of the local symmetry axes on which they are located, the possible types of capsomeres are restricted by the spectrum of local symmetry axes allowed by the theory. Based on this we determine the spectrum of all capsomere types that may occur in viral capsids and give explicit examples for the lower-order cases.

Uni-directional polymerization leading to homochirality in the RNA world

Abstract: The differences between uni-directional and bi-directional polymerization are considered. The uni-directional case is discussed in the framework of the RNA world. Similar to earlier models of this type, where polymerization was assumed to proceed in a bi-directional fashion (presumed to be relevant to peptide nucleic acids), left-handed and right-handed monomers are produced via an autocatalysis from an achiral substrate. The details of the bifurcation from a racemic solution to a homochiral state of either handedness is shown to be remarkably independent of whether the polymerization in uni-directional or bi-directional. Slightly larger differences are seen when dissociation is allowed and the dissociation fragments are being recycled into the achiral substrate.

The posterior-Viterbi: a new decoding algorithm for hidden Markov models

Abstract: Background: Hidden Markov models (HMM) are powerful machine learning tools successfully applied to problems of computational Molecular Biology. In a predictive task, the HMM is endowed with a decoding algorithm in order to assign the most probable state path, and in turn the class labeling, to an unknown sequence. The Viterbi and the posterior decoding algorithms are the most common. The former is very efficient when one path dominates, while the latter, even though does not guarantee to preserve the automaton grammar, is more effective when several concurring paths have similar probabilities. A third good alternative is 1-best, which was shown to perform equal or better than Viterbi. Results: In this paper we introduce the posterior-Viterbi (PV) a new decoding which combines the posterior and Viterbi algorithms. PV is a two step process: first the posterior probability of each state is computed and then the best posterior allowed path through the model is evaluated by a Viterbi algorithm. Conclusions: We show that PV decoding performs better than other algorithms first on toy models and then on the computational biological problem of the prediction of the topology of beta-barrel membrane proteins.

The effect of genome length on ejection forces in bacteriophage lambda

Abstract: A variety of viruses tightly pack their genetic material into protein capsids that are barely large enough to enclose the genome. In particular, in bacteriophages, forces as high as 60 pN are encountered during packaging and ejection, produced by DNA bending elasticity and self-interactions. The high forces are believed to be important for the ejection process, though the extent of their involvement is not yet clear. As a result, there is a need for quantitative models and experiments that reveal the nature of the forces relevant to DNA ejection. Here we report measurements of the ejection forces for two different mutants of bacteriophage lambda, lambda b221cI26 and lambda cI60, which differ in genome length by ~30%. As expected for a force-driven ejection mechanism, the osmotic pressure at which DNA release is completely inhibited varies with the genome length: we find inhibition pressures of 15 atm and 25 atm, respectively, values that are in agreement with our theoretical calculations.

Heisenberg uncertainty principle for thermal response of the microtubules excited by ultra-short laser pulses

Abstract: In this paper the heat signaling in microtubules (MT) is investigated It is argued that for the description of the heat signaling phenomena in MT, the hyperbolic heat transport (HHT) equation must be used It is shown that HHT is the Klein-Gordon (K-G) equation The general solution for the K-G equation for MT is obtained For the undistorted signal propagation in MT the Heisenberg uncertainty principle is formulated and discussed Key words: Microtubules; Heat signaling; Klein-Gordon equation; Heisenberg principle

Symmetry and Dynamics in living organisms: The self-similarity principle governs gene expression dynamics

Abstract: The ambitious and ultimate research purpose in Systems Biology is the understanding and modelling of the cell's system. Although a vast number of models have been developed in order to extract biological knowledge from complex systems composed of basic elements as proteins, genes and chemical compounds, a need remains for improving our understanding of dynamical features of the systems (i.e., temporal-dependence). In this article, we analyze the gene expression dynamics (i.e., how the genes expression fluctuates in time) by using a new constructive approach. This approach is based on only two fundamental ingredients: symmetry and the Markov property of dynamics. First, by using experimental data of human and yeast gene expression time series, we found a symmetry in short-time transition probability from time $t$ to time $t+1$. We call it self-similarity symmetry (i.e., surprisingly, the gene expression short-time fluctuations contain a repeating pattern of smaller and smaller parts that are like the whole, but different in size). Secondly, the Markov property of dynamics reflects that the short-time fluctuation governs the full-time behaviour of the system. Here, we succeed in reconstructing naturally the global behavior of the observed distribution of gene expression (i.e., scaling-law) and the local behaviour of the power-law tail of this distribution, by using only these two ingredients: symmetry and the Markov property of dynamics. This approach may represent a step forward toward an integrated image of the basic elements of the whole cell.

Design of amino acid sequences to fold into C_alpha-model proteins

Abstract: In order to extend the results obtained with minimal lattice models to more realistic systems, we study a model where proteins are described as a chain of 20 kinds of structureless amino acids moving in a continuum space and interacting through a contact potential controlled by a 20x20 quenched random matrix. The goal of the present work is to design and characterize amino acid sequences folding to the SH3 conformation, a 60-residues recognition domain common to many regulatory proteins. We show that a number of sequences can fold, starting from a random conformation, to within a distance root mean square deviation (dRMSD) of 2.6A from the native state. Good folders are those sequences displaying in the native conformation an energy lower than a sequence--independent threshold energy.

In silicio stretching of chromatin

Abstract: We present Monte-Carlo (MC) simulations of the stretching of a single 30 nm chromatin fiber The model approximates the DNA by a flexible polymer chain with Debye-H\"uckel electrostatics and uses a two-angle zig-zag model for the geometry of the linker DNA connecting the nucleosomes The latter are represented by flat disks interacting via an attractive Gay-Berne potential Our results show that the stiffness of the chromatin fiber strongly depends on the linker DNA length Furthermore, changing the twisting angle between nucleosomes from 90 deg to 130 deg increases the stiffness significantly An increase in the opening angle from 22 deg to 34 deg leads to softer fibers for small linker lengths We observe that fibers containing a linker histone at each nucleosome are stiffer compared to those without the linker histone The simulated persistence lengths and elastic moduli agree with experimental data Finally, we show that the chromatin fiber does not behave as an isotropic elastic rod, but its rigidity depends on the direction of deformation: chromatin is much more resistant to stretching than to bending

On the Way to Protein Alphabet: Informational Classification of Amino Acids in Comparison to Other Classifications

Abstract: What proteins are made from, as the working parts of the living cells protein machines? To answer this question, we need a technology to disassemble proteins onto elementary functional details and to prepare lumped description of such details. This lumped description might have a multiple material realization (in amino acids). Our hypothesis is that informational approach to this problem is possible. We propose a way of hierarchical classification that makes the primary structure of protein maximally non-random and compare them with other classifications. The first step of the suggested research program is realized: the analysis of protein binary alphabet in comparison with other amino acid classifications.

The energy cost of protein messages lead to a new protein information law

Abstract: By considering the energy cost of messages carried by proteins as proportional to their information content we found experimental proof that proteins from all living organisms tend to have their estimated semantic content of information per unit mass, statistically, close to a constant. Thus, in the message carried by proteins -to achieve minimum energy waste- the rate of information content per unit mass tends to be optimized in living organisms. The experimental evidence of this new information law resembles a marathon where highly optimized proteins correspond to advanced runners followed by a main bunch and the stragglers -lowly optimized proteins. Our results suggest the existence of a continuous optimization process that living organisms had to face, in which a compromise between biological functionality, economic feasibility and the survival requirements is established.

A New Form-Factor Method for the Analysis of Tissue Fluorescence

Abstract: The aim of this article is to present a developed method that decomposes the autofluorescence spectrum into the spectra of naturally occurring biochemical components of biotissue. It requires knowledge of detailed spectrum behaviour of different endogenous fluorophores. We have studied the main bio-markers in human tissue and proposed a simple modelling algorithm for their spectra shapes. The empirical method was tested theoretically by quantum-mechanical calculations of the spectra in the unharmonic Morse potential approach.