Giacomo Fiorin

Exploring Multidimensional Free Energy Landscapes Using Time-Dependent Biases on Collective Variables

A new implementation of the adaptive biasing force (ABF) method is described. This implementation supports a wide range of collective variables and can be applied to the computation of multidimensional energy profiles. It is provided to the community as part of a code that implements several analogous methods, including metadynamics. ABF and metadynamics have not previously been tested side by side on identical systems. Here, numerical tests are carried out on processes including conformational changes in model peptides and translocation of a halide ion across a lipid membrane through a peptide nanotube. On the basis of these examples, we discuss similarities and differences between the ABF and metadynamics schemes. Both approaches provide enhanced sampling and free energy profiles in quantitative agreement with each other in different applications. The method of choice depends on the dimension of the reaction coordinate space, the height of the barriers, and the relaxation times of degrees of freedom in the orthogonal space, which are not explicitly described by the chosen collective variables.

Using collective variables to drive molecular dynamics simulations

A software framework is introduced that facilitates the application of biasing algorithms to collective variables of the type commonly employed to drive massively parallel molecular dynamics (MD) simulations. The modular framework that is presented enables one to combine existing collective variables into new ones, and combine any chosen collective variable with available biasing methods. The latter include the classic time-dependent biases referred to as steered MD and targeted MD, the temperature-accelerated MD algorithm, as well as the adaptive free-energy biases called metadynamics and adaptive biasing force. The present modular software is extensible, and portable between commonly used MD simulation engines.

Structure and mechanism of proton transport through the transmembrane tetrameric M2 protein bundle of the influenza A virus

The M2 proton channel from influenza A virus is an essential protein that mediates transport of protons across the viral envelope. This protein has a single transmembrane helix, which tetramerizes into the active channel. At the heart of the conduction mechanism is the exchange of protons between the His37 imidazole moieties of M2 and waters confined to the M2 bundle interior. Protons are conducted as the total charge of the four His37 side chains passes through 2+ and 3+ with a pKa near 6. A 1.65 Å resolution X-ray structure of the transmembrane protein (residues 25–46), crystallized at pH 6.5, reveals a pore that is lined by alternating layers of sidechains and well-ordered water clusters, which offer a pathway for proton conduction. The His37 residues form a box-like structure, bounded on either side by water clusters with well-ordered oxygen atoms at close distance. The conformation of the protein, which is intermediate between structures previously solved at higher and lower pH, suggests a mechanism by which conformational changes might facilitate asymmetric diffusion through the channel in the presence of a proton gradient. Moreover, protons diffusing through the channel need not be localized to a single His37 imidazole, but instead may be delocalized over the entire His-box and associated water clusters. Thus, the new crystal structure provides a possible unification of the discrete site versus continuum conduction models.

Structure and inhibition of the drug-resistant S31N mutant of the M2 ion channel of influenza A virus

The influenza A virus M2 proton channel (A/M2) is the target of the antiviral drugs amantadine and rimantadine, whose use has been discontinued due to widespread drug resistance. Among the handful of drug-resistant mutants, S31N is found in more than 95% of the currently circulating viruses and shows greatly decreased inhibition by amantadine. The discovery of inhibitors of S31N has been hampered by the limited size, polarity, and dynamic nature of its amantadine-binding site. Nevertheless, we have discovered small-molecule drugs that inhibit S31N with potencies greater than amantadine’s potency against WT M2. Drug binding locks the protein into a well-defined conformation, and the NMR structure of the complex shows the drug bound in the homotetrameric channel, threaded between the side chains of Asn31. Unrestrained molecular dynamics simulations predicted the same binding site. This S31N inhibitor, like other potent M2 inhibitors, contains a charged ammonium group. The ammonium binds as a hydrate to one of three sites aligned along the central cavity that appear to be hotspots for inhibition. These sites might stabilize hydronium-like species formed as protons diffuse through the outer channel to the proton-shuttling residue His37 near the cytoplasmic end of the channel.

Functional studies and modeling of pore-lining residue mutants of the influenza A virus M2 ion channel

