The energy landscape in non-biological and biological molecules

Specificity of Monoterpene Interactions with Insect Octopamine and Tyramine Receptors: Insights from in Silico Sequence and Structure Comparison

Octopamine and tyramine receptors (OARs/TARs) are interesting targets for new insecticide development due to their unique roles in insects' physiological and cellular response and their specificity to invertebrates. Monoterpene compounds that bear resemblance to the natural ligands have been shown to bind to the OARs/TARs but elicit varied responses in different insect species. Using in silico methods, we attempt to investigate the molecular basis of monoterpene interactions and their specificity in different OARs and TARs of damaging or beneficial insects. Sequence and structure comparison revealed that the OARs/TARs studied generally have more similarities in terms of structure rather than sequence identity. Together with clustering and network analyses, we also revealed that the role of IL3 might be crucial in the identification of OAR and TAR and their unique function. Among the 35 monoterpenes subjected to ensemble docking, carvacrol had the most negative average binding energies with the target insect OARs and TARs. The differences in the key interacting residues of carvacrol with insect OARs and TARs could be the origin of variation in the responses of insect species to this monoterpene. Results suggest that carvacrol may be a potential natural-product-based insecticide, targeting multiple insect pests while being nonharmful to honeybees and Asian swallowtail butterflies. This work could provide insights into the development of effective species-specific natural-product-based insecticides that are more environmentally friendly than conventional insecticides.

Predictive landscapes hidden beneath biological cellular automata

To celebrate Hans Frauenfelder's achievements, we examine energy(-like) "landscapes" for complex living systems. Energy landscapes summarize all possible dynamics of some physical systems. Energy(-like) landscapes can explain some biomolecular processes, including gene expression and, as Frauenfelder showed, protein folding. But energy-like landscapes and existing frameworks like statistical mechanics seem impractical for describing many living systems. Difficulties stem from living systems being high dimensional, nonlinear, and governed by many, tightly coupled constituents that are noisy. The predominant modeling approach is devising differential equations that are tailored to each living system. This ad hoc approach faces the notorious "parameter problem": models have numerous nonlinear, mathematical functions with unknown parameter values, even for describing just a few intracellular processes. One cannot measure many intracellular parameters or can only measure them as snapshots in time. Another modeling approach uses cellular automata to represent living systems as discrete dynamical systems with binary variables. Quantitative (Hamiltonian-based) rules can dictate cellular automata (e.g., Cellular Potts Model). But numerous biological features, in current practice, are qualitatively described rather than quantitatively (e.g., gene is (highly) expressed or not (highly) expressed). Cellular automata governed by verbal rules are useful representations for living systems and can mitigate the parameter problem. However, they can yield complex dynamics that are difficult to understand because the automata-governing rules are not quantitative and much of the existing mathematical tools and theorems apply to continuous but not discrete dynamical systems. Recent studies found ways to overcome this challenge. These studies either discovered or suggest an existence of predictive "landscapes" whose shapes are described by Lyapunov functions and yield "equations of motion" for a "pseudo-particle." The pseudo-particle represents the entire cellular lattice and moves on the landscape, thereby giving a low-dimensional representation of the cellular automata dynamics. We outline this promising modeling strategy.

Fitness landscapes for coupled map lattices

Our goal is to match some dynamical aspects of biological systems with that of networks of coupled logistic maps. With these networks we generate sequences of iterates, convert them to symbol sequences by coarse-graining, and count the number of times combinations of symbols occur. Comparison of this with the number of times these combinations occur in experimental data-a sequence of interbeat intervals for example-is a measure of the fitness of each network to describe the target data. The most fit networks provide a cartoon that suggests a decomposition of the experimental data into a component that may be produced by a simple dynamical subsystem, and a residual component, the result of detailed, particular characteristics of the system that generated the target data. In the space of all network parameters, each point corresponds to a particular network. We construct a fitness landscape when we assign a fitness to each point. Because the parameters are distributed continuously over their ranges, and because fitnesses are estimated numerically, any plot of the landscape involves a finite sample of parameter values. We'll investigate how the local landscape geometry changes when the array of sample parameters is refined, and use the landscape geometry to explore complex relations between local fitness maxima.

Ligand-Induced Conformational Dynamics of A Tyramine Receptor from Sitophilus oryzae

Tyramine receptor (TyrR) is a biogenic amine G protein-coupled receptor (GPCR) associated with many important physiological functions in insect locomotion, reproduction, and pheromone response. Binding of specific ligands to the TyrR triggers conformational changes, relays the signal to G proteins, and initiates an appropriate cellular response. Here, we monitor the binding effect of agonist compounds, tyramine and amitraz, to a Sitophilus oryzae tyramine receptor (SoTyrR) homology model and their elicited conformational changes. All-atom molecular dynamics (MD) simulations of SoTyrR-ligand complexes have shown varying dynamic behavior, especially at the intracellular loop 3 (IL3) region. Moreover, in contrast to SoTyrR-tyramine, SoTyrR-amitraz and non-liganded SoTyrR shows greater flexibility at IL3 residues and were found to be coupled to the most dominant motion in the receptor. Our results suggest that the conformational changes induced by amitraz are different from the natural ligand tyramine, albeit being both agonists of SoTyrR. This is the first attempt to understand the biophysical implication of amitraz and tyramine binding to the intracellular domains of TyrR. Our data may provide insights into the early effects of ligand binding to the activation process of SoTyrR.

