The effects of solvent composition on the affinity of a peptide towards hair keratin: experimental and molecular dynamics data

The study of the interaction between hair filaments and formulations or peptides is of utmost importance in fields like cosmetic research. Keratin intermediate filaments structure is not fully described, limiting the molecular dynamics (MD) studies in this field although its high potential to improve the area. We developed a computational model of a truncated protofibril, simulated its behavior in alcoholic based formulations and with one peptide. The simulations showed a strong interaction between the benzyl alcohol molecules of the formulations and the model, leading to the disorganization of the keratin chains, which regress with the removal of the alcohol molecules. This behavior can explain the increase of peptide uptake in hair shafts evidenced in fluorescence microscopy pictures. The model developed is valid to computationally reproduce the interaction between hair and alcoholic formulations and provide a robust base for new MD studies about hair properties. It is shown that the MD simulations can improve hair cosmetic research, improving the uptake of a compound of interest.

Influence of the hip joint modeling approaches on the kinematics of human gait

The influence of the hip joint formulation on the kinematic response of the model of human gait is investigated throughout this work. To accomplish this goal, the fundamental issues of the modeling process of a planar hip joint under the framework of multibody systems are revisited. In particular, the formulations for the ideal, dry, and lubricated revolute joints are described and utilized for the interaction of femur head inside acetabulum or the hip bone. In this process, the main kinematic and dynamic aspects of hip joints are analyzed. In a simple manner, the forces that are generated during human gait, for both dry and lubricated hip joint models, are computed in terms of the system’s state variables and subsequently introduced into the dynamics equations of motion of the multibody system as external generalized forces. Moreover, a human multibody model is considered, which incorporates the different approaches for the hip articulation, namely ideal joint, dry, and lubricated models. Finally, several computational simulations based on different approaches are performed, and the main results presented and compared to identify differences among the methodologies and procedures adopted in this work. The input conditions to the models correspond to the experimental data capture from an adult male during normal gait. In general, the obtained results in terms of positions do not differ significantly when the different hip joint models are considered. In sharp contrast, the velocity and acceleration plotted vary significantly. The effect of the hip joint modeling approach is clearly measurable and visible in terms of peaks and oscillations of the velocities and accelerations. In general, with the dry hip model, intra-joint force peaks can be observed, which can be associated with the multiple impacts between the femur head and the cup. In turn, when the lubricant is present, the system’s response tends to be smoother due to the damping effects of the synovial fluid.

Fundamental concepts in multibody dynamics

In this chapter, the fundamental ingredients related to formulation of the equations of motion for multibody systems are described. In particular, aspects such as degrees of freedom, types of coordinates, basic kinematics joints and types of analysis in multibody systems are briefly characterized. Illustrative examples of application are also presented to better clarify the fundamental issues for spatial rigid multibody systems, which are of crucial importance in the formulation development of mathematical models of mechanical systems, as well as its computational implementation.

On the frictional contacts in multibody system dynamics

"Series title: Computational methods in applied sciences, ISSN1871-3033, vol. 42"

Dinámica y cinética de los procesos físicos y químicos del radical OH con compuestos orgánicos volátiles (COVs) de relevancia atmosférica. Dynamics and kinetics of the physical and chemical processes of OH radicals with atmospheric relevant volatile organic compounds (VOCs)

