PubMed Health⌕ Search

Biomedical subjects

K Schulten

Publications and source records attributed to K Schulten.

At least 19 recordsLinked to original sources
PubMed2006

Possible pathway for ubiquinone shuttling in Rhodospirillum rubrum revealed by molecular dynamics simulation.

In the last decade, the structures of many components of the photosynthetic apparatus of purple bacteria, as well as the mutual organization of these components within the purple membrane, were resolved. One key question that emerged concerned the assembly of the core complex consisting of the reaction center (RC) and the light-harvesting 1 (LH1) complex. In some species, like Rhodobacter sphaeroides, the ring-shaped LH1 complex was found to be open, whereas other species, like Rhodospirillum rubrum, have a closed ring surrounding the reaction center. This poses the question of how the ubiquinone molecule that transports electrons and protons from the RC to the cytochrome bc(1) complex overcomes the apparent barrier of the LH1 ring. In this study, we investigated how, in the case of a closed LH1 ring, the ubiquinone molecule diffuses through the LH1 ring. For this purpose, the LH1 structure of R. rubrum was modeled and the potential of mean force along the diffusion pathway through the LH1 was determined by steered molecular-dynamics simulations. The potential was reconstructed using the fluctuation theorem in combination with the stiff spring approximation. An upper limit for the mean first-passage time for diffusion of ubiquinone through the LH1 ring, based on a worst-case scenario potential, was calculated as approximately 8 x 10(-3) s, which is still in agreement with known turnover rates of RC and RC-LH1 complexes in the range of approximately 1000 Hz.

Biological Transport↗
PubMed2005

The electromechanics of DNA in a synthetic nanopore.

We have explored the electromechanical properties of DNA on a nanometer-length scale using an electric field to force single molecules through synthetic nanopores in ultrathin silicon nitride membranes. At low electric fields, E < 200 mV/10 nm, we observed that single-stranded DNA can permeate pores with a diameter >/=1.0 nm, whereas double-stranded DNA only permeates pores with a diameter >/=3 nm. For pores <3.0 nm diameter, we find a threshold for permeation of double-stranded DNA that depends on the electric field and pH. For a 2 nm diameter pore, the electric field threshold is approximately 3.1 V/10 nm at pH = 8.5; the threshold decreases as pH becomes more acidic or the diameter increases. Molecular dynamics indicates that the field threshold originates from a stretching transition in DNA that occurs under the force gradient in a nanopore. Lowering pH destabilizes the double helix, facilitating DNA translocation at lower fields.

Biophysics↗
PubMed2005

Approaches to developing biological H(2)-photoproducing organisms and processes.

The development of efficient biological systems for the direct photoproduction of H(2) gas from water faces several challenges, the more serious of which is the sensitivity of the H(2)-evolving enzymes (hydrogenases) to O(2), an obligatory by-product of photosynthesis. This high sensitivity is common to both FeFe and NiFe hydrogenases, and is caused by O(2) binding to their respective metallocatalytic sites. This overview describes approaches to (i) molecular engineering of algal FeFe-hydrogenase to prevent O(2) access to its catalytic site; (ii) transform a cyanobacterium with an O(2)-tolerant bacterial NiFe hydrogenase or (c) partially inactivate algal O(2)-evolution activity to create physiologically anaerobiosis and induce hydrogenase expression.

Animals↗
PubMed2005

Molecular dynamics and experimental investigation of H(2) and O(2) diffusion in [Fe]-hydrogenase.

The [Fe]-hydrogenase enzymes are highly efficient H(2) catalysts found in ecologically and phylogenetically diverse microorganisms, including the photosynthetic green alga, Chlamydomonas reinhardtii. Although these enzymes can occur in several forms, H(2) catalysis takes place at a unique [FeS] prosthetic group or H-cluster, located at the active site. Significant to the function of hydrogenases is how the surrounding protein structure facilitates substrate-product transfer, and protects the active site H-cluster from inactivation. To elucidate the role of protein structure in O(2) inactivation of [Fe]-hydrogenases, experimental and theoretical investigations have been performed. Molecular dynamics was used to comparatively investigate O(2) and H(2) diffusion in CpI ([Fe]-hydrogenase I from Clostridium pasteurianum). Our preliminary results suggest that H(2) diffuses more easily and freely than O(2), which is restricted to a small number of allowed pathways to and from the active site. These O(2) pathways are located in the conserved active site domain, shown experimentally to have an essential role in active site protection.

