A Quantum Dynamics Simulation of Molecular Vibrations by the Pechukas Method

Abstract

A quantum dynamics simulation of vibrations of molecules including transition from state to state is demonstrated based upon the Pechukas method. The method has been examined to clarify characteristics in relation to the simulation. It may present a lot of useful information of the vibrational relaxation and thermal excitation of the molecule.     

分子动力学模拟与分子生物力学

生物大分子的微观结构动力学决定其生物学功能,其力学-化学耦合规律是分子生物力学的重点关注方向。分子动力学模拟是耦合生物大分子力学-化学性质微观结构动力学基础的有效手段,其结果可用于预测结构-功能关系、指导实验设计和诠释实验结果。本文简要介绍了分子动力学模拟的方法学特点、基本工作原理及其在分子生物力学中的应用,并展望了未来可能的发展方向和应用前景。     

Multiplexed-Replica Exchange Molecular Dynamics Method for Protein Folding Simulation

Simulating protein folding thermodynamics starting purely from a protein sequence is a grand challenge of computational biology. Here, we present an algorithm to calculate a canonical distribution from molecular dynamics simulation of protein folding. This algorithm is based on the replica exchange method where the kinetic trapping problem is overcome by exchanging noninteracting replicas simulated at different temperatures. Our algorithm uses multiplexed-replicas with a number of independent molecular dynamics runs at each temperature. Exchanges of configurations between these multiplexed-replicas are also tried, rendering the algorithm applicable to large-scale distributed computing (i.e., highly heterogeneous parallel computers with processors having different computational power). We demonstrate the enhanced sampling of this algorithm by simulating the folding thermodynamics of a 23 amino acid miniprotein. We show that better convergence is achieved compared to constant temperature molecular dynamics simulation, with an efficient scaling to large number of computer processors. Indeed, this enhanced sampling results in (to our knowledge) the first example of a replica exchange algorithm that samples a folded structure starting from a completely unfolded state.     

A Method for Molecular Dynamics on Curved Surfaces

Dynamics simulations of constrained particles can greatly aid in understanding the temporal and spatial evolution of biological processes such as lateral transport along membranes and self-assembly of viruses. Most theoretical efforts in the field of diffusive transport have focused on solving the diffusion equation on curved surfaces, for which it is not tractable to incorporate particle interactions even though these play a crucial role in crowded systems. We show here that it is possible to take such interactions into account by combining standard constraint algorithms with the classical velocity Verlet scheme to perform molecular dynamics simulations of particles constrained to an arbitrarily curved surface. Furthermore, unlike Brownian dynamics schemes in local coordinates, our method is based on Cartesian coordinates, allowing for the reuse of many other standard tools without modifications, including parallelization through domain decomposition. We show that by applying the schemes to the Langevin equation for various surfaces, we obtain confined Brownian motion, which has direct applications to many biological and physical problems. Finally we present two practical examples that highlight the applicability of the method: 1) the influence of crowding and shape on the lateral diffusion of proteins in curved membranes; and 2) the self-assembly of a coarse-grained virus capsid protein model.     

A Novel Rabbit Monoclonal Antibody Platform To Dissect the Diverse Repertoire of Antibody Epitopes for HIV-1 Env Immunogen Design

The majority of available monoclonal antibodies (MAbs) in the current HIV vaccine field are generated from HIV-1-infected people. In contrast, preclinical immunogenicity studies have mainly focused on polyclonal antibody responses in experimental animals. Although rabbits have been widely used for antibody studies, there has been no report of using rabbit MAbs to dissect the specificity of antibody responses for AIDS vaccine development. Here we report on the production of a panel of 12 MAbs from a New Zealand White (NZW) rabbit that was immunized with an HIV-1 JR-FL gp120 DNA prime and protein boost vaccination regimen. These rabbit MAbs recognized a diverse repertoire of envelope (Env) epitopes ranging from the highly immunogenic V3 region to several previously underappreciated epitopes in the C1, C4, and C5 regions. Nine MAbs showed cross-reactivity to gp120s of clades other than clade B. Increased somatic mutation and extended CDR3 were observed with Ig genes of several molecularly cloned rabbit MAbs. Phylogenic tree analysis showed that the heavy chains of MAbs recognizing the same region on gp120 tend to segregate into an independent subtree. At least three rabbit MAbs showed neutralizing activities with various degrees of breadth and potency. The establishment of this rabbit MAb platform will significantly enhance our ability to test optimal designs of Env immunogens to gain a better understanding of the structural specificity and evolution process of Env-specific antibody responses elicited by candidate AIDS vaccines.     

From Molecular to Process Simulation: Novel Approaches to the Prediction of Phase Equilibria and PVT Behavior Based on Molecular/Quantum Mechanics and Molecular Dynamics Simulations

Abstract

Actually, in modern process simulators, more than 75% of the code implemented is dedicated to physical properties estimation, calculation and predictions. Data banks storing pure component parameters and binary interaction parameters for phase equilibrium calculations are extensively used and continuously implemented in actual process simulators. This gives an idea of the important role physical properties availability plays in process simulation.

