Exploring unfrequent events with accelerated molecular dynamics simulations: from photochemistry to drug design

La meva incorporació al grup de recerca del Prof. McCammon (University of California San Diego) en qualitat d’investigador post doctoral amb una beca Beatriu de Pinós, va tenir lloc el passat 1 de desembre de 2010; on vaig dur a terme les meves tasques de recerca fins al darrer 1 d’abril de 2012. El Prof. McCammon és un referent mundial en l’aplicació de simulacions de dinàmica molecular (MD) en sistemes biològics d’interès humà. La contribució més important del Prof. McCammon en la simulació de sistemes biològics és el desenvolupament del mètode de dinàmiques moleculars accelerades (AMD). Les simulacions MD convencionals, les quals estan limitades a l’escala de temps del nanosegon (~10-9s), no son adients per l’estudi de sistemes biològics rellevants a escales de temps mes llargues (μs, ms...). AMD permet explorar fenòmens moleculars poc freqüents però que son clau per l’enteniment de molts sistemes biològics; fenòmens que no podrien ser observats d’un altre manera. Durant la meva estada a la “University of California San Diego”, vaig treballar en diferent aplicacions de les simulacions AMD, incloent fotoquímica i disseny de fàrmacs per ordinador. Concretament, primer vaig desenvolupar amb èxit una combinació dels mètodes AMD i simulacions Car-Parrinello per millorar l’exploració de camins de desactivació (interseccions còniques) en reaccions químiques fotoactivades. En segon lloc, vaig aplicar tècniques estadístiques (Replica Exchange) amb AMD en la descripció d’interaccions proteïna-lligand. Finalment, vaig dur a terme un estudi de disseny de fàrmacs per ordinador en la proteïna-G Rho (involucrada en el desenvolupament de càncer humà) combinant anàlisis estructurals i simulacions AMD. Els projectes en els quals he participat han estat publicats (o estan encara en procés de revisió) en diferents revistes científiques, i han estat presentats en diferents congressos internacionals. La memòria inclosa a continuació conté més detalls de cada projecte esmentat.

The dynamics of research networks in Brite-Euram

Network formation within the BRITE--EURAM program is investigated.Wedescribe the role of the hub of the network, which is defined as the setofmain contractors that account for most of the participations. We studytheeffects that the conflict of objectives within European research fundingbetween pre-competitive research vs. European cohesion has on theformationof networks and on the relationship between different partnersof the network. \\A panel data set is constructed including the second and third frameworkof theBrite--Euram program. A model of joint production of research results isusedto test for changes in the behavior of partners within the twoframeworks. \\The main findings are that participations are very concentrated, that isasmall group of institutions account for most of the participations, butgoingfrom the second to the third framework the presence of subcontractorsand singleparticipants increases substantially. This result is reinforced by the factthat main contractors receive smaller spill-ins within networks, butspill-insincrease from the second to the third framework.

Renormalization-group improvement of the spectrum of hydrogenlike atoms with massless fermions

We obtain the next-to-next-to-leading-logarithmic renormalization-group improvement of the spectrum of hydrogenlike atoms with massless fermions by using potential NRQED. These results can also be applied to the computation of the muonic hydrogen spectrum where we are able to reproduce some known double logarithms at O(m¿s6). We compare with other formalisms dealing with logarithmic resummation available in the literature.

On the many-time formulation of classical particle dynamics

Starting from the standard one-time dynamics of n nonrelativistic particles, the n-time equations of motion are inferred, and a variational principle is formulated. A suitable generalization of the classical LieKnig theorem is demonstrated, which allows the determination of all the associated presymplectic structures. The conditions under which the action of an invariance group is canonical are studied, and a corresponding Noether theorem is deduced. A formulation of the theory in terms of n first-class constraints is recovered by means of coisotropic imbeddings. The proposed approach also provides for a better understanding of the relativistic particle dynamics, since it shows that the different roles of the physical positions and the canonical variables is not peculiar to special relativity, but rather to any n-time approach: indeed a nonrelativistic no-interaction theorem is deduced.

Correlations between negative-induced cerebral dynamics and idiopathic intellectual disability: an EEG study

