Computational investigations of ion conduction mechanisms using enhanced sampling methods

Proper ion channels’ functioning is a prerequisite for a normal cell and disorders involving ion channels, or channelopathies, underlie many human diseases. Long QT syndromes (LQTS) for example may arise from the malfunctioning of hERG channel, caused either by the binding of drugs or mutations in HERG gene. In the first part of this thesis I present a framework to investigate the mechanism of ion conduction through hERG channel. The free energy profile governing the elementary steps of ion translocation in the pore was computed by means of umbrella sampling simulations. Compared to previous studies, we detected a different dynamic behavior: according to our data hERG is more likely to mediate a conduction mechanism which has been referred to as “single-vacancy-like” by Roux and coworkers (2001), rather then a “knock-on” mechanism. The same protocol was applied to a model of hERG presenting the Gly628Ser mutation, found to be cause of congenital LQTS. The results provided interesting insights about the reason of the malfunctioning of the mutant channel. Since they have critical functions in viruses’ life cycle, viral ion channels, such as M2 proton channel, are considered attractive targets for antiviral therapy. A deep knowledge of the mechanisms that the virus employs to survive in the host cell is of primary importance in the identification of new antiviral strategies. In the second part of this thesis I shed light on the role that M2 plays in the control of electrical potential inside the virus, being the charge equilibration a condition required to allow proton influx. The ion conduction through M2 was simulated using metadynamics technique. Based on our results we suggest that a potential anion-mediated cation-proton exchange, as well as a direct anion-proton exchange could both contribute to explain the activity of the M2 channel.

Molecular interaction between perthiolated [beta]-cyclodextrin (CD) and the guests molecules adamantaneacetic acid (AD) and ferroceneacetic acid (FC); and the effect of the interaction on the electron transition of CD anchored particles

The molecular interactions between the host molecule, perthiolated beta-cyclodextrin (CD), and the guest molecules, adamantaneacetic acid (AD) and ferroceneacetic acid (FC), have been inestigated theoretically in both the gas and aqueous phases. The major computations have been carried out at the theoretical levels, RHF/6-31G and B3LYP/6- 31G. MP2 electronic energies were also computed based at the geometries optimized by both the RHF and B3LYP methods in the gas phase to establish a better estimate of the correlation effect. The solvent phase computations were completed at the RHF/6-31G and B3LYP/6-31G levels using the PCM model. The most stable structures optimized in gas phase by both the RHF and B3LYP methods were used for the computations in solution. A method to systematically manipulate the relative position and orientation between the interacting molecules is proposed. In the gas phase, six trials with different host-guest relative positions and orientations were completed successfully with the B3LYP method for both the CD-AD and CD-FC complexes. Only four trials were completed with RHF method. In the gas phase, the best results from the RHF method gives for the association Gibbs free energy (ΔG°) values equal to -32.21kj/mol for CD-AD and -25.73kj/mol for CD-FC. And the best results from the B3LYP method have ΔG° equal to -47.57kj/mol for CD-AD and -41.09kj/mol for CD-FC. The MP2 correction significantly lowers ΔG° based on the geometries from both methods. For the RHF structure, the MP2 computations lowered ΔG° to -60.64kj/mol for CD-AD and -54.10 for CD-FC. For the structure from the B3LYP method, it was reduced to -59.87 kj/mol for CD-AD and -54.84 kj/mol for CDFC. The RHF solvent phase calculations yielded following results: ΔG°(aq) equals 107.2kj/mol for CD-AD and 111.4kj/mol for CD-FC. Compared with the results from the RHF method, the B3LYP method provided clearly better solvent phase results with ΔG° (aq) equal to 38.64kj/mol for CD-AD and 39.61kj/mol for CD-FC. These results qualitatively explain the experimental observations. However quantitatively they are in poor agreement with the experimental values available in the literature and those recently published by Liu et al. And the reason is believed to be omission of hydrophobic contribution to the association. Determining the global geometrical minima for these very large systems was very difficult and computationally time consuming, but after a very thorough search, these were identified. A relevant result of this search is that when the complexes, CD-AD and CD-FC, are formed, the AD and FC molecules are only partially embedded inside the CD cavity. The totally embedded complexes were found to have significantly higher energies. The semiempirical method, ZINDO, was employed to investigate the effect of complexation on the first electronic excitation of CD anchored to a metal nano-particle. The computational results revealed that after complexation to FC, the transition intensity declines to about 25% of the original value, and after complexation with AD, the intensity drops almost 50%. The tighter binding and transition intensity of CD-AD qualitatively agrees with the experimental result that the addition of AD to a solution of CD and FC restores the fluorescence of CD that was quenched by the addition of FC. A method to evaluate the “hydrophobic force” effect is proposed for future work.

