Folding simulations of a three-stranded antiparallel β-sheet peptide

Protein folding is a grand challenge of the postgenomic era. In this paper, 58 folding events sampled during 47 molecular dynamics trajectories for a total simulation time of more than 4 μs provide an atomic detail picture of the folding of a 20-residue synthetic peptide with a stable three-stranded antiparallel β-sheet fold. The simulations successfully reproduce the NMR solution conformation, irrespective of the starting structure. The sampling of the conformational space is sufficient to determine the free energy surface and localize the minima and transition states. The statistically predominant folding pathway involves the formation of contacts between strands 2 and 3, starting with the side chains close to the turn, followed by association of the N-terminal strand onto the preformed 2–3 β-hairpin. The folding mechanism presented here, formation of a β-hairpin followed by consolidation, is in agreement with a computational study of the free energy surface of another synthetic three-stranded antiparallel β-sheet by Bursulaya and Brooks [(1999) J. Am. Chem. Soc. 121, 9947–9951]. Hence, it might hold in general for antiparallel β-sheets with short turns.

Top-down free-energy minimization on protein potential energy landscapes

The hierarchical properties of potential energy landscapes have been used to gain insight into thermodynamic and kinetic properties of protein ensembles. It also may be possible to use them to direct computational searches for thermodynamically stable macroscopic states, i.e., computational protein folding. To this end, we have developed a top-down search procedure in which conformation space is recursively dissected according to the intrinsic hierarchical structure of a landscape's effective-energy barriers. This procedure generates an inverted tree similar to the disconnectivity graphs generated by local minima-clustering methods, but it fundamentally differs in the manner in which the portion of the tree that is to be computationally explored is selected. A key ingredient is a branch-selection algorithm that takes advantage of statistically predictive properties of the landscape to guide searches down the tree branches that are most likely to lead to the physically relevant macroscopic states. Using the computational folding of a β-hairpin-forming peptide as an example, we show that such predictive properties indeed exist and can be used for structure prediction by free-energy global minimization.

Lifetimes of intermediates in the β-sheet to α-helix transition of β-lactoglobulin by using a diffusional IR mixer

The extremely slow α-helix/β-sheet transition of proteins is a crucial step in amylogenic diseases and represents an internal rearrangement of local contacts in an already folded protein. These internal structural rearrangements within an already folded protein are a critical aspect of biological action and are a product of conformational flow along unknown metastable local minima of the energy landscape of the compact protein. We use a diffusional IR mixer with time-resolved Fourier transform IR spectroscopy capable of 400-μs time resolution to show that the trifluoroethanol driven β-sheet to α-helix transition of β-lactoglobulin proceeds via a compact β-sheet intermediate with a lifetime of 7 ms, small compared with the overall folding time of β-lactoglobulin.

Circadian-Regulated Transcription of the psbD Light-Responsive Promoter in Wheat Chloroplasts1

The level of mRNAs derived from the plastid-encoded psbD light-responsive promoter (LRP) is controlled by a circadian clock(s) in wheat (Triticum aestivum). The circadian oscillations in the psbD LRP mRNA level persisted for at least three cycles in continuous light and for one cycle in continuous dark, with maxima in subjective morning and minima in subjective early night. In vitro transcription in chloroplast extracts revealed that the circadian cycles in the psbD LRP mRNA level were dominantly attributed to the circadian-regulated transcription of the psbD LRP. The effects of various mutations introduced into the promoter region on the psbD LRP activity in vitro suggest the existence of two positive elements located between −54 and −36, which generally enhance the transcription activity, and an anomalous core promoter structure lacking the functional “−35” element, which plays a crucial role in the circadian fluctuation and light dependency of psbD LRP transcription activity.

Accumulation of a Clock-Regulated Transcript during Flower-Inductive Darkness in Pharbitis nil1

To clarify the molecular basis of the photoperiodic induction of flowering in the short-day plant Pharbitis nil cv Violet, we examined changes in the level of mRNA in cotyledons during the flower-inductive photoperiod using the technique of differential display by the polymerase chain reaction. A transcript that accumulated during the inductive dark period was identified and a cDNA corresponding to the transcript, designated PnC401 (P. nil C401), was isolated. RNA-blot hybridization verified that levels of PnC401 mRNA fluctuated with a circadian rhythm, with maxima between 12 and 16 h after the beginning of the dark period) and minima of approximately 0. This oscillation continued even during an extended dark period but was damped under continuous light. Accumulation of PnC401 mRNA was reduced by a brief exposure to red light at the 8th h of the dark period (night-break treatment) or by exposure to far-red light at the end of the light period (end-of-day far-red treatment). These results suggest that fluctuations in levels of PnC401 mRNA are regulated by phytochrome(s) and a circadian clock and that they are associated with photoperiodic events that include induction of flowering.

