Molecular dynamics simulations of class C beta-lactamase from Citrobacter freundii: insights into the base catalyst for acylation

Herein, we present results from molecular dynamics (MD) simulations of the class C beta-lactamase from Citrobacter freundii and its Michaelis complex with aztreonam. Four different configurations of the active site were modeled in aqueous solution, and their relative stability was estimated by means of quantum mechanical energy calculations. For the free enzyme, the energetically most stable configurations present a neutral Lys67 residue or an anionic Tyr150 side chain. Our calculations predict that these two configurations are quite close in terms of free energy, the anionic Tyr150 state being favored by approximately 1 kcal/mol. In contrast, for the noncovalent complex formed between the C. freundii enzyme and aztreonam, the energetic analyses predict that the configuration with the neutral Lys67 residue is much more stable than the anionic Tyr150 one (approximately 20 kcal/mol). Moreover, the MD simulations reveal that the neutral Lys67 state results in a proper enzyme-aztreonam orientation for nucleophilic attack and in a very stable contact between the nucleophilic hydroxyl group of Ser64 and the neutral amino side chain of Lys67. Thus, both the computed free energies and the structural analyses support the assignation of Lys67 as the base catalyst for the acylation step in the native form of the C. freundii enzyme.     

Mutant bovine odorant-binding protein: Temperature affects the protein stability and dynamics as revealed by infrared spectroscopy and molecular dynamics simulations

Bovine odorant-binding protein (bOBP), a member of the lipocalin family, presents the so-called 3D "domain-swapped" protein structure. In fact, in solution, it appears as a dimer in which each monomer is composed by the classical lipocalin fold, with a central beta-barrel followed by a stretch of residues and the alpha-helix domain protruding out of the barrel and crossing the dimer interface. Recently, a deswapped mutant form of bOBP was obtained, in which a Gly residue was inserted after position 121 and the two residues in position 64 and 156 were replaced by Cys residues for restoring the disulfide bridge common to the lipocalin family. In this work, we used Fourier transform infrared spectroscopy and molecular dynamics simulations to investigate the effect of temperature on the structural stability and conformational dynamics of the mutant bOBP. The spectroscopic and molecular simulation data pointed out that the hydrophobic regions of the protein matrix appear to be an important factor for the protein stability and integrity. In addition, it was also found that the mutant bOBP is significantly stabilized by the binding of the ligand, which may have an impact on the biological function of bOBP. The obtained results will allow for a better use of this protein as probe for the design of advanced protein-based biosensors for the detection of compounds used in the fabrication of explosive powders.     

Static and dynamic structure of Pu3+ in molten LiCl-KCl eutectic revealed by first-principles molecular dynamics simulations

The local coordination structure, as well as its structural dynamics, of molten LiCl-KCl-PuCl₃ mixture is revealed by first-principles molecular dynamics (FPMD) simulations in terms of radial distribution functions, partial structural factors, and van Hove correlation functions. The self-diffusion coefficients and ionic conductivities are evaluated from the mean square displacements of ions, and the viscosity is then estimated from the self-diffusion coefficient of Pu³⁺. It is concluded that the FPMD simulations give satisfactory results compared with existing literature values. Therefore, the FPMD simulations can provide basic physicochemical information if such information is unavailable in literature and microscopic insight into the mechanism of the pyroprocessing of spent nuclear fuels.     

Excited-state dynamics of protochlorophyllide revealed by subpicosecond infrared spectroscopy

To gain a better understanding of the light-induced reduction of protochlorophyllide (PChlide) to chlorophyllide as a key regulatory step in chlorophyll synthesis, we performed transient infrared absorption measurements on PChlide in d4-methanol. Excitation in the Q-band at 630 nm initiates dynamics characterized by three time constants: τ₁ = 3.6 ± 0.2, τ₂ = 38 ± 2, and τ₃ = 215 ± 8 ps. As indicated by the C13′=O carbonyl stretching mode in the electronic ground state at 1686 cm⁻¹, showing partial ground-state recovery, and in the excited electronic state at 1625 cm⁻¹, showing excited-state decay, τ₂ describes the formation of a state with a strong change in electronic structure, and τ₃ represents the partial recovery of the PChlide electronic ground state. Furthermore, τ₁ corresponds with vibrational energy relaxation. The observed kinetics strongly suggest a branched reaction scheme with a branching ratio of 0.5 for the path leading to the PChlide ground state on the 200 ps timescale and the path leading to a long-lived state (>>700 ps). The results clearly support a branched reaction scheme, as proposed previously, featuring the formation of an intramolecular charge transfer state with ∼25 ps, its decay into the PChlide ground state with 200 ps, and a parallel reaction path to the long-lived PChlide triplet state.     

