953 resultados para Path Integral, Molecular Dynamics, Statistical Mechanics
Resumo:
The intensely studied MHC has become the paradigm for understanding the architectural evolution of vertebrate multigene families. The 4-Mb human MHC (also known as the HLA complex) encodes genes critically involved in the immune response, graft rejection, and disease susceptibility. Here we report the continuous 1,796,938-bp genomic sequence of the HLA class I region, linking genes between MICB and HLA-F. A total of 127 genes or potentially coding sequences were recognized within the analyzed sequence, establishing a high gene density of one per every 14.1 kb. The identification of 758 microsatellite provides tools for high-resolution mapping of HLA class I-associated disease genes. Most importantly, we establish that the repeated duplication and subsequent diversification of a minimal building block, MIC-HCGIX-3.8–1-P5-HCGIV-HLA class I-HCGII, engendered the present-day MHC. That the currently nonessential HLA-F and MICE genes have acted as progenitors to today’s immune-competent HLA-ABC and MICA/B genes provides experimental evidence for evolution by “birth and death,” which has general relevance to our understanding of the evolutionary forces driving vertebrate multigene families.
Resumo:
The vibrational energy relaxation of carbon monoxide in the heme pocket of sperm whale myoglobin was studied by using molecular dynamics simulation and normal mode analysis methods. Molecular dynamics trajectories of solvated myoglobin were run at 300 K for both the δ- and ɛ-tautomers of the distal His-64. Vibrational population relaxation times of 335 ± 115 ps for the δ-tautomer and 640 ± 185 ps for the ɛ-tautomer were estimated by using the Landau–Teller model. Normal mode analysis was used to identify those protein residues that act as the primary “doorway” modes in the vibrational relaxation of the oscillator. Although the CO relaxation rates in both the ɛ- and δ-tautomers are similar in magnitude, the simulations predict that the vibrational relaxation of the CO is faster in the δ-tautomer with the distal His playing an important role in the energy relaxation mechanism. Time-resolved mid-IR absorbance measurements were performed on photolyzed carbonmonoxy hemoglobin (Hb13CO). From these measurements, a T1 time of 600 ± 150 ps was determined. The simulation and experimental estimates are compared and discussed.
Resumo:
We investigated the relative free energies of hapten binding to the germ line and mature forms of the 48G7 antibody Fab fragments by applying a continuum model to structures sampled from molecular dynamics simulations in explicit solvent. Reasonable absolute and very good relative free energies were obtained. As a result of nine somatic mutations that do not contact the hapten, the affinity-matured antibody binds the hapten >104 tighter than the germ line antibody. Energetic analysis reveals that van der Waals interactions and nonpolar contributions to solvation are similar and drive the formations of both the germ line and mature antibody–hapten complexes. Affinity maturation of the 48G7 antibody therefore appears to occur through reorganization of the combining site geometry in a manner that optimizes the balance of gaining favorable electrostatic interactions with the hapten and losing those with solvent during the binding process. As reflected by lower rms fluctuations in the antibody–hapten complex, the mature complex undergoes more restricted fluctuations than the germ line complex. The dramatically increased affinity of the 48G7 antibody over its germ line precursor is thus made possible by electrostatic optimization.
Resumo:
NMR analysis and molecular dynamics simulations of d(GGTAATTACC)2 and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar–phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving ‘lock and key’ recognition.
Resumo:
Classical molecular dynamics is applied to the rotation of a dipolar molecular rotor mounted on a square grid and driven by rotating electric field E(ν) at T ≃ 150 K. The rotor is a complex of Re with two substituted o-phenanthrolines, one positively and one negatively charged, attached to an axial position of Rh\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{2}^{4+}}}\end{equation*}\end{document} in a [2]staffanedicarboxylate grid through 2-(3-cyanobicyclo[1.1.1]pent-1-yl)malonic dialdehyde. Four regimes are characterized by a, the average lag per turn: (i) synchronous (a < 1/e) at E(ν) = |E(ν)| > Ec(ν) [Ec(ν) is the critical field strength], (ii) asynchronous (1/e < a < 1) at Ec(ν) > E(ν) > Ebo(ν) > kT/μ, [Ebo(ν) is the break-off field strength], (iii) random driven (a ≃ 1) at Ebo(ν) > E(ν) > kT/μ, and (iv) random thermal (a ≃ 1) at kT/μ > E(ν). A fifth regime, (v) strongly hindered, W > kT, Eμ, (W is the rotational barrier), has not been examined. We find Ebo(ν)/kVcm−1 ≃ (kT/μ)/kVcm−1 + 0.13(ν/GHz)1.9 and Ec(ν)/kVcm−1 ≃ (2.3kT/μ)/kVcm−1 + 0.87(ν/GHz)1.6. For ν > 40 GHz, the rotor behaves as a macroscopic body with a friction constant proportional to frequency, η/eVps ≃ 1.14 ν/THz, and for ν < 20 GHz, it exhibits a uniquely molecular behavior.
Resumo:
Many elementary chemical and physical processes such as the breaking of a chemical bond or the vibrational motion of atoms within a molecule take place on a femtosecond (fs = 10−15 s) or picosecond (ps = 10−12 s) time scale. It is now possible to monitor these events as a function of time with temporal resolution well below 100 fs. This capability is based on the pump-probe technique where one optical pulse triggers a reaction and a second delayed optical pulse probes the changes that ensue. To illustrate this capability, the dynamics of ligand motion within a protein are presented. Moving beyond casual observation of a reaction to active control of its outcome requires additional experimental and theoretical effort. To illustrate the concept of control, the effect of optical pulse duration on the vibrational dynamics of a tri-atomic molecule are discussed. The experimental and theoretical resources currently available are poised to make the dream of reaction control a reality for certain molecular systems.
Resumo:
A model of interdependent decision making has been developed to understand group differences in socioeconomic behavior such as nonmarital fertility, school attendance, and drug use. The statistical mechanical structure of the model illustrates how the physical sciences contain useful tools for the study of socioeconomic phenomena.
Resumo:
Molecular dynamics simulations of the oligonucleotide duplex d(CGCGCG)2 in aqueous solution are used to investigate the glass transition phenomenon. The simulations were performed at temperatures in the 20 K to 340 K range. The mean square atomic fluctuations showed that the behavior of the oligonucleotide duplex was harmonic at low temperatures. A glass transition temperature at 223 K to 234 K was inferred for the oligonucleotide duplex, which is in agreement with experimental observations. The largest number of hydrogen bounds between the polar atoms of the oligonucleotide duplex and the water molecules was obtained at the glass transition temperature. With increasing temperature we observed a decrease in the average lifetime of the hydrogen bonds to water molecules.
Resumo:
Correlations in low-frequency atomic displacements predicted by molecular dynamics simulations on the order of 1 ns are undersampled for the time scales currently accessible by the technique. This is shown with three different representations of the fluctuations in a macromolecule: the reciprocal space of crystallography using diffuse x-ray scattering data, real three-dimensional Cartesian space using covariance matrices of the atomic displacements, and the 3N-dimensional configuration space of the protein using dimensionally reduced projections to visualize the extent to which phase space is sampled.
Resumo:
Gli acidi peptido nucleici sono potenti strumenti utilizzati in ambito biotecnologico per colpire DNA o RNA. PNA contenenti basi o backbone modificati sono attualmente studiati per migliorarne le proprietà in ambito biologico. Bersagliare i micro RNA (anti-miR) è particolarmente interessante nell’ottica di future applicazioni terapeutiche, ma strumenti computazionali che aiutino nel design di nuovi PNA anti-miR non sono stati ancora completamente sviluppati. Le proprietà conformazionali del singolo filamento di PNA (non modificato o recante modificazioni in γ) e dei duplex PNA:RNA e i processi di re-annealing e melting sono stati studiati tramite Dinamica Molecolare e Metadinamica. L’approccio computazionale consolidato, assieme a un programma modificato per la generazione delle strutture dei duplex contenenti PNA, è stato utilizzato per il virtual screening di PNA contenenti basi modificate. Sono state inoltre sintetizzate le unità per l’ottenimento del composto più promettente e una funzione idrolitica da legare al monomero finale.
Resumo:
The formation and rupture of atomic-sized contacts is modelled by means of molecular dynamics simulations. Such nano-contacts are realized in scanning tunnelling microscope and mechanically controlled break junction experiments. These instruments routinely measure the conductance across the nano-sized electrodes as they are brought into contact and separated, permitting conductance traces to be recorded that are plots of conductance versus the distance between the electrodes. One interesting feature of the conductance traces is that for some metals and geometric configurations a jump in the value of the conductance is observed right before contact between the electrodes, a phenomenon known as jump-to-contact. This paper considers, from a computational point of view, the dynamics of contact between two gold nano-electrodes. Repeated indentation of the two surfaces on each other is performed in two crystallographic orientations of face-centred cubic gold, namely (001) and (111). Ultimately, the intention is to identify the structures at the atomic level at the moment of first contact between the surfaces, since the value of the conductance is related to the minimum cross-section in the contact region. Conductance values obtained in this way are determined using first principles electronic transport calculations, with atomic configurations taken from the molecular dynamics simulations serving as input structures.
Resumo:
"August 20, 1965"
Resumo:
Thesis--Illinois.
Resumo:
Mode of access: Internet.