ALBERT

Description: We study four citrate synthase homodimeric proteins within a structure-based coarse-grained model. Two of these proteins come from thermophilic bacteria, one from a cryophilic bacterium and one from a mesophilic organism; three are in the closed and two in the open conformations. Even though the proteins belong to the same fold, the model distinguishes the properties of these proteins in a way which is consistent with experiments. For instance, the thermophilic proteins are more stable thermodynamically than their mesophilic and cryophilic homologues, which we observe both in the magnitude of thermal fluctuations near the native state and in the kinetics of thermal unfolding. The level of stability correlates with the average coordination number for amino acid contacts and with the degree of structural compactness. The pattern of positional fluctuations along the sequence in the closed conformation is different than in the open conformation, including within the active site. The modes of correlated and anticorrelated movements of pairs of amino acids forming the active site are very different in the open and closed conformations. Taken together, our results show that the precise location of amino acid contacts in the native structure appears to be a critical element in explaining the similarities and differences in the thermodynamic properties, local flexibility, and collective motions of the different forms of the enzyme.

Description: In this paper, we examine the behaviour of basic autocatalytic feedback modules involving a species catalyzing its own production, either directly or indirectly. We first perform a systematic study of the autocatalytic feedback module in isolation, examining the effect of different factors, showing how this module is capable of exhibiting monostable threshold and bistable switch-like behaviour. We then study the behaviour of this module embedded in different kinds of basic networks including (essentially) irreversible cycles, open and closed reversible chains, and networks with additional feedback. We study the behaviour of the networks deterministically and also stochastically, using simulations, analytical work, and bifurcation analysis. We find that (i) there are significant differences between the behaviour of this module in isolation and in a network: thresholds may be altered or destroyed and bistability may be destroyed or even induced, even when the ambient network is simple. The global characteristics and topology of this network and the position of the module in the ambient network can play important and unexpected roles. (ii) There can be important differences between the deterministic and stochastic dynamics of the module embedded in networks, which may be accentuated by the ambient network. This provides new insights into the functioning of such enzymatic modules individually and as part of networks, with relevance to other enzymatic signalling modules as well.

Description: Recent experiments indicate a connection between the structure of amyloid aggregates and their cytotoxicity as related to neurodegenerative diseases. Of particular interest is the Iowa Mutant, which causes early-onset of Alzheimer's disease. While wild-type Amyloid β -peptides form only parallel beta-sheet aggregates, the mutant also forms meta-stable antiparallel beta sheets. Since these structural variations may cause the difference in the pathological effects of the two A β -peptides, we have studied in silico the relative stability of the wild type and Iowa mutant in both parallel and antiparallel forms. We compare regular molecular dynamics simulations with such where the viscosity of the samples is reduced, which, we show, leads to higher sampling efficiency. By analyzing and comparing these four sets of all-atom molecular dynamics simulations, we probe the role of the various factors that could lead to the structural differences. Our analysis indicates that the parallel forms of both wild type and Iowa mutant aggregates are stable, while the antiparallel aggregates are meta-stable for the Iowa mutant and not stable for the wild type. The differences result from the direct alignment of hydrophobic interactions in the in-register parallel oligomers, making them more stable than the antiparallel aggregates. The slightly higher thermodynamic stability of the Iowa mutant fibril-like oligomers in its parallel organization over that in antiparallel form is supported by previous experimental measurements showing slow inter-conversion of antiparallel aggregates into parallel ones. Knowledge of the mechanism that selects between parallel and antiparallel conformations and determines their relative stability may open new avenues for the development of therapies targeting familial forms of early-onset Alzheimer's disease.

Description: Hemoprotein mimics, cobalt picket fence porphyrins have been prepared in the gas phase as neutral molecules for the first time. Their ligation properties have been studied with 1-methylimidazole and compared with those of other cobalt porphyrins, tetraphenyl porphyrin, and cobalt protoporphyrin IX chloride, in view of studying the sterical properties of the ligation. It is shown that the cobalt picket fence porphyrin can only accept one 1-methylimidazole ligand in contrast to less sterically crowded porphyrins like cobalt tetraphenylporphyrin that present two accessible ligation sites. The femtosecond dynamics of these ligated systems have been studied after excitation at 400 nm, in comparison with the unligated ones. The observed transients are formed in much shorter times, 30 fs for the ligated species, as compared to free species (100 fs), supporting the porphyrin to metal charge transfer nature of these transients. The similar decays of the ligated transients

Description: DNA-binding proteins locate and bind their target sequences positioned on DNA in crowded environments, but the molecular crowding effect on this search process is not clear. Using analytical techniques and Langevin dynamics simulations in two dimensions (2D), we find that the essential physics for facilitated diffusion in 2D search and 3D search is the same. We observe that the average search times have minima at the same optimal nonspecific binding energy for the cases with and without the crowding particle. Moreover, the molecular crowding increases the search time by increasing the average search rounds and the one-dimensional (1D) sliding time of a round, but almost not changing the average 2D diffusion time of a round. In addition, the fraction of 1D sliding time out of the total search time increases with increasing the concentration of crowders. For 2D diffusion, the molecular crowding decreases the jumping length and narrows its distribution due to the cage effect from crowders. These results shed light on the role of facilitated diffusion in DNA targeting kinetics in living cells.