UNDERSTANDING PEOPLE'S ADOPTION OF INDIVIDUAL-USE RENEWABLE ENERGY TECHNOLOGY

This thesis attempts to understand why people adopt or reject individual-use renewable energy technologies (IURET). I used factors from Everett Rogers' Diffusion of Innovation Theory to understand how people's perceptions towards the characteristics of a given IURET (such as price, compatibility, complexity, etc.), the characteristics of the individual adopter (such as innovativeness and environmental awareness), and the communication network (inter-personal communications and mass media) can influence adoption. An online questionnaire was sent to 101randomly selected Michigan households (using random digit dialing) to ask people whether or not they had adopted at least one IURET and to assess the above-mentioned factors from Rogers' theory. Data analysis was then conducted in SPSS using Chi-squared and binary logistic regression to determine the relationship between adoption behaviors (the dependent variable) and the factors from Rogers' theory (the independent variables) while controlling for education. The results show that Rogers' factors of price and observability and the control variable of education were all significant in explaining adoption but the other factors of Rogers' theory were not. For example, if individuals perceive the price of IURET to be reasonable or if they observe their neighbors using these technologies, then they are more likely to adopt. These results indicate that, if we want to promote greater adoption of IURET, we should focus our efforts on making the price of IURET more affordable through incentives and other mechanisms. Adopters should also be given some form of reward if they provide free demonstrations of their IURET in use to their neighbors to take advantage of the observability effects.

The free escape continuum of diffuse ions upstream of the Earth's quasi-parallel bow shock

The Earth's bow shock is very efficient in accelerating ions out of the incident solar wind distribution to high energies (≈ 200 keV/e). Fluxes of energetic ions accelerated at the quasi-parallel bow shock, also known as diffuse ions, are best represented by exponential spectra in energy/charge, which require additional assumptions to be incorporated into these model spectra. One of these assumptions is a so-called "free escape boundary" along the interplanetary magnetic field into the upstream direction. Locations along the IBEX orbit are ideally suited for in situ measurements to investigate the existence of an upstream free escape boundary for bow shock accelerated ions. In this study we use 2 years of ion measurements from the background monitor on the IBEX spacecraft, supported by ACE solar wind observations. The IBEX Background Monitor is sensitive to protons > 14 keV, which includes the energy of the maximum flux for diffuse ions. With increasing distance from the bow shock along the interplanetary magnetic field, the count rates for diffuse ions stay constant for ions streaming away from the bow shock, while count rates for diffuse ions streaming toward the shock gradually decrease from a maximum value to ~1/e at distances of about 10 RE to 14 RE. These observations of a gradual decrease support the transition to a free escape continuum for ions of energy >14 keV at distances from 10 RE to 14 RE from the bow shock.

Tracking individual membrane proteins and their biochemistry: The power of direct observation

The advent of single molecule fluorescence microscopy has allowed experimental molecular biophysics and biochemistry to transcend traditional ensemble measurements, where the behavior of individual proteins could not be precisely sampled. The recent explosion in popularity of new super-resolution and super-localization techniques coupled with technical advances in optical designs and fast highly sensitive cameras with single photon sensitivity and millisecond time resolution have made it possible to track key motions, reactions, and interactions of individual proteins with high temporal resolution and spatial resolution well beyond the diffraction limit. Within the purview of membrane proteins and ligand gated ion channels (LGICs), these outstanding advances in single molecule microscopy allow for the direct observation of discrete biochemical states and their fluctuation dynamics. Such observations are fundamentally important for understanding molecular-level mechanisms governing these systems. Examples reviewed here include the effects of allostery on the stoichiometry of ligand binding in the presence of fluorescent ligands; the observation of subdomain partitioning of membrane proteins due to microenvironment effects; and the use of single particle tracking experiments to elucidate characteristics of membrane protein diffusion and the direct measurement of thermodynamic properties, which govern the free energy landscape of protein dimerization. The review of such characteristic topics represents a snapshot of efforts to push the boundaries of fluorescence microscopy of membrane proteins to the absolute limit.