Intellectual disability has long been associated with deficits in socio-emotional processing. However, studies investigating brain dynamics of maladaptive socio-emotional skills associated with intellectual disability are scarce. Here, we compared differences in brain activity between low intelligence quotient (I.Q.<75, N=13) and normal controls (N=15) while evaluating their subjective emotions. Positive (P) and negative (N) valenced pictures were presented one at a time to participants of both groups, at a rate of ¾. The task required that each participant evaluate their subjective emotion and press a predefined push-button when done, alternatively P and N. Electroencephalographic (EEG) signals were continuously recorded, and the 1000ms time window following each picture was analyzed offline for power in frequency domain. Alpha low (8-10Hz) and upper (10-13Hz) frequency bands were then compared for both groups and for both P and N emotions in 12 distributed scalp electrodes. The qualitative evaluation of emotions was similar between both groups, with constant longer reaction times for the low IQ participants. The EEG signal comparison shows marked power decrease in upper alpha frequency range for N emotions in low intelligence group. Otherwise no significant difference was noticed between low and normal IQ. Main findings of the present study are (1) results do not support the hypothesis that impairment in developmental intelligence roots in maladaptive emotional processing; (2) the strong alpha power suppression during negative-induced emotions suggests the involvement of an extended neural network and more effortful inhibition processes than positive ones. We call for further studies with a larger sample.

Manual reaction times and brain dynamics after 'awake surgery' of slow-growing tumors invading the parietal area. A case report

Primary objectives: Awake surgeries of slow-growing tumours invading the brain and guided by direct electrical stimulation induce major brain reorganizations accompanied with slight impairments post-operatively. In most cases, these deficits are so slight after a few days that they are often not detectable on classical neuropsychological evaluations. Consequently, this study investigated whether simple visuo-manual reaction time paradigms would sign some level of functional asymmetries between both hemispheres. Importantly, the visual stimulus was located in the saggital plane in order to limit attentional biases and to focus mainly on the inter-hemispheric asymmetry. Methods and procedures: Three patients (aged 41, 59 and 59 years) after resections in parietal regions and a control group (age¼44, SD¼6.9) were compared during simple uni- and bimanual reaction times (RTs). Main outcomes and results: Longer RTs were observed for the contralesional compared to the ipsilesional hand in the unimanual condition. This asymmetry was reversed for the bimanual condition despite longer RTs. Conclusion and clinical implications: Reaction time paradigms are useful in these patients to monitor more precisely their functional deficits, especially their level of functional asymmetry, and to understand brain (re)organization following slowgrowing lesions.

Telomere length dynamics in normal individuals and in patients with hematopoietic stem cell-associated disorders.

The telomere length in nucleated peripheral blood (PB) cells indirectly reflects the mitotic history of their precursors: the hematopoietic stem cells (HSCs). The average length of telomeres in PB leukocytes can be measured using fluorescence in situ hybridization and flow cytometry (flow FISH). We previously used flow FISH to characterize the age-related turnover of HSCs in healthy individuals. In this review, we describe results of recent flow FISH studies in patients with selected hematopoietic stem cell-associated disorders: chronic myelogenous leukemia (CML) and several bone marrow failure syndromes. CML is characterized by a marked expansion of myeloid Philadelphia chromosome positive (Ph+) cells. Nevertheless, nonmalignant (Ph-) HSCs typically coexist in the bone marrow of CML patients. We analyzed the telomere length in > 150 peripheral blood leukocytes (PBLs) and bone marrow samples of patients with CML as well as samples of Ph- T-lymphocytes. Compared to normal controls, the overall telomere fluorescence in PBLs of patients with CML was significantly reduced. However, no telomere shortening was observed in Ph- T-lymphocytes. Patients in late chronic phase (CP) had significantly shorter telomeres than those assessed earlier in CP. Our data suggest that progressive telomere shortening is correlated with disease progression in CML. Within the group of patients with bone marrow failure syndromes, we only found significantly shortened telomeres (compared to age-adjusted controls) in granulocytes from patients with aplastic anemia (AA). Strikingly, the telomere length in granulocytes from AA patients who had recovered after immunosuppressive therapy (recAA) did not differ significantly from controls, whereas untreated patients and nonresponders with persistent severe pancytopenia (sAANR) showed marked and significant telomere shortening compared to healthy donors and patients with recAA. Furthermore, an inverse correlation between age-adjusted telomere length and peripheral blood counts was found in support of a model in which the degree of cytopenia and the amount of telomere shortening are correlated. These results support the concept of extensive proliferation of HSCs in subgroups of AA patients and suggest a potential use of telomere-length measurements as a prognostic tool in this group of disorders as well.

