Biomedical subjects

Klaus Schulten

Publications and source records attributed to Klaus Schulten.

At least 19 recordsLinked to original sources

PubMed2006

Association of nuclear pore FG-repeat domains to NTF2 import and export complexes.

Transport into and out of the nucleus is regulated by the nuclear pore complex. Vital to this regulation are nuclear pore proteins with FG sequence repeats, which have been shown to be crucial for cell viability and which interact with nuclear transport receptors. Here we use molecular dynamics simulations to investigate the binding of FG-repeat peptides to the surface of NTF2, the Ran importer. The simulations, covering over 254 ns, agree with previous X-ray, mutational, NMR, and computational data in identifying four binding spots. They also serve to provide an all-atom view of binding at each spot, whereas FG-repeat binding has been only directly observed at a single spot. Furthermore, the simulations identify two novel binding spots in addition to the four others. All six binding spots broadly form a stripe across the surface of NTF2. The resulting regularity and proximity of binding spots on the surface may be necessary for identification of the transport receptor by the FG-repeats in the nuclear pore complex and for the successful transit of NTF2 through the pore.

Active Transport, Cell Nucleus↗

PubMed2006

Ion conduction through MscS as determined by electrophysiology and simulation.

The mechanosensitive channel of small conductance (MscS) is a membrane protein thought to act as a safety valve in bacteria, regulating the release of ions and small solutes when gated by membrane tension under challenging osmotic conditions. The influence of voltage on channel activation and the functional state depicted by the available crystal structure of MscS remain debated. Therefore, in an effort to relate electrophysiological measurements on MscS and properties of the MscS crystal conformation, we report here MscS's response to voltage and pressure as determined by patch-clamp experiments, as well as MscS electrostatics and transport properties as determined through all-atom molecular dynamics simulations of the protein embedded in a lipid bilayer, a 224,000-atom system. The experiments reveal that MscS is a slightly anion-selective channel with a conductance of approximately 1 ns, activated by pressure and inactivated in a voltage-dependent manner. On the other hand, the simulations, covering over 200 ns and including biasing electrostatic potentials, show that MscS restrained to the crystal conformation exhibits low conductance; unrestrained it increases the channel radius upon application of a large electrostatic bias and exhibits then ion conduction that matches experimentally determined conductances. The simulated conductance stems mainly from Cl- ions.

Cell Membrane↗

PubMed2006

Water-silica force field for simulating nanodevices.

Amorphous silica is an inorganic material that is central for many nanotechnology applications, such as nanoelectronics, microfluidics, and nanopore sensors. To use molecular dynamics (MD) simulations to study the behavior of biomolecules interacting with silica, we developed a force field for amorphous silica surfaces based on their macroscopic wetting properties that is compatible with the CHARMM force field and TIP3P water model. The contact angle of a water droplet on a silica surface served as a criterion to tune the intermolecular interactions. The resulting force field was used to study the permeation of water through silica nanopores, illustrating the influence of the surface topography and the intermolecular parameters on permeation kinetics. We find that minute modeling of the amorphous surface is critical for MD studies, since the particular arrangement of surface atoms controls sensitively electrostatic interactions between silica and water.

Nanostructures↗

PubMed2006

How the headpiece hinge angle is opened: New insights into the dynamics of integrin activation.

How the integrin head transitions to the high-affinity conformation is debated. Although experiments link activation with the opening of the hinge angle between the betaA and hybrid domains in the ligand-binding headpiece, this hinge is closed in the liganded alpha(v)beta3 integrin crystal structure. We replaced the RGD peptide ligand of this structure with the 10th type III fibronectin module (FnIII10) and discovered through molecular dynamics (MD) equilibrations that when the conformational constraints of the leg domains are lifted, the betaA/hybrid hinge opens spontaneously. Together with additional equilibrations on the same nanosecond timescale in which small structural variations impeded hinge-angle opening, these simulations allowed us to identify the allosteric pathway along which ligand-induced strain propagates via elastic distortions of the alpha1 helix to the betaA/hybrid domain hinge. Finally, we show with steered MD how force accelerates hinge-angle opening along the same allosteric pathway. Together with available experimental data, these predictions provide a novel framework for understanding integrin activation.