Description: Protein conformation and orientation in the lipid membrane plays a key role in many cellular processes. Here we use molecular dynamics simulation to investigate the relaxation and C-terminus diffusion of a model helical peptide: beta-amyloid (Aβ) in a lipid membrane. We observed that after the helical peptide was initially half-embedded in the extracelluar leaflet of phosphatidylcholine (PC) or PC/cholesterol (PC/CHOL) membrane, the C-terminus diffused across the membrane and anchored to PC headgroups of the cytofacial lipid leaflet. In some cases, the membrane insertion domain of the Aβ was observed to partially unfold. Applying a sigmoidal fit to the process, we found that the characteristic velocity of the C-terminus, as it moved to its anchor site, scaled with θ u −4/3 , where θ u is the fraction of the original helix that was lost during a helix to coil transition. Comparing this scaling with that of bead-spring models of polymer relaxation suggests that the C-terminus velocity is highly regulated by the peptide helical content, but that it is independent of the amino acid type. The Aβ was stabilized by the attachment of the positive Lys28 side chain to the negative phosphate of PC or 3β oxygen of CHOL in the extracellular lipid leaflet and of the C-terminus to its anchor site in the cytofacial lipid leaflet.

Description: We present a method to obtain numerically accurate values of configurational free energies of semiflexible macromolecular systems, based on the technique of thermodynamic integration combined with normal-mode analysis of a reference system subject to harmonic constraints. Compared with previous free-energy calculations that depend on a reference state, our approach introduces two innovations, namely, the use of internal coordinates to constrain the reference states and the ability to freely select these reference states. As a consequence, it is possible to explore systems that undergo substantially larger fluctuations than those considered in previous calculations, including semiflexible biopolymers having arbitrary ratios of contour length L to persistence length P . To validate the method, high accuracy is demonstrated for free energies of prime DNA knots with L / P = 20 and L / P = 40, corresponding to DNA lengths of 3000 and 6000 base pairs, respectively. We then apply the method to study the free-energy landscape for a model of a synaptic nucleoprotein complex containing a pair of looped domains, revealing a bifurcation in the location of optimal synapse (crossover) sites. This transition is relevant to target-site selection by DNA-binding proteins that occupy multiple DNA sites separated by large linear distances along the genome, a problem that arises naturally in gene regulation, DNA recombination, and the action of type-II topoisomerases.

Description: Since the time of Kirkwood, observed deviations in magnitude of the dielectric constant of aqueous protein solution from that of neat water (∼80) and slower decay of polarization have been subjects of enormous interest, controversy, and debate. Most of the common proteins have large permanent dipole moments (often more than 100 D) that can influence structure and dynamics of even distant water molecules, thereby affecting collective polarization fluctuation of the solution, which in turn can significantly alter solution's dielectric constant. Therefore, distance dependence of polarization fluctuation can provide important insight into the nature of biological water. We explore these aspects by studying aqueous solutions of four different proteins of different characteristics and varying sizes, chicken villin headpiece subdomain (HP-36), immunoglobulin binding domain protein G (GB1), hen-egg white lysozyme (LYS), and Myoglobin (MYO). We simulate fairly large systems consisting of single protein molecule and 20000–30000 water molecules (varied according to the protein size), providing a concentration in the range of ∼2–3 mM. We find that the calculated dielectric constant of the system shows a noticeable increment in all the cases compared to that of neat water. Total dipole moment auto time correlation function of water ⟨δ M W (0)δ M W ( t )⟩ is found to be sensitive to the nature of the protein. Surprisingly, dipole moment of the protein and total dipole moment of the water molecules are found to be only weakly coupled. Shellwise decomposition of water molecules around protein reveals higher density of first layer compared to the succeeding ones. We also calculate heuristic effective dielectric constant of successive layers and find that the layer adjacent to protein has much lower value (∼50). However, progressive layers exhibit successive increment of dielectric constant, finally reaching a value close to that of bulk 4–5 layers away. We also calculate shellwise orientational correlation function and tetrahedral order parameter to understand the local dynamics and structural re-arrangement of water. Theoretical analysis providing simple method for calculation of shellwise local dielectric constant and implication of these findings are elaborately discussed in the present work.

Description: Hydrogenase enzymes are important because they can reversibly catalyze the production of molecular hydrogen. Proton transport mechanisms have been previously studied in residue pathways that lead to the active site of the enzyme via residues Cys299 and Ser319. The importance of this pathway and these residues has been previously exhibited through site-specific mutations, which were shown to interrupt the enzyme activity. It has been shown recently that a separate water channel (WC2) is coupled with electron transport to the active site of the [FeFe]-hydrogenase. The water-mediated proton transport mechanisms of the enzyme in different electronic states have been studied using the multistate empirical valence bond reactive molecular dynamics method, in order to understand any role WC2 may have in facilitating the residue pathway in bringing an additional proton to the enzyme active site. In a single electronic state A 2− , a water wire was formed through which protons can be transported with a low free energy barrier. The remaining electronic states were shown, however, to be highly unfavorable to proton transport in WC2. A double amino acid substitution is predicted to obstruct proton transport in electronic state A 2- by closing a cavity that could otherwise fill with water near the proximal Fe of the active site.

Description: We propose a discrete transition-based reweighting analysis method (dTRAM) for analyzing configuration-space-discretized simulation trajectories produced at different thermodynamic states (temperatures, Hamiltonians, etc.) dTRAM provides maximum-likelihood estimates of stationary quantities (probabilities, free energies, expectation values) at any thermodynamic state. In contrast to the weighted histogram analysis method (WHAM), dTRAM does not require data to be sampled from global equilibrium, and can thus produce superior estimates for enhanced sampling data such as parallel/simulated tempering, replica exchange, umbrella sampling, or metadynamics. In addition, dTRAM provides optimal estimates of Markov state models (MSMs) from the discretized state-space trajectories at all thermodynamic states. Under suitable conditions, these MSMs can be used to calculate kinetic quantities (e.g., rates, timescales). In the limit of a single thermodynamic state, dTRAM estimates a maximum likelihood reversible MSM, while in the limit of uncorrelated sampling data, dTRAM is identical to WHAM. dTRAM is thus a generalization to both estimators.