Temperature dependence of the isotope chemistry of the heavy elements.

The temperature coefficient of equilibrium isotope fractionation in the heavy elements is shown to be larger at high temperatures than that expected from the well-studied vibrational isotope effects. The difference in the isotopic behavior of the heavy elements as compared with the light elements is due to the large nuclear isotope field shifts in the heavy elements. The field shifts introduce new mechanisms for maxima, minima, crossovers, and large mass-independent isotope effects in the isotope chemistry of the heavy elements. The generalizations are illustrated by the temperature dependence of the isotopic fractionation in the redox reaction between U(VI) and U(IV) ions.

Betidamino acids: versatile and constrained scaffolds for drug discovery.

Betidamino acids (a contraction of "beta" position and "amide") are N'-monoacylated (optionally, N'-monoacylated and N-mono- or N,N'-dialkylated) aminoglycine derivatives in which each N'acyl/alkyl group may mimic naturally occurring amino acid side chains or introduce novel functionalities. Betidamino acids are most conveniently generated on solid supports used for the synthesis of peptides by selective acylation of one of the two amino functions of orthogonally protected aminoglycine(s) to generate the side chain either prior to or after the elongation of the main chain. We have used unresolved Nalpha-tert-butyloxycarbonyl-N'alpha-fluorenylmethoxycarbonyl++ + aminoglycine, and Nalpha-(Nalpha-methyl)-tert-butyloxycarbonyl-N'alpha-fluo renylmethoxycarbonyl aminoglycine as the templates for the introduction of betidamino acids in Acyline [Ac-D2Nal-D4Cpa-D3Pal-Ser-4Aph(Ac)-D4Aph(A c)-Leu-Ilys-Pro-DAla-NH2, where 2Nal is 2-naphthylalanine, 4Cpa is 4-chlorophenylalanine, 3Pal is 3-pyridylalanine, Aph is 4-aminophenylalanine, and Ilys is Nepsilon-isopropyllysine], a potent gonadotropin-releasing hormone antagonist, in order to test biocompatibility of these derivatives. Diasteremneric peptides could be separated in most cases by reverse-phase HPLC. Biological results indicated small differences in relative potencies (<5-fold) between the D and L nonalkylated betidamino acid-containing Acyline derivatives. Importantly, most betide diastereomers were equipotent with Acyline. In an attempt to correlate structure and observed potency, Ramachandran-type plots were calculated for a series of betidamino acids and their methylated homologs. According to these calculations, betidamino acids have access to a more limited and distinct number of conformational states (including those associated with alpha-helices, beta-sheets, or turn structures), with deeper minima than those observed for natural amino acids.

Unusual dynamics of extinction in a simple ecological model.

Studies on natural populations and harvesting biological resources have led to the view, commonly held, that (i) populations exhibiting chaotic oscillations run a high risk of extinction; and (ii) a decrease in emigration/exploitation may reduce the risk of extinction. Here we describe a simple ecological model with emigration/depletion that shows behavior in contrast to this. This model displays unusual dynamics of extinction and survival, where populations growing beyond a critical rate can persist within a band of high depletion rates, whereas extinction occurs for lower depletion rates. Though prior to extinction at lower depletion rates the population exhibits chaotic dynamics with large amplitudes of variation and very low minima, at higher depletion rates the population persists at chaos but with reduced variation and increased minima. For still higher values, within the band of persistence, the dynamics show period reversal leading to stability. These results illustrate that chaos does not necessarily lead to population extinction. In addition, the persistence of populations at high depletion rates has important implications in the considerations of strategies for the management of biological resources.

The old problems of glass and the glass transition, and the many new twists.

