PubMed Health⌕ Search

Biomedical subjects

Klaus Schulten

Publications and source records attributed to Klaus Schulten.

At least 19 recordsLinked to original sources
PubMed2004

The mechanism of proton exclusion in aquaporin channels.

The mechanism of proton exclusion in aquaporin channels is elucidated through free energy calculations of the pathway of proton transport. The second generation multistate empirical valence bond (MS-EVB2) model was applied to simulate the interaction of an excess proton with the channel environment. Jarzynski's equality was employed for rapid convergence of the free energy profile. A barrier sufficiently high to block proton transport is located near the channel center at the NPA motif-a site involved in bi-orientational ordering of the embedded water-wire in absence of the excess proton. A second and lower barrier is observed at the selectivity filter near the periplasmic outlet where the channel is narrowest. This secondary barrier may be essential in filtering other large solutes and cations.

Aquaporins↗
PubMed2004

Tuning the mechanical stability of fibronectin type III modules through sequence variations.

Cells can switch the functional states of extracellular matrix proteins by stretching them while exerting mechanical force. Using steered molecular dynamics, we investigated how the mechanical stability of FnIII modules from the cell adhesion protein fibronectin is affected by natural variations in their amino acid sequences. Despite remarkably similar tertiary structures, FnIII modules share low sequence homology. Conversely, the sequence homology for the same FnIII module across multiple species is notably higher, suggesting that sequence variability is functionally significant. Our studies find that the mechanical stability of FnIII modules can be tuned through substitutions of just a few key amino acids by altering access of water molecules to hydrogen bonds that break early in the unfolding pathway. Furthermore, the FnIII hierarchy of mechanical unfolding can be changed by environmental conditions, such as pH for FnIII10, or by forming complexes with other molecules, such as heparin binding to FnIII13.

Amino Acid Sequence↗
PubMed2004

Structural basis for cooperative DNA binding by CAP and lac repressor.

Catabolite gene activator protein (CAP) and lac repressor (LR) are celebrated transcription-regulating proteins that bind to DNA cooperatively forming a ternary complex with the promoter loop. Here we present a multiscale model of the ternary complex derived from crystal structures of the proteins and a continuous structure of the DNA loop built using the theory of elasticity. We predict that the loop is underwound in the binary complex with the LR, whereas in the ternary complex with the LR and CAP, the loop is overwound and extended due to an upstream relocation of a DNA binding hand of LR. The computed relocation distance matches the experimental observations and the energy balance of the system explains the cooperativity effect. Using the multiscale approach, we build an all-atom model of the ternary complex that suggests a series of further experimental investigations.

Binding, Competitive↗
PubMed2004

Theory and simulation of water permeation in aquaporin-1.

We discuss the difference between osmotic permeability pf and diffusion permeability pd of single-file water channels and demonstrate that the pf/pd ratio corresponds to the number of effective steps a water molecule needs to take to permeate a channel. While pd can be directly obtained from equilibrium molecular dynamics simulations, pf can be best determined from simulations in which a chemical potential difference of water has been established on the two sides of the channel. In light of this, we suggest a method to induce in molecular dynamics simulations a hydrostatic pressure difference across the membrane, from which pf can be measured. Simulations using this method are performed on aquaporin-1 channels in a lipid bilayer, resulting in a calculated pf of 7.1 x 10(-14) cm(3)/s, which is in close agreement with observation. Using a previously determined pd value, we conclude that pf/pd for aquaporin-1 measures approximately 12. This number is explained in terms of channel architecture and conduction mechanism.

Aquaporin 1↗
PubMed2004

Mechanism of anionic conduction across ClC.

ClC chloride channels are voltage-gated transmembrane proteins that have been associated with a wide range of regulatory roles in vertebrates. To accomplish their function, they allow small inorganic anions to efficiently pass through, while blocking the passage of all other particles. Understanding the conduction mechanism of ClC has been the subject of many experimental investigations, but until now, the detailed dynamic mechanism was not known despite the availability of crystallographic structures. We investigate Cl(-) conduction by means of an all-atom molecular dynamics simulation of the ClC channel in a membrane environment. Based on our simulation results, we propose a king-of-the-hill mechanism for permeation, in which a lone ion bound to the center of the ClC pore is pushed out by a second ion that enters the pore and takes its place. Although the energy required to extract the single central ion from the pore is enormous, by resorting to this two-ion process, the largest free energy barrier for conduction is reduced to 4 kcal/mol. At the narrowest part of the pore, residues Tyr-445 and Ser-107 stabilize the central ion. There, the bound ion blocks the pore, disrupting the formation of a continuous water file that could leak protons, possibly preventing the passage of uncharged solutes.