The Armc10/SVH gene: Genome context, regulation of mitochondrial dynamics and protection against Aβ-induced mitochondrial fragmentation

Mitochondrial function and dynamics are essential for neurotransmission, neural function and neuronal viability. Recently, we showed that the eutherian-specific Armcx gene cluster (Armcx1-6 genes), located in the X chromosome, encodes for a new family of proteins that localise to mitochondria, regulating mitochondrial trafficking. The Armcx gene cluster evolved by retrotransposition of the Armc10 gene mRNA, which is present in all vertebrates and is considered to be the ancestor gene. Here we investigate the genomic organisation, mitochondrial functions and putative neuroprotective role of the Armc10 ancestor gene. The genomic context of the Armc10 locus shows considerable syntenic conservation among vertebrates, and sequence comparisons and CHIP-data suggest the presence of at least three conserved enhancers. We also show that the Armc10 protein localises to mitochondria and that it is highly expressed in the brain. Furthermore, we show that Armc10 levels regulate mitochondrial trafficking in neurons, but not mitochondrial aggregation, by controlling the number of moving mitochondria. We further demonstrate that the Armc10 protein interacts with the KIF5/Miro1-2/Trak2 trafficking complex. Finally, we show that overexpression of Armc10 in neurons prevents A beta-induced mitochondrial fission and neuronal death. Our data suggest both conserved and differential roles of the Armc10/Armcx gene family in regulating mitochondrial dynamics in neurons, and underscore a protective effect of the Armc10 gene against A beta-induced toxicity. Overall, our findings support a further degree of regulation of mitochondrial dynamics in the brain of more evolved mammals.

The Armc10/SVH gene: Genome context, regulation of mitochondrial dynamics and protection against Aβ-induced mitochondrial fragmentation

Mitochondrial function and dynamics are essential for neurotransmission, neural function and neuronal viability. Recently, we showed that the eutherian-specific Armcx gene cluster (Armcx1-6 genes), located in the X chromosome, encodes for a new family of proteins that localise to mitochondria, regulating mitochondrial trafficking. The Armcx gene cluster evolved by retrotransposition of the Armc10 gene mRNA, which is present in all vertebrates and is considered to be the ancestor gene. Here we investigate the genomic organisation, mitochondrial functions and putative neuroprotective role of the Armc10 ancestor gene. The genomic context of the Armc10 locus shows considerable syntenic conservation among vertebrates, and sequence comparisons and CHIP-data suggest the presence of at least three conserved enhancers. We also show that the Armc10 protein localises to mitochondria and that it is highly expressed in the brain. Furthermore, we show that Armc10 levels regulate mitochondrial trafficking in neurons, but not mitochondrial aggregation, by controlling the number of moving mitochondria. We further demonstrate that the Armc10 protein interacts with the KIF5/Miro1-2/Trak2 trafficking complex. Finally, we show that overexpression of Armc10 in neurons prevents A beta-induced mitochondrial fission and neuronal death. Our data suggest both conserved and differential roles of the Armc10/Armcx gene family in regulating mitochondrial dynamics in neurons, and underscore a protective effect of the Armc10 gene against A beta-induced toxicity. Overall, our findings support a further degree of regulation of mitochondrial dynamics in the brain of more evolved mammals.

Nuclear magnetic resonance investigation of structure and dynamics of some electrolyte solutions

Electrolyte solutions are of importance in a wide range of scientific contexts and as such have attracted considerable theoretical and experimental effort over many years. Nuclear Magnetic resonance provides a precise and versatile tool for investigation of electrolyte solutions, both in water and in organic solvents. Many structural and dynamic properties can be obtained through NMR experiments. The solution of aluminum chloride in water was studied. Different concentrations were taken for investigation. Independence of maximum line shift from concentration and acidity was shown. Six-coordinated structure of solvation shell was confirmed by experiments on 'H and 27A1 nuclei. Diffusion coefficients were studied. The solution of nickel chloride in methanol was studied. Lines, corresponding to coordinated and bulk methanol were found. Four-, five- and six-coordinated structures were found in different temperatures. The line for coordinated -OD group of deuterated methanol was observed on 2H spectrum for the first time. Partial deuteration of CH3 group was detected. Inability to observe coordinated -OH group was explained.

