Hierarchical sampling for metastable conformers determines biomolecular recognition: the case of malectin and diglucosylated N-glycan interactions

Structural information over the entire course of binding interactions based on the analyses of energy landscapes is described, which provides a framework to understand the events involved during biomolecular recognition. Conformational dynamics of malectin's exquisite selectivity for diglucosylated N-glycan (Dig-N-glycan), a highly flexible oligosaccharide comprising of numerous dihedral torsion angles, are described as an example. For this purpose, a novel approach based on hierarchical sampling for acquiring metastable molecular conformations constituting low-energy minima for understanding the structural features involved in a biologic recognition is proposed. For this purpose, four variants of principal component analysis were employed recursively in both Cartesian space and dihedral angles space that are characterized by free energy landscapes to select the most stable conformational substates. Subsequently, k-means clustering algorithm was implemented for geometric separation of the major native state to acquire a final ensemble of metastable conformers. A comparison of malectin complexes was then performed to characterize their conformational properties. Analyses of stereochemical metrics and other concerted binding events revealed surface complementarity, cooperative and bidentate hydrogen bonds, water-mediated hydrogen bonds, carbohydrate-aromatic interactions including CH-pi and stacking interactions involved in this recognition. Additionally, a striking structural transition from loop to beta-strands in malectin CRD upon specific binding to Dig-N-glycan is observed. The interplay of the above-mentioned binding events in malectin and Dig-N-glycan supports an extended conformational selection model as the underlying binding mechanism.

Surface attachment of protein fibrils via covalent modification strategies.

Chemical control of surface functionality and topography is an essential requirement for many technological purposes. In particular, the covalent attachment of monomeric proteins to surfaces has been the object of intense studies in recent years, for applications as varied as electrochemistry, immuno-sensing, and the production of biocompatible coatings. Little is known, however, about the characteristics and requirements underlying surface attachment of supramolecular protein nanostructures. Amyloid fibrils formed by the self-assembly of peptide and protein molecules represent one important class of such structures. These highly organized beta-sheet-rich assemblies are a hallmark of a range of neurodegenerative disorders, including Alzheimer's disease and type II diabetes, but recent findings suggest that they have much broader significance, potentially representing the global free energy minima of the energy landscapes of proteins and having potential applications in material science. In this paper, we describe strategies for attaching amyloid fibrils formed from different proteins to gold surfaces under different solution conditions. Our methods involve the reaction of sulfur containing small molecules (cystamine and 2-iminothiolane) with the amyloid fibrils, enabling their covalent linkage to gold surfaces. We demonstrate that irreversible attachment using these approaches makes possible quantitative analysis of experiments using biosensor techniques, such as quartz crystal microbalance (QCM) assays that are revolutionizing our understanding of the mechanisms of amyloid growth and the factors that determine its kinetic behavior. Moreover, our results shed light on the nature and relative importance of covalent versus noncovalent forces acting on protein superstructures at metal surfaces.

III-V semiconductor nanowires for optoelectronic device applications

Semiconductor nanowires have recently emerged as a new class of materials with significant potential to reveal new fundamental physics and to propel new applications in quantum electronic and optoelectronic devices. Semiconductor nanowires show exceptional promise as nanostructured materials for exploring physics in reduced dimensions and in complex geometries, as well as in one-dimensional nanowire devices. They are compatible with existing semiconductor technologies and can be tailored into unique axial and radial heterostructures. In this contribution we review the recent efforts of our international collaboration which have resulted in significant advances in the growth of exceptionally high quality IIIV nanowires and nanowire heterostructures, and major developments in understanding the electronic energy landscapes of these nanowires and the dynamics of carriers in these nanowires using photoluminescence, time-resolved photoluminescence and terahertz conductivity spectroscopy. © 2011 Elsevier Ltd. All rights reserved.

Kinetics and Statistical Distributions of Single-Molecule Conformational Dynamics

