scholarly journals Thinking outside the Laboratory: Analyses of Antibody Structure and Dynamics within Different Solvent Environments in Molecular Dynamics (MD) Simulations

Antibodies ◽  
2018 ◽  
Vol 7 (3) ◽  
pp. 21 ◽  
Author(s):  
Mohammed Al Qaraghuli ◽  
Karina Kubiak-Ossowska ◽  
Paul Mulheran
Molecules ◽  
2018 ◽  
Vol 24 (1) ◽  
pp. 99 ◽  
Author(s):  
Siddharth Gautam ◽  
Tingting Liu ◽  
David Cole

Silicalite is an important nanoporous material that finds applications in several industries, including gas separation and catalysis. While the sorption, structure, and dynamics of several molecules confined in the pores of silicalite have been reported, most of these studies have been restricted to low pressures. Here we report a comparative study of sorption, structure, and dynamics of CO2 and ethane in silicalite at high pressures (up to 100 bar) using a combination of Monte Carlo (MC) and molecular dynamics (MD) simulations. The behavior of the two fluids is studied in terms of the simulated sorption isotherms, the positional and orientational distribution of sorbed molecules in silicalite, and their translational diffusion, vibrational spectra, and rotational motion. Both CO2 and ethane are found to exhibit orientational ordering in silicalite pores; however, at high pressures, while CO2 prefers to reside in the channel intersections, ethane molecules reside mostly in the sinusoidal channels. While CO2 exhibits a higher self-diffusion coefficient than ethane at low pressures, at high pressures, it becomes slower than ethane. Both CO2 and ethane exhibit rotational motion at two time scales. At both time scales, the rotational motion of ethane is faster. The differences observed here in the behavior of CO2 and ethane in silicalite pores can be seen as a consequence of an interplay of the kinetic diameter of the two molecules and the quadrupole moment of CO2.


2011 ◽  
Vol 2011 ◽  
pp. 1-9 ◽  
Author(s):  
Ronak Y. Patel ◽  
Petety V. Balaji

Glycolipids are important constituents of biological membranes, and understanding their structure and dynamics in lipid bilayers provides insights into their physiological and pathological roles. Experimental techniques have provided details into their behavior at model and biological membranes; however, computer simulations are needed to gain atomic level insights. This paper summarizes the insights obtained from MD simulations into the conformational and orientational dynamics of glycosphingolipids and their exposure, hydration, and hydrogen-bonding interactions in membrane environment. The organization of glycosphingolipids in raft-like membranes and their modulation of lipid membrane structure are also reviewed.


2016 ◽  
Vol 18 (15) ◽  
pp. 10507-10515 ◽  
Author(s):  
Vitaly V. Chaban

Thermodynamics, structure, and dynamics of diethyl sulfoxide (DESO) and ethyl methyl sulfoxide (EMSO) were investigated using ab initio calculations and non-polarizable potential based molecular dynamics (MD) simulations.


2018 ◽  
Author(s):  
Aoife M Harbison ◽  
Lorna P Brosnan ◽  
Keith Fenlon ◽  
Elisa Fadda

AbstractFc glycosylation of human immunoglobulins G (IgGs) is essential for their structural integrity and activity. Interestingly, the specific nature of the Fc glycoforms is known to modulate the IgG effector function. Indeed, while core-fucosylation of IgG Fc-glycans greatly affects the antibody-dependent cell-mediated cytotoxicity (ADCC) function, with obvious repercussions in the design of therapeutic antibodies, sialylation can reverse the antibody inflammatory response, and galactosylation levels have been linked to aging, to the onset of inflammation, and to the predisposition to rheumatoid arthritis. Within the framework of a structure-to-function relationship, we have studied the role of the N-glycan sequence on its intrinsic conformational propensity. Here we report the results of a systematic study, based on extensive molecular dynamics (MD) simulations in excess of 62 µs of cumulative simulation time, on the effect of sequence on the structure and dynamics of increasingly larger, complex biantennary N-glycoforms, i.e. from chitobiose to the larger N-glycan species commonly found in the Fc region of human IgGs. Our results show that while core fucosylation and sialylation do not affect the intrinsic dynamics of the isolated (unbound) N-glycans, galactosylation of the α(1-6) arm shifts dramatically its conformational equilibrium from an outstretched to a folded conformation. These findings are in agreement with and can help rationalize recent experimental evidence showing a differential recognition of positional isomers in glycan array data and also the preference of sialyltransferase for the more reachable, outstretched α(1-3) arm in both isolated and Fc-bound N-glycans.


2007 ◽  
Vol 5 (3) ◽  
pp. 635-671 ◽  
Author(s):  
Carmen Domene

AbstractDespite the complexity of ion-channels, MD simulations based on realistic all-atom models have become a powerful technique for providing accurate descriptions of the structure and dynamics of these systems, complementing and reinforcing experimental work. Successful multidisciplinary collaborations, progress in the experimental determination of three-dimensional structures of membrane proteins together with new algorithms for molecular simulations and the increasing speed and availability of supercomputers, have made possible a considerable progress in this area of biophysics. This review aims at highlighting some of the work in the area of potassium channels and molecular dynamics simulations where numerous fundamental questions about the structure, function, folding and dynamics of these systems remain as yet unresolved challenges.