Catalysis↗
PubMed2004

Modelling DNA loops using continuum and statistical mechanics.

The classical Kirchhoff elastic-rod model applied to DNA is extended to account for sequence-dependent intrinsic twist and curvature, anisotropic bending rigidity, electrostatic force interactions, and overdamped Brownian motion in a solvent. The zero-temperature equilibrium rod model is then applied to study the structural basis of the function of the lac repressor protein in the lac operon of Escherichia coli. The structure of a DNA loop induced by the clamping of two distant DNA operator sites by lac repressor is investigated and the optimal geometries for the loop of length 76 bp are predicted. Further, the mimicked binding of catabolite gene activator protein (CAP) inside the loop provides solutions that might explain the experimentally observed synergy in DNA binding between the two proteins. Finally, a combined Monte Carlo and Brownian dynamics solver for a worm-like chain model is described and a preliminary analysis of DNA loop-formation kinetics is presented.

Base Sequence↗
PubMed2004

Multiple probes reveal a native-like intermediate during low-temperature refolding of ubiquitin.

We investigate the refolding of ubiquitin Phe45Trp/Ile61Ala (Ub(*)I61A) in a low-temperature, high-viscosity buffer, where folding is slowed so that apparent two-state and three-state mechanisms are readily distinguishable. Ub(*)I61A forms a compact ensemble rapidly (as judged from stopped-flow, small-angle X-ray scattering) with a secondary structure signature similar to that of the native state (as judged from stopped-flow circular dichroism from 215 nm to 250 nm), but the fluorescence signature still resembles the guanidinium-denatured state. The compact ensemble forms over a range of solvent and temperature conditions. The native fluorescence signature, which requires the tryptophan residue to be packed tightly, is acquired at least 500 times more slowly. Molecular dynamics simulations at 495 K show no contraction of the backbone in ethylene glycol buffer compared to pure aqueous buffer, and no significant effect on the local backbone structure of the unfolded protein. Only at higher simulation temperature does a backbone contraction appear. Thus, it appears unlikely that the aqueous ethylene glycol buffer fundamentally changes the folding mechanism of ubiquitin. We suggest that ubiquitin forms a compact ensemble with native-like secondary structure, but without tight packing, long before the native state.

Circular Dichroism↗
PubMed2004

Forced dissociation of the strand dimer interface between C-cadherin ectodomains.

The force-induced dissociation of the strand dimer interface in C-cadherin has been studied using steered molecular dynamics simulations. The dissociation occurred, without domain unraveling, after the extraction of the conserved trypthophans (Trp2) from their respective hydrophobic pockets. The simulations revealed two stable positions for the Trp2 side chain inside the pocket. The most internal stable position involved a hydrogen bond between the ring Nepsilon of Trp2 and the backbone carbonyl of Glu90. In the second stable position, the aromatic ring is located at the pocket entrance. After extracting the two tryptophans from their pockets, the complex exists in an intermediate bound state that involves a close packing of the tryptophans with residues Asp1 and Asp27 from both domains. Dissociation occurred after this residue association was broken. Simulations carried out with a complex formed between W2A mutants showed that the mutant complex dissociates more easily than the wild type complex does. These results correlate closely with the role of the conserved tryptophans suggested previously by site directed mutagenesis.

Amino Acids↗
PubMed2003

Forced detachment of the CD2-CD58 complex.

The force-induced detachment of the adhesion protein complex CD2-CD58 was studied by steered molecular dynamics simulations. The forced detachment of CD2 and CD58 shows that the system can respond to an external force by two mechanisms, which depend on the loading rate. At the rapid loading rates of 70 and 35 pN/ps (pulling speeds of 1 and 0.5 A/ps) the two proteins unfold before they separate, whereas at slower loading rates of 7 and 3.5 pN/ps (pulling speeds of 0.1 and 0.05 A/ps), the proteins separate before the domains can unfold. When subjected to a constant force of 400 pN, the two proteins separated without significant structural distortion. These findings suggest that protein unfolding is not coupled to the adhesive function of CD2 and CD58. The simulations further confirm that salt bridges primarily determine the tensile strength of the protein-to-protein bond, and that the order of salt bridge rupture depends mainly on the position of the bond, relative to the line of action of the applied force. Salt bridges close to this line break first. The importance of each of the salt bridges for adhesion, determined from the simulations, correlates closely with their role in cell-to-cell adhesion and equilibrium binding determined by site-directed mutagenesis experiments.

