The evolutionary background and functional consequences of the rs2071307 polymorphism in human tropoelastin.

Elastin is a major polymeric protein of the extracellular matrix, providing critical properties of extensibility and elastic recoil. The rs2071307 genomic polymorphism, resulting in the substitution of a serine for a glycine residue in a VPG motif in tropoelastin, has an unusually high minor allele frequency in humans. A consequence of such allelic heterozygosity would be the presence of a heterogeneous elastin polymer in up to 50% of the population, a situation which appears to be unique to Homo sapiens. VPG motifs are extremely common in hydrophobic domains of tropoelastins and are the sites of transient ß-turns that are essential for maintaining the conformational flexibility required for its function as an entropic elastomer. Earlier data demonstrated that single amino acid substitutions in tropoelastin can have functional consequences for polymeric elastin, particularly when present in mixed polymers. Here, using NMR and molecular dynamics approaches, we show the rs2071307 polymorphism reduces local propensity for ß-turn formation, with a consequent increase in polypeptide hydration and an expansion of the conformational ensemble manifested as an increased hydrodynamic radius, radius of gyration and asphericity. Furthermore, this substitution affects functional properties of polymeric elastin, particularly in heterogeneous polymers mimicking allelic heterozygosity. We discuss whether such effects, together with the unusually high minor allele frequency of the polymorphism, could imply some some evolutionary advantage for the heterozygous state.

A sulfur-aromatic gate latch is essential for opening of the Orai1 channel pore.

Sulfur-aromatic interactions occur in the majority of protein structures, yet little is known about their functional roles in ion channels. Here, we describe a novel molecular motif, the M101 gate latch, which is essential for gating of human Orai1 channels via its sulfur-aromatic interactions with the F99 hydrophobic gate. Molecular dynamics simulations of different Orai variants reveal that the gate latch is mostly engaged in open but not closed channels. In experimental studies, we use metal-ion bridges to show that promoting an M101-F99 bond directly activates Orai1, whereas disrupting this interaction triggers channel closure. Mutational analysis demonstrates that the methionine residue at this position has a unique combination of length, flexibility, and chemistry to act as an effective latch for the phenylalanine gate. Because sulfur-aromatic interactions provide additional stabilization compared to purely hydrophobic interactions, we infer that the six M101-F99 pairs in the hexameric channel provide a substantial energetic contribution to Orai1 activation.

Molecular docking with Gaussian Boson Sampling.

Gaussian Boson Samplers are photonic quantum devices with the potential to perform intractable tasks for classical systems. As with other near-term quantum technologies, an outstanding challenge is to identify specific problems of practical interest where these devices can prove useful. Here, we show that Gaussian Boson Samplers can be used to predict molecular docking configurations, a central problem for pharmaceutical drug design. We develop an approach where the problem is reduced to finding the maximum weighted clique in a graph, and show that Gaussian Boson Samplers can be programmed to sample large-weight cliques, i.e., stable docking configurations, with high probability, even with photon losses. We also describe how outputs from the device can be used to enhance the performance of classical algorithms. To benchmark our approach, we predict the binding mode of a ligand to the tumor necrosis factor-α converting enzyme, a target linked to immune system diseases and cancer.

The basic residues in the Orai1 channel inner pore promote opening of the outer hydrophobic gate.

Store-operated Orai1 channels regulate a wide range of cellular functions from gene expression to cell proliferation. Previous studies have shown that gating of Orai1 channels is regulated by the outer pore residues V102 and F99, which together function as a hydrophobic gate to block ion conduction in resting channels. Opening of this gate occurs through a conformational change that moves F99 away from the permeation pathway, leading to pore hydration and ion conduction. In addition to this outer hydrophobic gate, several studies have postulated the presence of an inner gate formed by the basic residues R91, K87, and R83 in the inner pore. These positively charged residues were suggested to block ion conduction in closed channels via mechanisms involving either electrostatic repulsion or steric occlusion by a bound anion plug. However, in contrast to this model, here we find that neutralization of the basic residues dose-dependently abolishes both STIM1-mediated and STIM1-independent activation of Orai1 channels. Molecular dynamics simulations show that loss of the basic residues dehydrates the pore around the hydrophobic gate and stabilizes the pore in a closed configuration. Likewise, the severe combined immunodeficiency mutation, Orai1 R91W, closes the channel by dewetting the hydrophobic stretch of the pore and stabilizing F99 in a pore-facing configuration. Loss of STIM1-gating in R91W and in the other basic residue mutants is rescued by a V102A mutation, which restores pore hydration at the hydrophobic gate to repermit ion conduction. These results indicate that the inner pore basic residues facilitate opening of the principal outer hydrophobic gate through a long-range effect involving hydration of the outer pore.