Amino Acids↗
PubMed2004

Insights into the molecular mechanism of rotation in the Fo sector of ATP synthase.

F(1)F(o)-ATP synthase is a ubiquitous membrane protein complex that efficiently converts a cell's transmembrane proton gradient into chemical energy stored as ATP. The protein is made of two molecular motors, F(o) and F(1), which are coupled by a central stalk. The membrane unit, F(o), converts the transmembrane electrochemical potential into mechanical rotation of a rotor in F(o) and the physically connected central stalk. Based on available data of individual components, we have built an all-atom model of F(o) and investigated through molecular dynamics simulations and mathematical modeling the mechanism of torque generation in F(o). The mechanism that emerged generates the torque at the interface of the a- and c-subunits of F(o) through side groups aSer-206, aArg-210, and aAsn-214 of the a-subunit and side groups cAsp-61 of the c-subunits. The mechanism couples protonation/deprotonation of two cAsp-61 side groups, juxtaposed to the a-subunit at any moment in time, to rotations of individual c-subunit helices as well as rotation of the entire c-subunit. The aArg-210 side group orients the cAsp-61 side groups and, thereby, establishes proton transfer via aSer-206 and aAsn-214 to proton half-channels, while preventing direct proton transfer between the half-channels. A mathematical model proves the feasibility of torque generation by the stated mechanism against loads typical during ATP synthesis; the essential model characteristics, e.g., helix and subunit rotation and associated friction constants, have been tested and furnished by steered molecular dynamics simulations.

Binding Sites↗
PubMed2003

Structure and functional significance of mechanically unfolded fibronectin type III1 intermediates.

Fibronectin (FN) forms fibrillar networks coupling cells to the extracellular matrix. The formation of FN fibrils, fibrillogenesis, is a tightly regulated process involving the exposure of cryptic binding sites in individual FN type III (FN-III) repeats presumably exposed by mechanical tension. The FN-III1 module has been previously proposed to contain such cryptic sites that promote the assembly of extracellular matrix FN fibrils. We have combined NMR and steered molecular dynamics simulations to study the structure and mechanical unfolding pathway of FN-III1. This study finds that FN-III1 consists of a beta-sandwich structure that unfolds to a mechanically stable intermediate about four times the length of the native folded state. Considering previous experimental findings, our studies provide a structural model by which mechanical stretching of FN-III1 may induce fibrillogenesis through this partially unfolded intermediate.

Amino Acid Sequence↗
PubMed2003

Mechanisms of selectivity in channels and enzymes studied with interactive molecular dynamics.

Interactive molecular dynamics, a new modeling tool for rapid investigation of the physical mechanisms of biological processes at the atomic level, is applied to study selectivity and regulation of the membrane channel protein GlpF and the enzyme glycerol kinase. These proteins facilitate the first two steps of Escherichia coli glycerol metabolism. Despite their different function and architecture the proteins are found to employ common mechanisms for substrate selectivity: an induced geometrical fit by structurally homologous binding sites and an induced rapid dipole moment reversal. Competition for hydrogen bonding sites with water in both proteins is critical for substrate motion. In glycerol kinase, it is shown that the proposed domain motion prevents competition with water, in turn regulating the binding of glycerol.

Aquaporins↗
PubMed2003

Water and proton conduction through carbon nanotubes as models for biological channels.

Carbon nanotubes, unmodified (pristine) and modified through charged atoms, were simulated in water, and their water conduction rates determined. The conducted water inside the nanotubes was found to exhibit a strong ordering of its dipole moments. In pristine nanotubes the water dipoles adopt a single orientation along the tube axis with a low flipping rate between the two possible alignments. Modification can induce in nanotubes a bipolar ordering as previously observed in biological water channels. Network thermodynamics was applied to investigate proton conduction through the nanotubes.

Biomimetic Materials↗
PubMed2003

Molecular dynamics simulation of bacteriorhodopsin's photoisomerization using ab initio forces for the excited chromophore.