Binding Sites↗
PubMed2002

The role of intersection topography in bond selectivity of cis-trans photoisomerization.

Ab initio methods are used to characterize the ground and first excited state of the chromophore in the rhodopsin family of proteins: retinal protonated Schiff base. Retinal protonated Schiff base has five double bonds capable of undergoing isomerization. Upon absorption of light, the chromophore isomerizes and the character of the photoproducts (e.g., 13-cis and 11-cis) depends on the environment, protein vs. solution. Our ab initio calculations show that, in the absence of any specific interactions with the environment (e.g., discrete ordered charges in a protein), energetic considerations cannot explain the observed bond selectivity. We instead attribute the origin of bond selectivity to the shape (topography) of the potential energy surfaces in the vicinity of points of true degeneracy (conical intersections) between the ground and first excited electronic states. This provides a molecular example where a competition between two distinct but nearly isoenergetic photochemical reaction pathways is resolved by a topographical difference between two conical intersections.

Biophysical Phenomena↗
PubMed2001

Molecular dynamics study of aquaporin-1 water channel in a lipid bilayer.

The aquaporin-1 water channel was modeled in a palmitoyl-oleoyl-phosphatidyl-choline lipid bilayer, by means of molecular dynamics simulations. Interaction of the protein with the membrane and inter-monomer interactions were analyzed. Structural features of the channel important for its biological function, including the Asn-Pro-Ala (NPA) motifs, and the diffusion of water molecules into the channels, were investigated. Simulations revealed the formation of single file water inside the channels for certain relative positions of the NPA motifs.

Aquaporin 1↗
PubMed2001

Comparison of the early stages of forced unfolding for fibronectin type III modules.

The structural changes accompanying stretch-induced early unfolding events were investigated for the four type III fibronectin (FN-III) modules, FN-III(7), FN-III(8), FN-III(9), and FN-III(10) by using steered molecular dynamics. Simulations revealed that two main energy barriers, I and II, have to be overcome to initiate unraveling of FN-III's tertiary structure. In crossing the first barrier, the two opposing beta-sheets of FN-III are rotated against each other such that the beta-strands of both beta-sheets align parallel to the force vector (aligned state). All further events in the unfolding pathway proceed from this intermediate state. A second energy barrier has to be overcome to break the first major cluster of hydrogen bonds between adjacent beta-strands. Simulations revealed that the height of barrier I varied significantly among the four modules studied, being largest for FN-III(7) and lowest for FN-III(10), whereas the height of barrier II showed little variation. Key residues affecting the mechanical stability of FN-III modules were identified. These results suggest that FN-III modules can be prestretched into an intermediate state with only minor changes to their tertiary structures. FN-III(10), for example, extends 12 A from the native "twisted" to the intermediate aligned state, and an additional 10 A from the aligned state to further unfolding where the first beta-strand is peeled away. The implications of the existence of intermediate states regarding the elasticity of fibrillar fibers and the stretch-induced exposure of cryptic sites are discussed.

Amino Acid Sequence↗
PubMed2001

Phylogenetic analysis of metabolic pathways.

The information provided by completely sequenced genomes can yield insights into the multi-level organization of organisms and their evolution. At the lowest level of molecular organization individual enzymes are formed, often through assembly of multiple polypeptides. At a higher level, sets of enzymes group into metabolic networks. Much has been learned about the relationship of species from phylogenetic trees comparing individual enzymes. In this article we extend conventional phylogenetic analysis of individual enzymes in different organisms to the organisms' metabolic networks. For this purpose we suggest a method that combines sequence information with information about the underlying reaction networks. A distance between pathways is defined as incorporating distances between substrates and distances between corresponding enzymes. The new analysis is applied to electron-transfer and amino acid biosynthesis networks yielding a more comprehensive understanding of similarities and differences between organisms.

Animals↗
PubMed2001

Simulated refolding of stretched titin immunoglobulin domains.

Steered molecular dynamics (SMD) is used to investigate forced unfolding and spontaneous refolding of immunoglobulin I27, a domain of the muscle protein titin. Previous SMD simulations revealed the events leading to stretch-induced unfolding of I27, the rupture of hydrogen bonds bridging beta-strands A and B, and those bridging beta-strands A' and G, the latter rupture occurring at an extension of approximately 15 A and preceding the complete unfolding. Simulations are now used to study the refolding of partially unfolded I27 domains. The results reveal that stretched domains with ruptured interstrand hydrogen bonds shrink along the extension direction. Two types of refolding patterns are recognized: for separated beta-strands A' and G, in most simulations five of the six hydrogen bonds between A' and G stably reformed in 2 ns, whereas for separated beta-strands A and B hydrogen bonds seldom reformed in eight 2-ns simulations. The mechanical stability of the partially refolded intermediates has been tested by re-stretching.