Quasiharmonic analysis of protein energy landscapes from pressure-temperature molecular dynamics simulations

Positional fluctuations of an atom in a protein can be described as motion in an effective local energy minimum created by the surrounding protein atoms. The dependence of atomic fluctuations on both temperature (T) and pressure (P) has been used to probe the nature of these minima, which are generally described as harmonic in experiments such as x-ray crystallography and neutron scattering. Here, a quasiharmonic analysis method is presented in which the P-T dependence of atomic fluctuations is in terms of an intrinsic isobaric thermal expansivity α_P and an intrinsic isothermal compressibility κ_T. The method is tested on previously reported mean-square displacements from P-T molecular dynamics simulations of lysozyme, which were interpreted to have a pressure-independent dynamical transition T_g near 200 K and a change in the pressure dependence near 480 MPa. Our quasiharmonic analysis of the same data shows that the P-T dependence can be described in terms of α_P and κ_T where below T_g, the temperature dependence is frozen at the T_g value. In addition, the purported transition at 480 MPa is reinterpreted as a consequence of the pressure dependence of T_g and the quasiharmonic frequencies. The former also indicates that barrier heights between substates are pressure dependent in these data. Furthermore, the insights gained from this quasiharmonic analysis, which was of the energy landscape near the native state of a protein, suggest that similar analyses of other simulations may be useful in understanding such phenomena as pressure-induced protein unfolding.

Estimating the mean first passage time of protein misfolding

Most theoretical and experimental studies confirm that proteins fold in the time scale of microseconds to milliseconds, but the kinetics of the protein misfolding remains largely unexplored. The kinetics of unfolding-folding-misfolding equilibrium in proteins is formulated in the analytical framework of the Master equation. The folded, unfolded and the misfolded state are characterized in terms of their respective contacts. The Mean First Passage Time (MFPT) to acquire the misfolded conformation from the native or folded state is derived from this equation with different boundary conditions. The MFPT is found to be practically independent of the length of the protein, the number of native contacts and the rate constant for the misfolded to the folded state. The results obtained from the survival probability are directly correlated to the age of onset and appearance of misfolding diseases in humans.

Assessment of the Sampling Performance of Multiple-Copy Dynamics versus a Unique Trajectory

The goal of the present study was to ascertain the differential performance of a long molecular dynamics trajectory versus several shorter ones starting from different points in the phase space and covering the same sampling time. For this purpose, we selected the 16-mer peptide Bak₁₆BH3 as a model for study and carried out several samplings in explicit solvent. These samplings included an 8 μs trajectory (sampling S1); two 4 μs trajectories (sampling S2); four 2 μs trajectories (sampling S3); eight 1 μs trajectories (sampling S4); 16 0.5 μs trajectories (sampling S5), and 80 0.1 μs trajectories (sampling S6). Moreover, the 8 μs trajectory was further extended to 16 μs to have reference values of the diverse properties measured. The diverse samplings were compared qualitatively and quantitatively. Among the former, we carried out a comparison of the conformational profiles of the peptide using cluster analysis. Moreover, we also gained insight into the interchange among these structures along the sampling process. Among the latter, we computed the number of new conformational patterns sampled with time using strings defined from the conformations attained by each of the residues in the peptide. We also compared the locations and depths of the obtained minima on the free energy surface using principal component analysis. Finally, we also compared the helical profiles per residue at the end of the sampling process. The results suggest that a few short molecular dynamics trajectories may provide better sampling than one unique trajectory. Moreover, this procedure can also be advantageous to avoid getting trapped in a local minimum. However, caution should be exercised since short trajectories need to be long enough to overcome local barriers surrounding the starting point and the required sampling time depends on the number of degrees of freedom of the system under study. An effective way to gain insight into the minimum MD trajectory length is to monitor the convergence of different structural features, as shown in the present work.

Conformational Sub-states and Populations in Enzyme Catalysis

Enzyme function involves substrate and cofactor binding, precise positioning of reactants in the active site, chemical turnover, and release of products. In addition to formation of crucial structural interactions between enzyme and substrate(s), coordinated motions within the enzyme-substrate complex allow reaction to proceed at a much faster rate, compared to the reaction in solution and in the absence of enzyme. An increasing number of enzyme systems show the presence of conserved protein motions that are important for function. A wide variety of motions are naturally sampled (over femtosecond to millisecond time-scales) as the enzyme complex moves along the energetic landscape, driven by temperature and dynamical events from the surrounding environment. Areas of low energy along the landscape form conformational sub-states, which show higher conformational populations than surrounding areas. A small number of these protein conformational sub-states contain functionally important structural and dynamical features, which assist the enzyme mechanism along the catalytic cycle. Identification and characterization of these higher-energy (also called excited) sub-states and the associated populations are challenging, as these sub-states are very short-lived and therefore rarely populated. Specialized techniques based on computer simulations, theoretical modeling, and nuclear magnetic resonance have been developed for quantitative characterization of these sub-states and populations. This chapter discusses these techniques and provides examples of their applications to enzyme systems.

Perturbation of the Conformational Dynamics of an Active-Site Loop Alters Enzyme Activity