2016 ◽  
Vol 5 (3) ◽  
Author(s):  
Jamal Hassan ◽  
Georgios Diamantopoulos ◽  
Dirar Homouz ◽  
Georgios Papavassiliou

AbstractStudying the properties of water confined in carbon nanotubes (CNTs) have gained a lot of interest in recent years due to the vast potential applications of systems in nanoscale liquid transport as well as biology functions. This article presents a comprehensive review of recent experimental and theoretical results using nuclear magnetic resonance (NMR) and molecular dynamics (MD) simulations. Different NMR methods including


2019 ◽  
Author(s):  
Enrico M. A. Fassi ◽  
Jacopo Sgrignani ◽  
Gianluca D’Agostino ◽  
Valentina Cecchinato ◽  
Maura Garofalo ◽  
...  

AbstractHigh-mobility Group Box 1 (HMGB1) is an abundant protein present in all mammalian cells and involved in several processes. During inflammation or tissue damage, HMGB1 is released in the extracellular space and, depending on its redox state, can form a heterocomplex with CXCL12. The heterocomplex acts exclusively on the chemokine receptor CXCR4 enhancing leukocyte recruitment.Here, we used multi-microsecond molecular dynamics (MD) simulations to elucidate the effect of the disulfide bond on the structure and dynamics of HMGB1.The results of the MD simulations show that the presence or lack of the disulfide bond between Cys23 and Cys45 modulates the conformational space explored by HMGB1, making the reduced protein more suitable to form a complex with CXCL12.


2012 ◽  
Vol 506 ◽  
pp. 258-261
Author(s):  
W. Soodsawang ◽  
T. Benchawan ◽  
U. Wichai ◽  
Y. Tantirungrotechai

Universal base is a man-made residue that can be incorporated into the DNA double strands without any discrimination against natural bases (A, C, G, T). The MD simulations with AMBER99 force field were employed to investigate the structure and dynamics of the modified 15-mer DNA duplexes containing carbazole-derived universal bases: carbazole (CBZ), 3,6-dicyanocarbazole (DCC), 3,6-dinitrocarbazole (DNC), and 3-nitro-6-cyanocarbazole (NCC), where X = CBZ, DCC, DNC, or NCC, respectively. The RMSD and B-factor of the modified DNAs backbones around the universal base unit fluctuate more than the reference sequence in the same position. The thermodynamic parameter for duplex stability was estimated by using MM-PBSA method. The averaged duplex formation free energy (ΔG) of all modified DNAs exhibited that the stability order was approximately DNC>NCC>CBZ>DCC, which differed from the reference sequence exceptional DNC unit. The averaged ΔG value of the DNC unit is very close to that of the reference sequence. This calculation indicated that the DNC unit can be considered as a good candidate for using as a universal base.


2021 ◽  
Author(s):  
Thomas Tarenzi ◽  
Marta Rigoli ◽  
Raffaello Potestio

The affinity of an antibody for its antigen is primarily determined by the specific sequence and structural arrangement of the complementarity-determining regions (CDRs). Recently, however, evidence has accumulated that points toward a nontrivial relation between the CDR and distal sites on the antibody structure: variations in the binding strengths have been observed upon mutating amino acids separated from the paratope by several nanometers, thus suggesting the existence of a communication network within antibodies whose extension and relevance might be deeper than insofar expected. In this work, we test this hypothesis by means of molecular dynamics (MD) simulations of the IgG4 monoclonal antibody pembrolizumab, an approved drug that targets the programmed cell death protein 1 (PD-1). The molecule is simulated in both the apo and holo states, totalling 4 μs of MD trajectory. The analysis of these simulations shows that the bound antibody explores a restricted range of conformations with respect to the apo one, and that the global conformation of the molecule correlates with that of the CDR; a pivotal role in this relationship is played by the relatively short hinge, which mechanically couples Fab and Fc domains. These results support the hypothesis that pembrolizumab behaves as a complex machinery, with a multi-scale hierarchy of global and local conformational changes that communicate with one another. The analysis pipeline developed in this work is general, and it can help shed further light on the mechanistic aspects of antibody function.


Author(s):  
Natalia Lukasheva ◽  
Dmitry Tolmachev ◽  
Hector Martinez-Seara ◽  
Mikko Karttunen

Electrostatic interactions have a determining role in conformational and dynamic behavior of polyelectrolyte molecules [1]. In this study, anionic polyelectrolyte molecules, poly(glutamic acid) (PGA) and poly(aspartic acid) (PASA), in water solution with the most commonly used K+ or Na+ counterions were investigated using atomistic molecular dynamics (MD) simulations. Seven common force fields, AMBER99SB-ILDN, AMBER14SB, AMBER-FB15, CHARMM22*, CHARMM27, CHARMM36m and OPLS-AA/L, both with their native parameters and with the non-bonded fix (NBFIX) and electronic continuum corrections (ECC) to were studied. These corrections have bene introduced to correct for the problem of overbinding of ions to the charged groups of polyelectrolytes. Physical properties, such as molecular sizes, local structure and dynamics, were studied using two types of common counterions, potassium and sodium. The results show that in some cases, the macroion size and dynamics depend strongly on the models (parameters) for the counterions due to strong overbinding of ions and charged side chain groups. The local structures and dynamics are more sensitive on dihedral angle parameterization resulting in a preference for defined monomer conformations amd the type of correction used.


Sign in / Sign up

Export Citation Format

Share Document