Retinal proteins are photoreceptors found in many living organisms. They possess a common chromophore, retinal, that upon absorption of light isomerizes and thereby triggers biological functions ranging from light energy conversion to phototaxis and vision. The photoisomerization of retinal is extremely fast, highly selective inside the protein matrix, and characterized through optimal sensitivity to incoming light. This article describes the first report of an ab initio quantum mechanical description of the in situ isomerization dynamics of retinal in bacteriorhodopsin, a microbial retinal protein that functions as a light-driven proton pump. The description combines ab initio multi-electronic state molecular dynamics of a truncated retinal chromophore model (N-methyl-gamma-methylpenta-2,4-dieniminium cation fragment) with molecular mechanics of the protein motion and unveils in complete detail the photoisomerization process. The results illustrate the essential role of the protein for the characteristic kinetics and high selectivity of the photoisomerization: the protein arrests inhomogeneous photoisomerization paths and funnels them into a single path that initiates the functional process. Supported by comparison with dynamic spectral modulations observed in femtosecond spectroscopy, the results identify the principal molecular motion during photoisomerization.

Bacteriorhodopsins↗
PubMed2003

Gating of MscL studied by steered molecular dynamics.

Steered molecular dynamics simulations of the mechanosensitive channel of large conductance, MscL, were used to investigate how forces arising from membrane tension induce gating of the channel. A homology model of the closed form of MscL from Escherichia coli was subjected to external forces of 35-70 pN applied to residues near the membrane-water interface. The magnitude and location of these forces corresponded to those determined from the lateral pressure profile computed from a lipid bilayer simulation. A fully expanded state was obtained on the 10-ns timescale that revealed the mechanism for transducing membrane forces into channel opening. The expanded state agrees well with proposed models of MscL gating, in that it entails an irislike expansion of the pore accompanied by tilting of the transmembrane helices. The channel was most easily opened when force was applied predominantly on the cytoplasmic side of MscL. Comparison of simulations in which gating progressed to varying degrees identified residues that pose steric hindrance to channel opening.

Binding Sites↗
PubMed2003

On the mechanism of ATP hydrolysis in F1-ATPase.

Most of the cellular ATP in living organisms is synthesized by the enzyme F(1)F(o)-ATP synthase. The water soluble F(1) part of the enzyme can also work in reverse and utilize the chemical energy released during ATP hydrolysis to generate mechanical motion. Despite the availability of a large amount of biochemical data and several x-ray crystallographic structures of F(1), there still remains a considerable lack of understanding as to how this protein efficiently converts the chemical energy released during the reaction ATP + H(2)O --> ADP + P(i) into mechanical motion of the stalk. We report here an ab initio QM/MM study of ATP hydrolysis in the beta(TP) catalytic site of F(1). Our simulations provide an atomic level description of the reaction path, its energetics, and the interaction of the nucleotide with the protein environment during catalysis. The simulations suggest that the reaction path with the lowest potential energy barrier proceeds via nucleophilic attack on the gamma-phosphate involving two water molecules. Furthermore, the ATP hydrolysis reaction in beta(TP) is found to be endothermic, demonstrating that the catalytic site is able to support the synthesis of ATP and does not promote ATP hydrolysis in the particular conformation studied.

Adenosine Triphosphate↗
PubMed2003

Electrostatic tuning of permeation and selectivity in aquaporin water channels.

Water permeation and electrostatic interactions between water and channel are investigated in the Escherichia coli glycerol uptake facilitator GlpF, a member of the aquaporin water channel family, by molecular dynamics simulations. A tetrameric model of the channel embedded in a 16:0/18:1c9-palmitoyloleylphosphatidylethanolamine membrane was used for the simulations. During the simulations, water molecules pass through the channel in single file. The movement of the single file water molecules through the channel is concerted, and we show that it can be described by a continuous-time random-walk model. The integrity of the single file remains intact during the permeation, indicating that a disrupted water chain is unlikely to be the mechanism of proton exclusion in aquaporins. Specific hydrogen bonds between permeating water and protein at the channel center (at two conserved Asp-Pro-Ala "NPA" motifs), together with the protein electrostatic fields enforce a bipolar water configuration inside the channel with dipole inversion at the NPA motifs. At the NPA motifs water-protein electrostatic interactions facilitate this inversion. Furthermore, water-water electrostatic interactions are in all regions inside the channel stronger than water-protein interactions, except near a conserved, positively charged Arg residue. We find that variations of the protein electrostatic field through the channel, owing to preserved structural features, completely explain the bipolar orientation of water. This orientation persists despite water translocation in single file and blocks proton transport. Furthermore, we find that for permeation of a cation, ion-protein electrostatic interactions are more unfavorable at the conserved NPA motifs than at the conserved Arg, suggesting that the major barrier against proton transport in aquaporins is faced at the NPA motifs.