Allosteric Regulation↗

PubMed2006

Charge delocalization in proton channels, I: the aquaporin channels and proton blockage.

The explicit contribution to the free energy barrier and proton conductance from the delocalized nature of the excess proton is examined in aquaporin channels using an accurate all-atom molecular dynamics computer simulation model. In particular, the channel permeation free energy profiles are calculated and compared for both a delocalized (fully Grotthuss shuttling) proton and a classical (nonshuttling) hydronium ion along two aquaporin channels, Aqp1 and GlpF. To elucidate the effects of the bipolar field thought to arise from two alpha-helical macrodipoles on proton blockage, free energy profiles were also calculated for computational mutants of the two channels where the bipolar field was eliminated by artificially discharging the backbone atoms. Comparison of the free energy profiles between the proton and hydronium cases indicates that the magnitude of the free energy barrier and position of the barrier peak for the fully delocalized and shuttling proton are somewhat different from the case of the (localized) classical hydronium. The proton conductance through the two aquaporin channels is also estimated using Poisson-Nernst-Planck theory for both the Grotthuss shuttling excess proton and the classical hydronium cation.

Aquaporin 1↗

PubMed2006

Coarse-grained molecular dynamics simulations of a rotating bacterial flagellum.

Many types of bacteria propel themselves using elongated structures known as flagella. The bacterial flagellar filament is a relatively simple and well-studied macromolecular assembly, which assumes different helical shapes when rotated in different directions. This polymorphism enables a bacterium to switch between running and tumbling modes; however, the mechanism governing the filament polymorphism is not completely understood. Here we report a study of the bacterial flagellar filament using numerical simulations that employ a novel coarse-grained molecular dynamics method. The simulations reveal the dynamics of a half-micrometer-long flagellum segment on a timescale of tens of microseconds. Depending on the rotation direction, specific modes of filament coiling and arrangement of monomers are observed, in qualitative agreement with experimental observations of flagellar polymorphism. We find that solvent-protein interactions are likely to contribute to the polymorphic helical shapes of the filament.

Bacteria↗

PubMed2006

Sugar binding and protein conformational changes in lactose permease.

Lactose permease is an integral membrane protein that uses the cell membrane's proton gradient for import of lactose. Based on extensive biochemical data and a substrate-bound crystal structure, intermediates involved in lactose/H(+) co-transport have been suggested. Yet, the transport mechanism, especially the coupling of protonation states of essential residues and protein conformational changes involved in the transport, is not understood. Here we report molecular-dynamics simulations of membrane-embedded lactose permease in different protonation states, both in the presence and in the absence of lactose. The results analyzed in terms of pore diameter, salt-bridge formation, and substrate motion, strongly implicate Glu(269) as one of the main proton translocation sites, whose protonation state controls several key steps of the transport process. A critical ion pair (Glu(269) and Arg(144)) was found to keep the cytoplasmic entrance open, but via a different mechanism than the currently accepted model. After protonation of Glu(269), the salt bridge between Glu(269) and Arg(144) was found to break, and Arg(144) to move away from Glu(269), establishing a new salt bridge with Glu(126); furthermore, neutralization of Glu(269) and the displacement of Arg(144) and consequently of water molecules from the interdomain region was seen to initiate the closing of the cytoplasmic half channel (2.6-4.0 A reduction in diameter in the cytoplasmic constriction region in 10 ns) by allowing hydrophobic surfaces of the N- and C-domains to fuse. Charged Glu(269) was found to strongly bind the lactose permeant, indicating that proton transfer from water or another residue to Glu(269) is a prerequisite for unbinding of lactose from the binding pocket.

Arginine↗

PubMed2006

