Maria Pechlaner

Last Name

Pechlaner

First Name

Maria

ORCID

0000-0002-1615-5600

Organisational unit

02021 - Inf.zentrum Chemie Biologie Pharmazie

Show all metadata

Search Results

Publications 1 - 10 of 14

On the use of intra-molecular distance and angle constraints to lengthen the time step in molecular and stochastic dynamics simulations of proteins
Pechlaner, Maria; van Gunsteren, Wilfred F. (2022)
Proteins: Structure, Function and Bioinformatics
Computer simulation of proteins in aqueous solution at the atomic level of resolution is still limited in time span and system size due to limited computing power available and thus employs a variety of time-saving techniques that trade some accuracy against computational effort. An example of such a time-saving technique is the application of constraints to particular degrees of freedom when integrating Newton's or Langevin's equations of motion in molecular dynamics (MD) or stochastic dynamics (SD) simulations, respectively. The application of bond-length constraints is standard practice in protein simulations and allows for a lengthening of the time step by a factor of three. Applying recently proposed algorithms to constrain bond angles or dihedral angles, it is investigated, using the protein trypsin inhibitor as test molecule, whether bond angles and dihedral angles involving hydrogen atoms or even stiff proper (torsional) dihedral angles as well as improper ones (maintaining particular tetrahedral or planar geometries) may be constrained without generating too many artificial side effects. Constraining the relative positions of the hydrogen atoms in the protein allows for a lengthening of the time step by a factor of two. Additionally constraining the improper dihedral angles and the stiff proper (torsional) dihedral angles in the protein does not allow for an increase of the MD or SD time step.
Algorithms to apply dihedral-angle constraints in molecular or stochastic dynamics simulations
Pechlaner, Maria; van Gunsteren, Wilfred F. (2020)
The Journal of Chemical Physics
Molecular structure refinement based on residual dipolar couplings using magnetic-field rotational sampling
Pechlaner, Maria; van Gunsteren, Wilfred F.; Smith, Lorna J.; et al. (2024)
The Journal of Chemical Physics
A method for structure refinement of molecules based on residual dipolar coupling (RDC) data is proposed. It calculates RDC values using magnetic-field rotational sampling of the rotational degrees of freedom of a molecule in conjunction with molecule-internal configurational sampling. By applying rotational sampling, as is occurring in the experiment, leading to observable RDCs, the method stays close to the experiment. It avoids the use of an alignment tensor and, therefore, the assumptions that the overall rotation of the molecule is decoupled from its internal motions and that the molecule is rigid. Two simple molecules, a relatively rigid and a very flexible cyclo-octane molecule with eight aliphatic side chains containing 24 united atoms, serve as so-called "toy model" test systems. The method demonstrates the influence of molecular flexibility, force-field dominance, and the number of RDC restraints available on the outcome of structure refinement based on RDCs. Magnetic-field rotational sampling is basically equivalent but more efficient than explicitly sampling the rotational degrees of freedom of the molecule. In addition, the performance of the method is less dependent on the number NRDC of measured RDC-values available. The restraining forces bias the overall orientation distribution of the molecule correctly. This study suggests that the information content of RDCs with respect to molecular structure is limited.
Saturation Mutagenesis by Efficient Free-Energy Calculation
Jandova, Zuzana; Fast, Daniel; Setz, Martina; et al. (2018)
Journal of Chemical Theory and Computation
Single-point mutations in proteins can greatly influence protein stability, binding affinity, protein function or its expression per se. Here, we present accurate and efficient predictions of the free energy of mutation of amino acids. We divided the complete mutational free energy into an uncharging step, which we approximate by a third-power fitting (TPF) approach, and an annihilation step, which we approximate using the one-step perturbation (OSP) method. As a diverse set of test systems, we computed the solvation free energy of all amino acid side chain analogues and obtained an excellent agreement with thermodynamic integration (TI) data. Moreover, we calculated mutational free energies in model tripeptides and established an efficient protocol involving a single reference state. Again, the approximate methods agreed excellently with the TI references, with a root-mean-square error of only 3.6 kJ/mol over 17 mutations. Our combined TPF+OSP approach does show not only a very good agreement but also a 2-fold higher efficiency than full blown TI calculations.
On the Effect of the Various Assumptions and Approximations used in Molecular Simulations on the Properties of Bio‐Molecular Systems: Overview and Perspective on Issues
van Gunsteren, Wilfred F.; Daura, Xavier; Fuchs, Patrick F.J.; et al. (2021)
ChemPhysChem
Computer simulations of molecular systems enable structure‐energy‐function relationships of molecular processes to be described at the sub‐atomic, atomic, supra‐atomic or supra‐molecular level and plays an increasingly important role in chemistry, biology and physics. To interpret the results of such simulations appropriately, the degree of uncertainty and potential errors affecting the calculated properties must be considered. Uncertainty and errors arise from (1) assumptions underlying the molecular model, force field and simulation algorithms, (2) approximations implicit in the interatomic interaction function (force field), or when integrating the equations of motion, (3) the chosen values of the parameters that determine the accuracy of the approximations used, and (4) the nature of the system and the property of interest. In this overview, advantages and shortcomings of assumptions and approximations commonly used when simulating bio‐molecular systems are considered. What the developers of bio‐molecular force fields and simulation software can do to facilitate and broaden research involving bio‐molecular simulations is also discussed. © 2020 Wiley
Structure and Conformational Dynamics of the Domain 5 RNA Hairpin of a Bacterial Group II Intron Revealed by Solution Nuclear Magnetic Resonance and Molecular Dynamics Simulations
Pechlaner, Maria; Sigel, R. K. O.; van Gunsteren, Wilfred F.; et al. (2013)
Biochemistry
Molecular Structure Refinement Based on Residual Dipolar Couplings: A Comparison of the Molecular Rotational-Sampling Method with the Alignment-Tensor Approach
Pechlaner, Maria; van Gunsteren, Wilfred F.; Smith, Lorna J.; et al. (2024)
Journal of Chemical Information and Modeling
In NMR experiments, residual dipolar couplings (RDCs) in a molecule can be measured by averaging the dipolar couplings (DCs) over the rotational motion of a molecule in an environment that induces a slight anisotropic orientation distribution of the molecule. Since the shape of the anisotropic distribution cannot be measured, it is standard practice to use a particular orientation distribution of the molecule with respect to the magnetic field, in the form of a so-called alignment tensor (AT), to calculate RDC-values for the molecule. Since the same alignment tensor is commonly used to calculate the different RDCs of a molecule, this approach rests on the assumption that the rotational motion of the molecule is decoupled from its internal motions and that the molecule is rigid. The validity of these two assumptions is investigated for a small, simple molecule, using a relatively rigid atomic interaction function or force field and a more flexible one. By simulating the molecule using an orientation-biasing force an anisotropic rotational distribution can be generated, for which RDCs can be obtained. Using these RDCs as target RDCs when applying one of the two approaches of structure refinement based on RDCs, it can be investigated how well the target RDCs are approximated in the RDC restraining and whether the corresponding nonuniform orientation distribution is reproduced. For the relatively rigid version of the molecule, the AT approach reproduces the target RDC-values, although the nonuniform orientation distribution of the angle theta(ab,H) between the vector r(ab)(->) connecting two atoms a and b in the molecule and the vector representing the direction of the magnetic field H--> as generated in the orientation-biasing simulation cannot be reproduced in the AT RDC-restraining simulation. For the relatively flexible version of the molecule, the AT approach fails to reproduce both the target RDC values and the nonuniform orientation distribution. For biomolecules with flexible parts, the application of the AT approach is thus not recommended. Instead, a method based on sampling of the rotational and internal degrees of freedom of the molecule should be applied in molecular structure determination or refinement based on measured RDCs.
Molecular dynamics simulation or structure refinement of proteins: are solvent molecules required? A case study using hen lysozyme
Pechlaner, Maria; van Gunsteren, Wilfred F.; Hansen, Niels; et al. (2022)
European Biophysics Journal
In protein simulation or structure refinement based on values of observable quantities measured in (aqueous) solution, solvent (water) molecules may be explicitly treated, omitted, or represented by a potential of mean-solvation-force term, depending on protein coordinates only, in the force field used. These three approaches are compared for hen egg white lysozyme (HEWL). This 129-residue non-spherical protein contains a variety of secondary-structure elements, and ample experimental data are available: 1630 atom-atom Nuclear Overhauser Enhancement (NOE) upper distance bounds, 213 (3) J-couplings and 200 S-2 order parameters. These data are used to compare the performance of the three approaches. It is found that a molecular dynamics (MD) simulation in explicit water approximates the experimental data much better than stochastic dynamics (SD) simulation in vacuo without or with a solvent-accessible-surface-area (SASA) implicit-solvation term added to the force field. This is due to the missing energetic and entropic contributions and hydrogen-bonding capacities of the water molecules and the missing dielectric screening effect of this high-permittivity solvent. Omission of explicit water molecules leads to compaction of the protein, an increased internal strain, distortion of exposed loop and turn regions and excessive intra-protein hydrogen bonding. As a consequence, the conformation and dynamics of groups on the surface of the protein, which may play a key role in protein-protein interactions or ligand or substrate binding, may be incorrectly modelled. It is thus recommended to include water molecules explicitly in structure refinement of proteins in aqueous solution based on nuclear magnetic resonance (NMR) or other experimentally measured data.
A Method to Derive Structural Information on Molecules from Residual Dipolar Coupling NMR Data
van Gunsteren, Wilfred F.; Pechlaner, Maria; Smith, Lorna J.; et al. (2022)
The Journal of Physical Chemistry B
A method for structure refinement of molecules based on residual dipolar coupling (RDC) data is proposed. It calculates RDC values using rotational and molecule-internal configurational sampling instead of the common refinement procedure that is based on the approximation of the nonuniform rotational distribution of the molecule by a single alignment tensor representing the average nonuniformity of this distribution. Using rotational sampling, as is occurring in the experiment leading to observable RDCs, the method stays close to the experiment. It avoids the use of an alignment tensor and thus the assumption that the overall rotation of the molecule is decoupled from its internal motions and that the molecule be rigid. Two simple molecules, two-united-atomic ethane and a cyclooctane molecule with eight side chains, containing 24 united atoms, serve as the so-called "toy model" test systems. The method demonstrates the influence of molecular flexibility and force-field deficiencies on the outcome of structure refinement based on RDCs. For a molecule of a given size (number of atoms Nat), there must be a sufficiently large number NRDC of measured RDC values available to allow the restraining forces to bias the overall orientation distribution of the molecule. If the ratio NRDC/Nat gets too small, the RDC-restraining forces will either not be strong enough to change the overall rotational direction of the molecule such that the target RDC values are approximated well or will be so strong that they induce a local deformation of the molecule. In the latter case, the size or inertia of the molecule hinders a restraining-induced overall rotation and the internal structure of the molecule is not strong enough to avoid local deformation due to the restraining forces.
Lessons learned from merging wet lab experiments with molecular simulation to improve mAb humanization
Schwaigerlehner, Linda; Pechlaner, Maria; Mayrhofer, Patrick; et al. (2018)
Protein Engineering, Design & Selection
Humanized monoclonal antibodies (mAbs) are among the most promising modern therapeutics, but defined engineering strategies are still not available. Antibody humanization often leads to a loss of affinity, as it is the case for our model antibody Ab2/3H6 (PDB entry 3BQU). Identifying appropriate back-to-mouse mutations is needed to restore binding affinity, but highly challenging. In order to get more insight, we have applied molecular dynamics simulations and correlated them to antibody binding and expression in wet lab experiments. In this study, we discuss six mAb variants and investigate a tyrosine conglomeration, an isopolar substitution and the improvement of antibody binding towards wildtype affinity. In the 3D structure of the mouse wildtype, residue R94h is surrounded by three tyrosines which form a so-called ‘tyrosine cage’. We demonstrate that the tyrosine cage has a supporting function for the CDRh3 loop conformation. The isopolar substitution is not able to mimic the function appropriately. Finally, we show that additional light chain mutations can restore binding to wildtype-comparable level, and also improve the expression of the mAb significantly. We conclude that the variable light chain of Ab2/3H6 is of underestimated importance for the interaction with its antigen mAb 2F5.

Publications 1 - 10 of 14

Maria Pechlaner

Last Name

First Name

ORCID

Organisational unit

Filters

Search Results