Aquaporins↗
PubMed2003

Glycerol conductance and physical asymmetry of the Escherichia coli glycerol facilitator GlpF.

The aquaglyceroporin GlpF is a transmembrane channel of Escherichia coli that facilitates the uptake of glycerol by the cell. Its high glycerol uptake rate is crucial for the cell to survive in very low glycerol concentrations. Although GlpF allows both influx and outflux of glycerol, its structure, similar to the structure of maltoporin, exhibits a significant degree of asymmetry. The potential of mean force characterizing glycerol in the channel shows a corresponding asymmetry with an attractive vestibule only at the periplasmic side. In this study, we analyze the potential of mean force, showing that a simplified six-step model captures the kinetics and yields a glycerol conduction rate that agrees well with observation. The vestibule improves the conduction rate by 40% and 75% at 10- micro M and 10-mM periplasmic glycerol concentrations, respectively. In addition, neither the conduction rate nor the conduction probability for a single glycerol (efficiency) depends on the orientation of GlpF. GlpF appears to conduct equally well in both directions under physiological conditions.

Aquaporins↗
PubMed2003

Serum levels of insulin-like growth factor (IGF)-I and IGF binding protein (IGFBP)-1 to -6 and their relationship to bone metabolism in osteoporosis patients.

BACKGROUND: Insulin-like growth factor (IGF) system components are important regulators of bone formation. Alterations of individual IGF system components have been described in osteoporosis (OP) patients; however, no study has addressed changes in free IGF-I and in all six IGF binding proteins (IGFBPs). METHODS: A cross-sectional study was performed in 45 OP patients and 100 healthy matched controls. Serum levels of free and total insulin-like growth factor I (IGF-I), IGFBP-1 through -6, intact parathyroid hormone (PTH), 25-OH-vitamin D(3) (25OHD(3)), 1,25-(OH)(2)-vitamin D(3) (1,25-(OH)(2)D(3)), osteocalcin (OSC), bone alkaline phosphatase (B-ALP), and carboxyterminal propeptide of type-I procollagen (PICP) were measured with specific assays. Bone mineral density (BMD) of the lumbar spine was determined by dual-energy X-ray absorptiometry (DEXA). RESULTS: Compared with age- and sex-matched control subjects, OP patients showed a 73% decrease in free IGF-I, a 29% decrease in total IGF-I, a 10% decrease in IGFBP-3, and a 52% decrease in IGFBP-5 levels; they had higher levels of IGFBP-1 (4.1-fold), IGFBP-2 (1.8-fold), IGFBP-4 (1.3-fold), and IGFBP-6 (2.1-fold). Alterations in IGF system components were most evident in 13 OP patients with vertebral fractures in the past 4 years compared to patients without fractures. In OP patients with fractures, the ratio between IGFBP-4 and IGFBP-5 was increased whereas levels of OSC were decreased. CONCLUSIONS: Our data provide strong indirect evidence for a functional connection between circulating IGF system components and bone metabolism and the susceptibility to fractures in OP patients.

Journal Article↗
PubMed2003

Molecular biomimetics: nanotechnology through biology.

Proteins, through their unique and specific interactions with other macromolecules and inorganics, control structures and functions of all biological hard and soft tissues in organisms. Molecular biomimetics is an emerging field in which hybrid technologies are developed by using the tools of molecular biology and nanotechnology. Taking lessons from biology, polypeptides can now be genetically engineered to specifically bind to selected inorganic compounds for applications in nano- and biotechnology. This review discusses combinatorial biological protocols, that is, bacterial cell surface and phage-display technologies, in the selection of short sequences that have affinity to (noble) metals, semiconducting oxides and other technological compounds. These genetically engineered proteins for inorganics (GEPIs) can be used in the assembly of functional nanostructures. Based on the three fundamental principles of molecular recognition, self-assembly and DNA manipulation, we highlight successful uses of GEPI in nanotechnology.

Amino Acid Sequence↗