The role of internal dynamics in enzyme function is highly debated. Specifically, how small changes in structure far away from the reaction site alter protein dynamics and overall enzyme mechanisms is of wide interest in protein engineering. Using RNase A as a model, we demonstrate that elimination of a single methyl group located >10 Å away from the reaction site significantly alters conformational integrity and binding properties of the enzyme. This A109G mutation does not perturb structure or thermodynamic stability, both in the apo and ligand-bound states. However, significant enhancement in conformational dynamics was observed for the bound variant, as probed over nano- to millisecond timescales, resulting in major ligand repositioning. These results illustrate the large effects caused by small changes in structure on long-range conformational dynamics and ligand specificities within proteins, further supporting the importance of preserving wild-type dynamics in enzyme systems that rely on flexibility for function.

Effects of pressure on the dynamics of an oligomeric protein from deep-sea hyperthermophile

Inorganic pyrophosphatase (IPPase) from Thermococcus thioreducens is a large oligomeric protein derived from a hyperthermophilic microorganism that is found near hydrothermal vents deep under the sea, where the pressure is up to 100 MPa (1 kbar). It has attracted great interest in biophysical research because of its high activity under extreme conditions in the seabed. In this study, we use the quasielastic neutron scattering (QENS) technique to investigate the effects of pressure on the conformational flexibility and relaxation dynamics of IPPase over a wide temperature range. The β-relaxation dynamics of proteins was studied in the time ranges from 2 to 25 ps, and from 100 ps to 2 ns, using two spectrometers. Our results indicate that, under a pressure of 100 MPa, close to that of the native environment deep under the sea, IPPase displays much faster relaxation dynamics than a mesophilic model protein, hen egg white lysozyme (HEWL), at all measured temperatures, opposite to what we observed previously under ambient pressure. This contradictory observation provides evidence that the protein energy landscape is distorted by high pressure, which is significantly different for hyperthermophilic (IPPase) and mesophilic (HEWL) proteins. We further derive from our observations a schematic denaturation phase diagram together with energy landscapes for the two very different proteins, which can be used as a general picture to understand the dynamical properties of thermophilic proteins under pressure.

Single particle nanomechanics in operando batteries via lensless strain mapping

We reveal three-dimensional strain evolution in situ of a single LiNi0.5Mn1.5O4 nanoparticle in a coin cell battery under operando conditions during charge/discharge cycles with coherent X-ray diffractive imaging. We report direct observation of both stripe morphologies and coherency strain at the nanoscale. Our results suggest the critical size for stripe formation is 50 nm. Surprisingly, the single nanoparticle elastic energy landscape, which we map with femtojoule precision, depends on charge versus discharge, indicating hysteresis at the single particle level. This approach opens a powerful new avenue for studying battery nanomechanics, phase transformations, and capacity fade under operando conditions at the single particle level that will enable profound insight into the nanoscale mechanisms that govern electrochemical energy storage systems.

The importance of individual protein molecule dynamics in developing and assessing solid state protein preparations

Processing protein solutions into the solid state is a common approach for generating stable amorphous protein mixtures that are suitable for long-term storage. Great care is typically given to protecting the protein native structure during the various drying steps that render it into the amorphous solid state. However, many studies illustrate that chemical and physical degradations still occur in spite of this amorphous material having good glassy properties and it being stored at temperatures below its glass transition temperature (Tg). Because of these persistent issues and recent biophysical studies that have refined the debate ascribing meaning to the molecular dynamical transition temperature and Tg of protein molecules, we provide an updated discussion on the impact of assessing and managing localized, individual protein molecule nondiffusive motions in the context of proteins being prepared into bulk amorphous mixtures. Our aim is to bridge the pharmaceutical studies addressing bulk amorphous preparations and their glassy behavior, with the biophysical studies historically focused on the nondiffusive internal protein dynamics and a protein's activity, along with their combined efforts in assessing the impact of solvent hydrogen-bonding networks on local stability. We also provide recommendations for future research efforts in solid-state formulation approaches.

Two-barrier stability that allows low-power operation in current-induced domain-wall motion

Energy barriers in magnetization reversal dynamics have long been of interest because the barrier height determines the thermal stability of devices as well as the threshold force triggering their dynamics. Especially in memory and logic applications, there is a dilemma between the thermal stability of bit data and the operation power of devices, because larger energy barriers for higher thermal stability inevitably lead to larger magnetic fields (or currents) for operation. Here we show that this is not the case for current-induced magnetic domain-wall motion induced by adiabatic spin-transfer torque. By quantifying domain-wall depinning energy barriers by magnetic field and current, we find that there exist two different pinning barriers, extrinsic and intrinsic energy barriers, which govern the thermal stability and threshold current, respectively. This unique two-barrier system allows low-power operation with high thermal stability, which is impossible in conventional single-barrier systems.

Probabilistic determination of probe locations from distance data

Distance constraints, in principle, can be employed to determine information about the location of probes within a three-dimensional volume. Traditional methods for locating probes from distance constraints involve optimization of scoring functions that measure how well the probe location fits the distance data, exploring only a small subset of the scoring function landscape in the process. These methods are not guaranteed to find the global optimum and provide no means to relate the identified optimum to all other optima in scoring space. Here, we introduce a method for the location of probes from distance information that is based on probability calculus. This method allows exploration of the entire scoring space by directly combining probability functions representing the distance data and information about attachment sites. The approach is guaranteed to identify the global optimum and enables the derivation of confidence intervals for the probe location as well as statistical quantification of ambiguities. We apply the method to determine the location of a fluorescence probe using distances derived by FRET and show that the resulting location matches that independently derived by electron microscopy.