The solution structure of the Zα domain of the human RNA editing enzyme ADAR1 reveals a prepositioned binding surface for Z-DNA

Double-stranded RNA deaminase I (ADAR1) contains the Z-DNA binding domain Zα. Here we report the solution structure of free Zα and map the interaction surface with Z-DNA, confirming roles previously assigned to residues by mutagenesis. Comparison with the crystal structure of the (Zα)2/Z-DNA complex shows that most Z-DNA contacting residues in free Zα are prepositioned to bind Z-DNA, thus minimizing the entropic cost of binding. Comparison with homologous (α+β)helix–turn–helix/B-DNA complexes suggests that binding of Zα to B-DNA is disfavored by steric hindrance, but does not eliminate the possibility that related domains may bind to both B- and Z-DNA.

Multiple domains of repressor activator protein 1 contribute to facilitated binding of glycolysis regulatory protein 1

The function of repressor activator protein 1 (Rap1p) at glycolytic enzyme gene upstream activating sequence (UAS) elements in Saccharomyces cerevisiae is to facilitate binding of glycolysis regulatory protein 1 (Gcr1p) at adjacent sites. Rap1p has a modular domain structure. In its amino terminus there is an asymmetric DNA-bending domain, which is distinct from its DNA-binding domain, which resides in the middle of the protein. In the carboxyl terminus of Rap1p lie its silencing and putative activation domains. We carried out a molecular dissection of Rap1p to identify domains contributing to its ability to facilitate binding of Gcr1p. We prepared full-length and three truncated versions of Rap1p and tested their ability to facilitate binding of Gcr1p by gel shift assay. The ability to detect ternary complexes containing Rap1p⋅DNA⋅Gcr1p depended on the presence of binding sites for both proteins in the probe DNA. The DNA-binding domain of Rap1p, although competent to bind DNA, was unable to facilitate binding of Gcr1p. Full-length Rap1p and the amino- and carboxyl-truncated versions of Rap1p were each able to facilitate binding of Gcr1p at an appropriately spaced binding site. Under these conditions, Gcr1p displayed an approximately 4-fold greater affinity for Rap1p-bound DNA than for otherwise identical free DNA. When spacing between Rap1p- and Gcr1p-binding sites was altered by insertion of five nucleotides, the ability to form ternary Rap1p⋅DNA⋅Gcr1p complexes was inhibited by all but the DNA-binding domain of Rap1p itself; however, the ability of each individual protein to bind the DNA probe was unaffected.

A streptavidin mutant with altered ligand-binding specificity

The biotin-binding site of streptavidin was modified to alter its ligand-binding specificity. In natural streptavidin, the side chains of N23 and S27 make two of the three hydrogen bonds with the ureido oxygen of biotin. These two residues were mutated to severely weaken biotin binding while attempting to maintain the affinity for two biotin analogs, 2-iminobiotin and diaminobiotin. Redesigning of the biotin-binding site used the difference in local electrostatic charge distribution between biotin and these biotin analogs. Free energy calculations predicted that the introduction of a negative charge at the position of S27 plus the mutation N23A should disrupt two of the three hydrogen bonds between natural streptavidin and the ureido oxygen of biotin. In contrast, the imino hydrogen of 2-iminobiotin should form a hydrogen bond with the side chain of an acidic amino acid at position 27. This should reduce the biotin-binding affinity by approximately eight orders of magnitude, while leaving the affinities for these biotin analogs virtually unaffected. In good agreement with these predictions, a streptavidin mutant with the N23A and S27D substitutions binds 2-iminobiotin with an affinity (Ka) of 1 × 106 M−1, two orders of magnitude higher than that for biotin (1 × 104 M−1). In contrast, the binding affinity of this streptavidin mutant for diaminobiotin (2.7 × 104 M−1) was lower than predicted (2.9 × 105 M−1), suggesting the position of the diaminobiotin in the biotin-binding site was not accurately determined by modeling.

A signature of the T → R transition in human hemoglobin

