Search

To become specialized binders, antibodies undergo a process called affinity maturation to maximize their binding affinity. Despite this process, some antibodies retain low-affinity binding to diverse epitopes in a phenomenon called polyreactivity. Here we seek to understand the molecular basis of this polyreactivity in antibodies. Our results highlight that polyreactive antigen-binding fragments (Fabs) bind their targets with low affinities, comparable to T cell receptor recognition of autologous classical major histocompatibility complex. Extensive mutagenic studies find no singular amino acid residue or biochemical property responsible for polyreactive interaction, suggesting that polyreactive antibodies use multiple strategies for engagement. Finally, our crystal structures and all-atom molecular dynamics simulations of polyreactive Fabs show increased rigidity compared to their monoreactive relatives, forming a neutral and accessible platform for diverse antigens to bind. Together, these data support a cooperative strategy of rigid neutrality in establishing the polyreactive status of an antibody molecule., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 The Author(s). Published by Elsevier Inc. All rights reserved.)

In plants, host response to pathogenic microbes is driven both by microbial perception and detection of modified-self. The Xanthomonas campestris effector protein AvrAC/XopAC uridylylates the Arabidopsis BIK1 kinase to dampen basal resistance and thereby promotes bacterial virulence. Here we show that PBL2, a paralog of BIK1, is similarly uridylylated by AvrAC. However, in contrast to BIK1, PBL2 uridylylation is specifically required for host recognition of AvrAC to trigger immunity, but not AvrAC virulence. PBL2 thus acts as a decoy and enables AvrAC detection. AvrAC recognition also requires the RKS1 pseudokinase of the ZRK family and the NOD-like receptor ZAR1, which is known to recognize the Pseudomonas syringae effector HopZ1a. ZAR1 forms a stable complex with RKS1, which specifically recruits PBL2 when the latter is uridylylated by AvrAC, triggering ZAR1-mediated immunity. The results illustrate how decoy substrates and pseudokinases can specify and expand the capacity of the plant immune system., (Copyright © 2015 Elsevier Inc. All rights reserved.)