El objetivo del presente proyecto es estudiar los procesos físicos y químicos del radical OH con compuestos orgánicos volátiles (COVs), con los cuales sea factible la formación de agregados de van der Waals (vdW) responsables de la curvatura en los gráficos de Arrhenius, empleando técnicas modernas, complementarias entre si y novedosas en el país. El problema será abordado desde tres perspectivas complementarias: 1) estudios cinéticos, 2) estudios mecanísticos y de distribución de productos y 3) estudios de la dinámica de los procesos físicos y químicos. La finalidad es alcanzar una mejor comprensión de los mecanismos que intervienen en el comportamiento químico de especies presentes en la atmósfera y obtener datos cinéticos de alta calidad que puedan alimentar modelos computacionales capaces de describir la composición de la atmósfera, presente y futura. Los objetivos son estudiar: 1) mediante fotólisis láser pulsada con detección por fluorescencia inducida por láser (PLP-LIF), en reactores de flujo, la cinética de reacción del radical OH(v”=0) con COVs que presentan gráficos de Arrhenius curvos con energías de activación negativas, tales como alcoholes insaturados, alquenos halogenados, éteres halogenados, ésteres alifáticos; 2) en una cámara de simulación de condiciones atmosféricas de gran volumen (4500 L), la identidad y el rendimiento de productos de las reacciones mencionadas, a fines de evaluar su impacto atmosférico y dilucidar los mecanismos de reacción; 3) mediante haces moleculares y espectroscopía láser, la estructura y reactividad de complejos de vdW entre alcoholes insaturados o aromáticos (cresoles) y el radical OH, como modelo de los aductos propuestos como responsables de la desviación al comportamiento de Arrhenius de las reacciones mencionadas; 4) mediante PLP-LIF y expansiones supersónicas, las constantes específicas estado a estado (ksts) de relajación/reacción del radical OH(v”=1-4) vibracionalmente excitado con los COVs mencionados. Los resultados experimentales obtenidos serán contrastados con cálculos ab-initio de estructura electrónica, los cuales apoyarán las interpretaciones, permitirán proponer estructuras de estados de transición y aductos colisionales, como así también calcular las frecuencias de vibración de los complejos de vdW para su posterior asignación en los espectros LIF y REMPI. Asimismo, los mecanismos de reacción propuestos y los parámetros cinéticos medidos experimentalmente serán comparados con aquellos obtenidos por cálculos teóricos. The aim of this project is to study the physical and chemical processes of OH radicals with volatile organic compounds (VOCs) with which the formation of van der Waals (vdW) clusters, responsible for the observed curvature in the Arrhenius plots, might be feasible. The problem will be addressed as follow : 1) kinetic studies; 2) products distribution and mechanistic studies and 3) dynamical studies of the physical and chemical processes. The purpose is to obtain a better understanding of the mechanisms that govern the chemical behavior of species present in the atmosphere and to obtain high quality kinetic data to be used as input to computational models. We will study: 1) the reaction kinetics of OH (v”=0) radicals with VOCs such as unsaturated alcohols, halogenated alkenes, halogenated ethers, aliphatic esters, which show curved Arrhenius plots and negative activation energies, by PLP-LIF, in flow systems; 2) in a large volume (4500 L) atmospheric simulation chamber, reaction products yields in order to evaluate their atmospheric impact and reaction mechanisms; 3) using molecular beams and laser spectroscopy, the structure and reactivity of the vdW complexes formed between the unsaturated or aromatic alcohols and the OH radicals as a model of the adducts proposed as responsible for the non-Arrhenius behavior; 4) the specific state-to-state relaxation/reaction rate constants (ksts) of the vibrationally excited OH (v”=1-4) radical with the VOCs by PLP-LIF and supersonic expansions. Ab-initio calculations will be carried out to support the interpretation of the experimental results, to obtain the transition state and collisional adducts structures, as well as to calculate the vibrational frequencies of the vdW complexes to assign to the LIF and REMPI spectra. Also, the proposed reaction mechanisms and the experimentally measured kinetic parameters will be compared with those obtained from theoretical calculations.

A Model-to-model analysis of the repeated Prisoners' Dilemma: genetic algorithms vs. evolutionary dynamics

We study the properties of the well known Replicator Dynamics when applied to a finitely repeated version of the Prisoners' Dilemma game. We characterize the behavior of such dynamics under strongly simplifying assumptions (i.e. only 3 strategies are available) and show that the basin of attraction of defection shrinks as the number of repetitions increases. After discussing the difficulties involved in trying to relax the 'strongly simplifying assumptions' above, we approach the same model by means of simulations based on genetic algorithms. The resulting simulations describe a behavior of the system very close to the one predicted by the replicator dynamics without imposing any of the assumptions of the analytical model. Our main conclusion is that analytical and computational models are good complements for research in social sciences. Indeed, while on the one hand computational models are extremely useful to extend the scope of the analysis to complex scenar

Complementarities between mathematical and computational modeling: the case of the repeated Prisoners' Dilemma

We study the properties of the well known Replicator Dynamics when applied to a finitely repeated version of the Prisoners' Dilemma game. We characterize the behavior of such dynamics under strongly simplifying assumptions (i.e. only 3 strategies are available) and show that the basin of attraction of defection shrinks as the number of repetitions increases. After discussing the difficulties involved in trying to relax the 'strongly simplifying assumptions' above, we approach the same model by means of simulations based on genetic algorithms. The resulting simulations describe a behavior of the system very close to the one predicted by the replicator dynamics without imposing any of the assumptions of the mathematical model. Our main conclusion is that mathematical and computational models are good complements for research in social sciences. Indeed, while computational models are extremely useful to extend the scope of the analysis to complex scenarios hard to analyze mathematically, formal models can be useful to verify and to explain the outcomes of computational models.

Computational drug design strategies applied to the modelling of human immunodeficiency virus-1 reverse transcriptase inhibitors

Reverse transcriptase (RT) is a multifunctional enzyme in the human immunodeficiency virus (HIV)-1 life cycle and represents a primary target for drug discovery efforts against HIV-1 infection. Two classes of RT inhibitors, the nucleoside RT inhibitors (NRTIs) and the nonnucleoside transcriptase inhibitors are prominently used in the highly active antiretroviral therapy in combination with other anti-HIV drugs. However, the rapid emergence of drug-resistant viral strains has limited the successful rate of the anti-HIV agents. Computational methods are a significant part of the drug design process and indispensable to study drug resistance. In this review, recent advances in computer-aided drug design for the rational design of new compounds against HIV-1 RT using methods such as molecular docking, molecular dynamics, free energy calculations, quantitative structure-activity relationships, pharmacophore modelling and absorption, distribution, metabolism, excretion and toxicity prediction are discussed. Successful applications of these methodologies are also highlighted.

The dynamics of Brazilian protozoology over the past century