Dynamic switching mechanisms in LOV1 and LOV2 domains of plant phototropins.

LOV domains are the light-sensitive portion of plant phototropins. They absorb light through a flavin cofactor, photochemically form a covalent bond between the chromophore and a cysteine residue in the protein, and proceed to mediate activation of an attached kinase domain. Although the photoreaction itself is now well-characterized experimentally and computationally, it is still unclear how the formation of the adduct leads to kinase activation. We have performed molecular dynamics simulations on the LOV1 domain of Chlamydomonas reinhardtii and the LOV2 domain of Avena sativa, both before and after the photoreaction, to answer this question. The extensive simulations, over 240 ns in duration, reveal significant differences in how the LOV1 and LOV2 domains respond to photoactivation. The simulations indicate that LOV1 activation is likely caused by a change in hydrogen bonding between protein and ligand that destabilizes a highly conserved salt bridge, whereas LOV2 activation seems to result from a change in the flexibility of a set of protein loops. Results of electrostatics calculations, principal component analysis, sequence alignments, and root mean-square deviation analysis corroborate the above findings.

Computer Simulation↗

PubMed2006

Assembly of lipoprotein particles revealed by coarse-grained molecular dynamics simulations.

High-density lipoproteins (HDL) function as cholesterol transporters, facilitating the removal of excess cholesterol from the body. Due to the heterogeneity of native HDL particles (both in size and shape), the details on how these protein-lipid particles form and the structure they assume in their lipid-associated states are not well characterized. We report here a study of the self-assembly of discoidal HDL particles using coarse-grained (CG) molecular dynamics. The microsecond simulations reveal the self-assembly of HDL particles from disordered protein-lipid complexes to form structures containing many of the features of the generally accepted double-belt model for discoidal HDL particles. HDL assembly is found to proceed in two broad steps, aggregation of proteins and lipids driven by the hydrophobic effect which occurs on a approximately 1 micros time scale, followed by the optimization of the protein structure driven by increasingly specific protein-protein interactions.

Lipoproteins, HDL↗

PubMed2006

Conformational equilibrium in alanine-rich peptides probed by reversible stretching simulations.

Reversible stretching of the alanine-rich peptide 3K (Proc. Natl. Acad. Sci. USA 1989, 86, 5286-5290) and its analogue MW (Nature 1992, 359, 653-655) is examined using molecular dynamics simulations in explicit water. In both cases, sampling of the extension pathway is obtained on the 10 ns time scale by applying an adaptive biasing force. The free energy profile reveals a single minimum associated with a contiguous alpha-helix. Short 3(10)-helical motifs are observed in folded as well as extended conformations, in accordance with their proposed role as folding intermediates. The native 3(10)-helical content of both peptides is found, however, to be no higher than a few percent. Difficulties in both the definition and the detection of secondary structure motifs, most notably in relation to bifurcated hydrogen bonds, are proposed to account for the discrepancy between 3(10)-helical propensities reported by several authors, based on experimental and computational results.

Alanine↗

PubMed2006

Color tuning in rhodopsins: the mechanism for the spectral shift between bacteriorhodopsin and sensory rhodopsin II.

The mechanism of color tuning in the rhodopsin family of proteins has been studied by comparing the optical properties of the light-driven proton pump bacteriorhodopsin (bR) and the light detector sensory rhodopsin II (sRII). Despite a high structural similarity, the maximal absorption is blue-shifted from 568 nm in bR to 497 nm in sRII. The molecular mechanism of this shift is still a matter of debate, and its clarification sheds light onto the general mechanisms of color tuning in retinal proteins. The calculations employ a combined quantum mechanical/molecular mechanical (QM/MM) technique, using a DFT-based method for ground state properties and the semiempirical OM2/MRCI method and ab initio SORCI method for excited state calculations. The high efficiency of the methodology has allowed us to study a wide variety of aspects including dynamical effects. The absorption shift as well as various mutation experiments and vibrational properties have been successfully reproduced. Our results indicate that several sources contribute to the spectral shift between bR and sRII. The main factors are the counterion region at the extracellular side of retinal and the amino acid composition of the binding pocket. Our analysis allows a distinction and identification of the different effects in detail and leads to a clear picture of the mechanism of color tuning, which is in good agreement with available experimental data.