A simplified model of local structure in aqueous proline amino acid revealed by first-principles molecular dynamics simulations

Aqueous proline solutions are deceptively simple as they can take on complex roles such as protein chaperones, cryoprotectants, and hydrotropic agents in biological processes. Here, a molecular level picture of proline/water mixtures is developed. Car-Parrinello ab initio molecular dynamics (CPAIMD) simulations of aqueous proline amino acid at the B-LYP level of theory, performed using IBM's Blue Gene/L supercomputer and massively parallel software, reveal hydrogen-bonding propensities that are at odds with the predictions of the CHARMM22 empirical force field but are in better agreement with results of recent neutron diffraction experiments. In general, the CPAIMD (B-LYP) simulations predict a simplified structural model of proline/water mixtures consisting of fewer distinct local motifs. Comparisons of simulation results to experiment are made by direct evaluation of the neutron static structure factor S(Q) from CPAIMD (B-LYP) trajectories as well as to the results of the empirical potential structure refinement reverse Monte Carlo procedure applied to the neutron data.     

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.     

Solution structure of the Lewis x oligosaccharide determined by NMR spectroscopy and molecular dynamics simulations.

The Lewis x (Lex) determinant is a trisaccharide fragment that has been implicated as a specific differentiation antigen, as a tumor antigen, and as a key component of the ligand for the endothelial leukocyte adhesion molecule, ELAM-1. High-resolution nuclear magnetic resonance spectroscopy shows it to have a relatively rigid structure. Only a small range of glycosidic dihedral angles in the trisaccharide produce simulated nuclear Overhauser effect spectra agreeing with data measured for the human milk pentasaccharide, lacto-N-fucopentaose-3, which contains the Lex determinant. Independently, the same average structure for the Lex determinant arises from in vacuo molecular dynamics simulations. The proposed conformation of the Lex trisaccharide is very similar to that recently determined for the closely related Lea trisaccharide. In agreement with the recent finding that both sialylated Lea and Lex react with ELAM-1, the results presented here show that the Lea and Lex determinants contain very similar carbohydrate domains.     

Beta-hairpin conformation of fibrillogenic peptides: structure and alpha-beta transition mechanism revealed by molecular dynamics simulations

Understanding the conformational transitions that trigger the aggregation and amyloidogenesis of otherwise soluble peptides at atomic resolution is of fundamental relevance for the design of effective therapeutic agents against amyloid-related disorders. In the present study the transition from ideal alpha-helical to beta-hairpin conformations is revealed by long timescale molecular dynamics simulations in explicit water solvent, for two well-known amyloidogenic peptides: the H1 peptide from prion protein and the Abeta(12-28) fragment from the Abeta(1-42) peptide responsible for Alzheimer's disease. The simulations highlight the unfolding of alpha-helices, followed by the formation of bent conformations and a final convergence to ordered in register beta-hairpin conformations. The beta-hairpins observed, despite different sequences, exhibit a common dynamic behavior and the presence of a peculiar pattern of the hydrophobic side-chains, in particular in the region of the turns. These observations hint at a possible common aggregation mechanism for the onset of different amyloid diseases and a common mechanism in the transition to the beta-hairpin structures. Furthermore the simulations presented herein evidence the stabilization of the alpha-helical conformations induced by the presence of an organic fluorinated cosolvent. The results of MD simulation in 2,2,2-trifluoroethanol (TFE)/water mixture provide further evidence that the peptide coating effect of TFE molecules is responsible for the stabilization of the soluble helical conformation.     

Studies on somatostatin with time-resolved spectroscopy and molecular dynamics simulations

Somatostatin was studied by time-resolved fluorescence spectroscopy and by molecular dynamics simulations. The results obtained indicate the existence of mainly one conformation in DMSO, and the existence of several conformations in other solvents. Molecular dynamics simulations showed that the most important region for activity (residue 7-9) is also the most flexible region in the peptide.     

Structure and interactions of RecA: plasticity revealed by molecular dynamics simulations