Postnatal mammalian retinal development: quantitative data and general rules

This article is aimed at providing comparative quantitative data about postnatal mammalian retina development, and at searching for some general rules at both the descriptive and the mechanistic level. In mammals the eye continues to grow, and the retina continues to expand, much after the end of retinal cytogenesis. Thus, although the total number of retinal cells remains constant after cessation of mitotic activity (and the end of 'physiological cell death'), the retinal surface area increases by a factor of two or more. In most mammals, ocular growth exceeds retinal expansion: the neural retina lines 70-80% of the inner ocular surface at the beginning but only about 40-60% in adults. Differential local expansion of the retina (the peripheral area increases more than the central one) can be explained by 'passive stretching' of the retinal tissue by the growing eyeball; it depends on the different biomechanical properties of the peripheral vs. central retinal tissue. The increasing retinal surface area allows for a re-distribution of cells such that the thickness of the (particularly, outer) nuclear layer(s) decreases proportional to the areal expansion. This causes a considerable developmental reduction of the number of cell nuclei 'stacked above each other' by a factor of more than two, and requires a translocation of the somata against their neighbors. We provide a physico-mathematical model of these oblique 'down-sliding' movements of the photoreceptor cell somata along the Müller cell process in the center of their columnar cell unit.

End-plate acetylcholine receptor: structure, mechanism, pharmacology, and disease

The synapse is a localized neurohumoral contact between a neuron and an effector cell and may be considered the quantum of fast intercellular communication. Analogously, the postsynaptic neurotransmitter receptor may be considered the quantum of fast chemical to electrical transduction. Our understanding of postsynaptic receptors began to develop about a hundred years ago with the demonstration that electrical stimulation of the vagus nerve released acetylcholine and slowed the heart beat. During the past 50 years, advances in understanding postsynaptic receptors increased at a rapid pace, owing largely to studies of the acetylcholine receptor (AChR) at the motor endplate. The endplate AChR belongs to a large superfamily of neurotransmitter receptors, called Cys-loop receptors, and has served as an exemplar receptor for probing fundamental structures and mechanisms that underlie fast synaptic transmission in the central and peripheral nervous systems. Recent studies provide an increasingly detailed picture of the structure of the AChR and the symphony of molecular motions that underpin its remarkably fast and efficient chemoelectrical transduction.

DNA energy landscapes via calorimetric detection of microstate ensembles of metastable macrostates and triplet repeat diseases

Biopolymers exhibit rough energy landscapes, thereby allowing biological processes to access a broad range of kinetic and thermodynamic states. In contrast to proteins, the energy landscapes of nucleic acids have been the subject of relatively few experimental investigations. In this study, we use calorimetric and spectroscopic observables to detect, resolve, and selectively enrich energetically discrete ensembles of microstates within metastable DNA structures. Our results are consistent with metastable, "native" DNA states being composed of an ensemble of discrete and kinetically stable microstates of differential stabilities, rather than exclusively being a single, discrete thermodynamic species. This conceptual construct is important for understanding the linkage between biopolymer conformational/configurational space and biological function, such as in protein folding, allosteric control of enzyme activity, RNA and DNA folding and function, DNA structure and biological regulation, etc. For the specific DNA sequences and structures studied here, the demonstration of discrete, kinetically stable microstates potentially has biological consequences for understanding the development and onset of DNA expansion and triplet repeat diseases.

Protonation and conformational dynamics of GFP mutants by two-photon excitation fluorescence correlation spectroscopy

GFP mutants are known to display fluorescence flickering, a process that occurs in a wide time range. Because serine 65, threonine 203, glutamate 222, and histidine 148 have been indicated as key residues in determining the GFP fluorescence photodynamics, we have focused here on the role of histidine 148 and glutamate 222 by studying the fluorescence dynamics of GFPmut2 (S65A, V68L, and S72A GFP) and its H148G (Mut2G) and E222Q (Mut2Q) mutants. Two relaxation components are found in the fluorescence autocorrelation functions of GFPmut2: a 10-100 micros pH-dependent component and a 100-500 micros laser-power-dependent component. The comparison of these three mutants shows that the mutation of histidine 148 to glycine induces a 3-fold increase in the protonation rate, thereby indicating that the protonation-deprotonation of the chromophore occurs via a proton exchange with the solution mediated by the histidine 148 residue. The power-dependent but pH-independent relaxation mode, which is not affected by the E222Q and H148G mutations, is due to an excited-state process that is probably related to conformational rearrangements of the chromophore after the photoexcitation, more than to the chromophore excited-state proton transfer.

AID-initiated purposeful mutations in immunoglobulin genes

