Membrane remodeling from N-BAR domain interactions: insights from multi-scale simulation

Liposome remodeling processes (e.g., vesiculation and tubulation) due to N-BAR domain interactions with the lipid bilayer are explored with a multi-scale simulation approach. Results from atomistic-level molecular dynamics simulations of membrane binding to the concave face of N-BAR domains are used along with discretized mesoscopic field-theoretic simulations to examine how the spontaneous curvature fields generated by N-BAR domains result in membrane remodeling. It is found that tubulation can be generated by anisotropic N-BAR spontaneous curvature fields, whereas vesiculation is only observed with isotropic N-BAR spontaneous curvature fields at high density. The results of the multi-scale simulations provide insight into recent experimental observations.     

The thermodynamics of protein-ligand interaction and solvation: insights for ligand design

Isothermal titration calorimetry is able to provide accurate information on the thermodynamic contributions of enthalpy and entropy changes to free energies of binding. The Structure/Calorimetry of Reported Protein Interactions Online database of published isothermal titration calorimetry studies and structural information on the interactions between proteins and small-molecule ligands is used here to reveal general thermodynamic properties of protein-ligand interactions and to investigate correlations with changes in solvation. The overwhelming majority of interactions are found to be enthalpically favoured. Synthetic inhibitors and biological ligands form two distinct subpopulations in the data, with the former having greater average affinity due to more favourable entropy changes on binding. The greatest correlation is found between the binding free energy and apolar surface burial upon complex formation. However, the free-energy contribution per unit area buried is only 30-50% of that expected from earlier studies of transfer free energies of small molecules. A simple probability-based estimator for the maximal affinity of a binding site in terms of its apolar surface area is proposed. Polar surface area burial also contributes substantially to affinity but is difficult to express in terms of unit area due to the small variation in the amount of polar surface buried and a tendency for cancellation of its enthalpic and entropic contributions. Conventionally, the contribution of apolar desolvation to affinity is attributed to gain of entropy due to solvent release. Although data presented here are supportive of this notion, because the correlation of entropy change with apolar surface burial is relatively weak, it cannot, on present evidence, be confidently considered to be correct. Further, thermodynamic changes arising from small differences between ligands binding to individual proteins are relatively large and, in general, uncorrelated with changes in solvation, suggesting that trends identified across widely differing proteins are of limited use in explaining or predicting the effects of ligand modifications.     

Dynamics of neuromuscular system recovery as a consequence of various massage techniques as determined from electroneuromyography parameters

By the method of stimulation electromyography, functional changes in the segmental system and skeletal muscles of the lower limbs of 16 persons (18-to 22-year-old men) were determined for 15 min after various massage techniques. Different massage techniques were clearly shown to differ in their effects on neuromuscular system parameters. Among the massage aftereffects were a rapid effect and a delayed influence of individual techniques.     

Dynamics driving function: new insights from electron transferring flavoproteins and partner complexes

Electron transferring flavoproteins (ETFs) are soluble heterodimeric FAD-containing proteins that function primarily as soluble electron carriers between various flavoprotein dehydrogenases. ETF is positioned at a key metabolic branch point, responsible for transferring electrons from up to 10 primary dehydrogenases to the membrane-bound respiratory chain. Clinical mutations of ETF result in the often fatal disease glutaric aciduria type II. Structural and biophysical studies of ETF in complex with partner proteins have shown that ETF partitions the functions of partner binding and electron transfer between (a) a 'recognition loop', which acts as a static anchor at the ETF-partner interface, and (b) a highly mobile redox-active FAD domain. Together, this enables the FAD domain of ETF to sample a range of conformations, some compatible with fast interprotein electron transfer. This 'conformational sampling' enables ETF to recognize structurally distinct partners, whilst also maintaining a degree of specificity. Complex formation triggers mobility of the FAD domain, an 'induced disorder' mechanism contrasting with the more generally accepted models of protein-protein interaction by induced fit mechanisms. We discuss the implications of the highly dynamic nature of ETFs in biological interprotein electron transfer. ETF complexes point to mechanisms of electron transfer in which 'dynamics drive function', a feature that is probably widespread in biology given the modular assembly and flexible nature of biological electron transfer systems.     

Structure and assembly of designed beta-hairpin peptides in crystals as models for beta-sheet aggregation