The dynamics of biofilm bacterial communities is driven by flow wax and wane in a temporary stream

Biofilm communities are exposed to long periods of desiccation in temporary streams. We investigated how water flow intermittency affected the bacterial community structure colonizing three different streambed compartments in a Mediterranean stream. Massive parallel sequencing revealed different bacterial communities in biofilms from sand sediments and cobbles. Bacterial communities were similar (62% of shared operational taxonomic units) in the epipsammic and hyporheic biofilms, and more diverse than those in the epilithic biofilms. The non-flow phase caused a decrease of bacterial diversity in the biofilms, when communities included only bacterial taxa assumed to be adapted to water stress. The most sensitive bacterial communities to flow intermittency were in the epilithic, where the exposure to physical stress was the highest. In sand sediments a wide group of bacterial taxa was tolerant to desiccation. During non-flow the proliferation of opportunistic taxa in the superficial compartments evidenced the biological link with the terrestrial environment. Bacterial communities better tolerate rewetting than desiccation, since a major number of taxa tolerant to rewetting occurred in all biofilms. Overall, bacterial communities in sandy compartments showed higher resistance to flow intermittency than those in epilithic biofilms

Trypanosoma vivax infection dynamics in a cattle herd maintained in a transition area between Pantanal lowlands and highlands of Mato Grosso do Sul, Brazil

Trypanosoma vivax outbreaks in beef cattle in the Pantanal region of Mato Grosso do Sul state, Brazil, causes relevant economical impact due to weight loss, abortion and mortality. Cattle moved from the Pantanal to adjacent areas of this ecosystem for breeding and fattening is a common feature. Therefore an epidemiological study on breeding cows in the transition area between Pantanal lowland and adjacent highlands of Mato Grosso do Sul was performed to determine the T. vivax infection dynamics and outbreak risk. Three experimental groups were formed: Group 1 consisted of cows parasitologically negative by the Woo test and in the enzyme-linked immunosorbent assay for T. vivax antibody detection (Tv-ELISA-Ab); Group 2 parasitologically negative and positive in the Tv-ELISA-Ab; and in Group 3 cows were parasitologically positive and with positive reactions in the Tv-ELISA-Ab. During 24 months, the cows' dislodgment between the above established groups was monitored by Woo test and Tv-ELISA-Ab exams. The tabanid population was also monitored and the highest number occurred during the rainy season. Although parasitemias were detected only in the first four samplings of the experimental period, the cows could be considered as trypanotolerant, because no clinical signs were observed. Despite the higher T. vivax incidence during the dry season, no disease symptoms were seen. Even though T. vivax epidemiological situation in the herd was characterized as endemic with seasonal variation, the probability of outbreaks was null within the conditions of the study.

On the numerical integration of rigid body nonlinear dynamics in presence of parameters singularities

One of the main complexities in the simulation of the nonlinear dynamics of rigid bodies consists in describing properly the finite rotations that they may undergo. It is well known that, to avoid singularities in the representation of the SO(3) rotation group, at least four parameters must be used. However, it is computationally expensive to use a four-parameters representation since, as only three of the parameters are independent, one needs to introduce constraint equations in the model, leading to differential-algebraic equations instead of ordinary differential ones. Three-parameter representations are numerically more efficient. Therefore, the objective of this paper is to evaluate numerically the influence of the parametrization and its singularities on the simulation of the dynamics of a rigid body. This is done through the analysis of a heavy top with a fixed point, using two three-parameter systems, Euler's angles and rotation vector. Theoretical results were used to guide the numerical simulation and to assure that all possible cases were analyzed. The two parametrizations were compared using several integrators. The results show that Euler's angles lead to faster integration compared to the rotation vector. An Euler's angles singular case, where representation approaches a theoretical singular point, was analyzed in detail. It is shown that on the contrary of what may be expected, 1) the numerical integration is very efficient, even more than for any other case, and 2) in spite of the uncertainty on the Euler's angles themselves, the body motion is well represented.