Effects of oscillatory behavior of the dipole function on the dissociation dynamics of the classical driven Morse oscillator

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

Role of the range of the dipole function in the classical dynamics of molecular dissociation

The dissociation dynamics of heteronuclear diatomic molecules induced by infrared laser pulses is investigated within the framework of the classical driven Morse oscillator. The interaction between the molecule and the laser field described in the dipole formulation is given by the product of a time-dependent external field with a position-dependent permanent dipole function. The effects of changing the spatial range of the dipole function in the classical dissociation dynamics of large ensembles of trajectories are studied. Numerical calculations have been performed for distinct amplitudes and carrier frequencies of the external pulses and also for ensembles with different initial energies. It is found that there exist a set of values of the dipole range for which the dissociation probability can be completely suppressed. The dependence of the dissociation on the dipole range is explained through the examination of the Fourier series coefficients of the dipole function in the angle variable of the free system. In particular, the suppression of dissociation corresponds to dipole ranges for which the Fourier coefficients associated with nonlinear resonances are null and the chaotic region in the phase space is reduced to thin layers. In this context, it is shown that the suppression of dissociation of heteronuclear molecules for certain frequencies of the external field is a consequence of the finite range of the corresponding permanent dipole. © 2013 American Physical Society.

Quantum-classical correspondence and the role of the dipole function in molecular dissociation

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Controle da dinâmica caótica do oscilador de morse forçado

In this work, the dissociation dynamics of heteronuclear diatomic molecules is investigated by means of the classical driven Morse oscillator. The interaction of the molecule and the laser field is represented through the product of the molecule dipole function and the electric field of the laser. This interaction may lead to the breaking of the chemical bound, that is, to the dissociation of the molecule. The work was developed in two parts. In the first part, we studied the dissociation as a function of the range of the permanent dipole. In the second part, we maximized the dissociation probability manipulating the parameters of the external field. We have observed that the dissociation can be controlled by means of variations of parameters associated with the range of the permanent dipole

Shape resonance spectra of uracil, 5-fluorouracil, and 5-chlorouracil.

We report on the shape resonance spectra of uracil, 5-fluorouracil, and 5-chlorouracil, as obtained from fixed-nuclei elastic scattering calculations performed with the Schwinger multichannel method with pseudopotentials. Our results are in good agreement with the available electron transmission spectroscopy data, and support the existence of three π* resonances in uracil and 5-fluorouracil. As expected, the anion states are more stable in the substituted molecules than in uracil. Since the stabilization is stronger in 5-chlorouracil, the lowest π* resonance in this system becomes a bound anion state. The present results also support the existence of a low-lying σ ∗ CCl shape resonance in 5- chlorouracil. Exploratory calculations performed at selected C–Cl bond lengths suggest that the σ ∗ CCl resonance could couple to the two lowest π* states, giving rise to a very rich dissociation dynamics. These facts would be compatible with the complex branching of the dissociative electron attachment cross sections, even though we cannot discuss any details of the vibration dynamics based only on the present fixed-nuclei results.

Ligand dissociation from estrogen receptor is mediated by receptor dimerization: Evidence from molecular dynamics simulations