De novo designed beta-hairpin peptides have generally been recalcitrant to crystallization. The crystal structures of four synthetic peptide beta-hairpins, Boc-Leu-Val-Val-DPro-Gly-Leu-Phe-Val-OMe (1), Boc-Leu-Phe-Val-DPro-Ala-Leu-Phe-Val-OMe (2), Boc-Leu-Val-Val-DPro-Aib-Leu-Val-Val-OMe (3), and Boc-Met-Leu-Phe-Val-DPro-Ala-Leu-Val-Val-Phe-OMe (4), are described. The centrally positioned DPro-Xxx segment promotes prime beta-turn formation, thereby nucleating beta-hairpin structures. In all four peptides well-defined beta-hairpins nucleated by central type II' DPro-Xxx beta-turns have been characterized by X-ray diffraction, providing a view of eight crystallographically independent hairpins. In peptides 1-3 three intramolecular cross-strand hydrogen bonds stabilized the observed beta-hairpin, with some fraying of the structures at the termini. In peptide 4, four intramolecular cross-strand hydrogen bonds stabilized the hairpin. Peptides 1-4 reveal common features of packing of beta-hairpins into crystals. Two-dimensional sheet formation mediated by intermolecular hydrogen bonds formed between antiparallel strands of adjacent molecule is a recurrent theme. The packing of two-dimensional sheets into the crystals is mediated in the third dimension by bridging solvents and interactions of projecting side chains, which are oriented on either face of the sheet. In all cases, solvation of the central DPro-Xxx peptide unit beta-turn is observed. The hairpins formed in the octapeptides are significantly buckled as compared to the larger hairpin in peptide 4, which is much flatter. The crystal structures provide insights into the possible modes of beta-sheet packing in regular crystalline arrays, which may provide a starting point for understanding beta-sandwich and cross-beta-structures in amyloid fibrils.     

How dopamine transporter interacts with dopamine: insights from molecular modeling and simulation

By performing homology modeling, molecular docking, and molecular dynamics simulations, we have developed three-dimensional (3D) structural models of both dopamine transporter and dopamine transporter-dopamine complex in the environment of lipid bilayer and solvent water. According to the simulated structure of dopamine transporter-dopamine complex, dopamine was orientated in a hydrophobic pocket at the midpoint of the membrane. The modeled 3D structures provide some detailed structural and mechanistic insights concerning how dopamine transporter (DAT) interacts with dopamine at atomic level, extending our mechanistic understanding of the dopamine reuptake with the help of Na(+) ions. The general features of the modeled 3D structures are consistent with available experimental data. Based on the modeled structures, our calculated binding free energy (DeltaG(bind) = -6.4 kcal/mol) for dopamine binding with DAT is also reasonably close to the experimentally derived DeltaG(bind) value of -7.4 kcal/mol. Finally, a possible dopamine-entry pathway, which involves formation and breaking of the salt bridge between side chains of Arg(85) and Asp(476), is proposed based on the results obtained from the modeling and molecular dynamics simulation. The new structural and mechanistic insights obtained from this computational study are expected to stimulate future, further biochemical and pharmacological studies on the detailed structures and mechanisms of DAT and other homologous transporters.     

Electrostatic contributions to the kinetics and thermodynamics of protein assembly

The role of electrostatic interactions in the assembly of a native protein structure was studied using fragment complementation. Contributions of salt, pH, or surface charges to the kinetics and equilibrium of calbindin D(9k) reconstitution was measured in the presence of Ca(2+) using surface plasmon resonance and isothermal titration calorimetry. Whereas surface charge substitutions primarily affect the dissociation rate constant, the association rates are correlated with subdomain net charge in a way expected for Coulomb interactions. The affinity is reduced in all mutants, with the largest effect (260-fold) observed for the double mutant K25E+K29E. At low net charge, detailed charge distribution is important, and charges remote from the partner EF-hand have less influence than close ones. The effects of salt and pH on the reconstitution are smaller than mutational effects. The interaction between the wild-type EF-hands occurs with high affinity (K(A) = 1.3 x 10(10) M(-1); K(D) = 80 pM). The enthalpy of association is overall favorable and there appears to be a very large favorable entropic contribution from the desolvation of hydrophobic surfaces that become buried in the complex. Electrostatic interactions contribute significantly to the affinity between the subdomains, but other factors, such as hydrophobic interactions, dominate.     

beta-hairpin stability and folding: molecular dynamics studies of the first beta-hairpin of tendamistat

The stability and (un)folding of the 19-residue peptide, SCVTLYQSWRYSQADNGCA, corresponding to the first beta-hairpin (residues 10 to 28) of the alpha-amylase inhibitor tendamistat (PDB entry 3AIT) has been studied by molecular dynamics simulations in explicit water under periodic boundary conditions at several temperatures (300 K, 360 K and 400 K), starting from various conformations for simulation lengths, ranging from 10 to 30 ns. Comparison of trajectories of the reduced and oxidized native peptides reveals the importance of the disulphide bridge closing the beta-hairpin in maintaining a proper turn conformation, thereby insuring a proper side-chain arrangement of the conserved turn residues. This allows rationalization of the conservation of those cysteine residues among the family of alpha-amylase inhibitors. High temperature simulations starting from widely different initial configurations (native beta-hairpin, alpha and left-handed helical and extended conformations) begin sampling similar regions of the conformational space within tens of nanoseconds, and both native and non-native beta-hairpin conformations are recovered. Transitions between conformational clusters are accompanied by an increase in energy fluctuations, which is consistent with the increase in heat capacity measured experimentally upon protein folding. The folding events observed in the various simulations support a model for beta-hairpin formation in which the turn is formed first, followed by hydrogen bond formation closing the hairpin, and subsequent stabilization by side-chain hydrophobic interactions.     

