47 research outputs found

    Extending fragment-based free energy calculations with library Monte Carlo simulation: Annealing in interaction space

    Get PDF
    Pre-calculated libraries of molecular fragment configurations have previously been used as a basis for both equilibrium sampling (via "library-based Monte Carlo") and for obtaining absolute free energies using a polymer-growth formalism. Here, we combine the two approaches to extend the size of systems for which free energies can be calculated. We study a series of all-atom poly-alanine systems in a simple dielectric "solvent" and find that precise free energies can be obtained rapidly. For instance, for 12 residues, less than an hour of single-processor is required. The combined approach is formally equivalent to the "annealed importance sampling" algorithm; instead of annealing by decreasing temperature, however, interactions among fragments are gradually added as the molecule is "grown." We discuss implications for future binding affinity calculations in which a ligand is grown into a binding site

    Simultaneous computation of dynamical and equilibrium information using a weighted ensemble of trajectories

    Get PDF
    Equilibrium formally can be represented as an ensemble of uncoupled systems undergoing unbiased dynamics in which detailed balance is maintained. Many non-equilibrium processes can be described by suitable subsets of the equilibrium ensemble. Here, we employ the "weighted ensemble" (WE) simulation protocol [Huber and Kim, Biophys. J., 1996] to generate equilibrium trajectory ensembles and extract non-equilibrium subsets for computing kinetic quantities. States do not need to be chosen in advance. The procedure formally allows estimation of kinetic rates between arbitrary states chosen after the simulation, along with their equilibrium populations. We also describe a related history-dependent matrix procedure for estimating equilibrium and non-equilibrium observables when phase space has been divided into arbitrary non-Markovian regions, whether in WE or ordinary simulation. In this proof-of-principle study, these methods are successfully applied and validated on two molecular systems: explicitly solvated methane association and the implicitly solvated Ala4 peptide. We comment on challenges remaining in WE calculations

    The HLA class II allele DRB1*1501 is over-represented in patients with idiopathic pulmonary fibrosis

    Get PDF
    Background: Idiopathic pulmonary fibrosis (IPF) is a progressive and medically refractory lung disease with a grim prognosis. Although the etiology of IPF remains perplexing, abnormal adaptive immune responses are evident in many afflicted patients. We hypothesized that perturbations of human leukocyte antigen (HLA) allele frequencies, which are often seen among patients with immunologic diseases, may also be present in IPF patients. Methods/Principal Findings: HLA alleles were determined in subpopulations of IPF and normal subjects using molecular typing methods. HLA-DRB1*15 was over-represented in a discovery cohort of 79 Caucasian IPF subjects who had lung transplantations at the University of Pittsburgh (36.7%) compared to normal reference populations. These findings were prospectively replicated in a validation cohort of 196 additional IPF subjects from four other U.S. medical centers that included both ambulatory patients and lung transplantation recipients. High-resolution typing was used to further define specific HLA-DRB1*15 alleles. DRB1*1501 prevalence in IPF subjects was similar among the 143 ambulatory patients and 132 transplant recipients (31.5% and 34.8%, respectively, p = 0.55). The aggregate prevalence of DRB1*1501 in IPF patients was significantly greater than among 285 healthy controls (33.1% vs. 20.0%, respectively, OR 2.0; 95%CI 1.3-2.9, p = 0.0004). IPF patients with DRB1*1501 (n = 91) tended to have decreased diffusing capacities for carbon monoxide (DLCO) compared to the 184 disease subjects who lacked this allele (37.8±1.7% vs. 42.8±1.4%, p = 0.036). Conclusions/Significance: DRB1*1501 is more prevalent among IPF patients than normal subjects, and may be associated with greater impairment of gas exchange. These data are novel evidence that immunogenetic processes can play a role in the susceptibility to and/or manifestations of IPF. Findings here of a disease association at the HLA-DR locus have broad pathogenic implications, illustrate a specific chromosomal area for incremental, targeted genomic study, and may identify a distinct clinical phenotype among patients with this enigmatic, morbid lung disease

    LIM Homeodomain Factors Lhx3 and Lhx4 Assign Subtype Identities for Motor Neurons

    Get PDF
    AbstractThe circuits that control movement are comprised of discrete subtypes of motor neurons. How motor neuron subclasses develop and extend axons to their correct targets is still poorly understood. We show that LIM homeodomain factors Lhx3 and Lhx4 are expressed transiently in motor neurons whose axons emerge ventrally from the neural tube (v-MN). Motor neurons develop in embryos deficient in both Lhx3 and Lhx4, but v-MN cells switch their subclass identity to become motor neurons that extend axons dorsally from the neural tube (d-MN). Conversely, the misexpression of Lhx3 in dorsal-exiting motor neurons is sufficient to reorient their axonal projections ventrally. Thus, Lhx3 and Lhx4 act in a binary fashion during a brief period in development to specify the trajectory of motor axons from the neural tube
    corecore