Exposure brings risk to all living organisms. Using a remarkably effective strategy, higher vertebrates mitigate risk by mounting a complex and sophisticated immune response to counter the potentially toxic invasion by a virtually limitless army of chemical and biological antagonists. Mutations are almost always deleterious, but in the case of antibody diversification there are mutations occurring at hugely elevated rates within the variable (V) and switch regions (SR) of the immunoglobulin (Ig) genes that are responsible for binding to and neutralizing foreign antigens throughout the body. These mutations are truly purposeful. This chapter is centered on activation-induced cytidine deaminase (AID). AID is required for initiating somatic hypermutation (SHM) in the V regions and class switch recombination (CSR) in the SR portions of Ig genes. By converting C --> U, while transcription takes place, AID instigates a cascade of mutational events involving error-prone DNA polymerases, base excision and mismatch repair enzymes, and recombination pathways. Together, these processes culminate in highly mutated antibody genes and the B cells expressing antibodies that have achieved optimal antigenic binding undergo positive selection in germinal centers. We will discuss the biological role of AID in this complex process, primarily in terms of its biochemical properties in relation to SHM in vivo. The chapter also discusses recent advances in experimental methods to characterize antibody dynamics as a function of SHM to help elucidate the role that the AID-induced mutations play in tailoring molecular recognition. The emerging experimental techniques help to address long-standing conundrums concerning evolution-imposed constraints on antibody structure and function.

Protein phase diagrams: the physics behind their elliptic shape

We relate the condition for the elliptic shape of the phase diagram of proteins to the degree of correlation in the fluctuations of the changes of enthalpy and volume at the denaturing-refolding transition. Since this degree cannot be larger than 1, hyperbolically shaped diagrams are not likely to exist. Experiments show that the correlation factor is actually quite low for proteins implying that one-order parameter is not enough to describe the folding-denaturing transition. These findings seem to be the thermodynamic manifestation of the glasslike properties of proteins despite the fact that the transition itself is of first order.

Conformer generation under restraints

Conformational sampling by direct optimization of an all-atom energy function is ineffective and inefficient because of the ruggedness of the energy landscape. Discrete sampling schemes represent an attractive alternative for generating ensembles of conformers consistent with spatial restraints derived from empirical data. Conformational sampling is becoming increasingly important for structure prediction as the bottleneck in accurate prediction shifts from energy functions to the methods used to find low-energy conformers. Experimental structure determination remains a perennial challenge as investigators tackle larger macromolecular systems, and begin to incorporate more complete descriptions of uncertainty, heterogeneity and dynamics into their models. Computational approaches that combine dense, discrete sampling with all-atom energy evaluation and refinement may help to overcome the remaining barriers to solving these problems.

Flexibility and Conformational Entropy in Protein-Protein Binding

To better understand the interplay between protein-protein binding and protein dynamics, we analyzed molecular dynamics simulations of 17 protein-protein complexes and their unbound components. Complex formation does not restrict the conformational freedom of the partner proteins as a whole, but, rather, it leads to a redistribution of dynamics. We calculate the change in conformational entropy for seven complexes with quasiharmonic analysis. We see significant loss, but also increased or unchanged conformational entropy. Where comparison is possible, the results are consistent with experimental data. However, stringent error estimates based on multiple independent simulations reveal large uncertainties that are usually overlooked. We observe substantial gains of pseudo entropy in individual partner proteins, and we observe that all complexes retain residual stabilizing intermolecular motions. Consequently, protein flexibility has an important influence on the thermodynamics of binding and may disfavor as well as favor association. These results support a recently proposed unified model for flexible protein-protein association.

Functional subconformations in protein folding: evidence from single-channel experiments

We study fluctuations in ion conductance and enzymatic rates of the sugar-specific channel-forming membrane protein, Maltoporin, at the single-molecule level. Specifically, we analyze time-persistent deviations in the transport parameters of individual channels from the multichannel averages and discuss our findings in the context of static disorder in protein folding. We show that the disorder responsible for variations in ion conductance does not affect sugar binding, suggesting that Maltoporin can exist in a wide set of fully functional, yet distinctly different, subconformations.

Acid destabilization of the solution conformation of Bcl-xL does not drive its pH-dependent insertion into membranes

Regulation of programmed cell death by Bcl-xL is dependent on both its solution and integral membrane conformations. A conformational change from solution to membrane is also important in this regulation. This conformational change shows a pH-dependence similar to the translocation domain of diphtheria toxin, where an acid-induced molten globule conformation in the absence of lipid vesicles mediates the change from solution to membrane conformations. By contrast, Bcl-xL deltaTM in the absence of lipid vesicles exhibits no gross conformational changes upon acidification as observed by near- and far-UV circular dichroism spectropolarimetry. Additionally, no significant local conformational changes upon acidification were observed by heteronuclear NMR spectroscopy of Bcl-xL deltaTM. Under conditions that favor the solution conformation (pH 7.4), the free energy of folding for Bcl-xL deltaTM (deltaG(o)) was determined to be 15.8 kcal x mol(-1). Surprisingly, under conditions that favor a membrane conformation (pH 4.9), deltaG(o) was 14.6 kcal x mol(-1). These results differ from those obtained with many other membrane-insertable proteins where acid-induced destabilization is important. Therefore, other contributions must be necessary to destabilize the solution conformation Bcl-xL and favor the membrane conformation at pH 4.9. Such contributions might include the presence of a negatively charged membrane or an electrostatic potential across the membrane. Thus, for proteins that adopt both solution and membrane conformations, an obligatory molten globule intermediate may not be necessary. The absence of a molten globule intermediate might have evolved to protect Bcl-xL from intracellular proteases as it undergoes this conformational change essential for its activity.

Thermodynamic and dynamic factors involved in the stability of native protein structure in amorphous solids in relation to levels of hydration