In this paper we propose a new way for coupling molecular and process simulation. The basic machinery is to resort to molecular/quantum mechanics and molecular dynamics simulation techniques for generating the parameters of some equations of state that will subsequently be used for the prediction of phase equilibria and PVT behavior of small and polymeric molecules as well. This information, in turn, will be used as input in the process simulator, thus creating a final and well-defined bridge between molecular and process simulations in chemical engineering.     

基于快速沃尔什变换的分子子序列识别

提出了一种基于快速沃尔什变换的分子子序列识别的方法。这种方法不仅能快速识别出子序列并确定子序列的位置,而且极大地降低了CPU运行时间和计算复杂度,结合实例对这种方法进行了分析。     

Molecular Dynamics Simulation on the Connection Machine

An algorithm for the efficient calculation of macromolecular force fields on the Connection Machine is described. The full force field is separated into bond interactions and non-bonding interactions. Only the latter are implemented on the Connection Machine, the former, less computationally intensive tasks are performed by an existing, conventional molecular dynamics code on the front end. Parallelization of the evaluation of non-bonding interactions is achieved by the Replicated Systolic Loop algorithm introduced in this paper. The algorithm is a variant of the Systolic Loop scheme often used for the computation of 2-particle forces for the classical N-particle problem.     

Model-free Analysis of Protein Dynamics: Assessment of Accuracy and Model Selection Protocols Based on Molecular Dynamics Simulation

The popular model-free approach to analyze NMR relaxation measurements has been examined using artificial amide (15)N relaxation data sets generated from a 10 nanosecond molecular dynamics trajectory of a dihydrofolate reductase ternary complex in explicit water. With access to a detailed picture of the underlying internal motions, the efficacy of model-free analysis and impact of model selection protocols on the interpretation of NMR data can be studied. In the limit of uncorrelated global tumbling and internal motions, fitting the relaxation data to the model-free models can recover a significant amount of quantitative information on the internal dynamics. Despite a slight overestimation, the generalized order parameter is quite accurately determined. However, the model-free analysis appears to be insensitive to the presence of nanosecond time scale motions with relatively small magnitude. For such cases, the effective correlation time can be significantly underestimated. As a result, proteins appear to be more rigid than they really are. The model selection protocols have a major impact on the information one can reliably obtain. The commonly employed protocol based on step-up hypothesis testing has severe drawbacks of oversimplification and underfitting. The consequences are that the order parameter is more severely overestimated and the correlation time more severely underestimated. Instead, model selection based on Bayesian Information Criteria (BIC), recently introduced to the model-free analysis by d'Auvergne and Gooley (2003), provides a better balance between bias and variance. More appropriate models can be selected, leading to improved estimate of both the order parameter and correlation time. In addition, the computational cost is significantly reduced and subjective parameters such as the significance level are unnecessary.     

Molecular Dynamics Simulation of Alkali Silicates Based on the Quantum Mechanical Potential Surfaces

Abstract

Two potential parameter sets for alkali silicates were derived on the basis of ab-initio MO calculations. One is a model containing completely ionic alkali (model I), and another is that derived from cluster calculations (model II). These sets were tested against the crystal, glass, and liquid of metasilicates. The model II can reproduce these structures well under constant pressure conditions, and is found to be better than model I as a whole.     

A Molecular Dynamics Simulation Study of Rigid and Non-rigid Hard Dumb-bells

Abstract

We present a comparative study, using molecular dynamics, of systems of diatomic, hard dumb-bell, molecules in which the interatomic distance is either constrained to a fixed value or is allowed to vary freely between preset limits. A significant improvement in simulation effciency can be attained by allowing the bond length to vary. We find that thermodynamic properties, and some time correlation functions, are only slightly affected by the removal of the rigid bond-length constraint. The atomic velocity correlation function responds dramatically at short times to changes in the degree of non-rigidity, but at long times these differences are much less important.     

A Molecular Dynamics Simulation of an Aqueous Beryllium Chloride Solution

Abstract

A Molecular Dynamics simulation of a 1.1 molal aqueous BeCl₂ solution was performed with the flexible BJH model for water and a newly developed three-body potential for Be²⁺ -H₂O interactions derived from ab-initio calculations. The properties of the potential are discussed and radial distribution functions, angular distributions and dynamic properties of the solution like vibrational modes and hindered rotations are analyzed.     

鼠抗人纤维蛋白抗体单链Fv片段的人源化分子设计

产生免疫原性的残基都是位于蛋白表面的暴露残基,为了消除鼠抗体对人的免疫原性,利用表面再塑方法对本室克隆的鼠抗人纤维蛋白抗体单链Fv片断进行了人源化分子设计。首先确定了鼠及人Fv表面残基,在此基础上分析了鼠与人Fv间表面残基的差异,将有差异的鼠表面残基换成人的。提出了残基最高频率人源化及最相似链人源化两种人源化方案。人源化后鼠抗人纤维蛋白抗体单链Fv的结构经Profile-3D验证是合理的,置换的表面残基溶液可及性未变,且未影响CDRs的结构,应不会影响与纤维蛋白的亲和力,为鼠抗体人源化实验研究奠定了基础。     