Allosteric effects in hemoglobin arise from the equilibrium between at least two energetic states of the molecule: a tense state, T, and a relaxed state, R. The two states differ from each other in the number and energy of the interactions between hemoglobin subunits. In the T state, constraints between subunits oppose the structural changes resulting from ligand binding. In the R state, these constraints are released, thus enhancing ligand-binding affinity. In the present work, we report the presence of four sites in hemoglobin that are structurally stabilized in the R relative to the T state. These sites are Hisα103(G10) and Hisα122(H5) in each α subunit of hemoglobin. They are located at the α1β1 and α2β2 interfaces of the hemoglobin tetramer, where the histidine side chains form hydrogen bonds with specific residues from the β chains. We have measured the solvent exchange rates of side chain protons of Hisα103(G10) and Hisα122(H5) in both deoxygenated and ligated hemoglobin by NMR spectroscopy. The exchange rates were found to be higher in the deoxygenated-T than in ligated-R state. Analysis of exchange rates in terms of the local unfolding model revealed that the structural stabilization free energy at each of these two histidines is larger by ≈1.5 kcal/(mol tetramer) in the R relative to the T state. The location of these histidines at the intradimeric α1β1 and α2β2 interfaces also suggests a role for these interfaces in the allosteric equilibrium of hemoglobin.

Mechanism for nucleic acid chaperone activity of HIV-1 nucleocapsid protein revealed by single molecule stretching

The nucleocapsid protein (NC) of HIV type 1 is a nucleic acid chaperone that facilitates the rearrangement of nucleic acids into conformations containing the maximum number of complementary base pairs. We use an optical tweezers instrument to stretch single DNA molecules from the helix to coil state at room temperature in the presence of NC and a mutant form (SSHS NC) that lacks the two zinc finger structures present in NC. Although both NC and SSHS NC facilitate annealing of complementary strands through electrostatic attraction, only NC destabilizes the helical form of DNA and reduces the cooperativity of the helix-coil transition. In particular, we find that the helix-coil transition free energy at room temperature is significantly reduced in the presence of NC. Thus, upon NC binding, it is likely that thermodynamic fluctuations cause continuous melting and reannealing of base pairs so that DNA strands are able to rapidly sample configurations to find the lowest energy state. The reduced cooperativity allows these fluctuations to occur in the middle of complex double-stranded structures. The reduced stability and cooperativity, coupled with the electrostatic attraction generated by the high charge density of NC, is responsible for the nucleic acid chaperone activity of this protein.

Mimicry of the calcium-induced conformational state of troponin C by low temperature under pressure.

Calcium binding to the N-domain of troponin C initiates a series of conformational changes that lead to muscle contraction. Calcium binding provides the free energy for a hydrophobic region in the core of N-domain to assume a more open configuration. Fluorescence measurements on a tryptophan mutant (F29W) show that a similar conformational change occurs in the absence of Ca2+ when the temperature is lowered under pressure. The conformation induced by subzero temperatures binds the hydrophobic probe bis-aminonaphthalene sulfonate, and the tryptophan has the same fluorescence lifetime (7 ns) as in the Ca2+-bound form. The decrease in volume (delta V = -25.4 ml/mol) corresponds to an increase in surface area. Thermodynamic measurements suggest an enthalpy-driven conformational change that leads to an intermediate with an exposed N-domain core and a high affinity for Ca2+.

Unraveling principles of lead discovery: from unfrustrated energy landscapes to novel molecular anchors.

The search for novel leads is a critical step in the drug discovery process. Computational approaches to identify new lead molecules have focused on discovering complete ligands by evaluating the binding affinity of a large number of candidates, a task of considerable complexity. A new computational method is introduced in this work based on the premise that the primary molecular recognition event in the protein binding site may be accomplished by small core fragments that serve as molecular anchors, providing a structurally stable platform that can be subsequently tailored into complete ligands. To fulfill its role, we show that an effective molecular anchor must meet both the thermodynamic requirement of relative energetic stability of a single binding mode and its consistent kinetic accessibility, which may be measured by the structural consensus of multiple docking simulations. From a large number of candidates, this technique is able to identify known core fragments responsible for primary recognition by the FK506 binding protein (FKBP-12), along with a diverse repertoire of novel molecular cores. By contrast, absolute energetic criteria for selecting molecular anchors are found to be promiscuous. A relationship between a minimum frustration principle of binding energy landscapes and receptor-specific molecular anchors in their role as "recognition nuclei" is established, thereby unraveling a mechanism of lead discovery and providing a practical route to receptor-biased computational combinatorial chemistry.