Bacteriorhodopsins↗

PubMed2006

Dynamics of K+ ion conduction through Kv1.2.

The crystallographic structure of a potassium channel, Kv1.2, in an open state makes it feasible to simulate entire K(+) ion permeation events driven by a voltage bias and, thereby, elucidate the mechanism underlying ion conduction and selectivity of this type of channel. This Letter demonstrates that molecular dynamics simulations can provide movies of the overall conduction of K(+) ions through Kv1.2. As suggested earlier, the conduction is concerted in the selectivity filter, involving 2-3 ions residing mainly at sites identified previously by crystallography and modeling. The simulations reveal, however, the jumps of ions between these sites and identify the sequence of multi-ion configurations involved in permeation.

Computer Simulation↗

PubMed2006

Molecular mechanisms of cellular mechanics.

Mechanical forces play an essential role in cellular processes as input, output, and signals. Various protein complexes in the cell are designed to handle, transform and use such forces. For instance, proteins of muscle and the extracellular matrix can withstand considerable stretching forces, hearing-related and mechanosensory proteins can transform weak mechanical stimuli into electrical signals, and regulatory proteins are suited to forcing DNA into loops to control gene expression. Here we review the structure-function relationship of four protein complexes with well defined and representative mechanical functions. The first example is titin, a protein that confers passive elasticity on muscle. The second system is the elastic extracellular matrix protein, fibronectin, and its cellular receptor integrin. The third protein system is the transduction apparatus in hearing and other mechanical senses, likely containing cadherin and ankyrin repeats. The last system is the lac repressor protein, which regulates gene expression by looping DNA. This review focuses on atomic level descriptions of the physical mechanisms underlying the various mechanical functions of the stated proteins.

Biomechanical Phenomena↗

PubMed2006

Structure-based model of the stepping motor of PcrA helicase.

DNA helicases are ubiquitous molecular motors involved in cellular DNA metabolism. They move along single-stranded DNA (ssDNA) and separate duplex DNA into its component strands, utilizing the free energy from ATP hydrolysis. The PcrA helicase from Bacillus stearothermophilus translocates as a monomer progressively from the 3' end to the 5' end of ssDNA and is one of the smallest motor proteins structurally known in full atomic detail. Using high-resolution crystal structures of the PcrA-DNA complex, we performed nanosecond molecular dynamics simulations and derived potential energy profiles governing individual domain movement of the PcrA helicase along ssDNA. Based on these profiles, the millisecond translocation of the helicase along ssDNA was described through Langevin dynamics. The calculations support a domain stepping mechanism of PcrA helicase, in which, during one ATP hydrolysis cycle, the pulling together and pushing apart of domains 2A and 1A are synchronized with alternating mobilities of the individual domains in such a fashion that PcrA moves unidirectionally along ssDNA. By combining short timescale (nanoseconds) molecular dynamics and long timescale (milliseconds) stochastic-dynamics descriptions, our study suggests a structure-based mechanism of the ATP-powered unidirectional movement of PcrA helicase.

Adenosine Triphosphate↗

PubMed2006

Molecular dynamics study of MscL interactions with a curved lipid bilayer.