Connectin↗
PubMed2001

Structural determinants of MscL gating studied by molecular dynamics simulations.

The mechanosensitive channel of large conductance (MscL) in prokaryotes plays a crucial role in exocytosis as well as in the response to osmotic downshock. The channel can be gated by tension in the membrane bilayer. The determination of functionally important residues in MscL, patch-clamp studies of pressure-conductance relationships, and the recently elucidated crystal structure of MscL from Mycobacterium tuberculosis have guided the search for the mechanism of MscL gating. Here, we present a molecular dynamics study of the MscL protein embedded in a fully hydrated POPC bilayer. Simulations totaling 3 ns in length were carried out under conditions of constant temperature and pressure using periodic boundary conditions and full electrostatics. The protein remained in the closed state corresponding to the crystal structure, as evidenced by its impermeability to water. Analysis of equilibrium fluctuations showed that the protein was least mobile in the narrowest part of the channel. The gating process was investigated through simulations of the bare protein under conditions of constant surface tension. Under a range of conditions, the transmembrane helices flattened as the pore widened. Implications for the gating mechanism in light of these and experimental results are discussed.

Computer Simulation↗
PubMed2001

Steered molecular dynamics and mechanical functions of proteins.

Atomic force microscopy of single molecules, steered molecular dynamics and the theory of stochastic processes have established a new field that investigates mechanical functions of proteins, such as ligand-receptor binding/unbinding and elasticity of muscle proteins during stretching. The combination of these methods yields information on the energy landscape that controls mechanical function and on the force-bearing components of proteins, as well as on the underlying physical mechanisms.

Avidin↗
PubMed2001

The mechanism of glycerol conduction in aquaglyceroporins.

BACKGROUND: The E. coli glycerol facilitator, GlpF, selectively conducts glycerol and water, excluding ions and charged solutes. The detailed mechanism of the glycerol conduction and its relationship to the characteristic secondary structure of aquaporins and to the NPA motifs in the center of the channel are unknown. RESULTS: Molecular dynamics simulations of GlpF reveal spontaneous glycerol and water conduction driven, on a nanosecond timescale, by thermal fluctuations. The bidirectional conduction, guided and facilitated by the secondary structure, is characterized by breakage and formation of hydrogen bonds for which water and glycerol compete. The conduction involves only very minor changes in the protein structure, and cooperativity between the GlpF monomers is not evident. The two conserved NPA motifs are strictly linked together by several stable hydrogen bonds and their asparagine side chains form hydrogen bonds with the substrates passing the channel in single file. CONCLUSIONS: A complete conduction of glycerol through the GlpF was deduced from molecular dynamics simulations, and key residues facilitating the conduction were identified. The nonhelical parts of the two half-membrane-spanning segments expose carbonyl groups towards the channel interior, establishing a curve-linear pathway. The conformational stability of the NPA motifs is important in the conduction and critical for selectivity. Water and glycerol compete in a random manner for hydrogen bonding sites in the protein, and their translocations in single file are correlated. The suggested conduction mechanism should apply to the whole family.

Amino Acid Motifs↗
PubMed2001

Steered molecular dynamics investigations of protein function.

Molecular recognition and mechanical properties of proteins govern molecular processes in the cell that can cause disease and can be targeted for drug design. Single molecule measurement techniques have greatly advanced knowledge but cannot resolve enough detail to be interpreted in terms of protein structure. We seek to complement the observations through so-called Steered Molecular Dynamics (SMD) simulations that link directly to experiments and provide atomic-level descriptions of the underlying events. Such a research program has been initiated in our group and has involved, for example, studies of elastic properties of immunoglobulin and fibronectin domains as well as the binding of biotin and avidin. In this article we explain the SMD method and suggest how it can be applied to the function of three systems that are the focus of modern molecular biology research: force transduction by the muscle protein titin and extracellular matrix protein fibronectin, recognition of antibody-antigene pairs, and ion selective conductivity of the K+ channel.

Antigen-Antibody Reactions↗