The internal, dynamical fluctuations of protein molecules exhibit many of the features typical of polymeric and bulk small molecule glass forming systems. The response of a protein's internal molecular mobility to temperature changes is similar to that of other amorphous systems, in that different types of motions freeze out at different temperatures, suggesting they exhibit the alpha-beta-modes of motion typical of polymeric glass formers. These modes of motion are attributed to the dynamic regimes that afford proteins the flexibility for function but that also develop into the large-scale collective motions that lead to unfolding. The protein dynamical transition, T(d), which has the same meaning as the T(g) value of other amorphous systems, is attributed to the temperature where protein activity is lost and the unfolding process is inhibited. This review describes how modulation of T(d) by hydration and lyoprotectants can determine the stability of protein molecules that have been processed as bulk, amorphous materials. It also examines the thermodynamic, dynamic, and molecular factors involved in stabilizing folded proteins, and the effects typical pharmaceutical processes can have on native protein structure in going from the solution state to the solid state.

Complementarity of structure ensembles in protein-protein binding

Protein-protein association is often accompanied by changes in receptor and ligand structure. This interplay between protein flexibility and protein-protein recognition is currently the largest obstacle both to our understanding of and to the reliable prediction of protein complexes. We performed two sets of molecular dynamics simulations for the unbound receptor and ligand structures of 17 protein complexes and applied shape-driven rigid body docking to all combinations of representative snapshots. The crossdocking of structure ensembles increased the likelihood of finding near-native solutions. The free ensembles appeared to contain multiple complementary conformations. These were in general not related to the bound structure. We suggest that protein-protein binding follows a three-step mechanism of diffusion, free conformer selection, and refolding. This model combines previously conflicting ideas and is in better agreement with the current data on interaction forces, time scales, and kinetics.

Computational analysis of ligand binding dynamics at the intermolecular hot spots with the aid of simulated tempering and binding free energy calculations

Equilibrium binding dynamics is studied for a panel of benzimidazole-containing compounds at the remodeled interface between human growth hormone (hGH) and the extracellular domain of its receptor (hGHbp), engineered by mutating to glycine hot spot residues T175 from the hormone and W104 from the receptor. Binding energetics is predicted in a good agreement with the experimental data for a panel of these small molecules that complement the engineered defect and restore the binding affinity of the wild-type hGH-hGHbp complex. The results of simulated tempering ligand dynamics at the protein-protein interface reveals a diversity of ligand binding modes that is consistent with the structural orientation of the benzimidazole ring which closely mimics the position of the mutated W104 hot spot residue in the wild-type hGH-hGHbp complex. This structural positioning of the benzimidazole core motif is shown to be a critical feature of the low-energy ligand conformations binding in the engineered cavity. The binding free energy analysis provides a plausible rationale behind the experimental dissociation constants and suggests a key role of ligand-protein van der Waals interactions in restoring binding affinity.

Conformational sampling for the impatient

Several new methods for sampling conformations of biomolecules have appeared recently. A brief review thereof is presented, with particular emphasis on applications that have been published, and suitability for different kinds of systems. Four methods (namely: RESPA, replica-exchange molecular dynamics, CONCOORD and Gaussian network method) are readily applicable for biomolecular systems.

Simultaneous time resolution of the emission spectra of fluorescent proteins and zooxanthellar chlorophyll in reef-building corals

Light is absorbed by photosynthetic algal symbionts (i.e. zooxanthellae) and by chromophoric fluorescent proteins (FP) in reef-building coral tissue. We used a streak-camera spectrograph equipped with a pulsed, blue laser diode (50 ps, 405 nm) to simultaneously resolve the fluorescence spectra and kinetics for both the FP and the zooxanthellae. Shallow water (<9 m)-dwelling Acropora spp. and Plesiastrea versipora specimens were collected from Okinawa, Japan, and Sydney, Australia, respectively. The main FP emitted light in the blue, blue-green and green emission regions with each species exhibiting distinct color morphs and spectra. All corals showed rapidly decaying species and reciprocal rises in greener emission components indicating Förster resonance energy transfer (FRET) between FP populations. The energy transfer modes were around 250 ps, and the main decay modes of the acceptor FP were typically 1900-2800 ps. All zooxanthellae emitted similar spectra and kinetics with peak emission (approximately 683 nm) mainly from photosystem II (PSII) chlorophyll (chl) a. Compared with the FP, the PSII emission exhibited similar rise times but much faster decay times, typically around 640-760 ps. The fluorescence kinetics and excitation versus emission mapping indicated that the FP emission played only a minor role, if any, in chl excitation. We thus suggest the FP could only indirectly act to absorb, screen and scatter light to protect PSII and underlying and surrounding animal tissue from excess visible and UV light. We conclude that our time-resolved spectral analysis and simulation revealed new FP emission components that would not be easily resolved at steady state because of their relatively rapid decays due to efficient FRET. We believe the methods show promise for future studies of coral bleaching and for potentially identifying FP species for use as genetic markers and FRET partners, like the related green FP from Aequorea spp.

Congenital myasthenic syndromes: progress over the past decade