Mechanosensitivity is a ubiquitous sensory mechanism found in living organisms. The simplest known mechanotransducing mechanism is found in bacteria in the form of the mechanosensitive membrane channel of large conductance, MscL. This channel has been studied extensively using a variety of methods at a functional and structural level. The channel is gated by membrane tension in the lipid bilayer alone. It serves as a safety valve protecting bacterial cells against hypoosmotic shock. MscL of Escherichia coli embedded in bilayers composed of asymmetric amounts of single-tailed and double-tailed lipids has been shown to gate spontaneously, even in the absence of membrane tension. To gain insight into the effect of the lipid membrane composition and geometry on MscL structure, a fully solvated, all-atom model of MscL in a stress-free curved bilayer composed of double- and single-tailed lipids was studied using a 9.5-ns molecular dynamics simulation. The bilayer was modeled as a domed structure accommodating the asymmetric composition of the monolayers. During the course of the simulation a spontaneous restructuring of the periplasmic loops occurred, leading to interactions between one of the loops and phospholipid headgroups. Previous experimental studies of the role of the loops agree with the observation that opening starts with a restructuring of the periplasmic loop, suggesting an effect of the curved bilayer. Because of limited resources, only one simulation of the large system was performed. However, the results obtained suggest that through the geometry and composition of the bilayer the protein structure can be affected even on short timescales.

Computer Simulation↗

PubMed2006

Imaging the migration pathways for O2, CO, NO, and Xe inside myoglobin.

Myoglobin (Mb) is perhaps the most studied protein, experimentally and theoretically. Despite the wealth of known details regarding the gas migration processes inside Mb, there exists no fully conclusive picture of these pathways. We address this deficiency by presenting a complete map of all the gas migration pathways inside Mb for small gas ligands (O2, NO, CO, and Xe). To accomplish this, we introduce a computational approach for studying gas migration, which we call implicit ligand sampling. Rather than simulating actual gas migration events, we infer the location of gas migration pathways based on a free-energy perturbation approach applied to simulations of Mb's dynamical fluctuations at equilibrium in the absence of ligand. The method provides complete three-dimensional maps of the potential of mean force of gas ligand placement anywhere inside a protein-solvent system. From such free-energy maps we identify each gas docking site, the pathways between these sites, to the heme and to the external solution. Our maps match previously known features of these pathways in Mb, but also point to the existence of additional exits from the protein matrix in regions that are not easily probed by experiment. We also compare the pathway maps of Mb for different gas ligands and for different animal species.

Carbon Monoxide↗

PubMed2006

Modeling DNA loops using the theory of elasticity.

An elastic rod model of a protein-bound DNA loop is adapted for application in multi-scale simulations of protein-DNA complexes. The classical Kirchhoff system of equations which describes the equilibrium structure of the elastic loop is modified to account for the intrinsic twist and curvature, anisotropic bending properties, and electrostatic charge of DNA. The effects of bending anisotropy and electrostatics are studied for the DNA loop clamped by the lac repressor protein. For two possible lengths of the loop, several topologically different conformations are predicted and extensively analyzed over the broad range of model parameters describing DNA bending and electrostatic properties. The scope and applications of the model in already accomplished and in future multi-scale studies of protein-DNA complexes are discussed.

Computer Simulation↗

PubMed2006

Coarse grained protein-lipid model with application to lipoprotein particles.

A coarse-grained model for molecular dynamics simulations is extended from lipids to proteins. In the framework of such models pioneered by Klein, atoms are described group-wise by beads, with the interactions between beads governed by effective potentials. The extension developed here is based on a coarse-grained lipid model developed previously by Marrink et al., although future versions will reconcile the approach taken with the systematic approach of Klein and other authors. Each amino acid of the protein is represented by two coarse-grained beads, one for the backbone (identical for all residues) and one for the side-chain (which differs depending on the residue type). The coarse-graining reduces the system size about 10-fold and allows integration time steps of 25-50 fs. The model is applied to simulations of discoidal high-density lipoprotein particles involving water, lipids, and two primarily helical proteins. These particles are an ideal test system for the extension of coarse-grained models. Our model proved to be reliable in maintaining the shape of preassembled particles and in reproducing the overall structural features of high-density lipoproteins accurately. Microsecond simulations of lipoprotein assembly revealed the formation of a protein-lipid complex in which two proteins are attached to either side of a discoidal lipid bilayer.

Algorithms↗