The fundamental biophysics underlying the selective movement of ions through ion channels was launched by George Eisenman in the 1960s, using glass electrodes. This minireview examines the insights from these early studies and the explosive progress made since then., (Copyright © 2014 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Actin, a highly conserved cytoskeletal protein found in all eukaryotic cells, facilitates cell motility and membrane remodeling via a directional polymerization cycle referred to as treadmilling. The nucleotide bound at the core of each actin subunit regulates this process. Although the biochemical kinetics of treadmilling has been well characterized, the atomistic details of how the nucleotide affects polymerization remain to be definitively determined. There is increasing evidence that the nucleotide regulation (and other characteristics) of actin cannot be fully described from the minimum energy structure, but rather depends on a dynamic equilibrium between conformations. In this work we explore the conformational mobility of the actin monomer (G-actin) in a coarse-grained subspace using umbrella sampling to bias all-atom molecular-dynamics simulations along the variables of interest. The results reveal that ADP-bound actin subunits are more conformationally mobile than ATP-bound subunits. We used a multiscale analysis method involving coarse-grained and atomistic representations of these simulations to characterize how the nucleotide affects the low-energy states of these systems. The interface between subdomains SD2-SD4, which is important for polymerization, is stabilized in an actin filament-like (F-actin) conformation in ATP-bound G-actin. Additionally, the nucleotide modulates the conformation of the SD1-SD3 interface, a region involved in the binding of several actin-binding proteins., (Copyright © 2014 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

The 2013 Nobel Prize in Chemistry has been awarded to Martin Karplus, Michael Levitt, and Arieh Warshel for their work on developing computational methods to study complex chemical systems. Their work has led to mechanistic critical insights into chemical systems both large and small and has enabled progress in a number of different fields, including structural biology., (Copyright © 2013 Elsevier Ltd. All rights reserved.)

Ionotropic glutamate receptors (iGluRs) transduce the chemical signal of neurotransmitter release into membrane depolarization at excitatory synapses in the brain. The opening of the transmembrane ion channel of these ligand-gated receptors is driven by conformational transitions that are induced by the association of glutamate molecules to the ligand-binding domains (LBDs). Here, we describe the crystal structure of a GluA2 LBD tetramer in a configuration that involves an ∼30° rotation of the LBD dimers relative to the crystal structure of the full-length receptor. The configuration is stabilized by an engineered disulfide crosslink. Biochemical and electrophysiological studies on full-length receptors incorporating either this crosslink or an engineered metal bridge show that this LBD configuration corresponds to an intermediate state of receptor activation. GluA2 activation therefore involves a combination of both intra-LBD (cleft closure) and inter-LBD dimer conformational transitions. Overall, these results provide a comprehensive structural characterization of an iGluR intermediate state., (Copyright © 2013 Elsevier Inc. All rights reserved.)

The structure of OmpF porin in complex with three common antibiotics (zwitterionic ampicillin, anionic ertapenem, and di-anionic carbenicillin) was determined using X-ray crystallography. The three antibiotics are found to bind within the extracellular and periplasmic pore vestibules, away from the narrow OmpF constriction zone. Using the X-ray structures as a starting point, nonequilibrium molecular dynamics simulations with an applied membrane voltage show that ionic current through the OmpF channel is blocked with bound ampicillin, but not with bound carbenicillin. The susceptibility of Escherichia coli expressing OmpF mutants to ampicillin and carbenicillin was also experimentally characterized using microbiologic assays. These results show that general diffusion by OmpF porins allows for transfer of molecules with varied charged states and give insights into the design of more efficient antibiotics. A better understanding of this mechanism will shed light on nature's way of devising channels able to enhance the transport of molecules through membranes., (Copyright © 2013 Elsevier Ltd. All rights reserved.)

Proteins may undergo multiple conformational changes required for their function. One strategy used to estimate target-site positions in unknown structural conformations involves single-pair resonance energy transfer (RET) distance measurements. However, interpretation of inter-residue distances is difficult when applied to three-dimensional structural rearrangements, especially in homomeric systems. We developed a positioning method using inverse trilateration/triangulation to map target sites within a homomeric protein in all defined states, with simultaneous functional recordings. The procedure accounts for probe diffusion to accurately determine the three-dimensional position and confidence region of lanthanide LRET donors attached to a target site (one per subunit), relative to a single fluorescent acceptor placed in a static site. As first application, the method is used to determine the position of a functional voltage-gated potassium channel's voltage sensor. Our results verify the crystal structure relaxed conformation and report on the resting and active conformations for which crystal structures are not available., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

In K+ channels, rearrangements of the pore outer vestibule have been associated with C-type inactivation gating. Paradoxically, the crystal structure of Open/C-type inactivated KcsA suggests these movements to be modest in magnitude. In this study, we show that under physiological conditions, the KcsA outer vestibule undergoes relatively large dynamic rearrangements upon inactivation. External Cd2+ enhances the rate of C-type inactivation in an cysteine mutant (Y82C) via metal-bridge formation. This effect is not present in a non-inactivating mutant (E71A/Y82C). Tandem dimer and tandem tetramer constructs of equivalent cysteine mutants in KcsA and Shaker K+ channels demonstrate that these Cd2+ metal bridges are formed only between adjacent subunits. This is well supported by molecular dynamics simulations. Based on the crystal structure of Cd2+ -bound Y82C-KcsA in the closed state, together with electron paramagnetic resonance distance measurements in the KcsA outer vestibule, we suggest that subunits must dynamically come in close proximity as the channels undergo inactivation., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

Inward rectifier potassium (Kir) channels act as cellular diodes, allowing unrestricted flow of potassium (K(+)) into the cell while preventing currents of large magnitude in the outward direction. The rectification mechanism by which this occurs involves a coupling between K(+) and intracellular blockers-magnesium (Mg(2+)) or polyamines-that simultaneously occupy the permeation pathway. In addition to the transmembrane pore, Kirs possess a large cytoplasmic domain (CD) that provides a favorable electronegative environment for cations. Electrophysiological experiments have shown that the CD is a key regulator of both conductance and rectification. In this study, we calculate and compare averaged equilibrium probability densities of K(+) and Cl(-) in open-pore models of the CDs of a weak (Kir1.1-ROMK) and a strong (Kir2.1-IRK) rectifier through explicit-solvent molecular-dynamics simulations in ~1 M KCl. The CD of both channels concentrates K(+) ions greater than threefold inside the cytoplasmic pore while IRK shows an additional K(+) accumulation region near the cytoplasmic entrance. Simulations carried out with Mg(2+) or spermine (SPM(4+)) show that these ions interact with pore-lining residues, shielding the surface charge and reducing K(+) in both channels. The results also show that SPM(4+) behaves differently inside these two channels. Although SPM(4+) remains inside the CD of ROMK, it diffuses around the entire volume of the pore. In contrast, this polyatomic cation finds long-lived conformational states inside the IRK pore, interacting with residues E224, D259, and E299. The strong rectifier CD is also capable of sequestering an additional SPM(4+) at the cytoplasmic entrance near a cluster of negative residues D249, D274, E275, and D276. Although understanding the actual mechanism of rectification blockade will require high-resolution structural information of the blocked state, these simulations provide insight into how sequence variation in the CD can affect the multi-ion distributions that underlie the mechanisms of conduction, rectification affinity, and kinetics., (Copyright © 2012 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Loop-gating is one of two voltage-dependent mechanisms that regulate the open probability of connexin channels. The loop-gate permeability barrier is formed by a segment of the first extracellular loop (E1) (the parahelix) and appears to be accompanied by straightening of the bend angle between E1 and the first transmembrane domain (TM1). Here, all-atom molecular dynamics simulations are used to identify and characterize interacting van der Waals and electrostatic networks that stabilize the parahelices and TM1/E1 bend angles of the open Cx26 hemichannel. Dynamic fluctuations in an electrostatic network in each subunit are directly linked to the stability of parahelix structure and TM1/E1 bend angle in adjacent subunits. The electrostatic network includes charged residues that are pore-lining and thus positioned to be voltage sensors. We propose that the transition to the closed state is initiated by voltage-driven disruption of the networks that stabilize the open-state parahelix configuration, allowing the parahelix to protrude into the channel pore to form the loop-gate barrier. Straightening of the TM1/E1 bend appears to be a consequence of the reorganization of the interacting networks that accompany the conformational change of the parahelix. The electrostatic network extends across subunit boundaries, suggesting a concerted gating mechanism., (Copyright © 2012 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

The accurate prediction of membrane-insertion probability for arbitrary protein sequences is a critical challenge to identifying membrane proteins and determining their folded structures. Although algorithms based on sequence statistics have had moderate success, a complete understanding of the energetic factors that drive the insertion of membrane proteins is essential to thoroughly meeting this challenge. In the last few years, numerous attempts to define a free-energy scale for amino-acid insertion have been made, yet disagreement between most experimental and theoretical scales persists. However, for a recently resolved water-to-bilayer scale, it is found that molecular dynamics simulations that carefully mimic the conditions of the experiment can reproduce experimental free energies, even when using the same force field as previous computational studies that were cited as evidence of this disagreement. Therefore, it is suggested that experimental and simulation-based scales can both be accurate and that discrepancies stem from disparities in the microscopic processes being considered rather than methodological errors. Furthermore, these disparities make the development of a single universally applicable membrane-insertion free energy scale difficult., (Copyright © 2012 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Voltage sensor domains (VSD) are transmembrane proteins that respond to changes in membrane voltage and modulate the activity of ion channels, enzymes, or in the case of proton channels allow permeation of protons across the cell membrane. VSDs consist of four transmembrane segments, S1-S4, forming an antiparallel helical bundle. The S4 segment contains several positively charged residues, mainly arginines, located at every third position along the helix. In the voltage-gated Shaker K(+) channel, the mutation of the first arginine of S4 to a smaller uncharged amino acid allows permeation of cations through the VSD. These currents, known as ω-currents, pass through the VSD and are distinct from K(+) currents passing through the main ion conduction pore. Here we report molecular dynamics simulations of the ω-current in the resting-state conformation for Kv1.2 and for four of its mutants. The four tested mutants exhibit various degrees of conductivity for K(+) and Cl(-) ions, with a slight selectivity for K(+) over Cl(-). Analysis of the ion permeation pathway, in the case of a highly conductive mutant, reveals a negatively charged constriction region near the center of the membrane that might act as a selectivity filter to prevent permeation of anions through the pore. The residues R1 in S4 and E1 in S2 are located at the narrowest region of the ω-pore for the resting state conformation of the VSD, in agreement with experiments showing that the largest increase in current is produced by the double mutation E1D and R1S., (Copyright © 2012 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Voltage-sensing domains (VSDs) undergo conformational changes in response to the membrane potential and are the critical structural modules responsible for the activation of voltage-gated channels. Structural information about the key conformational states underlying voltage activation is currently incomplete. Through the use of experimentally determined residue-residue interactions as structural constraints, we determine and refine a model of the Kv channel VSD in the resting conformation. The resulting structural model is in broad agreement with results that originate from various labs using different techniques, indicating the emergence of a consensus for the structural basis of voltage sensing., (Copyright © 2011 Elsevier Inc. All rights reserved.)

Revealing the three-dimensional organization of large dynamic protein complexes in solution is challenging. To tackle this problem, Różycki and colleagues (2011) design a method combining small angle X-ray scattering (SAXS) data with the results of computer simulations. Their study offers new insights into the conformational transition induced by salt that occurs in an endosome-associated ESCRT-III CHMP3 domain., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

For ion channels, the transmembrane potential plays a critical role by acting as a driving force for permeant ions. At the microscopic level, the transmembrane potential is thought to decay nonlinearly across the ion permeation pathway because of the irregular three-dimensional shape of the channel's pore. By taking advantage of the current structural and functional understanding of cyclic nucleotide-gated channels, in this study we experimentally explore the transmembrane potential's distribution across the open pore. As a readout for the voltage drop, we engineered cysteine residues along the selectivity filter and scanned the sensitivity of their modification rates by Ag(+) to the transmembrane potential. The experimental data, which indicate that the majority of the electric field drops across the selectivity filter, are in good agreement with continuum electrostatic calculations using a homology model of an open CNG channel. By focusing the transmembrane potential across the selectivity filter, the electromotive driving force is coupled with the movement of permeant ions in the filter, maximizing the efficiency of this process., (Copyright © 2010 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

The transmembrane conformation of Thermotoga maritima CorA, a magnesium transport system, has been studied in its Mg(2+)-bound form by site-directed spin labeling and electron paramagnetic resonance spectroscopy. Probe mobility together with accessibility data were used to evaluate the overall dynamics and relative arrangement of individual transmembrane segments TM1 and TM2. TM1 extends toward the cytoplasmic side creating a water-filled cavity, while TM2 is located in the periphery of the oligomer, contacting the lipid bilayer. A structural model for the conserved extracellular loop was generated based on EPR data and MD simulations, in which residue E316 is located toward the five-fold symmetry axis in position to electrostatically influence divalent ion translocation. Electrostatic analysis of our model suggest that, in agreement with the crystal structure, Mg(2+) -bound CorA is in a closed conformation. The present results suggest that long-range structural rearrangements are necessary to allow Mg(2+) translocation., (Copyright 2010 Elsevier Ltd. All rights reserved.)

The ability to discriminate between different cations efficiently is essential for the proper physiological functioning of many membrane transport proteins. One obvious mechanism of ion selectivity is when a binding site is structurally constrained by the protein architecture and its geometry is precisely adapted to fit an ion of a given size. This mechanism is not effective in the case of flexible protein binding sites that are able to deform structurally or to adapt to a bound ion. In this study, the concept of nontrivial ion selectivity arising in a highly flexible protein binding site conceptually represented as a microdroplet of ligands confined to a small volume is explored. The environment imposed by the spatial confinement is a critical feature of the reduced models. A large number of reduced binding site models (1077) comprising typical ion-coordinating ligands (carbonyl, hydroxyl, carboxylate, water) are constructed and characterized for Na(+)/K(+) and Ca(2+)/Ba(2+) size selectivity using free energy perturbation molecular dynamics simulations. Free energies are highly correlated with the sum of ion-ligand and ligand-ligand mean interactions, but the relative balance of those two contributions is different for K(+)-selective and Na(+)-selective binding sites. The analysis indicates that both the number and the type of ligands are important factors in ion selectivity., ((c) 2010 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

The atomic models of the Kv1.2 potassium channel in the active and resting state, originally presented elsewhere, are here refined using molecular dynamics simulations in an explicit membrane-solvent environment. With a minor adjustment of the orientation of the first arginine along the S4 segment, the total gating charge of the channel determined from >0.5 mus of molecular dynamics simulation is approximately 12-12.7 e, in good accord with experimental estimates for the Shaker potassium channel, indicating that the final models offer a realistic depiction of voltage-gating. In the resting state of Kv1.2, the S4 segment in the voltage-sensing domain (VSD) spontaneously converts into a 3(10) helix over a stretch of 10 residues. The 3(10) helical conformation orients the gating arginines on S4 toward a water-filled crevice within the VSD and allows salt-bridge interactions with negatively charged residues along S2 and S3. Free energy calculations of the fractional transmembrane potential, acting upon key charged residues of the VSD, reveals that the applied field varies rapidly over a narrow region of 10-15 A corresponding to the outer leaflet of the bilayer. The focused field allows the transfer of a large gating charge without translocation of S4 across the membrane., (Copyright 2010 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

We present a coarse residue-based computational method to rapidly compute the solution scattering profile from a protein with dynamical fluctuations. The method is built upon a coarse-grained (CG) representation of the protein. This CG representation takes advantage of the intrinsic low-resolution and CG nature of solution scattering data. It allows rapid scattering determination from a large number of conformations that can be extracted from CG simulations to obtain scattering characterization of protein conformations. The method includes several important elements, effective residue structure factors derived from the Protein Data Bank, explicit treatment of water molecules in the hydration layer at the surface of the protein, and an ensemble average of scattering from a variety of appropriate conformations to account for macromolecular flexibility. This simplified method is calibrated and illustrated to accurately reproduce the experimental scattering curve of Hen egg white lysozyme. We then illustrated the applications of this CG method by computing the solution scattering patterns of several representative protein folds and multiple conformational states. The results suggest that solution scattering data, when combined with the reliable computational method that we developed, show great potential for a better structural description of multidomain complexes in different functional states, and for recognizing structural folds when sequence similarity to a protein of known structure is low.

Long timescale (>1 micros) molecular dynamics simulations of protein folding offer a powerful tool for understanding the atomic-scale interactions that determine a protein's folding pathway and stabilize its native state. Unfortunately, when the simulated protein fails to fold, it is often unclear whether the failure is due to a deficiency in the underlying force fields or simply a lack of sufficient simulation time. We examine one such case, the human Pin1 WW domain, using the recently developed deactivated morphing method to calculate free energy differences between misfolded and folded states. We find that the force field we used favors the misfolded states, explaining the failure of the folding simulations. Possible further applications of deactivated morphing and implications for force field development are discussed.

We present an approach for incorporating solvent accessibility data from electron paramagnetic resonance experiments in the structural refinement of membrane proteins through restrained molecular dynamics simulations. The restraints have been parameterized from oxygen (PiO(2)) and nickel-ethylenediaminediacetic acid (PiNiEdda) collision frequencies, as indicators of lipid or aqueous exposed spin-label sites. These are enforced through interactions between a pseudoatom representation of the covalently attached Nitroxide spin-label and virtual "solvent" particles corresponding to O(2) and NiEdda in the surrounding environment. Interactions were computed using an empirical potential function, where the parameters have been optimized to account for the different accessibilities of the spin-label pseudoatoms to the surrounding environment. This approach, "pseudoatom-driven solvent accessibility refinement", was validated by refolding distorted conformations of the Streptomyces lividans potassium channel (KcsA), corresponding to a range of 2-30 A root mean-square deviations away from the native structure. Molecular dynamics simulations based on up to 58 electron paramagnetic resonance restraints derived from spin-label mutants were able to converge toward the native structure within 1-3 A root mean-square deviations with minimal computational cost. The use of energy-based ranking and structure similarity clustering as selection criteria helped in the convergence and identification of correctly folded structures from a large number of simulations. This approach can be applied to a variety of integral membrane protein systems, regardless of oligomeric state, and should be particularly useful in calculating conformational changes from a known reference crystal structure.

A theoretical framework is elaborated to account for the effect of a transmembrane potential in computer simulations. It is shown that a simulation with a constant external electric field applied in the direction normal to the membrane is equivalent to the influence of surrounding infinite baths maintained to a voltage difference via ion-exchanging electrodes connected to an electromotive force. It is also shown that the linearly-weighted displacement charge within the simulation system tracks the net flow of charge through the external circuit comprising the electromotive force and the electrodes. Using a statistical mechanical reduction of the degrees of freedom of the external system, three distinct theoretical routes are formulated and examined for the purpose of characterizing the free energy of a protein embedded in a membrane that is submitted to a voltage difference. The W-route is constructed from the variations in the voltage-dependent potential of mean force along a reaction path connecting two conformations of the protein. The Q-route is based on the average displacement charge as a function of the conformation of the protein. Finally, the G-route considers the relative charging free energy of specific residues, with and without applied membrane potentials. The theoretical formulation is illustrated with a simple model of an ion crossing a vacuum slab surrounded by two aqueous bulk phases and with a fragment of the voltage-sensor of the KvAP potassium channel.

All-atom molecular dynamics simulations are used to better understand the dynamic environment experienced by the Kv1.2 channel in a lipid membrane. The structure of the channel is stable during the trajectories. The pore domain keeps a well-defined conformation, whereas the voltage-sensing domains undergo important lateral fluctuations, consistent with their modular nature. A channel-like region at the center of the S1-S4 helical bundle fills rapidly with water, reminiscent of the concept of high-dielectric aqueous crevices. The first two arginines along S4 (R294 and R297) adopt an interfacial position where they interact favorably with water and the lipid headgroups. The following two arginines (R300 and R303) interact predominantly with water and E226 in S2. Despite the absence of a structurally permanent gating pore formed by protein residues and surrounding the S4 helix, as traditionally pictured, the charged residues are located in a favorable environment and are not extensively exposed to the membrane nonpolar region. Continuum electrostatic computations indicate that the transmembrane potential sensed by the charged residues in the voltage sensor varies abruptly over the outer half of the membrane in the arginine-rich region of S4; thus, the voltage gradient or membrane electric field is "focused". Interactions of basic residues with the lipid headgroups at the intracellular membrane-solution interface reduce the membrane thickness near the channel, resulting in an increased transmembrane field.

Membrane depolarization causes voltage-gated ion channels to transition from a resting/closed conformation to an activated/open conformation. We used voltage-clamp fluorometry to measure protein motion at specific regions of the Shaker Kv channel. This enabled us to construct new structural models of the resting/closed and activated/open states based on the Kv1.2 crystal structure using the Rosetta-Membrane method and molecular dynamics simulations. Our models account for the measured gating charge displacement and suggest a molecular mechanism of activation in which the primary voltage sensors, S4s, rotate by approximately 180 degrees as they move "outward" by 6-8 A. A subsequent tilting motion of the S4s and the pore domain helices, S5s, of all four subunits induces a concerted movement of the channel's S4-S5 linkers and S6 helices, allowing ion conduction. Our models are compatible with a wide body of data and resolve apparent contradictions that previously led to several distinct models of voltage sensing.

Ionotropic glutamate receptors are ligand-gated transmembrane ion channels activated by the binding of glutamate. The free energy landscapes governing the opening/closing of the GluR2 S1S2 ligand-binding domain in the apo, DNQX-, and glutamate-bound forms are computed by using all-atom molecular dynamics simulations with explicit solvent, in conjunction with an umbrella sampling strategy. The apo S1S2 easily accesses low-energy conformations that are more open than observed in X-ray crystal structures. A free energy of 9-12 kcal/mol becomes available upon glutamate binding for driving conformational changes in S1S2 associated with receptor activation. Small-angle X-ray scattering profiles calculated from computed ensemble averages agree better with experimental results than profiles calculated from static X-ray crystal structures. Water molecules in the cleft may contribute to stabilizing the apo S1S2 in open conformations. Free energy landscapes were also computed for the glutamate-bound T686A and T686S S1S2 mutants, and the results elaborate on findings from experimental functional studies.

The gating mechanism of K(v) channels is not known. In this issue of Neuron, Soler-Llavina et al. present fascinating results that support the concept of relatively independent voltage-sensing modules. However, they also find that its interactions with the pore domain are rather complex, with specific S4-S5 intersubunit contacts underlying the concerted transition leading to the channel opening.

The absolute (standard) binding free energy of eight FK506-related ligands to FKBP12 is calculated using free energy perturbation molecular dynamics (FEP/MD) simulations with explicit solvent. A number of features are implemented to improve the accuracy and enhance the convergence of the calculations. First, the absolute binding free energy is decomposed into sequential steps during which the ligand-surrounding interactions as well as various biasing potentials restraining the translation, orientation, and conformation of the ligand are turned "on" and "off." Second, sampling of the ligand conformation is enforced by a restraining potential based on the root mean-square deviation relative to the bound state conformation. The effect of all the restraining potentials is rigorously unbiased, and it is shown explicitly that the final results are independent of all artificial restraints. Third, the repulsive and dispersive free energy contribution arising from the Lennard-Jones interactions of the ligand with its surrounding (protein and solvent) is calculated using the Weeks-Chandler-Andersen separation. This separation also improves convergence of the FEP/MD calculations. Fourth, to decrease the computational cost, only a small number of atoms in the vicinity of the binding site are simulated explicitly, while all the influence of the remaining atoms is incorporated implicitly using the generalized solvent boundary potential (GSBP) method. With GSBP, the size of the simulated FKBP12/ligand systems is significantly reduced, from approximately 25,000 to 2500. The computations are very efficient and the statistical error is small ( approximately 1 kcal/mol). The calculated binding free energies are generally in good agreement with available experimental data and previous calculations (within approximately 2 kcal/mol). The present results indicate that a strategy based on FEP/MD simulations of a reduced GSBP atomic model sampled with conformational, translational, and orientational restraining potentials can be computationally inexpensive and accurate.

We investigate methods for extracting the potential of mean force (PMF) governing ion permeation from molecular dynamics simulations (MD) using gramicidin A as a prototypical narrow ion channel. It is possible to obtain well-converged meaningful PMFs using all-atom MD, which predict experimental observables within order-of-magnitude agreement with experimental results. This was possible by careful attention to issues of statistical convergence of the PMF, finite size effects, and lipid hydrocarbon chain polarizability. When comparing the modern all-atom force fields of CHARMM27 and AMBER94, we found that a fairly consistent picture emerges, and that both AMBER94 and CHARMM27 predict observables that are in semiquantitative agreement with both the experimental conductance and dissociation coefficient. Even small changes in the force field, however, result in significant changes in permeation energetics. Furthermore, the full two-dimensional free-energy surface describing permeation reveals the location and magnitude of the central barrier and the location of two binding sites for K(+) ion permeation near the channel entrance--i.e., an inner site on-axis and an outer site off-axis. We conclude that the MD-PMF approach is a powerful tool for understanding and predicting the function of narrow ion channels in a manner that is consistent with the atomic and thermally fluctuating nature of proteins.

Signal transduction in cell growth and proliferation involves regulation of kinases through long-range allostery between remote protein regions. Molecular dynamics free energy calculations are used to clarify the coupling between the catalytic domain of Src kinase Hck and its N-terminal end connecting to the regulatory SH2 and SH3 modules. The N-terminal end is stable in the orientation required for the regulatory modules to remain properly bound only in the inactive catalytic domain. In the active catalytic domain, the N-terminal end prefers a different conformation consistent with dissociation of the regulatory modules. The free energy surface shows that the N-terminal end acts as a reversible two-state conformational switch coupling the catalytic domain to the regulatory modules. Structural analogy with insulin receptor kinase and c-Src suggests that such reversible conformational switching in a critical hinge region could be a common mechanism in long-range allosteric regulation of protein kinase activity.

The essential function of potassium channels is to selectively facilitate the passage of K(+) ions across the cell membrane by lowering the energy barrier that opposes this process. In this issue of Neuron, Chatelain et al. report the results from an elegant study that sheds new light on the microscopic factors that modulate the function of these important membrane channels.

The selectivity filter of potassium channels is the structural element directly responsible for the selective and rapid conduction of K+, whereas other parts of the protein are thought to function as a molecular gate that either permits or blocks the passage of ions. However, whether the selectivity filter itself also possesses the ability to play the role of a gate is an unresolved question. Using free energy molecular dynamics simulations, it is shown that the reorientation of two peptide linkages in the selectivity filter of the KcsA K+ channel can lead to a stable nonconducting conformational state. Two microscopic factors influence the transition toward such a conformational state: the occupancy of one specific cation binding site in the selectivity filter (S2), and the strength of intersubunit interactions involving the GYG signature sequence. These results suggest that such conformational transitions occurring in the selectivity filter might be related to different K+ channel gating events, including C-type (slow) inactivation.

Intracellular tetraethylammonium (TEA) inhibition was studied at the single-channel level in the KcsA potassium channel reconstituted in planar lipid bilayers. TEA acts as a fast blocker (resulting in decreased current amplitude) with an affinity in the 75 mM range even at high bandwidth. Studies over a wide voltage range reveal that TEA block has a complex voltage-dependence that also depends on the ionic conditions. These observations are examined in the context of permeation models to extend our understanding of the coupling between permeant ions and TEA blockade.

Identification of the molecular interaction governing ion conduction through biological pores is one of the most important goals of modern electrophysiology. Grand canonical Monte Carlo Brownian dynamics (GCMC/BD) and three-dimensional Poisson-Nernst-Plank (3d-PNP) electrodiffusion algorithms offer powerful and general approaches to study of ion permeation through wide molecular pores. A detailed analysis of ion flows through the staphylococcal alpha-hemolysin channel based on series of simulations at different concentrations and transmembrane potentials is presented. The position-dependent diffusion coefficient is approximated on the basis of a hydrodynamic model. The channel conductance calculated by GCMC/BD is approximately 10% higher than (electrophysiologically measured) experimental values, whereas results from 3d-PNP are always 30-50% larger. Both methods are able to capture all important electrostatic interactions in equilibrium conditions. The asymmetric conductance upon the polarity of the transmembrane potential observed experimentally is reproduced by GCMC/BD and 3d-PNP. The separation of geometrical and energetic influence of the channel on ion conduction reveals that such asymmetries arise from the permanent charge distribution inside the pore. The major determinant of the asymmetry is unbalanced charge in the triad of polar residues D127, D128, and K131. The GCMC/BD or 3d-PNP calculations reproduce also experimental reversal potentials and permeability rations in asymmetric ionic solutions. The weak anionic selectivity of the channel results from the presence of the salt bridge between E111 and K147 in the constriction zone. The calculations also reproduce the experimentally derived dependence of the reversible potential to the direction of the salt gradient. The origin of such effect arises from the asymmetrical distribution of energetic barriers along the channel axis, which modulates the preferential ion passage in different directions., (Copyright 2004 Biophysical Society)

The determination of the structure of several members of the K+ channel and aquaporin family represents a unique opportunity to explain the mechanism of these biomolecular systems. With their ability to go beyond static structures, molecular dynamics simulations offer a unique route for relating functional properties to membrane channel structure. The recent progress in this area is reviewed.

Water transport channels in membrane proteins of the aquaporin superfamily are impermeable to ions, including H+ and OH-. We examine the molecular basis for the blockage of proton translocation through the single-file water chain in the pore of a bacterial aquaporin, GlpF. We compute the reversible thermodynamic work for the two complementary steps of the Grotthuss "hop-and-turn" relay mechanism: consecutive transfers of H+ along the hydrogen-bonded chain (hop) and conformational reorganization of the chain (turn). In the absence of H+, the strong preference for the bipolar orientation of water around the two Asn-Pro-Ala (NPA) motifs lining the pore over both unidirectional polarization states of the chain precludes the reorganization of the hydrogen-bonded network. Inversely, translocation of an excess proton in either direction is opposed by a free-energy barrier centered at the NPA region. Both hop and turn steps of proton translocation are opposed by the electrostatic field of the channel.

A recently proposed model for voltage-dependent activation in K+ channels, largely influenced by the KvAP X-ray structure, suggests that S4 is located at the periphery of the channel and moves through the lipid bilayer upon depolarization. To investigate the physical distance between S4 and the pore domain in functional channels in a native membrane environment, we engineered pairs of cysteines, one each in S4 and the pore of Shaker channels, and identified two instances of spontaneous intersubunit disulfide bond formation, between R362C/A419C and R362C/F416C. After reduction, these cysteine pairs bound Cd2+ with high affinity, verifying that the residues are in atomic proximity. Molecular modeling based on the MthK structure revealed a single position for S4 that was consistent with our results and many other experimental constraints. The model predicts that S4 is located in the groove between pore domains from different subunits, rather than at the periphery of the protein.

Computational methods are used to determine the three-dimensional structure of the Agitoxin (AgTx2)-Shaker complex. In a first stage, a large number of models of the complex are generated using high temperature molecular dynamics, accounting for side chain flexibility with distance restraints deduced from thermodynamic analysis of double mutant cycles. Four plausible binding mode candidates are found using this procedure. In a second stage, the quality and validity of the resulting complexes is assessed by examining the stability of the binding modes during molecular dynamics simulations with explicit water molecules and by calculating the binding free energies of mutant proteins using a continuum solvent representation and comparing with experimental data. The docking protocol and the continuum solvent model are validated using the Barstar-Barnase and the lysozyme-antibody D1.2 complexes, for which there are high-resolution structures as well as double mutant data. This combination of computational methods permits the identification of two possible structural models of AgTx2 in complex with the Shaker K+ channel, additional structural analysis providing further evidence in favor of a single model. In this final complex, the toxin is bound to the extracellular entrance of the channel along the pore axis via a combination of hydrophobic, hydrogen bonding, and electrostatic interactions. The magnitude of the buried solvent accessible area corresponding to the protein-protein contact is on the order of 1000 A with roughly similar contributions from each of the four subunits. Some side chains of the toxin adopt different conformation than in the experimental solution structure, indicating the importance of an induced-fit upon the formation of the complex. In particular, the side chain of Lys-27, a residue highly conserved among scorpion toxins, points deep into the pore with its positively charge amino group positioned at the outer binding site for K+. Specific site-directed mutagenesis experiments are suggested to verify and confirm the structure of the toxin-channel complex.

The conduction of protons in the hydrogen-bonded chain of water molecules (or "proton wire") embedded in the lumen of gramicidin A is studied with molecular dynamics free energy simulations. The process may be described as a "hop-and-turn" or Grotthuss mechanism involving the chemical exchange (hop) of hydrogen nuclei between hydrogen-bonded water molecules arranged in single file in the lumen of the pore, and the subsequent reorganization (turn) of the hydrogen-bonded network. Accordingly, the conduction cycle is modeled by two complementary steps corresponding respectively to the translocation 1) of an ionic defect (H+) and 2) of a bonding defect along the hydrogen-bonded chain of water molecules in the pore interior. The molecular mechanism and the potential of mean force are analyzed for each of these two translocation steps. It is found that the mobility of protons in gramicidin A is essentially determined by the fine structure and the dynamic fluctuations of the hydrogen-bonded network. The translocation of H+ is mediated by spontaneous (thermal) fluctuations in the relative positions of oxygen atoms in the wire. In this diffusive mechanism, a shallow free-energy well slightly favors the presence of the excess proton near the middle of the channel. In the absence of H+, the water chain adopts either one of two polarized configurations, each of which corresponds to an oriented donor-acceptor hydrogen-bond pattern along the channel axis. Interconversion between these two conformations is an activated process that occurs through the sequential and directional reorientation of water molecules of the wire. The effect of hydrogen-bonding interactions between channel and water on proton translocation is analyzed from a comparison to the results obtained previously in a study of model nonpolar channels, in which such interactions were missing. Hydrogen-bond donation from water to the backbone carbonyl oxygen atoms lining the pore interior has a dual effect: it provides a coordination of water molecules well suited both to proton hydration and to high proton mobility, and it facilitates the slower reorientation or turn step of the Grotthuss mechanism by stabilizing intermediate configurations of the hydrogen-bonded network in which water molecules are in the process of flipping between their two preferred, polarized states. This mechanism offers a detailed molecular model for the rapid transport of protons in channels, in energy-transducing membrane proteins, and in enzymes.

The atomic force microscope (AFM) was used to image native OmpF porin and to detect the electrostatic potential generated by the protein. To this end the OmpF porin trimers from Escherichia coli was reproducibly imaged at a lateral resolution of approximately 0.5 nm and a vertical resolution of approximately 0.1 nm at variable electrolyte concentrations of the buffer solution. At low electrolyte concentrations the charged AFM probe not only contoured structural details of the membrane protein surface but also interacted with local electrostatic potentials. Differences measured between topographs recorded at variable ionic strength allowed mapping of the electrostatic potential of OmpF porin. The potential map acquired by AFM showed qualitative agreement with continuum electrostatic calculations based on the atomic OmpF porin embedded in a lipid bilayer at the same electrolyte concentrations. Numerical simulations of the experimental conditions showed the measurements to be reproduced quantitatively when the AFM probe was included in the calculations. This method opens a novel avenue to determine the electrostatic potential of native protein surfaces at a lateral resolution better than 1 nm and a vertical resolution of approximately 0.1 nm.

The side chain of Glu-71 of the KcsA K(+) channel, an important residue in the vicinity of the selectivity filter, was not resolved in the crystallographic structure of Doyle et al. (Doyle, D. A., J. M. Cabral, R. A. Pfuetzner, A. Kuo, J. M. Gulbis, S. L. Cohen, B. T. Chait, and R. MacKinnon. 1998. Science. 280:69-77). Its atomic coordinates are undetermined and its ionization state is unknown. For meaningful theoretical and computational studies of the KcsA K(+) channel, it is essential to address questions about the conformation and the ionization state of this residue in detail. In previous MD simulations in which the side chain of Glu-71 is protonated and forming a strong hydrogen bond with Asp-80 it was observed that the channel did not deviate significantly from the crystallographic structure (Bernèche, S., and B. Roux. 2000. Biophys. J. 78:2900-2917). In contrast, we show here that the structure of the selectivity filter of the KcsA channel is significantly disrupted when these side chains are fully ionized on each of the four monomers. To further resolve questions about the ionization state of Glu-71 we calculated the pK(a) value of this residue using molecular dynamics free energy simulations (MD/FES) with a fully flexible system including explicit solvent and membrane and finite-difference Poisson-Boltzmann (PB) continuum electrostatics. It is found that the pK(a) of Glu-71 is shifted by approximately +10 pK(a) units. These results strongly suggest that Glu-71 is protonated under normal conditions.

The present study is an application of an approach recently developed by the authors for describing the structure of the hydrocarbon chains of lipid-bilayer membranes (LBMs) around embedded protein inclusions ( Biophys. J. 79:2867-2879). The approach is based on statistical mechanical integral equation theories developed for the study of dense liquids. First, the configurations extracted from molecular dynamics simulations of pure LBMs are used to extract the lateral density-density response function. Different pure LBMs composed of different lipid molecules were considered: dioleoyl phosphatidylcholine (DOPC), palmitoyl-oleoyl phosphatidylcholine (POPC), dipalmitoyl phosphatidylcholine (DPPC), and dimyristoyl phosphatidylcholine (DMPC). The results for the lateral density-density response function was then used as input in the integral equation theory. Numerical calculations were performed for protein inclusions of three different sizes. For the sake of simplicity, protein inclusions are represented as hard smooth cylinders excluding the lipid hydrocarbon core from a small cylinder of 2.5 A radius, corresponding roughly to one aliphatic chain, a medium cylinder of 5 A radius, corresponding to one alpha-helix, and a larger cylinder of 9 A radius, representing a small protein such as the gramicidin channel. The lipid-mediated interaction between protein inclusions was calculated using a closed-form expression for the configuration-dependent free energy. This interaction was found to be repulsive at intermediate range and attractive at short range for two small cylinders in POPC, DPPC, and DMPC bilayers, whereas it oscillates between attractive and repulsive values in DOPC bilayers. For medium size cylinders, it is again repulsive at intermediate range and attractive at short range, but for every model LBM considered here. In the case of a large cylinder, the lipid-mediated interaction was shown to be repulsive for both short and long ranges for the DOPC, POPC, and DPPC bilayers, whereas it is again repulsive and attractive for DMPC bilayers. The results indicate that the packing of the hydrocarbon chains around protein inclusions in LBMs gives rise to a generic (i.e., nonspecific) lipid-mediated interaction which favors the association of two alpha-helices and depends on the lipid composition of the membrane.

The effect of C-terminal tyrosine phosphorylation on molecular motions in the Src kinases Hck and c-Src is investigated by molecular dynamics simulations. The SH2 and SH3 domains of the inactive kinases are seen to be tightly coupled by the connector between them, impeding activation. Dephosphorylation of the tail reduces the coupling between the SH2 and SH3 domains in the simulations, as does replacement of connector residues with glycine. A mutational analysis of c-Src expressed in Schizosaccharomyces pombe demonstrates that replacement of residues in the SH2-SH3 connector with glycine activates c-Src. The SH2-SH3 connector appears to be an inducible "snap lock" that clamps the SH2 and SH3 domains upon tail phosphorylation, but which allows flexibility when the tail is released.

This paper describes a framework model for proton conduction through gramicidin; a model designed to incorporate information from molecular dynamics and use this to predict conductance properties. The state diagram describes both motion of an excess proton within the pore as well as the reorientation of waters within the pore in the absence of an excess proton. The model is constructed as the diffusion limit of a random walk, allowing control over the boundary behavior of trajectories. Simple assumptions about the boundary behavior are made, which allow an analytical solution for the proton current and conductance. This is compared with corresponding expressions from statistical mechanics. The random walk construction allows diffusing trajectories underlying the model to be simulated in a simple way. Details of the numerical algorithm are described.