Substrate-Based Allosteric Regulation of a Homodimeric Enzyme.

Many enzymes operate through half-of-the sites reactivity wherein a single protomer is catalytically engaged at one time. In the case of the homodimeric enzyme, fluoroacetate dehalogenase, substrate binding triggers closing of a regulatory cap domain in the empty protomer, preventing substrate access to the remaining active site. However, the empty protomer serves a critical role by acquiring more disorder upon substrate binding, thereby entropically favoring the forward reaction. Empty protomer dynamics are also allosterically coupled to the bound protomer, driving conformational exchange at the active site and progress along the reaction coordinate. Here, we show that at high concentrations, a second substrate binds along the substrate-access channel of the occupied protomer, thereby dampening interprotomer dynamics and inhibiting catalysis. While a mutation (K152I) abrogates second site binding and removes inhibitory effects, it also precipitously lowers the maximum catalytic rate, implying a role for the allosteric pocket at low substrate concentrations, where only a single substrate engages the enzyme at one time. We show that this outer pocket first desolvates the substrate, whereupon it is deposited in the active site. Substrate binding to the active site then triggers the empty outer pocket to serve as an interprotomer allosteric conduit, enabling enhanced dynamics and sampling of activation states needed for catalysis. These allosteric networks and the ensuing changes resulting from second substrate binding are delineated using rigidity-based allosteric transmission theory and validated by nuclear magnetic resonance and functional studies. The results illustrate the role of dynamics along allosteric networks in facilitating function.

Mapping the functional anatomy of Orai1 transmembrane domains for CRAC channel gating.

Store-operated Orai1 channels are activated through a unique inside-out mechanism involving binding of the endoplasmic reticulum Ca2+ sensor STIM1 to cytoplasmic sites on Orai1. Although atomic-level details of Orai structure, including the pore and putative ligand binding domains, are resolved, how the gating signal is communicated to the pore and opens the gate is unknown. To address this issue, we used scanning mutagenesis to identify 15 residues in transmembrane domains (TMs) 1-4 whose perturbation activates Orai1 channels independently of STIM1. Cysteine accessibility analysis and molecular-dynamics simulations indicated that constitutive activation of the most robust variant, H134S, arises from a pore conformational change that opens a hydrophobic gate to augment pore hydration, similar to gating evoked by STIM1. Mutational analysis of this locus suggests that H134 acts as steric brake to stabilize the closed state of the channel. In addition, atomic packing analysis revealed distinct functional contacts between the TM1 pore helix and the surrounding TM2/3 helices, including one set mediated by a cluster of interdigitating hydrophobic residues and another by alternative ridges of polar and hydrophobic residues. Perturbing these contacts via mutagenesis destabilizes STIM1-mediated Orai1 channel gating, indicating that these bridges between TM1 and the surrounding TM2/3 ring are critical for conveying the gating signal to the pore. These findings help develop a framework for understanding the global conformational changes and allosteric interactions between topologically distinct domains that are essential for activation of Orai1 channels.

Structural basis for gating pore current in periodic paralysis.

Potassium-sensitive hypokalaemic and normokalaemic periodic paralysis are inherited skeletal muscle diseases characterized by episodes of flaccid muscle weakness1,2. They are caused by single mutations in positively charged residues ('gating charges') in the S4 transmembrane segment of the voltage sensor of the voltage-gated sodium channel Nav1.4 or the calcium channel Cav1.11,2. Mutations of the outermost gating charges (R1 and R2) cause hypokalaemic periodic paralysis1,2 by creating a pathogenic gating pore in the voltage sensor through which cations leak in the resting state3,4. Mutations of the third gating charge (R3) cause normokalaemic periodic paralysis 5 owing to cation leak in both activated and inactivated states 6 . Here we present high-resolution structures of the model bacterial sodium channel NavAb with the analogous gating-charge mutations7,8, which have similar functional effects as in the human channels. The R2G and R3G mutations have no effect on the backbone structures of the voltage sensor, but they create an aqueous cavity near the hydrophobic constriction site that controls gating charge movement through the voltage sensor. The R3G mutation extends the extracellular aqueous cleft through the entire length of the activated voltage sensor, creating an aqueous path through the membrane. Conversely, molecular modelling shows that the R2G mutation creates a continuous aqueous path through the membrane only in the resting state. Crystal structures of NavAb(R2G) in complex with guanidinium define a potential drug target site. Molecular dynamics simulations illustrate the mechanism of Na+ permeation through the mutant gating pore in concert with conformational fluctuations of the gating charge R4. Our results reveal pathogenic mechanisms of periodic paralysis at the atomic level and suggest designs of drugs that may prevent ionic leak and provide symptomatic relief from hypokalaemic and normokalaemic periodic paralysis.