We developed a coarse-grained yet microscopic detailed model to study the statistical fluctuations of single-molecule protein conformational dynamics of adenylate kinase. We explored the underlying conformational energy landscape and found that the system has two basins of attractions, open and closed conformations connected by two separate pathways. The kinetics is found to be nonexponential, consistent with single-molecule conformational dynamics experiments. Furthermore, we found that the statistical distribution of the kinetic times for the conformational transition has a long power law tail, reflecting the exponential density of state of the underlying landscape. We also studied the joint distribution of the two pathways and found memory effects.

Quantifying intrinsic specificity: A potential complement to affinity in drug screening

We report here the investigation of a novel description of specificity in protein-ligand binding based on energy landscape theory. We define a new term, intrinsic specificity ratio (ISR), which describes the level of discrimination in binding free energies of the native basin for a protein-ligand complex from the weaker binding states of the same ligand. We discuss the relationship between the intrinsic specificity we defined here and the conventional definition of specificity. In a docking study of molecules with the enzyme COX-2, we demonstrate a statistical correspondence between ISR value and geometrical shapes of the small molecules binding to COX-2. We further observe that the known selective (nonselective) inhibitors of COX-2 have higher (lower) ISR values. We suggest that intrinsic specificity ratio may be a useful new criterion and a complement to affinity in drug screening and in searching for potential drug lead compounds.

Using sketch-map coordinates to analyze and bias molecular dynamics simulations

When examining complex problems, such as the folding of proteins, coarse grained descriptions of the system drive our investigation and help us to rationalize the results. Oftentimes collective variables (CVs), derived through some chemical intuition about the process of interest, serve this purpose. Because finding these CVs is the most difficult part of any investigation, we recently developed a dimensionality reduction algorithm, sketch-map, that can be used to build a low-dimensional map of a phase space of high-dimensionality. In this paper we discuss how these machine-generated CVs can be used to accelerate the exploration of phase space and to reconstruct free-energy landscapes. To do so, we develop a formalism in which high-dimensional configurations are no longer represented by low-dimensional position vectors. Instead, for each configuration we calculate a probability distribution, which has a domain that encompasses the entirety of the low-dimensional space. To construct a biasing potential, we exploit an analogy with metadynamics and use the trajectory to adaptively construct a repulsive, history-dependent bias from the distributions that correspond to the previously visited configurations. This potential forces the system to explore more of phase space by making it desirable to adopt configurations whose distributions do not overlap with the bias. We apply this algorithm to a small model protein and succeed in reproducing the free-energy surface that we obtain from a parallel tempering calculation.

Probing the Unfolded Configurations of a β-Hairpin Using Sketch-Map

This work examines the conformational ensemble involved in β-hairpin folding by means of advanced molecular dynamics simulations and dimensionality reduction. A fully atomistic description of the protein and the surrounding solvent molecules is used, and this complex energy landscape is sampled by means of parallel tempering metadynamics simulations. The ensemble of configurations explored is analyzed using the recently proposed sketch-map algorithm. Further simulations allow us to probe how mutations affect the structures adopted by this protein. We find that many of the configurations adopted by a mutant are the same as those adopted by the wild-type protein. Furthermore, certain mutations destabilize secondary-structure-containing configurations by preventing the formation of hydrogen bonds or by promoting the formation of new intramolecular contacts. Our analysis demonstrates that machine-learning techniques can be used to study the energy landscapes of complex molecules and that the visualizations that are generated in this way provide a natural basis for examining how the stabilities of particular configurations of the molecule are affected by factors such as temperature or structural mutations.

Visualization of Protein Folding Funnels in Lattice Models

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Protein folding kinetics exhibit an Arrhenius temperature dependence when corrected for the temperature dependence of protein stability