PR-39, a proline-rich antibacterial peptide that inhibits phagocyte NADPH oxidase activity by binding to Src homology 3 domains of p47 phox.

Reactive oxygen intermediates generated by the phagocyte NADPH oxidase are critically important components of host defense. However, these highly toxic oxidants can cause significant tissue injury during inflammation; thus, it is essential that their generation and inactivation are tightly regulated. We show here that an endogenous proline-arginine (PR)-rich antibacterial peptide, PR-39, inhibits NADPH oxidase activity by blocking assembly of this enzyme through interactions with Src homology 3 domains of a cytosolic component. This neutrophil-derived peptide inhibited oxygen-dependent microbicidal activity of neutrophils in whole cells and in a cell-free assay of NADPH oxidase. Both oxidase inhibitory and direct antimicrobial activities were defined within the amino-terminal 26 residues of PR-39. Oxidase inhibition was attributed to binding of PR-39 to the p47phox cytosolic oxidase component. Its effects involve both a polybasic amino-terminal segment and a proline-rich core region of PR-39 that binds to the p47phox Src homology 3 domains and, thereby, inhibits interaction with the small subunit of cytochrome b558, p22phox. These findings suggest that PR-39, which has been shown to be involved in tissue repair processes, is a multifunctional peptide that can regulate NADPH oxidase production of superoxide anion O2-. thus limiting excessive tissue damage during inflammation.

Identification, purification, and molecular cloning of autonomously replicating sequence-binding protein 1 from fission yeast Schizosaccharomyces pombe.

Autonomously replicating sequence (ARS) elements of the fission yeast Schizosaccharomyces pombe contain multiple imperfect copies of the consensus sequence reported by Maundrell et al. [Maundrell K., Hutchison, A. & Shall, S. (1988) EMBO J. 7, 2203-2209]. When cell free extracts of S. pombe were incubated with a dimer or tetramer of an oligonucleotide containing the ARS consensus sequence, several complexes were detected using a gel mobility-shift assay. The proteins forming these complexes also bind ars3002, which is the most active origin in the ura4 region of chromosome III of S. pombe. One protein, partly responsible for the binding activity observed with crude extracts, was purified to near homogeneity. It is a 60-kDa protein and was named ARS-binding protein 1 (Abp1). Abp1 preferentially binds to multiple sites in ARS 3002 and to the DNA polymer poly[d(A.T)]. The cloning and sequence of the gene coding for Abp1 revealed that it encodes a protein of 59.8 kDa (522 amino acids). Abp1 has significant homology (25% identity, 50% similarity) to the N-terminal region (approximately 300 amino acids) of the human and mouse centromere DNA-binding protein CENP-B. Because centromeres of S. pombe contain a high density of ARS elements, Abp1 may play a role connecting DNA replication and chromosome segregation.

Structure, dynamics, and energetics of siRNA-cationic vector complexation:a molecular dynamics study

The design and synthesis of safe and efficient nonviral vectors for gene delivery has attracted significant attention in recent years. Previous experiments have revealed that the charge density of a polycation (the carrier) plays a crucial role in complexation and the release of the gene from the complex in the cytosol. In this work, we adopt an atomistic molecular dynamics simulation approach to study the complexation of short strand duplex RNA with six cationic carrier systems of varying charge and surface topology. The simulations reveal detailed molecular-level pictures of the structures and dynamics of the RNA-polycation complexes. Estimates for the binding free energy indicate that electrostatic contributions are dominant followed by van der Waals interactions. The binding free energy between the 8(+)polymers and the RNA is found to be larger than that of the 4(+)polymers, in general agreement with previously published data. Because reliable binding free energies provide an effective index of the ability of the polycationic carrier to bind the nucleic acid and also carry implications for the process of gene release within the cytosol, these novel simulations have the potential to provide us with a much better understanding of key mechanistic aspects of gene-polycation complexation and thereby advance the rational design of nonviral gene delivery systems.