Structures of closed and open states of a voltage-gated sodium channel.

Bacterial voltage-gated sodium channels (BacNavs) serve as models of their vertebrate counterparts. BacNavs contain conserved voltage-sensing and pore-forming domains, but they are homotetramers of four identical subunits, rather than pseudotetramers of four homologous domains. Here, we present structures of two NaVAb mutants that capture tightly closed and open states at a resolution of 2.8-3.2 Å. Introduction of two humanizing mutations in the S6 segment (NaVAb/FY: T206F and V213Y) generates a persistently closed form of the activation gate in which the intracellular ends of the four S6 segments are drawn tightly together to block ion permeation completely. This construct also revealed the complete structure of the four-helix bundle that forms the C-terminal domain. In contrast, truncation of the C-terminal 40 residues in NavAb/1-226 captures the activation gate in an open conformation, revealing the open state of a BacNav with intact voltage sensors. Comparing these structures illustrates the full range of motion of the activation gate, from closed with its orifice fully occluded to open with an orifice of ∼10 Å. Molecular dynamics and free-energy simulations confirm designation of NaVAb/1-226 as an open state that allows permeation of hydrated Na+, and these results also support a hydrophobic gating mechanism for control of ion permeation. These two structures allow completion of a closed-open-inactivated conformational cycle in a single voltage-gated sodium channel and give insight into the structural basis for state-dependent binding of sodium channel-blocking drugs.

STIM1 activates CRAC channels through rotation of the pore helix to open a hydrophobic gate.

Store-operated Ca2+ release-activated Ca2+ (CRAC) channels constitute a major pathway for Ca2+ influx and mediate many essential signalling functions in animal cells, yet how they open remains elusive. Here, we investigate the gating mechanism of the human CRAC channel Orai1 by its activator, stromal interacting molecule 1 (STIM1). We find that two rings of pore-lining residues, V102 and F99, work together to form a hydrophobic gate. Mutations of these residues to polar amino acids produce channels with leaky gates that conduct ions in the resting state. STIM1-mediated channel activation occurs through rotation of the pore helix, which displaces the F99 residues away from the pore axis to increase pore hydration, allowing ions to flow through the V102-F99 hydrophobic band. Pore helix rotation by STIM1 also explains the dynamic coupling between CRAC channel gating and ion selectivity. This hydrophobic gating mechanism has implications for CRAC channel function, pharmacology and disease-causing mutations.

The role of dimer asymmetry and protomer dynamics in enzyme catalysis.

Freeze-trapping x-ray crystallography, nuclear magnetic resonance, and computational techniques reveal the distribution of states and their interconversion rates along the reaction pathway of a bacterial homodimeric enzyme, fluoroacetate dehalogenase (FAcD). The crystal structure of apo-FAcD exhibits asymmetry around the dimer interface and cap domain, priming one protomer for substrate binding. This asymmetry is dynamically averaged through conformational exchange on a millisecond time scale. During catalysis, the protomer conformational exchange rate becomes enhanced, the empty protomer exhibits increased local disorder, and water egresses. Computational studies identify allosteric pathways between protomers. Water release and enhanced dynamics associated with catalysis compensate for entropic losses from substrate binding while facilitating sampling of the transition state. The studies provide insights into how substrate-coupled allosteric modulation of structure and dynamics facilitates catalysis in a homodimeric enzyme.

Structure and Dynamics of Extracellular Loops in Human Aquaporin-1 from Solid-State NMR and Molecular Dynamics.

Multiple moderate-resolution crystal structures of human aquaporin-1 have provided a foundation for understanding the molecular mechanism of selective water translocation in human cells. To gain insight into the interfacial structure and dynamics of human aquaporin-1 in a lipid environment, we performed nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics simulations. Using magic angle spinning solid-state NMR, we report a near complete resonance assignment of the human aquaporin-1. Chemical shift analysis of the secondary structure identified pronounced deviations from crystallographic structures in extracellular loops A and C, including the cis Y37-P38 bond in loop A, as well as ordering and immobilization of loop C. Site-specific H/D exchange measurements identify a number of protected nitrogen-bearing side chains and backbone amide groups, involved in stabilizing the loops. A combination of molecular dynamics simulations with NMR-derived restraints and filtering based on solvent accessibility allowed for the determination of a structural model of extracellular loops largely consistent with NMR results. The simulations reveal loop stabilizing interactions that alter the extracellular surface of human AQP1, with possible implications for water transport regulation through the channel. Modulation of water permeation may occur as a result of rearrangement of side chains from loop C in the extracellular vestibule of hAQP1, affecting the aromatic arginine selectivity filter.