The anomalous temperature dependence of protein folding has received considerable attention. Here we show that the temperature dependence of the folding of protein L becomes extremely simple when the effects of temperature on protein stability are corrected for; the logarithm of the folding rate is a linear function of 1/T on constant stability contours in the temperature–denaturant plane. This convincingly demonstrates that the anomalous temperature dependence of folding derives from the temperature dependence of the interactions that stabilize proteins, rather than from the super Arrhenius temperature dependence predicted for the configurational diffusion constant on a rough energy landscape. However, because of the limited temperature range accessible to experiment, the results do not rule out models with higher order temperature dependences. The significance of the slope of the stability-corrected Arrhenius plots is discussed.

Evidence for nonrandom hydrophobicity structures in protein chains.

The question of whether proteins originate from random sequences of amino acids is addressed. A statistical analysis is performed in terms of blocked and random walk values formed by binary hydrophobic assignments of the amino acids along the protein chains. Theoretical expectations of these variables from random distributions of hydrophobicities are compared with those obtained from functional proteins. The results, which are based upon proteins in the SWISS-PROT data base, convincingly show that the amino acid sequences in proteins differ from what is expected from random sequences in a statistically significant way. By performing Fourier transforms on the random walks, one obtains additional evidence for nonrandomness of the distributions. We have also analyzed results from a synthetic model containing only two amino acid types, hydrophobic and hydrophilic. With reasonable criteria on good folding properties in terms of thermodynamical and kinetic behavior, sequences that fold well are isolated. Performing the same statistical analysis on the sequences that fold well indicates similar deviations from randomness as for the functional proteins. The deviations from randomness can be interpreted as originating from anticorrelations in terms of an Ising spin model for the hydrophobicities. Our results, which differ from some previous investigations using other methods, might have impact on how permissive with respect to sequence specificity protein folding process is-only sequences with nonrandom hydrophobicity distributions fold well. Other distributions give rise to energy landscapes with poor folding properties and hence did not survive the evolution.

Chaperonin-facilitated protein folding: optimization of rate and yield by an iterative annealing mechanism.

We develop a heuristic model for chaperonin-facilitated protein folding, the iterative annealing mechanism, based on theoretical descriptions of "rugged" conformational free energy landscapes for protein folding, and on experimental evidence that (i) folding proceeds by a nucleation mechanism whereby correct and incorrect nucleation lead to fast and slow folding kinetics, respectively, and (ii) chaperonins optimize the rate and yield of protein folding by an active ATP-dependent process. The chaperonins GroEL and GroES catalyze the folding of ribulose bisphosphate carboxylase at a rate proportional to the GroEL concentration. Kinetically trapped folding-incompetent conformers of ribulose bisphosphate carboxylase are converted to the native state in a reaction involving multiple rounds of quantized ATP hydrolysis by GroEL. We propose that chaperonins optimize protein folding by an iterative annealing mechanism; they repeatedly bind kinetically trapped conformers, randomly disrupt their structure, and release them in less folded states, allowing substrate proteins multiple opportunities to find pathways leading to the most thermodynamically stable state. By this mechanism, chaperonins greatly expand the range of environmental conditions in which folding to the native state is possible. We suggest that the development of this device for optimizing protein folding was an early and significant evolutionary event.

The distribution of macrobenthos in surface sediments at Fladen Ground stations, North Sea