Eleven independent simulations, each involving three consecutive molecules in the RecA filament, carried out on the protein from Mycobacterium tuberculosis, Mycobacterium smegmatis and Escherichia coli and their Adenosine triphosphate (ATP) complexes, provide valuable information which is complementary to that obtained from crystal structures, in addition to confirming the robust common structural framework within which RecA molecules from different eubacteria function. Functionally important loops, which are largely disordered in crystal structures, appear to adopt in each simulation subsets of conformations from larger ensembles. The simulations indicate the possibility of additional interactions involving the P-loop which remains largely invariant. The phosphate tail of the ATP is firmly anchored on the loop while the nucleoside moiety exhibits substantial structural variability. The most important consequence of ATP binding is the movement of the ‘switch’ residue. The relevant simulations indicate the feasibility of a second nucleotide binding site, but the pathway between adjacent molecules in the filament involving the two nucleotide binding sites appears to be possible only in the mycobacterial proteins.     

A noncanonical binding site of chloramphenicol revealed via molecular dynamics simulations

Chloramphenicol, an antibiotic belonging to the family of amphenicols, is an inhibitor of translation. On the basis of X–ray structural analysis of the binding of chloramphenicol to free bacterial ribosomes, the chloramphenicol action mechanism that consists in preventing the binding of aminoacyl-tRNA to the A–site of the large subunit of the ribosome was adopted. However, the known structures of chloramphenicol complexes with bacterial ribosomes poorly explain the results of the experiments on the chemical modification of 23S rRNA, the resistance to chloramphenicol caused by mutations in 23S rRNA and, which is particularly important, the selectivity of chloramphenicol in suppression of translation, depending on the amino acid sequence of the nascent peptide. In the present study the putative structure of the chloramphenicol complex with a bacterial ribosome in the A,A/P,P–state has been obtained by molecular dynamics simulations methods. The proposed structure of the complex allows us to explain the results of biochemical studies of the interaction of chloramphenicol with the bacterial ribosome.     

Structured pathway across the transition state for peptide folding revealed by molecular dynamics simulations

Small globular proteins and peptides commonly exhibit two-state folding kinetics in which the rate limiting step of folding is the surmounting of a single free energy barrier at the transition state (TS) separating the folded and the unfolded states. An intriguing question is whether the polypeptide chain reaches, and leaves, the TS by completely random fluctuations, or whether there is a directed, stepwise process. Here, the folding TS of a 15-residue β-hairpin peptide, Peptide 1, is characterized using independent 2.5 μs-long unbiased atomistic molecular dynamics (MD) simulations (a total of 15 μs). The trajectories were started from fully unfolded structures. Multiple (spontaneous) folding events to the NMR-derived conformation are observed, allowing both structural and dynamical characterization of the folding TS. A common loop-like topology is observed in all the TS structures with native end-to-end and turn contacts, while the central segments of the strands are not in contact. Non-native sidechain contacts are present in the TS between the only tryptophan (W11) and the turn region (P7-G9). Prior to the TS the turn is found to be already locked by the W11 sidechain, while the ends are apart. Once the ends have also come into contact, the TS is reached. Finally, along the reactive folding paths the cooperative loss of the W11 non-native contacts and the formation of the central inter-strand native contacts lead to the peptide rapidly proceeding from the TS to the native state. The present results indicate a directed stepwise process to folding the peptide.     

Biochemical and molecular characterization of the oxidative branch of glycerol utilization by Citrobacter freundii.

Glycerol dehydrogenase (EC 1.1.1.6) and dihydroxyacetone kinase (EC 2.7.1.29) were purified from Citrobacter freundii. The dehydrogenase is a hexamer of a polypeptide of 43,000 Da. The enzyme exhibited a rather broad substrate specificity, but glycerol was the preferred substrate in the physiological direction. The apparent Kms of the enzyme for glycerol and NAD+ were 1.27 mM and 57 microM, respectively. The kinase is a dimer of a polypeptide of 57,000 Da. The enzyme was highly specific for the substrates dihydroxyacetone and ATP; the apparent Kms were 30 and 70 microM, respectively. The DNA region which contained the genes encoding glycerol dehydrogenase (dhaD) and dihydroxyacetone kinase (dhaK) was cloned and sequenced. Both genes were identified by N-terminal sequence comparison. The deduced dhaD gene product (365 amino acids) exhibited high degrees of homology to glycerol dehydrogenases from other organisms and less homology to type III alcohol dehydrogenases, whereas the dhaK gene product (552 amino acids) revealed no significant homology to any other protein in the databases. A large gene (dhaR) of 1,929 bp was found downstream from dhaD. The deduced gene product (641 amino acids) showed significant similarities to members of the sigma 54 bacterial enhancer-binding protein family.     

A highly strained nuclear conformation of the exportin Cse1p revealed by molecular dynamics simulations