Congenital myasthenic syndromes (CMS) stem from defects in presynaptic, synaptic basal lamina, and postsynaptic proteins. The presynaptic CMS are associated with defects that curtail the evoked release of acetylcholine (ACh) quanta or ACh resynthesis. Defects in ACh resynthesis have now been traced to mutations in choline acetyltransferase. A basal lamina CMS is caused by mutations in the collagenic tail subunit (ColQ) of the endplate species of acetylcholinesterase that prevent the tail subunit from associating with catalytic subunits or from becoming inserted into the synaptic basal lamina. Most postsynaptic CMS are caused by mutations in subunits of the acetylcholine receptor (AChR) that alter the kinetic properties or decrease the expression of AChR. The kinetic mutations increase or decrease the synaptic response to ACh and result in slow- and fast-channel syndromes, respectively. Most low-expressor mutations reside in the AChR epsilon subunit and are partially compensated by residual expression of the fetal type gamma subunit. In a subset of CMS patients, endplate AChR deficiency is caused by mutations in rapsyn, a molecule that plays a critical role in concentrating AChR in the postsynaptic membrane.

Monte Carlo simulations of the peptide recognition at the consensus binding site of the constant fragment of human immunoglobulin G: the energy landscape analysis of a hot spot at the intermolecular interface

Monte Carlo simulations of molecular recognition at the consensus binding site of the constant fragment (Fc) of human immunoglobulin G (Ig) protein have been performed to analyze structural and thermodynamic aspects of binding for the 13-residue cyclic peptide DCAWHLGELVWCT. The energy landscape analysis of a hot spot at the intermolecular interface using alanine scanning and equilibrium-simulated tempering dynamics with the simplified, knowledge-based energy function has enabled the role of the protein hot spot residues in providing the thermodynamic stability of the native structure to be determined. We have found that hydrophobic interactions between the peptide and the Met-252, Ile-253, His-433, and His-435 protein residues are critical to guarantee the thermodynamic stability of the crystallographic binding mode of the complex. Binding free energy calculations, using a molecular mechanics force field and a solvation energy model, combined with alanine scanning have been conducted to determine the energetic contribution of the protein hot spot residues in binding affinity. The conserved Asn-434, Ser-254, and Tyr-436 protein residues contribute significantly to the binding affinity of the peptide-protein complex, serving as an energetic hot spot at the intermolecular interface. The results suggest that evolutionary conserved hot spot protein residues at the intermolecular interface may be partitioned in fulfilling thermodynamic stability of the native binding mode and contributing to the binding affinity of the complex.

Solution structure and dynamics of the lipoic acid-bearing domain of human mitochondrial branched-chain alpha-keto acid dehydrogenase complex

The lipoyl-bearing domain (LBD) of the transacylase (E2) subunit of the branched-chain alpha-keto acid dehydrogenase complex plays a central role in substrate channeling in this mitochondrial multienzyme complex. We have employed multidimensional heteronuclear NMR techniques to determine the structure and dynamics of the LBD of the human branched-chain alpha-keto acid dehydrogenase complex (hbLBD). Similar to LBD from other members of the alpha-keto acid dehydrogenase family, the solution structure of hbLBD is a flattened beta-barrel formed by two four-stranded antiparallel beta-sheets. The lipoyl Lys(44) residue resides at the tip of a beta-hairpin comprising a sharp type I beta-turn and the two connecting beta-strands 4 and 5. A prominent V-shaped groove formed by a surface loop, L1, connecting beta 1- and beta 2-strands and the lipoyl lysine beta-hairpin constitutes the functional pocket. We further applied reduced spectral density functions formalism to extract dynamic information of hbLBD from (15)N-T(1), (15)N-T(2), and ((1)H-(15)N) nuclear Overhauser effect data obtained at 600 MHz. The results showed that residues surrounding the lipoyl lysine region comprising the L1 loop and the Lys(44) beta-turn are highly flexible, whereas beta-sheet S1 appears to display a slow conformational exchange process.

Complexity and simplicity of ligand-macromolecule interactions: the energy landscape perspective

The energy landscape approach has contributed to recent progress in understanding the complexity and simplicity of ligand-macromolecule interactions. Significant advances in computational structure prediction of ligand-protein complexes have been made using approaches that include the effects of protein flexibility and incorporate a hierarchy of energy functions. The results suggest that the complexity of structure prediction in molecular recognition may be determined by low-resolution properties of the underlying binding energy landscapes and by the nature of the energy funnels near the native structures of the complexes.

Kinetic transition in model proteins with a denatured native spinodal

The relaxation kinetics of three-dimensional lattice model proteins with Gō potential is studied. A kinetic transition from an exponential behavior to a nonexponential one with a denatured native spinodal is characterized. The transition temperatures T(k), obtained from simulations and a semiquantitative estimation, are found to be the same. The change in free energy landscape during the transition is discussed microscopically by studying the detailed folding processes of various paths. The connection of T(k) with the foldability is also studied by a Z-score-like quantity T(f)/T(k).

The red-edge effects: 30 years of exploration

In 1970, three laboratories independently made a discovery that, for aromatic fluorophores embedded into different rigid and highly viscous media, the spectroscopic properties do not conform to classical rules. The fluorescence spectra can depend on excitation wavelength, and the excited-state energy transfer, if present, fails at the "red" excitation edge. These red-edge effects were related to the existence of excited-state distribution of fluorophores on their interaction energy with the environment and the slow rate of dielectric relaxation of this environment. In these conditions the site-selection can be provided by variation of the energy of illuminating light quanta, and the behaviour of selected species can be followed as a function of time and other variables. These observations found extensive application in different areas of research: colloid and polymer science, molecular biophysics, photochemistry and photobiology. In particular, they led to the development of very productive methods of studying the dynamics of dielectric relaxations in protein and membranes, using the tryptophan emission and the emission of a variety of probes. These studies were extended to the time domain with the observation of new site-selective effects in emission intensity and anisotropy decays. They stimulated the emergence and development of cryogenic energy-selective and single-molecular techniques that became valuable tools in their own right in chemistry and biophysics research. Site-selection effects were discovered for electron-transfer and proton-transfer reactions if they depended on the dynamics of the environment. This review is focused on the progress in the field of red-edge effects, their applications and prospects.