A general study of structure, biomass estimates and dynamics on the macrofauna was carried out in August 1975 and March 1976 during PREFLEX (1975) and FLEX (1976), the Fladen Ground Experiment. On the basis of these data an attempt was made to estimate macrobenthic production expressed as minimum production (MP). The macrobenthic production is discussed together with meiobenthic annual production and with indirectly estimated microbenthic production in relation to an energy input from the water column of about 25 g C m**-2 year**-1. From the production estimates of the three benthic components a rough energy budget is proposed. Sampling was performed at five stations for endofauna twice during the time of investigation and for epifauna once. At each station two replicate box core samples (30 X 20 cm) were taken for endofauna. Epifauna was sampled with an Agassiz trawl once at each station. The total numbers of endofauna increased from station 1 to 5. This was valid as well for August 1975 (4,233-12,166 individuals per m**2 and 10 cm sediment depth) as for March 1976 (1,008-2,925 individuals). The polychaetes were the dominant organisms with a share of 33 to 62 %. The densities for the endofauna decreased from August 1975 to March 1976 by a mean factor of 2.8. Abundances of epifauna amounted to values between 11 and 102 individuals per 1000 m**2. The biomass dry weights (DWT) for macrobenthic endofauna varied between 0.97 g DWT m**-2 and 6.42 g DWT m**-2 in August 1975 and between 0.27 g DWT m**-2 and 2.64 g DWT m**-2 in March 1976. The mean amounted to 1.74 g DWT m**-2. Dry weights of epifauna biomass gave values between 4.9 and 83.1 g DWT * 1000 m**-2. The minimum production for the total macro-endofauna at Fladen Ground amounted to 1.43 g DWT m**-2 yr**-1 or 0.82 g C m**-2 yr**-1. This resulted in a minimum turnover rate (P/B) of 0.8. The share produced by the polychaetes amounted to 1.06g DWT m**-2 yr**-1 or 74 %.

Quantifying protein folding pathways

A deep understanding of the proteins folding dynamics can be get quantifying folding landscape by calculating how the number of microscopic configurations (entropy) varies with the energy of the chain, Ω=Ω(E). Because of the incredibly large number of microstates available to a protein, direct enumeration of Ω(E) is not possible on realistic computer simulations. An estimate of Ω(E) can be obtained by use of a combination of statistical mechanics and thermodynamics. By combining different definitions of entropy that are valid for a system whose probability for occupying a state is given by the canonical Boltzmann probability, computers allow the determination of Ω(E). ^ The energy landscapes of two similar, but not identical model proteins were studied. One protein contains no kinetic tracks. Results show a smooth funnel for the folding landscape. That allows the contour determination of the folding funnel. Also it was presented results for the folding landscape for a modified protein with kinetic traps. Final results show that the computational approach is able to distinguish and explore regions of the folding landscape that are due to kinetic traps from the native state folding funnel.^

Practical engineering of hard spin-glass instances

Recent technological developments in the field of experimental quantum annealing have made prototypical annealing optimizers with hundreds of qubits commercially available. The experimental demonstration of a quantum speedup for optimization problems has since then become a coveted, albeit elusive goal. Recent studies have shown that the so far inconclusive results, regarding a quantum enhancement, may have been partly due to the benchmark problems used being unsuitable. In particular, these problems had inherently too simple a structure, allowing for both traditional resources and quantum annealers to solve them with no special efforts. The need therefore has arisen for the generation of harder benchmarks which would hopefully possess the discriminative power to separate classical scaling of performance with size from quantum. We introduce here a practical technique for the engineering of extremely hard spin-glass Ising-type problem instances that does not require "cherry picking" from large ensembles of randomly generated instances. We accomplish this by treating the generation of hard optimization problems itself as an optimization problem, for which we offer a heuristic algorithm that solves it. We demonstrate the genuine thermal hardness of our generated instances by examining them thermodynamically and analyzing their energy landscapes, as well as by testing the performance of various state-of-the-art algorithms on them. We argue that a proper characterization of the generated instances offers a practical, efficient way to properly benchmark experimental quantum annealers, as well as any other optimization algorithm.

Folding funnels and frustration in off-lattice minimalist protein landscapes

A full quantitative understanding of the protein folding problem is now becoming possible with the help of the energy landscape theory and the protein folding funnel concept. Good folding sequences have a landscape that resembles a rough funnel where the energy bias towards the native state is larger than its ruggedness. Such a landscape leads not only to fast folding and stable native conformations but, more importantly, to sequences that are robust to variations in the protein environment and to sequence mutations. In this paper, an off-lattice model of sequences that fold into a β-barrel native structure is used to describe a framework that can quantitatively distinguish good and bad folders. The two sequences analyzed have the same native structure, but one of them is minimally frustrated whereas the other one exhibits a high degree of frustration.