To investigate the stability of the open nuclear state of the exportin Cse1p and its closing mechanism at the atomic level, we have performed multiple molecular dynamics simulations. The simulations revealed a strikingly fast transition of Cse1p from the open conformation to the closed cytoplasmic form, consistent with the proposal that Cse1p represents a "spring-loaded molecule." The structure of the ring-shaped state obtained in the simulations is remarkably close to the crystal structure of the cytoplasmic state, though the open nuclear structure was used as the only input. The conformational change is initially driven by release of strain due to RanGTP/importin-alpha binding. Subsequently, a stable closed state is formed, driven by attraction of electrostatically complementary interfaces. These results are consistent with and extend previous proposals. Reverse-charge and neutral mutants remained in an open state. The simulations predict a detailed reaction pathway and resolve the role of suggested hinge regions.     

Unfolding and extraction of a transmembrane alpha-helical peptide: dynamic force spectroscopy and molecular dynamics simulations

An atomic force microscope (AFM) was used to visualize CWALP(19)23 peptides ((+)H(3)N-ACAGAWWLALALALALALALWWA-COO(-)) inserted in gel-phase DPPC and DSPC bilayers. The peptides assemble in stable linear structures and domains. A model for the organization of the peptides is given from AFM images and a 20 ns molecular dynamics (MD) simulation. Gold-coated AFM cantilevers were used to extract single peptides from the bilayer through covalent bonding to the cystein residue. Experimental and simulated force curves show two distinct force maxima. In the simulations these two maxima correspond to the extraction of the two pairs of tryptophan residues from the membrane. Unfolding of the peptide precedes extraction of the second distal set of tryptophans. To probe the energies involved, AFM force curves were obtained from 10 to 10(4) nm/s and MD force curves were simulated with 10(8)-10(11) nm/s pulling velocities (V). The velocity relationship with the force, F, was fitted to two fluctuation adhesive potential models. The first assumes the pulling produces a constant bias in the potential and predicts an F approximately ln (V) relationship. The second takes into account the ramped bias that the linker feels as it is being driven out of the adhesion complex and scales as F approximately (ln V)2/3.     

DNA tension-modulated translocation and loop extrusion by SMC complexes revealed by molecular dynamics simulations

Structural Maintenance of Chromosomes (SMC) complexes play essential roles in genome organization across all domains of life. To determine how the activities of these large (≈50 nm) complexes are controlled by ATP binding and hydrolysis, we developed a molecular dynamics model that accounts for conformational motions of the SMC and DNA. The model combines DNA loop capture with an ATP-induced ‘power stroke’ to translocate the SMC complex along DNA. This process is sensitive to DNA tension: at low tension (0.1 pN), the model makes loop-capture steps of average 60 nm and up to 200 nm along DNA (larger than the complex itself), while at higher tension, a distinct inchworm-like translocation mode appears. By tethering DNA to an experimentally-observed additional binding site (‘safety belt’), the model SMC complex can perform loop extrusion (LE). The dependence of LE on DNA tension is distinct for fixed DNA tension vs. fixed DNA end points: LE reversal occurs above 0.5 pN for fixed tension, while LE stalling without reversal occurs at about 2 pN for fixed end points. Our model matches recent experimental results for condensin and cohesin, and makes testable predictions for how specific structural variations affect SMC function.     

Dynamics of the base of ribosomal A-site finger revealed by molecular dynamics simulations and Cryo-EM

Helix 38 (H38) of the large ribosomal subunit, with a length of 110 Å, reaches the small subunit through intersubunit bridge B1a. Previous cryo-EM studies revealed that the tip of H38 moves by more than 10 Å from the non-ratcheted to the ratcheted state of the ribosome while mutational studies implicated a key role of flexible H38 in attenuation of translocation and in dynamical signaling between ribosomal functional centers. We investigate a region including the elbow-shaped kink-turn (Kt-38) in the Haloarcula marismortui archaeal ribosome, and equivalently positioned elbows in three eubacterial species, located at the H38 base. We performed explicit solvent molecular dynamics simulations on the H38 elbows in all four species. They are formed by at first sight unrelated sequences resulting in diverse base interactions but built with the same overall topology, as shown by X-ray crystallography. The elbows display similar fluctuations and intrinsic flexibilities in simulations indicating that the eubacterial H38 elbows are structural and dynamical analogs of archaeal Kt-38. We suggest that this structural element plays a pivotal role in the large motions of H38 and may act as fulcrum for the abovementioned tip motion. The directional flexibility inferred from simulations correlates well with the cryo-EM results.     