Backbone dynamics of Escherichia coli thioesterase/protease I: evidence of a flexible active-site environment for a serine protease

Escherichia coli thioesterase/protease I (TEP-I) is a member of a novel subclass of the lipolytic enzymes with a distinctive GDSLS motif. In addition to possessing thioesterase and protease activities, TEP-I also exhibits arylesterase activity. We have determined the (15)N nuclear magnetic spin relaxation rates, R(1) and R(2), and the steady state (1)H-(15)N heteronuclear Overhauser effect, measured at both 11.74 T and 14.09 T, of (u-(15)N) TEP-I. These data were analyzed using model-free formalism (with axially symmetric rotational diffusion anisotropy) to extract the backbone dynamics of TEP-I. The results reveal that the core structure of the central beta-sheet and the long alpha-helices are rigid, while the binding pocket appears to be rather flexible. The rigid core serves as a scaffold to anchor the essential loops, which form the binding pocket. The most flexible residues display large amplitude fast (ps/ns time-scale) motion and lie on one stripe whose orientation is presumed to be the ligand-binding orientation. We also detected the presence of several residues displaying slow (microseconds/ms time-scale) conformational exchanging processes. These residues lie around the binding pocket and are oriented perpendicularly to the orientation of the flexible stripe. Two of the putative catalytic triads, Ser10 and His157, and their neighbors show motion on the microseconds/ms time-scale, suggesting that their slow motion may have a role in catalysis, in addition to their possible roles in ligand binding. The presence of a flexible substrate-binding pocket may also facilitate binding to a wide range of substrates and confer the versatile functional property of this protein.

Trehalose effect on low temperature protein dynamics: fluctuation and relaxation phenomena

We performed spectral diffusion experiments in trehalose-enriched glycerol/buffer-glass on horseradish peroxidase where the heme was replaced by metal-free mesoporphyrin IX, and compared them with the respective behavior in a pure glycerol/buffer-glass (Schlichter et al., J. Chem. Phys. 2000, 112:3045-3050). Trehalose has a significant influence: spectral diffusion broadening speeds up compared to the trehalose-free glass. This speeding up is attributed to a shortening of the correlation time of the frequency fluctuations most probably by preventing water molecules from leaving the protein interior. Superimposed to the frequency fluctuation dynamics is a relaxation dynamics that manifests itself as an aging process in the spectral diffusion broadening. Although the trehalose environment speeds up the fluctuations, it does not have any influence on the relaxation. Both relaxation and fluctuations are governed by power laws in time. The respective exponents do not seem to change with the protein environment. From the spectral dynamics, the mean square displacement in conformation space can be determined. It is governed by anomalous diffusion. The associated frequency correlation time is incredibly long, demonstrating that proteins at low temperatures are truly nonergodic systems.

Enhanced internal dynamics of a membrane transport protein during substrate translocation

Conformational changes are essential for the activity of many proteins. If, or how fast, internal fluctuations are related to slow conformational changes that mediate protein function is not understood. In this study, we measure internal fluctuations of the transport protein lactose permease in the presence and absence of substrate by tryptophan fluorescence spectroscopy. We demonstrate that nanosecond fluctuations of alpha-helices are enhanced when the enzyme transports substrate. This correlates with previously published kinetic data from transport measurements showing that millisecond conformational transitions of the substrate-loaded carrier are faster than those in the absence of substrate. These findings corroborate the hypothesis of the hierarchical model of protein dynamics that predicts that slow conformational transitions are based on fast, thermally activated internal motions.

Fundamental gating mechanism of nicotinic receptor channel revealed by mutation causing a congenital myasthenic syndrome

We describe the genetic and kinetic defects in a congenital myasthenic syndrome due to the mutation epsilonA411P in the amphipathic helix of the acetylcholine receptor (AChR) epsilon subunit. Myasthenic patients from three unrelated families are either homozygous for epsilonA411P or are heterozygous and harbor a null mutation in the second epsilon allele, indicating that epsilonA411P is recessive. We expressed human AChRs containing wild-type or A411P epsilon subunits in 293HEK cells, recorded single channel currents at high bandwidth, and determined microscopic rate constants for individual channels using hidden Markov modeling. For individual wild-type and mutant channels, each rate constant distributes as a Gaussian function, but the spread in the distributions for channel opening and closing rate constants is greatly expanded by epsilonA411P. Prolines engineered into positions flanking residue 411 of the epsilon subunit greatly increase the range of activation kinetics similar to epsilonA411P, whereas prolines engineered into positions equivalent to epsilonA411 in beta and delta subunits are without effect. Thus, the amphipathic helix of the epsilon subunit stabilizes the channel, minimizing the number and range of kinetic modes accessible to individual AChRs. The findings suggest that analogous stabilizing structures are present in other ion channels, and possibly allosteric proteins in general, and that they evolved to maintain uniformity of activation episodes. The findings further suggest that the fundamental gating mechanism of the AChR channel can be explained by a corrugated energy landscape superimposed on a steeply sloped energy well.