Mechanism of Amyloidogenesis of a Bacterial AAA+ Chaperone.

Amyloids are fibrillar protein superstructures that are commonly associated with diseases in humans and with physiological functions in various organisms. The precise mechanisms of amyloid formation remain to be elucidated. Surprisingly, we discovered that a bacterial Escherichia coli chaperone-like ATPase, regulatory ATPase variant A (RavA), and specifically the LARA domain in RavA, forms amyloids under acidic conditions at elevated temperatures. RavA is involved in modulating the proper assembly of membrane respiratory complexes. LARA contains an N-terminal loop region followed by a ß-sandwich-like folded core. Several approaches, including nuclear magnetic resonance spectroscopy and molecular dynamics simulations, were used to determine the mechanism by which LARA switches to an amyloid state. These studies revealed that the folded core of LARA is amyloidogenic and is protected by its N-terminal loop. At low pH and high temperatures, the interaction of the N-terminal loop with the folded core is disrupted, leading to amyloid formation.

The molecular mechanism of Zinc acquisition by the neisserial outer-membrane transporter ZnuD.

Invading bacteria from the Neisseriaceae, Acinetobacteriaceae, Bordetellaceae and Moraxellaceae families express the conserved outer-membrane zinc transporter zinc-uptake component D (ZnuD) to overcome nutritional restriction imposed by the host organism during infection. Here we demonstrate that ZnuD is required for efficient systemic infections by the causative agent of bacterial meningitis, Neisseria meningitidis, in a mouse model. We also combine X-ray crystallography and molecular dynamics simulations to gain insight into the mechanism of zinc recognition and transport across the bacterial outer-membrane by ZnuD. Because ZnuD is also considered a promising vaccine candidate against N. meningitidis, we use several ZnuD structural intermediates to map potential antigenic epitopes, and propose a mechanism by which ZnuD can maintain high sequence conservation yet avoid immune recognition by altering the conformation of surface-exposed loops.

Modification and periplasmic translocation of the biofilm exopolysaccharide poly-ß-1,6-N-acetyl-D-glucosamine.

Poly-ß-1,6-N-acetyl-D-glucosamine (PNAG) is an exopolysaccharide produced by a wide variety of medically important bacteria. Polyglucosamine subunit B (PgaB) is responsible for the de-N-acetylation of PNAG, a process required for polymer export and biofilm formation. PgaB is located in the periplasm and likely bridges the inner membrane synthesis and outer membrane export machinery. Here, we present structural, functional, and molecular simulation data that suggest PgaB associates with PNAG continuously during periplasmic transport. We show that the association of PgaB's N- and C-terminal domains forms a cleft required for the binding and de-N-acetylation of PNAG. Molecular dynamics (MD) simulations of PgaB show a binding preference for N-acetylglucosamine (GlcNAc) to the N-terminal domain and glucosammonium to the C-terminal domain. Continuous ligand binding density is observed that extends around PgaB from the N-terminal domain active site to an electronegative groove on the C-terminal domain that would allow for a processive mechanism. PgaB's C-terminal domain (PgaB310-672) directly binds PNAG oligomers with dissociation constants of ∼1-3 mM, and the structures of PgaB310-672 in complex with ß-1,6-(GlcNAc)6, GlcNAc, and glucosamine reveal a unique binding mode suitable for interaction with de-N-acetylated PNAG (dPNAG). Furthermore, PgaB310-672 contains a ß-hairpin loop (ßHL) important for binding PNAG that was disordered in previous PgaB42-655 structures and is highly dynamic in the MD simulations. We propose that conformational changes in PgaB310-672 mediated by the ßHL on binding of PNAG/dPNAG play an important role in the targeting of the polymer for export and its release.

Inclusion of trial functions in the Langevin equation path integral ground state method: application to parahydrogen clusters and their isotopologues.