The A/M2 protein of influenza A virus forms a tetrameric proton-selective pH-gated ion channel. The H(37)xxxW(41) motif located in the channel pore is responsible for its gating and proton selectivity. Channel activation most likely involves protonation of the H37 residues, while the conductive state of the channel is characterized by two or three charged His residues in a tetrad. A/M2 channel activity is inhibited by the antiviral drug amantadine. Although a large number of functional amantadine-resistant mutants of A/M2 have been observed in vitro, only a few are observed in highly transmissible viruses in the presence or absence of amantadine. We therefore examined 49 point mutants of the pore-lining residues, representing both natural and nonnatural variants. Their ion selectivity, amantadine sensitivity, specific activity, and pH-dependent conductance were measured in Xenopus oocytes. These measurements showed how variations in the sequence lead to variations in the proton conduction. The results are consistent with a multistep mechanism that allows the protein to fine-tune its pH-rate profile over a wide range of proton concentrations, hypothesized to arise from different protonation states of the H37 tetrad. Mutations that give native-like conductance at low pH as well as minimal leakage current at pH 7.0 were surprisingly rare. Moreover, the results are consistent with a location of the amantadine-binding site inside the channel pore. These findings have helped to define the set of functionally fit mutants that should be targeted when considering the design of novel drugs that inhibit amantadine-resistant strains of influenza A virus.

Discovery of Novel Dual Inhibitors of the Wild-Type and the Most Prevalent Drug-Resistant Mutant, S31N, of the M2 Proton Channel from Influenza A Virus

Anti-influenza drugs, amantadine and rimantadine, targeting the M2 channel from influenza A virus are no longer effective because of widespread drug resistance. S31N is the predominant and amantadine-resistant M2 mutant, present in almost all of the circulating influenza A strains as well as in the pandemic 2009 H1N1 and the highly pathogenic H5N1 flu strains. Thus, there is an urgent need to develop second-generation M2 inhibitors targeting the S31N mutant. However, the S31N mutant presents a huge challenge to drug discovery, and it has been considered undruggable for several decades. Using structural information, classical medicinal chemistry approaches, and M2-specific biological testing, we discovered benzyl-substituted amantadine derivatives with activity against both S31N and WT, among which 4-(adamantan-1-ylaminomethyl)-benzene-1,3-diol (44) is the most potent dual inhibitor. These inhibitors demonstrate that S31N is a druggable target and provide a new starting point to design novel M2 inhibitors that address the problem of drug-resistant influenza A infections.

Transmembrane orientation and possible role of the fusogenic peptide from parainfluenza virus 5 (PIV5) in promoting fusion.

Membrane fusion is required for diverse biological functions ranging from viral infection to neurotransmitter release. Fusogenic proteins increase the intrinsically slow rate of fusion by coupling energetically downhill conformational changes of the protein to kinetically unfavorable fusion of the membrane–phospholipid bilayers. Class I viral fusogenic proteins have an N-terminal hydrophobic fusion peptide (FP) domain, important for interaction with the target membrane, plus a C-terminal transmembrane (C-term-TM) helical membrane anchor. The role of the water-soluble regions of fusogenic proteins has been extensively studied, but the contributions of the membrane-interacting FP and C-term-TM peptides are less well characterized. Typically, FPs are thought to bind to membranes at an angle that allows helix penetration but not traversal of the lipid bilayer. Here, we show that the FP from the paramyxovirus parainfluenza virus 5 fusogenic protein, F, forms an N-terminal TM helix, which self-associates into a hexameric bundle. This FP also interacts strongly with the C-term-TM helix. Thus, the fusogenic F protein resembles SNARE proteins involved in vesicle fusion by having water-soluble coiled coils that zipper during fusion and TM helices in both membranes. By analogy to mechanosensitive channels, the force associated with zippering of the water-soluble coiled-coil domain is expected to lead to tilting of the FP helices, promoting interaction with the C-term-TM helices. The energetically unfavorable dehydration of lipid headgroups of opposing bilayers is compensated by thermodynamically favorable interactions between the FP and C-term-TM helices as the coiled coils zipper into the membrane phase, leading to a pore lined by both lipid and protein.

Flipping in the Pore: Discovery of Dual Inhibitors That Bind in Different Orientations to the Wild-Type versus the Amantadine-Resistant S31N Mutant of the Influenza A Virus M2 Proton Channel.