The mechanism of HIV-1 core assembly: insights from three-dimensional reconstructions of authentic virions

Infectious HIV particles contain a characteristic cone-shaped core encasing the viral RNA and replication proteins. The core exhibits significant heterogeneity in size and shape, yet consistently forms a well-defined structure. The mechanism by which the core is assembled in the maturing virion remains poorly understood. Using cryo-electron tomography, we have produced three-dimensional reconstructions of authentic, unstained HIV-1. These reveal the viral morphology with unprecedented clarity and suggest the following mechanism for core formation inside the extracellular virion: core growth initiates at the narrow end of the cone and proceeds toward the distal side of the virion until limited by the viral membrane. Curvature and closure of the broad end of the core are then directed by the inner surface of the viral membrane. This mechanism accommodates significant flexibility in lattice growth while ensuring the closure of cores of variable size and shape.     

Transmembrane signaling of chemotaxis receptor tar: insights from molecular dynamics simulation studies

Transmembrane signaling of chemotaxis receptors has long been studied, but how the conformational change induced by ligand binding is transmitted across the bilayer membrane is still elusive at the molecular level. To tackle this problem, we carried out a total of 600-ns comparative molecular dynamics simulations (including model-building simulations) of the chemotaxis aspartate receptor Tar (a part of the periplasmic domain/transmembrane domain/HAMP domain) in explicit lipid bilayers. These simulations reveal valuable insights into the mechanistic picture of Tar transmembrane signaling. The piston-like movement of a transmembrane helix induced by ligand binding on the periplasmic side is transformed into a combination of both longitudinal and transversal movements of the helix on the cytoplasmic side as a result of different protein-lipid interactions in the ligand-off and ligand-on states of the receptor. This conformational change alters the dynamics and conformation of the HAMP domain, which is presumably a mechanism to deliver the signal from the transmembrane domain to the cytoplasmic domain. The current results are consistent with the previously suggested dynamic bundle model in which the HAMP dynamics change is a key to the signaling. The simulations provide further insights into the conformational changes relevant to the HAMP dynamics changes in atomic detail.     

Microtubule assembly dynamics: new insights at the nanoscale

Although the dynamic self-assembly behavior of microtubule ends has been well characterized at the spatial resolution of light microscopy (~200 nm), the single-molecule events that lead to these dynamics are less clear. Recently, a number of in vitro studies used novel approaches combining laser tweezers, microfabricated chambers, and high-resolution tracking of microtubule-bound beads to characterize mechanochemical aspects of MT dynamics at nanometer scale resolution. In addition, computational modeling is providing a framework for integrating these experimental results into physically plausible models of molecular scale microtubule dynamics. These nanoscale studies are providing new fundamental insights about microtubule assembly, and will be important for advancing our understanding of how microtubule dynamic instability is regulated in vivo via microtubule-associated proteins, therapeutic agents, and mechanical forces.     

Protein structure and dynamics in nonaqueous solvents: insights from molecular dynamics simulation studies

Protein structure and dynamics in nonaqueous solvents are here investigated using molecular dynamics simulation studies, by considering two model proteins (ubiquitin and cutinase) in hexane, under varying hydration conditions. Ionization of the protein groups is treated assuming "pH memory," i.e., using the ionization states characteristic of aqueous solution. Neutralization of charged groups by counterions is done by considering a counterion for each charged group that cannot be made neutral by establishing a salt bridge with another charged group; this treatment is more physically reasonable for the nonaqueous situation, contrasting with the usual procedures. Our studies show that hydration has a profound effect on protein stability and flexibility in nonaqueous solvents. The structure becomes more nativelike with increasing values of hydration, up to a certain point, when further increases render it unstable and unfolding starts to occur. There is an optimal amount of water, approximately 10% (w/w), where the protein structure and flexibility are closer to the ones found in aqueous solution. This behavior can explain the experimentally known bell-shaped dependence of enzyme catalysis on hydration, and the molecular reasons for it are examined here. Water and counterions play a fundamental and dynamic role on protein stabilization, but they also seem to be important for protein unfolding at high percentages of bound water.     