In this paper I review the ways in which the glassy state is obtained both in nature and in materials science and highlight a "new twist"--the recent recognition of polymorphism within the glassy state. The formation of glass by continuous cooling (viscous slowdown) is then examined, the strong/fragile liquids classification is reviewed, and a new twist-the possibility that the slowdown is a result of an avoided critical point-is noted. The three canonical characteristics of relaxing liquids are correlated through the fragility. As a further new twist, the conversion of strong liquids to fragile liquids by pressure-induced coordination number increases is demonstrated. It is then shown that, for comparable systems, it is possible to have the same conversion accomplished via a first-order transition within the liquid state during quenching. This occurs in the systems in which "polyamorphism" (polymorphism in the glassy state) is observed, and the whole phenomenology is accounted for by Poole's bond-modified van der Waals model. The sudden loss of some liquid degrees of freedom through such weak first-order transitions is then related to the polyamorphic transition between native and denatured hydrated proteins, since the latter are also glass-forming systems--water-plasticized, hydrogen bond-cross-linked chain polymers (and single molecule glass formers). The circle is closed with a final new twist by noting that a short time scale phenomenon much studied by protein physicists-namely, the onset of a sharp change in d/dT ( is the Debye-Waller factor)--is general for glass-forming liquids, including computer-simulated strong and fragile ionic liquids, and is closely correlated with the experimental glass transition temperature. The latter thus originates in strong anharmonicity in certain components of the vibrational density of states, which permits the system to access the multiple minima of its configuration space. The connection between the anharmonicity in these modes, vibrational localization, the Kauzmann temperature, and the fragility of the liquid is proposed as the key problem in glass science.

The impact of a future solar minimum on climate change projections in the Northern Hemisphere

Solar variability represents a source of uncertainty in the future forcings used in climate model simulations. Current knowledge indicates that a descent of solar activity into an extended minimum state is a possible scenario. With aid of experiments from a state-of-the-art Earth system model, we investigate the impact of a future solar minimum on Northern Hemisphere climate change projections. This scenario is constructed from recent 11 year solar-cycle minima of the solar spectral irradiance, and is therefore more conservative than the 'grand' minima employed in some previous modeling studies. Despite the small reduction in total solar irradiance (0.36 W m^-2), relatively large responses emerge in the winter Northern Hemisphere, with a reduction in regional-scale projected warming by up to 40%. To identify the origin of the enhanced regional signals, we assess the role of the different mechanisms by performing additional experiments forced only by irradiance changes at different wavelengths of the solar spectrum. We find that a reduction in visible irradiance drives changes in the stationary wave pattern of the North Pacific and sea-ice cover. A decrease in UV irradiance leads to smaller surface signals, although its regional effects are not negligible. These results point to a distinct but additive role of UV and visible irradiance in the Earth's climate, and stress the need to account for solar forcing as a source of uncertainty in regional scale projections.

Hybrid convex combinations for IIR system identification.

The low complexity of IIR adaptive filters (AFs) is specially appealing to realtime applications but some drawbacks have been preventing their widespread use so far. For gradient based IIR AFs, adverse operational conditions cause convergence problems in system identification scenarios: underdamped and clustered poles, undermodelling or non-white input signals lead to error surfaces where the adaptation nearly stops on large plateaus or get stuck at sub-optimal local minima that can not be identified as such a priori. Furthermore, the non-stationarity in the input regressor brought by the filter recursivity and the approximations made by the update rules of the stochastic gradient algorithms constrain the learning step size to small values, causing slow convergence. In this work, we propose IIR performance enhancement strategies based on hybrid combinations of AFs that achieve higher convergence rates than ordinary IIR AFs while keeping the stability.

Detección automática con snakes y representación 3D sobre imágenes cerebrales