Brazilian scientists have been contributing to the protozoology field for more than 100 years with important discoveries of new species such asTrypanosoma cruzi and Leishmania spp. In this work, we used a Brazilian thesis database (Coordination for the Improvement of Higher Education Personnel) covering the period from 1987-2011 to identify researchers who contributed substantially to protozoology. We selected 248 advisors by filtering to obtain researchers who supervised at least 10 theses. Based on a computational analysis of the thesis databases, we found students who were supervised by these scientists. A computational procedure was developed to determine the advisors’ scientific ancestors using the Lattes Platform. These analyses provided a list of 1,997 researchers who were inspected through Lattes CV examination and allowed the identification of the pioneers of Brazilian protozoology. Moreover, we investigated the areas in which researchers who earned PhDs in protozoology are now working. We found that 68.4% of them are still in protozoology, while 16.7% have migrated to other fields. We observed that support for protozoology by national or international agencies is clearly correlated with the increase of scientists in the field. Finally, we described the academic genealogy of Brazilian protozoology by formalising the “forest” of Brazilian scientists involved in the study of protozoa and their vectors over the past century.

A Negotiation Style Recommender Based on Computational Ecology in Open Negotiation Environments

The system described herein represents the first example of a recommender system in digital ecosystems where agents negotiate services on behalf of small companies. The small companies compete not only with price or quality, but with a wider service-by-service composition by subcontracting with other companies. The final result of these offerings depends on negotiations at the scale of millions of small companies. This scale requires new platforms for supporting digital business ecosystems, as well as related services like open-id, trust management, monitors and recommenders. This is done in the Open Negotiation Environment (ONE), which is an open-source platform that allows agents, on behalf of small companies, to negotiate and use the ecosystem services, and enables the development of new agent technologies. The methods and tools of cyber engineering are necessary to build up Open Negotiation Environments that are stable, a basic condition for predictable business and reliable business environments. Aiming to build stable digital business ecosystems by means of improved collective intelligence, we introduce a model of negotiation style dynamics from the point of view of computational ecology. This model inspires an ecosystem monitor as well as a novel negotiation style recommender. The ecosystem monitor provides hints to the negotiation style recommender to achieve greater stability of an open negotiation environment in a digital business ecosystem. The greater stability provides the small companies with higher predictability, and therefore better business results. The negotiation style recommender is implemented with a simulated annealing algorithm at a constant temperature, and its impact is shown by applying it to a real case of an open negotiation environment populated by Italian companies

Ab initio absorption spectrum of NiO combining molecular dynamics with the embedded cluster approach in a discrete reaction field

We developed a procedure that combines three complementary computational methodologies to improve the theoretical description of the electronic structure of nickel oxide. The starting point is a Car-Parrinello molecular dynamics simulation to incorporate vibrorotational degrees of freedom into the material model. By means ofcomplete active space self-consistent field second-order perturbation theory (CASPT2) calculations on embedded clusters extracted from the resulting trajectory, we describe localized spectroscopic phenomena on NiO with an efficient treatment of electron correlation. The inclusion of thermal motion into the theoretical description allowsus to study electronic transitions that, otherwise, would be dipole forbidden in the ideal structure and results in a natural reproduction of the band broadening. Moreover, we improved the embedded cluster model by incorporating self-consistently at the complete active space self-consistent field (CASSCF) level a discrete (or direct) reaction field (DRF) in the cluster surroundings. The DRF approach offers an efficient treatment ofelectric response effects of the crystalline embedding to the electronic transitions localized in the cluster. We offer accurate theoretical estimates of the absorption spectrum and the density of states around the Fermi level of NiO, and a comprehensive explanation of the source of the broadening and the relaxation of the charge transferstates due to the adaptation of the environment

Patient-specific computational hemodynamics of intracranial aneurysms from 3D rotational angiography and CT angiography: an in vivo reproducibility study

Patient-specific simulations of the hemodynamics in intracranial aneurysms can be constructed by using image-based vascular models and CFD techniques. This work evaluates the impact of the choice of imaging technique on these simulations

Dynamics of driven three-dimensional thin films: from hdrophilic to superhydrophobic substrates

We study the forced displacement of a thin film of fluid in contact with vertical and inclined substrates of different wetting properties, that range from hydrophilic to hydrophobic, using the lattice-Boltzmann method. We study the stability and pattern formation of the contact line in the hydrophilic and superhydrophobic regimes, which correspond to wedge-shaped and nose-shaped fronts, respectively. We find that contact lines are considerably more stable for hydrophilic substrates and small inclination angles. The qualitative behavior of the front in the linear regime remains independent of the wetting properties of the substrate as a single dispersion relation describes the stability of both wedges and noses. Nonlinear patterns show a clear dependence on wetting properties and substrate inclination angle. The effect is quantified in terms of the pattern growth rate, which vanishes for the sawtooth pattern and is finite for the finger pattern. Sawtooth shaped patterns are observed for hydrophilic substrates and low inclination angles, while finger-shaped patterns arise for hydrophobic substrates and large inclination angles. Finger dynamics show a transient in which neighboring fingers interact, followed by a steady state where each finger grows independently.