Actin in striated muscle: recent insights into assembly and maintenance

Striated muscle cells are characterised by a para-crystalline arrangement of their contractile proteins actin and myosin in sarcomeres, the basic unit of the myofibrils. A multitude of proteins is required to build and maintain the structure of this regular arrangement as well as to ensure regulation of contraction and to respond to alterations in demand. This review focuses on the actin filaments (also called thin filaments) of the sarcomere and will discuss how they are assembled during myofibrillogenesis and in hypertrophy and how their integrity is maintained in the working myocardium.     

Suppression mechanisms of COX assembly defects in yeast and human: insights into the COX assembly process

Eukaryotic cytochrome c oxidase (COX) is the terminal enzyme of the mitochondrial respiratory chain. COX is a multimeric enzyme formed by subunits of dual genetic origin whose assembly is intricate and highly regulated. In addition to the structural subunits, a large number of accessory factors are required to build the holoenzyme. The function of these factors is required in all stages of the assembly process. They are relevant to human health because devastating human disorders have been associated with mutations in nuclear genes encoding conserved COX assembly factors. The study of yeast strains and human cell lines from patients carrying mutations in structural subunits and COX assembly factors has been invaluable to attain the current state of knowledge, even if still fragmentary, of the COX assembly process. After the identification of the genes involved, the isolation and characterization of genetic and metabolic suppressors of COX assembly defects, reviewed here, have become a profitable strategy to gain insight into their functions and the pathways in which they operate. Additionally, they have the potential to provide useful information for devising therapeutic approaches to combat human disorders associated with COX deficiency.     

Redox-coupled proton pumping in cytochrome c oxidase: further insights from computer simulation

The membrane-bound enzyme cytochrome c oxidase, the terminal member in the respiratory chain, converts oxygen into water and generates an electrochemical gradient by coupling the electron transfer to proton pumping across the membrane. Here we have investigated the dynamics of an excess proton and the surrounding protein environment near the active sites. The multi-state empirical valence bond (MS-EVB) molecular dynamics method was used to simulate the explicit dynamics of proton transfer through the critically important hydrophobic channel between Glu242 (bovine notation) and the D-propionate of heme a3 (PRDa3) for the first time. The results from these molecular dynamics simulations indicate that the PRDa3 can indeed re-orientate and dissociate from Arg438, despite the high stability of such an ion pair, and has the ability to accept protons via bound water molecules. Any large conformational change of the adjacent heme a D-propionate group is, however, sterically blocked directly by the protein. Free energy calculations of the PRDa3 side chain isomerization and the proton translocation between Glu242 and the PRDa3 site have also been performed. The results exhibit a redox state-dependent dynamical behavior and indicate that reduction of the low-spin heme a may initiate internal transfer of the pumped proton from Glu242 to the PRDa3 site.     

New insights into receptor ligand binding domains from a novel assembly assay

Previous studies have demonstrated that hormone binding stabilizes the ligand binding domain (LBD) of the nuclear hormone receptors against proteolysis. We have confirmed and extended this observation using a newly developed assembly assay. In this assay, the LBD is divided into two parts, of which one includes the first helix of this domain and the other corresponds to the remainder of the LBD. Several independent criteria demonstrate that these two fragments can assemble into a functional LBD in the presence of a ligand, but not in its absence, and that this is a reflection of the stabilizing effect of ligand. We have also used this assay to demonstrate that binding of the nuclear receptor corepressor NCoR can directly stabilize the LBD. Overall, these results highlight the dynamic nature of the LBD and suggest that current models for activation based solely on allosteric effects on the C-terminal helix may be too limited.     

Toward quantitative simulation of germinal center dynamics: biological and modeling insights from experimental validation.

As models of immune system dynamics are developed, it is important to validate them with specific experimental data in order to understand their shortcomings and guide them toward becoming predictive. In this paper, we examine whether a particular mathematical model of germinal center dynamics, proposed by Oprea and Perelson, can reproduce experimental data from two specific primary responses, namely those directed against the haptens 2-phenyl-5-oxazolone and (4-hydroxy-3-nitrophenyl)acetyl. We develop formulas for estimating response-specific model parameters, as well as constraints for validating the model. In addition, we outline a general methodology for translating a continuous/deterministic model, expressed as a set of ordinary differential equations, into a discrete/stochastic framework. This methodology is used to create a new implementation of the Oprea and Perelson model that enables comparison with data on individual germinal centers. We conclude that while the model can reproduce the average dynamics of splenic germinal centers, it is at best incomplete and does not reproduce the distribution of individual germinal center behaviors. In addition to suggesting possible extensions to the model which can reconcile the dynamics with some aspects of the experimental data, we make a number of specific predictions that can be tested by in vivo experiments to obtain further insights and validation.