Comunicación presentada en la VI Conferencia de la Asociación Española para la Inteligencia Artificial (CAEPIA'95), Alicante, 15-17 noviembre 1995.

Registro de excentricidad de alta frecuencia en los sedimentos continentales de la Cuenca de Guadix (Cordillera Bética, España)

Los sedimentos continentales (Plioceno-Cuaternario) que afloran en el sector central de la Cuenca de Guadix (Sur de España) muestran una ciclicidad de 100 ka consistente en la alternancia de depósitos de abanicos aluviales y sedimentos fluvio-lacustres. Durante el Plioceno y el Pleistoceno la Cuenca de Guadix era endorreica, y se caracterizaba por la existencia de un sistema axial fluvial y una orla marginal de abanicos aluviales transversales. En la zona de estudio, estos sistemas estaban relacionados lateralmente, ocupando de forma alterna el valle axial en el sector central de la Cuenca de Guadix. La edad estimada para la alternancia, ca. 100 ka, cae en la banda de excentricidad de alta frecuencia de Milankovitch. Estas fases podrían interpretarse como el resultado de máximos de excentricidad (inviernos más largos y fríos, con mayor volumen de precipitaciones que favorecerían las progradaciones de los abanicos) o de excentricidad mínima (períodos más secos y fríos, con una cubierta vegetal más escasa en las áreas fuente y, por tanto, un mayor aporte de sedimento por precipitaciones muy concentradas en el tiempo al sistema aluvial). Se muestra cómo los datos paleomagnéticos no son lo suficientemente precisos para proporcionar una buena correlación de las fases de progradación con la curva de excentricidad de Laskar et al. (2004), por lo que se pone en duda su precisión a la hora de determinar el significado climático de la ciclicidad.

The Use of the Relative Interior in Integration with Nonconvex Data

Introducing an appropriate inclusion between approximate minima associated with two nonconvex functions, we derive explicit relations between the closed convex hulls of these functions. The formula we obtain goes beyond the so-called epi-pointed property of functions which is usually concerned with such a topic.

(Appendix) Oxygen isotope record and carbonate mineralogy of the top part of ODP Hole 101-633A

Hole 633A was drilled in the southern part of Exuma Sound on the toe-of-slope of the southeastern part of Great Bahama Bank during ODP Leg 101. The top 55 m, collected as a suite of six approximately 9.5-m-long hydraulic piston cores, represents a Pliocene-Pleistocene sequence of periplatform carbonate ooze, a mixture of pelagic calcite (foraminifer and coccolith tests), some pelagic aragonite (pteropod tests), and bank-derived fine aragonite and magnesian calcite. A 1.6-m.y.-long hiatus was identified at 43.75 mbsf using calcareous nannofossil biostratigraphy and magnetostratigraphy. The 43.75-m-thick periplatform sequence above the hiatus is a complete late Pliocene-Quaternary record of the past 2.15 m.y. The d18O curve, primarily based on Globigerinoides sacculifera, clearly displays high-frequency/low-amplitude cycles during the early Pleistocene and low-frequency/high-amplitude cycles during the middle and late Pleistocene. Variations in aragonite content in the fine fraction of the periplatform ooze show a cyclic pattern throughout the Pleistocene, as previously observed in piston cores of the upper Pleistocene. These variations correlate well with the d18O record: high aragonite corresponds to light interglacial d18O values, and vice versa. Comparison of the d18O record and the aragonite curve helps to identify 23 interglacial and glacial oxygen-isotope stages, corresponding to 10.5 aragonite cycles (labeled A to K) commonly established during the middle and late Pleistocene (0.9 Ma-present). Strictly based on the aragonite curve, another 11 aragonite cycles, labeled L to V, were identified for the early Pleistocene (0.9 to 1.6 Ma). Mismatches between the d18O record and the aragonite curve occur mainly at some of the glacial-to-interglacial transitions, where aragonite increases usually lag behind d18O depletion. When one visually connects the minima on the Pleistocene aragonite curve, low-frequency (0.4 to 0.5 m.y.) supercycles seem to be superimposed on the high-frequency cycles. The timing of this supercycle roughly matches the timing of the Pleistocene carbonate preservation supercycles described in the Pacific, Indian, and Atlantic oceans. Mismatches between aragonite and d18O cycles are even more obvious for the late Pliocene (1.6 to 2.15 Ma). Irregular aragonite variations are observed for the late Pliocene, although after the onset of late Pleistocene-like glaciations in the North Atlantic Ocean 2.4 m.y. ago the d18O record has shown a mode of high-frequency/low-amplitude cycles. Initiation of climatically induced aragonite cycles occurs only at the Pliocene-Pleistocene transition, 1.6 m.y. ago. After that time, aragonite cycles are fully developed throughout the Quaternary. The 11-m-thick periplatform sequence below the hiatus represents a lower Pliocene interval between 3.75 and 4.45 Ma. The bottom half (4.25-4.45 Ma) has a fairly constant, high aragonite content (averaging 60%) and high sedimentation rates (28 m/m.y.) and corresponds to the end of the prolonged early Pliocene interglacial interval (4.1-5.0 Ma), established as a worldwide high sea-level stand. The second half (3.75-4.25 Ma), in which aragonite content decreases by successive steps, paralleled by a gradual 5180 enrichment in Globigerinoides sacculifera and low sedimentation rates (10 m/m.y), corresponds to the climatic deterioration established worldwide between 4.1 and 3.8 Ma, to a decrease of carbonate preservation observed in the equatorial Pacific Ocean, and to a global sea-level decline. Dolomite, a ubiquitous secondary component in the lower Pliocene, is interpreted as being authigenic and possibly related to diagenetic transformation of primary bank-derived fine magnesian calcite. Transformation of the primary mineralogical composition of the periplatform ooze was evidently minor, as the sediments have retained a detailed record of the Pliocene-Pleistocene climatic evolution. Clear evidence of diagenetic transformations in the periplatform ooze includes (1) the disappearance of magnesian calcite in the upper 20 m of Hole 633A, (2) the occurrence of calcite overgrowths on foraminiferal tests and microclasts at intermittent chalky core levels, and (3) the ubiquitous presence of authigenic dolomite in the lower Pliocene.