Estrogen Receptor (ER) is an important target for pharmaceutical design. Like other ligand-dependent transcription factors, hormone binding regulates ER transcriptional activity. Nevertheless, the mechanisms by which ligands enter and leave ERs and other nuclear receptors remain poorly understood. Here, we report results of locally enhanced sampling molecular dynamics simulations to identify dissociation pathways of two ER ligands [the natural hormone 17 beta-estradiol (E-2) and the selective ER modulator raloxifene (RAL)] from the human ER alpha ligand-binding domain in monomeric and dimeric forms. E-2 dissociation occurs via three different pathways in ER monomers. One resembles the mousetrap mechanism (Path I), involving repositioning of helix 12 (H12), others involve the separation of H8 and H11 (Path II), and a variant of this pathway at the bottom of the ligand-binding domain (Path II`). RAL leaves the receptor through Path I and a Path I variant in which the ligand leaves the receptor through the loop region between H11 and H12 (Path I`). Remarkably, ER dimerization strongly suppresses Paths II and II` for E-2 dissociation and modifies RAL escape routes. We propose that differences in ligand release pathways detected in the simulations for ER monomers and dimers provide an explanation for previously observed effects of ER quaternary state on ligand dissociation rates and suggest that dimerization may play an important, and hitherto unexpected, role in regulation of ligand dissociation rates throughout the nuclear receptor family.

Vibration and dissociation coupling behind strong shock waves. Paper presented at the Symposium on Dynamics of Manned Lifting Planetary Entry,

Cover title.

Dynamics of trapped Bose-Einstein condensates: Beyond the mean-field

Since dilute Bose gas condensates were first experimentally produced, the Gross-Pitaevskii equation has been successfully used as a descriptive tool. As a mean-field equation, it cannot by definition predict anything about the many-body quantum statistics of condensate. We show here that there are a class of dynamical systems where it cannot even make successful predictions about the mean-field behavior, starting with the process of evaporative cooling by which condensates are formed. Among others are parametric processes, such as photoassociation and dissociation of atomic and molecular condensates.

NMR structure and backbone dynamics of a concatemer of epidermal growth factor homology modules of the human low-density lipoprotein receptor

The ligand-binding region of the low-density lipoprotein (LDL) receptor is formed by seven N-terminal, imperfect, cysteine-rich (LB) modules. This segment is followed by an epidermal growth factor precursor homology domain with two N-terminal, tandem, EGF-like modules that are thought to participate in LDL binding and recycling of the endocytosed receptor to the cell surface. EGF-A and the concatemer, EGF-AB, of these modules were expressed in Escherichia coli. Correct protein folding of EGF-A and the concatemer EGF-AB was achieved in the presence or absence of calcium ions, in contrast to the LB modules, which require them for correct folding. Homonuclear and heteronuclear H-1-N-15 NMR spectroscopy at 17.6 T was used to determine the three-dimensional structure of the concatemer. Both modules are formed by two pairs of short, anti-parallel beta -strands. In the concatemer, these modules have a fixed relative orientation, stabilized by calcium ion-binding and hydrophobic interactions at the interface. N-15 longitudinal and transverse relaxation rates, and {H-1}-N-15 heteronuclear NOEs were used to derive a model-free description of the backbone dynamics of the molecule. The concatemer appears relatively rigid, particularly near the calcium ion-binding site at the module interface, with an average generalized order parameter of 0.85 +/- 0.11. Some mutations causing familial hypercholesterolemia may now be rationalized. Mutations of D41, D43 and E44 in the EGF-B calcium ion-binding region may affect the stability of the linker and thus the orientation of the tandem modules. The diminutive core also provides little structural stabilization, necessitating the presence of disulfide bonds. The structure and dynamics of EGF-AB contrast with the N-terminal LB modules, which require calcium ions both for folding to form the correct disulfide connectivities and for maintenance of the folded structure, and are connected by highly mobile linking peptides. (C) 2001 Academic Press.

Calculation of resonances and product state distributions for the unimolecular dissociation of H2S

Resonance phenomena associated with the unimolecular dissociation of H2S --> SH + H have been investigated quantum mechanically by the Lanczos homogeneous filter diagonalization method using a newly developed potential energy surface (J. Chem. Phys. 2001, 114, 320). Resonance energies, widths (rates), and product state distributions have been obtained. Both dissociation rates and product state distributions of SH show, strong fluctuations, indicating that the dissociation of H2S is essentially irregular. Statistical analysis of neighboring level spacing and width distributions also confirms this behavior. The dissociation rates and product state distributions are compared to the predictions of quantum phase space theory.

Autoinhibition of TBCB regulates EB1-mediated microtubule dynamics

Tubulin cofactors (TBCs) participate in the folding, dimerization, and dissociation pathways of the tubulin dimer. Among them, TBCB and TBCE are two CAP-Gly domain-containing proteins that together efficiently interact with and dissociate the tubulin dimer. In the study reported here we showed that TBCB localizes at spindle and midzone microtubules during mitosis. Furthermore, the motif DEI/M-COO− present in TBCB, which is similar to the EEY/F-COO− element characteristic of EB proteins, CLIP-170, and α-tubulin, is required for TBCE–TBCB heterodimer formation and thus for tubulin dimer dissociation. This motif is responsible for TBCB autoinhibition, and our analysis suggests that TBCB is a monomer in solution. Mutants of TBCB lacking this motif are derepressed and induce microtubule depolymerization through an interaction with EB1 associated with microtubule tips. TBCB is also able to bind to the chaperonin complex CCT containing α-tubulin, suggesting that it could escort tubulin to facilitate its folding and dimerization, recycling or degradation.

Protein-protein interaction investigated by steered molecular dynamics: the TCR-pMHC complex.

We present a novel steered molecular dynamics scheme to induce the dissociation of large protein-protein complexes. We apply this scheme to study the interaction of a T cell receptor (TCR) with a major histocompatibility complex (MHC) presenting a peptide (p). Two TCR-pMHC complexes are considered, which only differ by the mutation of a single amino acid on the peptide; one is a strong agonist that produces T cell activation in vivo, while the other is an antagonist. We investigate the interaction mechanism from a large number of unbinding trajectories by analyzing van der Waals and electrostatic interactions and by computing energy changes in proteins and solvent. In addition, dissociation potentials of mean force are calculated with the Jarzynski identity, using an averaging method developed for our steering scheme. We analyze the convergence of the Jarzynski exponential average, which is hampered by the large amount of dissipative work involved and the complexity of the system. The resulting dissociation free energies largely underestimate experimental values, but the simulations are able to clearly differentiate between wild-type and mutated TCR-pMHC and give insights into the dissociation mechanism.

How T cell receptors interact with peptide-MHCs: A multiple steered molecular dynamics study.

The T-cell receptor (TCR) interaction with antigenic peptides (p) presented by the major histocompatibility complex (MHC) molecule is a key determinant of immune response. In addition, TCR-pMHC interactions offer examples of features more generally pertaining to protein-protein recognition: subtle specificity and cross-reactivity. Despite their importance, molecular details determining the TCR-pMHC binding remain unsolved. However, molecular simulation provides the opportunity to investigate some of these aspects. In this study, we perform extensive equilibrium and steered molecular dynamics simulations to study the unbinding of three TCR-pMHC complexes. As a function of the dissociation reaction coordinate, we are able to obtain converged H-bond counts and energy decompositions at different levels of detail, ranging from the full proteins, to separate residues and water molecules, down to single atoms at the interface. Many observed features do not support a previously proposed two-step model for TCR recognition. Our results also provide keys to interpret experimental point-mutation results. We highlight the role of water both in terms of interface resolvation and of water molecules trapped in the bound complex. Importantly, we illustrate how two TCRs with similar reactivity and structures can have essentially different binding strategies. Proteins 2011; © 2011 Wiley-Liss, Inc.

Brain dynamics sustaining rapid rule extraction from speech

Language acquisition is a complex process that requires the synergic involvement of different cognitive functions, which include extracting and storing the words of the language and their embedded rules for progressive acquisition of grammatical information. As has been shown in other fields that study learning processes, synchronization mechanisms between neuronal assemblies might have a key role during language learning. In particular, studying these dynamics may help uncover whether different oscillatory patterns sustain more item-based learning of words and rule-based learning from speech input. Therefore, we tracked the modulation of oscillatory neural activity during the initial exposure to an artificial language, which contained embedded rules. We analyzed both spectral power variations, as a measure of local neuronal ensemble synchronization, as well as phase coherence patterns, as an index of the long-range coordination of these local groups of neurons. Synchronized activity in the gamma band (2040 Hz), previously reported to be related to the engagement of selective attention, showed a clear dissociation of local power and phase coherence between distant regions. In this frequency range, local synchrony characterized the subjects who were focused on word identification and was accompanied by increased coherence in the theta band (48 Hz). Only those subjects who were able to learn the embedded rules showed increased gamma band phase coherence between frontal, temporal, and parietal regions.

Elementary steps in aqueous proton transfer reactions : a first principles molecular dynamics study

La nature des acides dans un environnement aqueux est primordiale dans de nombreux aspects de la chimie et de la biologie. La caractéristique principale d'un acide est sa capacité à transférer un proton vers une molécule d'eau ou vers n'importe quelle base, mais ce procédé n'est pas aussi simple qu'il y paraît. Il peut au contraire être extrêmement complexe et dépendre de manière cruciale de la solvatation des différents intermédiaires de réaction impliqués. Cette thèse décrit les études computationnelles basées sur des simulations de dynamique moléculaire ab initio qui ont pour but d'obtenir une description à l'échelle moléculaire des divers procédés de transferts de proton entre acide et bases dans un milieu aqueux. Pour cela, nous avons étudié une serie de système, dont l'acide hydrofluorique aqueux, l'acide trifluoroacétique aqueux, et un système modèle constitué d'un phénol et d'une entité carboxylate reliés entre eux par une molécule d'eau en solution aqueuse. Deux états intermédiaires ont été identifiés pour le transfert d'un proton depuis un acide. Ces intermédiaires apparaissent stabilisés par un motif local de solvatation via des ponts H. Leurs signatures spectroscopiques ont été caractérisées au moyen de la spectroscopie infrarouge, en utilisant le formalisme de la dynamique moléculaire ab initio, qui inclut l'effet quantique nucléaire de manière explicite. Cette étude a aussi identifié trois chemins de réaction élémentaire, qui sont responsable pour le transfert d'un proton d'un acide à une base, ainsi que leurs échelles de temps caractéristiques. Les conclusions tirées de ces études sont discutées dans les détails, au niveau moléculaire, avec une emphase sur les comparaisons entre les résultats théoriques et les mesures expérimentales obtenues dans a littérature ou via des collaborateurs.