A Molecular Dynamics Simulation Study of Coaxial Stacking in RNA

Abstract

We report on unrestrained molecular dynamics simulations of an RNA tetramer binding to a tetra-nucleotide overhang at the 5′-end of an RNA hairpin (nicked structure) and of the corresponding continuous hairpin with Na⁺ as counterions. The simulations lead to stable structures and in this way a structural model for the coaxially stacked RNA hairpin is generated. The stacking interface in the coaxially stacked nicked hairpin structure is characterized by a reduced twist and shift and a slightly increased propeller twist as compared to the continuous system. This leads to an increased overlap between C22 and G23 in the stacking interface of the nicked structure. In the simulations the continuous RNA hairpin has an almost straight helical axis. On the other hand, the corresponding axis for the nicked structure exhibits a marked kink of 39°. The stacking interface exhibits no increased flexibility as compared to the corresponding base pair step in the continuous structure.     

Grand Molecular Dynamics: A Method for Open Systems

We present a new molecular dynamics method for studying the dynamics of open systems. The method couples a classical system to a chemical potential reservior. In the formulation, following the extended system dynamics approach, we introduce a variable, v to represent the coupling to the chemical potential reservoir. The new variable governs the dynamics of the variation of number of particles in the system. The number of particles is determined by taking the integer part of v. The fractional part of the new variable is used to scale the potential energy and the kinetic energy of an additional particle: i.e., we introduce a fractional particle. We give the ansatz Lagrangians and equations of motion for both the isothermal and the adiabatic forms of grand molecular dynamics. The averages calculated over the trajectories generated by these equations of motion represent the classical grand canonical ensemble (μVT) and the constant chemical potential adiabatic ensemble (μVL) averages, respectively. The microcanonical phase space densities of the adiabatic and isothermal forms the molecular dynamics method are shown to be equivalent to adiabatic constant chemical potential ensemble, and grand canonical ensemble partition functions. We also discuss the extension to multi-component systems, molecular fluids, ionic solutions and the problems and solutions associated with the implementation of the method. The statistical expressions for thermodynamic functions such as specific heat; adiabatic bulk modulus, Grüneissen parameter and number fluctuations are derived. These expressions are used to analyse trajectories of constant chemical potential systems.     

Molecular Dynamics Simulation of Hexamine and Suberic Acid

In order to perform a molecular dynamics (MD) simulation of the incommensurate crystalline structure hexamethylenetetramine suberate (C 6 H 12 N 4 )(HOOC-(CH 2 ) 6 -COOH), we present in a first step the separate simulations of the crystalline structure of each of the two pure components, hexamethylenetetramine (HMT) and suberic acid. The domain decomposition parallel MD program ddgmq is used for this purpose. A second-generation consistent force field (CFF91) is employed to describe the interactions between atoms. Starting from experimental crystal structures, both pure components were heated from low to high temperatures. Our MD results show that the HMT system can be well represented by CFF91. In the case of suberic acid the layered structure of the crystal is largely preserved although deviations in the unit cell lengths from the experimental values are ～10%. Rather than attempt a complex re-parametrisation of CFF91 we chose to impose a fixed compensating external pressure tensor to correct for the deficiencies of the chosen force field. After optimising this compensating external pressure tensor at one temperature we find that experimental lattice constants and angles can be well reproduced over a range of temperatures.     

Humanization and Characterization of an Anti-Ricin Neutralization Monoclonal Antibody

Ricin is regarded as a high terrorist risk for the public due to its high toxicity and ease of production. Currently, there is no therapeutic or vaccine available against ricin. D9, a murine monoclonal antibody developed previously in our laboratory, can strongly neutralize ricin and is therefore a good candidate for humanization. Humanization of D9 variable regions was achieved by a complementarity-determining region grafting approach. The humanized D9 (hD9) variable regions were further grafted onto human heavy and light chain constant regions to assemble the complete antibody gene. A foot-and-mouth-disease virus-derived 2A self-processing sequence was introduced between heavy and light chain DNA sequences to cleave the recombinant protein into a functional full-length antibody molecule from a single open reading frame driven by a single promoter in an adenoviral vector. After expression in mammalian cells and purification, the hD9 was demonstrated to have equimolar expression of the full-length antibody heavy and light chains. More importantly, the hD9 exhibited high affinity to ricin with K_D of 1.63 nM, comparable to its parental murine D9 (2.55 nM). In a mouse model, intraperitoneal (i.p.) administration of hD9, at a low dose of 5 µg per mouse, 4 hours after the i.p. challenge with 5×LD50 ricin was found to rescue 100% of the mice. In addition, administered 6 hours post-challenge, hD9 could still rescue 50% of the mice. The hD9 has the potential to be used for prophylactic or therapeutic purposes against ricin poisoning.