Influenza virus infections lead to numerous deaths and millions of hospitalizations each year. One challenge facing anti-influenza drug development is the heterogeneity of the circulating influenza viruses, which comprise several strains with variable susceptibility to antiviral drugs. For example, the wild-type (WT) influenza A viruses, such as the seasonal H1N1, tend to be sensitive to antiviral drugs, amantadine and rimantadine, while the S31N mutant viruses, such as the pandemic 2009 H1N1 (H1N1pdm09) and seasonal H3N2, are resistant to this class of drugs. Thus, drugs targeting both WT and the S31N mutant are highly desired. We report our design of a novel class of dual inhibitors along with their ion channel blockage and antiviral activities. The potency of the most active compound 11 in inhibiting WT and the S31N mutant influenza viruses is comparable with that of amantadine in inhibiting WT influenza virus. Solution NMR studies and molecular dynamics (MD) simulations of drug-M2 interactions supported our design hypothesis: namely, the dual inhibitor binds in the WT M2 channel with an aromatic group facing down toward the C-terminus, while the same drug binds in the S31N M2 channel with its aromatic group facing up toward the N-terminus. The flip-flop mode of drug binding correlates with the structure–activity relationship (SAR) and has paved the way for the next round of rational design of broad-spectrum antiviral drugs.

Supramolecular polymerization of benzene-1,3,5-tricarboxamide: a molecular dynamics simulation study.

Supramolecular polymerization in the family of benzene-1,3,5-tricarboxamide (BTA) has been investigated using atomistic molecular dynamics (MD) simulations. Gas phase calculations using a nonpolarizable force field reproduce the cooperativity in binding energy and intermolecular structure seen in quantum chemical calculations. Both quantum chemical and force field based calculations suggest that the ground state structure of the BTA dimer contains two donor hydrogen bonds and one acceptor hydrogen bond rather than the conjectured three-donor and zero-acceptor hydrogen-bonded state. MD simulations of BTA molecules in a realistic solvent, n-nonane, demonstrate the self-assembly process. The free energy (FE) of dimerization and of solvation has been determined. The solvated dimer of BTA with hexyl tails is more stable than two monomers by about 13 kcal/mol. Furthermore, the FE of association of a BTA molecule to an oligomer exhibits a dependence on the oligomer size, which is a robust signature of cooperative self-assembly.

Pore waters regulate ion permeation in a calcium release-activated calcium channel

Significance Calcium channel dysfunction is implicated in cardiac arrhythmias and immunodeficiency disorders. The recent publication of the crystal structure of a calcium ion channel denoted CRAC has allowed us to use atomic-level computer simulation to begin to understand the mechanism of ion permeation through this ubiquitous and vital cell entry pathway. Our computations identify pore water molecules as a crucial contributor in governing the channel conductance characteristics not only of the wild-type protein but also a disease-related mutant. The present study should be helpful for research designed to improve treatment of immune diseases caused by disrupted calcium signaling. The recent crystal structure of Orai, the pore unit of a calcium release-activated calcium (CRAC) channel, is used as the starting point for molecular dynamics and free-energy calculations designed to probe this channel’s conduction properties. In free molecular dynamics simulations, cations localize preferentially at the extracellular channel entrance near the ring of Glu residues identified in the crystal structure, whereas anions localize in the basic intracellular half of the pore. To begin to understand ion permeation, the potential of mean force (PMF) was calculated for displacing a single Na+ ion along the pore of the CRAC channel. The computed PMF indicates that the central hydrophobic region provides the major hindrance for ion diffusion along the permeation pathway, thereby illustrating the nonconducting nature of the crystal structure conformation. Strikingly, further PMF calculations demonstrate that the mutation V174A decreases the free energy barrier for conduction, rendering the channel effectively open. This seemingly dramatic effect of mutating a nonpolar residue for a smaller nonpolar residue in the pore hydrophobic region suggests an important role for the latter in conduction. Indeed, our computations show that even without significant channel-gating motions, a subtle change in the number of pore waters is sufficient to reshape the local electrostatic field and modulate the energetics of conduction, a result that rationalizes recent experimental findings. The present work suggests the activation mechanism for the wild-type CRAC channel is likely regulated by the number of pore waters and hence pore hydration governs the conductance.