The gating mechanism of the bacterial mechanosensitive channel MscL revealed by molecular dynamics simulations

One of the ultimate goals of the study on mechanosensitive (MS) channels is to understand the biophysical mechanisms of how the MS channel protein senses forces and how the sensed force induces channel gating. The bacterial MS channel MscL is an ideal subject to reach this goal owing to its resolved 3D protein structure in the closed state on the atomic scale and large amounts of electrophysiological data on its gating kinetics. However, the structural basis of the dynamic process from the closed to open states in MscL is not fully understood. In this study, we performed molecular dynamics (MD) simulations on the initial process of MscL opening in response to a tension increase in the lipid bilayer. To identify the tension-sensing site(s) in the channel protein, we calculated interaction energy between membrane lipids and candidate amino acids (AAs) facing the lipids. We found that Phe78 has a conspicuous interaction with the lipids, suggesting that Phe78 is the primary tension sensor of MscL. Increased membrane tension by membrane stretch dragged radially the inner (TM1) and outer (TM2) helices of MscL at Phe78, and the force was transmitted to the pentagon-shaped gate that is formed by the crossing of the neighboring TM1 helices in the inner leaflet of the bilayer. The radial dragging force induced radial sliding of the crossing portions, leading to a gate expansion. Calculated energy for this expansion is comparable to an experimentally estimated energy difference between the closed and the first subconductance state, suggesting that our model simulates the initial step toward the full opening of MscL. The model also successfully mimicked the behaviors of a gain of function mutant (G22N) and a loss of function mutant (F78N), strongly supporting that our MD model did simulate some essential biophysical aspects of the mechano-gating in MscL.     

Membrane-binding mechanism of the EEA1 FYVE domain revealed by multi-scale molecular dynamics simulations

Early Endosomal Antigen 1 (EEA1) is a key protein in endosomal trafficking and is implicated in both autoimmune and neurological diseases. The C-terminal FYVE domain of EEA1 binds endosomal membranes, which contain phosphatidylinositol-3-phosphate (PI(3)P). Although it is known that FYVE binds PI(3)P specifically, it has not previously been described of how FYVE attaches and binds to endosomal membranes. In this study, we employed both coarse-grained (CG) and atomistic (AT) molecular dynamics (MD) simulations to determine how FYVE binds to PI(3)P-containing membranes. CG-MD showed that the dominant membrane binding mode resembles the crystal structure of EEA1 FYVE domain in complex with inositol-1,3-diphospate (PDB ID 1JOC). FYVE, which is a homodimer, binds the membrane via a hinge mechanism, where the C-terminus of one monomer first attaches to the membrane, followed by the C-terminus of the other monomer. The estimated total binding energy is ~70 kJ/mol, of which 50–60 kJ/mol stems from specific PI(3)P-interactions. By AT-MD, we could partition the binding mode into two types: (i) adhesion by electrostatic FYVE-PI(3)P interaction, and (ii) insertion of amphipathic loops. The AT simulations also demonstrated flexibility within the FYVE homodimer between the C-terminal heads and coiled-coil stem. This leads to a dynamic model whereby the 200 nm long coiled coil attached to the FYVE domain dimer can amplify local hinge-bending motions such that the Rab5-binding domain at the other end of the coiled coil can explore an area of 0.1 μm² in the search for a second endosome with which to interact.     

Multi-step formation of a hemifusion diaphragm for vesicle fusion revealed by all-atom molecular dynamics simulations

Membrane fusion is essential for intracellular trafficking and virus infection, but the molecular mechanisms underlying the fusion process remain poorly understood. In this study, we employed all-atom molecular dynamics simulations to investigate the membrane fusion mechanism using vesicle models which were pre-bound by inter-vesicle Ca^2 +-lipid clusters to approximate Ca^2 +-catalyzed fusion. Our results show that the formation of the hemifusion diaphragm for vesicle fusion is a multi-step event. This result contrasts with the assumptions made in most continuum models. The neighboring hemifused states are separated by an energy barrier on the energy landscape. The hemifusion diaphragm is much thinner than the planar lipid bilayers. The thinning of the hemifusion diaphragm during its formation results in the opening of a fusion pore for vesicle fusion. This work provides new insights into the formation of the hemifusion diaphragm and thus increases understanding of the molecular mechanism of membrane fusion. This article is part of a Special Issue entitled: Membrane Structure and Function: Relevance in the Cell's Physiology, Pathology and Therapy.