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.     

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

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

基于B细胞表位制备抗肝细胞生成素的抗体

目的：基于B细胞表位制备抗肝细胞生成素（HPO）的抗体。方法：根据HPO的空间结构选择了2个候选B细胞表位,展示在T7噬菌体的表面,将提取的重组噬菌体免疫动物,采用ELISA法检测抗血清的效价,通过杂交瘤技术制备针对HPOC端表位的单克隆抗体。结果：2个候选B细胞表位KDGSCD和DGWKDGSC均能诱导抗相应表位多肽的多克隆抗体的产生,免疫6周后血清中抗体效价均达到1∶103,产生的抗体还能够特异识别HPO全蛋白;针对HPOC端表位KDGSCD的单克隆抗体也能识别HPO全蛋白,且具有良好的特异性。结论：基于T7噬菌体展示的B细胞表位可作为免疫原用于制备识别该B细胞表位来源的全蛋白质的抗体。     

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.     

一种水蛭素类融合蛋白与凝血酶作用的动力学模拟

通过将凝血酶抑制剂水蛭素的C端20肽片段嫁接到血小板结合蛋白AnnexinⅤ上,可以期望获得既具有抗凝血活性,同时又具有导向性的新型工程蛋白质分子.利用计算机辅助分子设计手段模拟了该融合蛋白的分子结构,并对该融合蛋白对凝血酶的抑制作用进行了分子动力学模拟,得到了支持上述想法的结果.     

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

提出了一种基于快速沃尔什变换的分子子序列识别的方法。这种方法不仅能快速识别出子序列并确定子序列的位置，而且极大地降低了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.     

Molecular Dynamics Simulation on a Parallel Computer

For the purpose of molecular dynamics simulations of large biopolymers we have built a parallel computer with a systolic loop architecture, based on Transputers as computational units, and have programmed it in occam II. The computational nodes of the computer are linked together in a systolic ring. The program based on this topology for large biopolymers increases its computational throughput nearly linearly with the number of computational nodes. The program developed is closely related to the simulation programs CHARMM and XPLOR, the input files required (force field, protein structure file, coordinates) and output files generated (sets of atomic coordinates representing dynamic trajectories and energies) are compatible with the corresponding files of these programs. Benchmark results of simulations of biopolymers comprising 66, 568, 3 634, 5 797 and 12 637 atoms are compared with XPLOR simulations on conventional computers (Cray, Convex, Vax). These results demonstrate that the software and hardware developed provide extremely cost effective biopolymer simulations. We present also a simulation (equilibrium of X-ray structure) of the complete photosynthetic reaction center of Rhodopseudomonas viridis (12 637 atoms). The simulation accounts for the Coulomb forces exactly, i.e. no cut-off had been assumed.     

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.     

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

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

Coarse-grained Modeling and Simulation of Actin Filament Behavior Based on Brownian Dynamics Method

The actin filament, which is the most abundant component of the cytoskeleton, plays important roles in fundamental cellular activities such as shape determination, cell motility, and mechanosensing. In each activity, the actin filament dynamically changes its structure by polymerization, depolymerization, and severing. These phenomena occur on the scales ranging from the dynamics of actin molecules to filament structural changes with its deformation due to the various forces, for example, by the membrane and solvent. To better understand the actin filament dynamics, it is important to focus on these scales and develop its mathematical model. Thus, the objectives of this study were to model and simulate actin filament polymerization, depolymerization, and severing based on the Brownian dynamics method. In the model, the actin monomers and the solvent were considered as globular particles and a continuum, respectively. The motion of the actin molecules was assumed to follow the Langevin equation. The polymerization, which increases the filament length, was determined by the distance between the center of the actin particle at the barbed end and actin particles in the solvent. The depolymerization, which decreases the filament length, was modeled such that the number of dissociation particles from the filament end per unit time was constant. In addition, the filament severing, in which one filament divides into two, was modeled to occur at an equal rate along the filament. Then, we simulated the actin filament dynamics using the developed model, and analyzed the filament elongation rate, its turnover, and the effects of filament severing on the polymerization and depolymerization. Results indicated that the model reproduced the linear dependence of the filament elongation on time, filament turnover process by polymerization and depolymerization, and acceleration of the polymerization and depolymerization by severing, which qualitatively agreed with those observed in experiments.     

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.     

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.