We developed and studied the implementation of trial wavefunctions in the newly proposed Langevin equation Path Integral Ground State (LePIGS) method [S. Constable, M. Schmidt, C. Ing, T. Zeng, and P.-N. Roy, J. Phys. Chem. A 117, 7461 (2013)]. The LePIGS method is based on the Path Integral Ground State (PIGS) formalism combined with Path Integral Molecular Dynamics sampling using a Langevin equation based sampling of the canonical distribution. This LePIGS method originally incorporated a trivial trial wavefunction, ψT, equal to unity. The present paper assesses the effectiveness of three different trial wavefunctions on three isotopes of hydrogen for cluster sizes N = 4, 8, and 13. The trial wavefunctions of interest are the unity trial wavefunction used in the original LePIGS work, a Jastrow trial wavefunction that includes correlations due to hard-core repulsions, and a normal mode trial wavefunction that includes information on the equilibrium geometry. Based on this analysis, we opt for the Jastrow wavefunction to calculate energetic and structural properties for parahydrogen, orthodeuterium, and paratritium clusters of size N = 4 - 19, 33. Energetic and structural properties are obtained and compared to earlier work based on Monte Carlo PIGS simulations to study the accuracy of the proposed approach. The new results for paratritium clusters will serve as benchmark for future studies. This paper provides a detailed, yet general method for optimizing the necessary parameters required for the study of the ground state of a large variety of systems.

Catalysis of Na+ permeation in the bacterial sodium channel Na(V)Ab.

Determination of a high-resolution 3D structure of voltage-gated sodium channel Na(V)Ab opens the way to elucidating the mechanism of ion conductance and selectivity. To examine permeation of Na(+) through the selectivity filter of the channel, we performed large-scale molecular dynamics simulations of Na(V)Ab in an explicit, hydrated lipid bilayer at 0 mV in 150 mM NaCl, for a total simulation time of 21.6 µs. Although the cytoplasmic end of the pore is closed, reversible influx and efflux of Na(+) through the selectivity filter occurred spontaneously during simulations, leading to equilibrium movement of Na(+) between the extracellular medium and the central cavity of the channel. Analysis of Na(+) dynamics reveals a knock-on mechanism of ion permeation characterized by alternating occupancy of the channel by 2 and 3 Na(+) ions, with a computed rate of translocation of (6 ± 1) × 10(6) ions⋅s(-1) that is consistent with expectations from electrophysiological studies. The binding of Na(+) is intimately coupled to conformational isomerization of the four E177 side chains lining the extracellular end of the selectivity filter. The reciprocal coordination of variable numbers of Na(+) ions and carboxylate groups leads to their condensation into ionic clusters of variable charge and spatial arrangement. Structural fluctuations of these ionic clusters result in a myriad of ion binding modes and foster a highly degenerate, liquid-like energy landscape propitious to Na(+) diffusion. By stabilizing multiple ionic occupancy states while helping Na(+) ions diffuse within the selectivity filter, the conformational flexibility of E177 side chains underpins the knock-on mechanism of Na(+) permeation.

Langevin equation path integral ground state.

We propose a Langevin equation path integral ground state (LePIGS) approach for the calculation of ground state (zero temperature) properties of molecular systems. The approach is based on a modification of the finite temperature path integral Langevin equation (PILE) method (J. Chem. Phys. 2010, 133, 124104) to the case of open Feynman paths. Such open paths are necessary for a ground state formulation. We illustrate the applicability of the method using model systems and the weakly bound water-parahydrogen dimer. We show that the method can lead to converged zero point energies and structural properties.

A path-integral Langevin equation treatment of low-temperature doped helium clusters.

We present an implementation of path integral molecular dynamics for sampling low temperature properties of doped helium clusters using Langevin dynamics. The robustness of the path integral Langevin equation and white-noise Langevin equation [M. Ceriotti, M. Parrinello, T. E. Markland, and D. E. Manolopoulos, J. Chem. Phys. 133, 124104 (2010)] sampling methods are considered for those weakly bound systems with comparison to path integral Monte Carlo (PIMC) in terms of efficiency and accuracy. Using these techniques, convergence studies are performed to confirm the systematic error reduction introduced by increasing the number of discretization steps of the path integral. We comment on the structural and energetic evolution of He(N)-CO(2) clusters from N = 1 to 20. To quantify the importance of both rotations and exchange in our simulations, we present a chemical potential and calculated band origin shifts as a function of cluster size utilizing PIMC sampling that includes these effects. This work also serves to showcase the implementation of path integral simulation techniques within the molecular modelling toolkit [K. Hinsen, J. Comp. Chem. 21, 79 (2000)], an open-source molecular simulation package.

Adhesion of cylindrical colloids to the surface of a membrane.

We study the system of cylindrical colloids adhering to an originally flat fluid membrane, on the basis of a full treatment of the Helfrich model. Our approach allows for numerical calculation of the free energy in both shallow- and deep-wrapping conformations. We show that the free energy of two cylinders adhering to the same side of a membrane has two branches corresponding to shallow and deep wrapping and that the system of two cylinders adhering to opposite sides of a membrane can undergo a first-order phase transition between two membrane-mediated attractive states.